Forum studentów IV roku EiT Politechniki Śląskiej

Gość

Witam Studentów 1 roku EiT. Każdy musi mieć kiedyś pierwszy temat i pierwszego posta więc oto mój. Wydaje mi się ze pewne grupy miały dzisiaj sprawdzian z Matlaba więc prosiłbym o podanie zadań (wraz z rozwiązaniami jeśli można), które się na nim znalazły. Z góry dzięki... :wink:

Gość

1. układ 3 równań z 3 niewiadomymi - wyliczyc te neiwiadome (np. metodą wyznaczikową)
2. podana funkcja (coś jak: a(x^3)+b(x^2)-x+7, obliczyć pierwiastki, miejsca zerwowe i wykres zrobić wszystkiego.

Rozwiązań chwilowo nie posiadam.

Gość

isohunt(dot)com/download/26498635/matlab <<= torrent do matlaba 2007b, może komuś się przyda [w miejsce (dot) wstawiamy kropke]

Zack de la Rocha · Dołączył: 10 Lis 2007 Posty: 187 Przeczytał: 0 tematów

Moglby ktos wrzucic rozwiazania na forum?

blue-code

ja tez bym bardzo prosił

mariuszl2

na pierwsze to chyba tak:
załóżmy że układ równań jest taki:
x1+2*x2+4*x3=2
3*x1+5*x2+7*x3=4
2*x1+8*x2+3*x3=3
to w matlabie piszemy:
a=[ 1 2 4; 3 5 7; 2 8 3]
b=[2; 4; 3]
a1=inv(a)
x=a1*b
( głowy nie daje że jest dobrze, trzeba by ręcznie przeliczyć i sprawdzić, ale to ponad me siły :wink: )

Tajok

Także dorzucam sie do prośby o rozwiąznia. Próbowałem w matlabie zrobić funkcję (zad2) i jak wprowadzam zewnętrznie a i b to wyskakuje błąd, w jaki sposób należy to zrobić?

mariuszl2

kawałek drugiego:
tworzymy m-plik z zawartością:
function y=funkcja1(a,b,x)
y=a*(x.^3)+b*(x.^2)-x+7
zapisujemy go jako funkcja1.m
w głównym oknie matlaba piszemy:
a=3
b=2
x=[-10:1:10]
funkcja1(a,b,x)
mamy teraz już obliczone igreki w przedziale -10 do 10
liczymy pierwiastki
pierw=roots([a b -1 7])
(w tym przypadku akurat wyszły też pierwiastki zespolone, tak jak poprzednio nie daje gwarancji że jest dobrze :wink: na razie nie daje wykresów bo mi cos nie wychodzi)[/b]

slomek

Moje rozwiazanie drugiego... nei wiem czy dobrze...

m-plik fun1
function [y]=fun1(a,b,x);
y=a*(x.^3)+b*(x.^2)-x+7;

m-plik zer1
function [y]=zer1 (x);
y=0;

Dalej w Command Window:
dla a=1, b=1:

fplot('[fun1(1,1,x),zer1(x)]',[-10;10])

pierw1=roots([1 1 -1 7])

hold on

plot(pierw1(1),0,'ro',pierw1(2),0,'bo',pierw1(3),0,'ro')

Rysuje wykres z zaznaczonymi pierwiastkami i "osią" OX, a przynajmniej powinno :wink:

GREY

function [y] = GRUPAE4_PONIKIERSKI_SPR1ZAD2(a,b,c,d)
x = -10:0.1:10;
y = a*x.*x.*x + b*x.*x + c*x + d;
plot(x,y);
pierwiastki = roots([a b c d]);
save Output.txt pierwiastki -ASCII

Gość

grey to nie dziala

GREY

Trudno XD Znaczy sie z tego co pamietam to dzialalo. Oczywiscie to musisz odpowiednio wywowalac jako funkcja z odpowiednimi parametrami.