Legend('X1-0.1','x1-0.0001'),title('Ур-е Ван-дер-Поля')

% М-сценарий сохранить под именем ret71.

% Формат записи функции odeset включает строковый служебный символ RelTol и задаваемую %величину относительной погрешности (0.1 для d1 и 0.2 для d2).

% d1 и d2 подставляются в соответствующие решатели (ode23).

Задание к примеру 1:

— выполнить с ode45, ode113;

— изменить коэффициент в М-функции van33: 0.3, 2.3;

— изменить для d1 относительную погрешность: 1, 0.5, 0.1, 0.01, 0.001, 0.0001 (для ode23).

Пример 2. Анализ нелинейной тест-системы с задаваемой относительной погрешностью вычислений в решателях MATLAB.

где — постоянные действительные числа.

2.1. Проинтегрируем тест-систему с коэффициентами и с относительной погрешностью, равной 0.1, 0.01, 0.001, 0.4, 0.3, 0.2

% Создадим М-функцию под именем dif21

function f22=dif21(t,x);

f22=[2.333*x(2)*x(3);-0.7789*(x(1)-0.1)*x(3);-0.7888*x(2)*x(1)];

% Создадим М-сценарий под именем syst21

t33=[0 20];

x0=[0;1;-1];

[t,x]=ode23('dif21',t33,x0);

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