plot(t,x),grid,title('Система с параметрами'),legend('x1','x2','x3')

% В записи функции ode23 символ [ ] указывает на опции (RelTol, AbsTol и др.), принятые по %умолчанию

Задание к примеру 4.

— Изменить параметр k1: 2, 1, 0.5, 0.2; (остальные исходные);

— Изменить параметр k2: -0.111,-0.333,-0.555,0.7789 (остальные исходные);

— Изменить параметр k3: -1.333, -0.333, -0.666, 0.7888 (остальные исходные);

— Изменить параметр a: -1, 1, 0, 0.5 (остальные исходные);

Произвести интегрирование системы при:

x0=[0;1;-1];

tt=[0,150];

k1=1

k2=-0.51

k3=-0.51

a=0

Пример 5. Интегрирование систем с циклическим изменением параметра.

% Используем М-функцию dif55

function f55=dif55(t,x,flag,k1,k2,k3,a);

Switch flag

case ''

f66=[k1*x(2)*x(3);k2*(x(1)+a)*x(3);k3*x(2)*x(1)];

End

% Видоизменим М-сценарий syst66 так, чтобы программно изменялся какой-либо параметр (например, k3) и происходило наложение графиков решения по заданной координате (например, по ). Новый М-сценарий будет с именем syst77.

x0=[0;1;-1];

t55=[0,100];

k1=1;

k2=-0.51;

a=0;

for i=-(0.508:0.001:0.511)

if i==-0.51 % Двойное равенство соответствует логической истине

[t,x]=ode23('dif55',t55,x0,[],k1,k2,-0.51,a);

Plot(t,x(:,2),'r'),title('Система с циклическим параметром')

Else

k3=i; % Без точки с запятой (;) выводятся значения k3

[t,x]=ode23('dif55',t55,x0,[],k1,k2,i,a);

Plot(t,x(:,2)),grid on,hold on,legend('','','x2-red','')

End

End

Задание к примеру 5.

— Расширить диапазон изменения параметра k3;

— Написать программу с изменением параметра k2;

— Написать программу с графическим выводом всех координат системы при k3=0.51 и только координаты при остальных значениях параметра k3.

Поиск по сайту: