Plot(t,x(:,1),'r'),grid,hold on,plot(t1,x1(:,1),'k'),

legend('x1','x1-d1'),title('Тест-система.d1=0.1')

% Проанализировать тест-систему по 2-й и 3-й координатам.

% Для графического вывода результатов по 2-й координате следует ввести plot в виде

% plot(t,x(:,2),'r'),grid,hold on,plot(t1,x1(:,2),'k'),

%legend('x2','x2-d1'),title('Тест-система.d1=0.1')

%% Для графического вывода результатов по 3-й координате следует ввести plot в виде

% plot(t,x(:,3),'r'),grid,hold on,plot(t1,x1(:,3),'k'),

%legend('x3','x3-d1'),title('Тест-система.d1=0.1')

Задание:

% Проанализировать тест-систему с относительной погрешностью, равной 0.01, 0.4, 0.3, 0.2

% Проанализировать тест-систему с помощью решателей ode45, ode113, сравнить результаты.

2. Задание абсолютной погрешности — AbsTol.

Абсолютная погрешность AbsTol контролирует разницу между ожидаемым решением и его действительным значением. AbsTol начинает проявляться, когда компонента (координата) решения становится неожидаемо большой. Влияние AbsTol зависит также от интервала и масштаба интегрирования. По умолчанию решатели дифференциальных уравнений Matlabустанавливают величину абсолютной погрешности, равную .

Пример3. Уравнение Ван-дер-Поля с заданной абсолютной погрешностью — 'AbsTol'.

% Используем имеющую М-функцию описания уравнения Ван-дер-Поля — van33

function f2=van33(t,X);

f2=[X(2);1*(1-X(1)^2)*X(2)-X(1)];

% Создадим М-сценарий (присвоить имя) с задаваемой абсолютной погрешностью по всем координатам и по %одной из возможных

T=[0 60];

X0=[3;0];

a1=odeset('AbsTol',0.5); %Скаляр 0.1 по всем координатам

[t1,x1]=ode23('van33',T,X0);

[t2,x2]=ode23('van33',T,X0,a1);

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