С заданной точностью и с параметрами
Решатели дифференциальных уравнений по умолчанию производят численное интегрирование с относительной погрешностью и вектором абсолютной погрешности со значением по всем компонентам вектора.
1. Установление заданной относительной погрешности RelTol— ODESET(odeset).
С помощью установления относительной погрешности RelTolконтролируется количество "правильных" цифр в решении дифференциального уравнения в соответствии с общей записью , где показатель степени P есть число контролируемых цифр в решении.
Пример1. Уравнение Ван-дер-Поля с заданной относительной погрешностью.
С помощью функции ODESET задаются опции решателя дифуравнений с помощью соответствующих строковых символов, которых всего может быть 18. Перечень строковых символов функции ODESET можно просмотреть из командной строки, набрав в ней ODESET
% Сформируем М-функцию для описания правых частей дифференциальных уравнений
function f2=van33(t,X);
f2=[X(2);1*(1-X(1)^2)*X(2)-X(1)];
% Сохранить под именем van33
% Создадим М-сценарий решения на основе ode23 и с относительной погрешностью, равной
% d1=0.1 и d2=0.0001
T=[0 60];
X0=[3;0];
d1=odeset('RelTol',0.1);
[t1,X]=ode23('van33',T,X0,d1);
d2=odeset('RelTol',0.0001);
[t2,x]=ode23('van33',T,X0,d2); 1 | 2 | 3 | 4 | 5 | 6 | 7 | Поиск по сайту:
|