IMAG0150 (11)

'pr jest umteszczom blok inicjujący zmienne takie jak: przedział czasowy, warunki dokładność i krok obliczeń numerycznych. W kolejnym kroku jest rozwiązywane -y różnicriaowe crr pomocy standardowej procedury MATLABA - ode23. Jest to procedura jjjjjjr.cmega obliczania p lerwiastków równania różniczkowego II i HI-rzędu metodą Runge-Kutty. ^^ptctus: wejściowymi frakcji są zainicjowane wcześniej zmienne i nazwa m-pliku l^cjjDeąo. w którym zdefiniowane jest rozwiązywane równanie (w naszym przypadku Obliczone pierwiastki są graficznie przedstawiane w postaci wykresu fazowego za koci funkcji pioi słsźępm i: tworzenia wykresów płaskich. Funkcje title, xlabel, ylabel służą do ^-orowadzama kocnectaizy do wykresów. Zaznaczona na wydruku, jako komentarz, opcja efpmi powoduje proporcjonalne wyskalowanie osi, zaś opcja grid rysuje pomocniczą siatkę wykresu. Dokładny opis sMaim poszczególnych funkcji można uzyskać w programie MATLAB, wpisjąc poprompzcii komendę help i nazwę interesującej nas funkcji, np. help ode23.

M-piik rohmwum ma następuj ącą postać:

fc rofcars.m

fe T. Salacinski 4 3. Sycz, Warszawa 1997

% roharss.- fimkMjia opisującą równanie oscylatora harmonicznego i§- reharm.r. jest: wdane przez funkcje ODE23 w pliku oharm.m

Lp=zer©s(2,1j    % Zerowy wektor dwuelementowy

Łp(l) = y:2)    S y(l)=y y(2)=dy/dt

W pliku Tokarzm występuje definicja funkcji różniczkowej, zdefiniowana w sposób jęto zumiały dla procedur- >yde23 (patrz równania w rozdziale 2). W pierwszej linii (nie uwzględniając V - — »ntgrra i m-pliku funkcyjnego musi występować nagłówek definiujący funkcję:

Fuocdon pisra argumemów wejściowych] = nazwa funkcji [lista argumentów wyjściowych]

■falezy zwrocie uwagę, ze nazwa funkcji ma być taka sama jak nazwa m-pliku (bez rozszerzenia), w którym tę funkcję zapisano W kolejnych wierszach zdefiniowano zerowy wektor dwuelementowy. p'ement .p    opissgi pierwsza pochodną element yp(2) opisuje drugą pochodną, poprzez

>zw iktame względem rut; równania różniczkowego (patrz rozdział 2).

ro wy konaniu ord, _ ©trzy mano następującą trajektorię fazową (rys.4):

6