MATLAB16

17

17

x

całka

pochodne

tO

tmax

x0

tol

ślad

-    wiersz wektorów kolumnowych reprezentujących wartości rozwiązania kolejnych zmiennych stanu odpowiadające wartościom zmiennej niezależnej

-    nazwa procedury całkowania (ode23, cde45),

-    nazwa funkcji wyznaczania prawych stron równań stanu,

-    początek przedziału całkowania,

-    koniec przedziału całkowania,

wektor kolumnowy określający warunki początkowe zmiennych stanu,

-    opcjonalny parametr określający tolerancję porównywaną z lokalnym błędem obcięcia; domyślne wartości 1 e—3 dla ode23 lub le-6 dla cde45,

-    opcjonalny parametr, domyślna wartość 0; dla wartości rożnych od zera powoduje wypisywanie kolejnych kroków działania metody w oknie poleceń.

Jako przykład wykorzystania funkcji całkowania rozpatrzmy zadanie modelowania, polegające na wyznaczeniu rodziny odpowiedzi *,(/),/ e[0, 5] układu dynamicznego opisanego równaniami stanu (1), dla parametru cc= 0.1, 0.2 i 0.3 oraz zerowych warunków początkowych. Rozwiązanie zadania uzyskamy wprowadzając następujący ciąg poleceń: x0=zeros(2, 1) ; global ALFA.

ALFA=0.1;

[t,x]=ode45('rstanu', 0, 5, x0) plot(t,x (:, 1))

hołd on % ustawienie stanu dopisywania przebiegów ALFAfO.2;

[t,x]=ode45('rstanu', 0, 5, x0) plot(t,x(:, 1))

ALFA=0.3;

[t,x]=ode45('rstanu', 0, 5, x0) plot(t,x(:, 1))

Uwzględniając powtarzalność trzech ostatnich poleceń możemy przygotować M-plik funkcyjny (plik eksp.m), zawierający polecenia:

function eksp(a,x0)

%eksp - wykonanie eksperyment dla parametru ALFA=a %    i warunków początkowych x0;

%    rozwiązanie x(l) wykreślone jest w oknie graficznym

global ALFA

ALFA=a;