Применение DDE23 к простому дифференциальному уравнению с задержкой - PullRequest
0 голосов
/ 11 ноября 2018

Здравствуйте, я пытаюсь применить dde23 к дифференциальному уравнению с задержкой

y '(t) = 3y (t - 2) с историей h (t) = 1, когда enter image description here на интервале enter image description here.

Я уже решил это уравнение, используя метод шагов для получения кусочного решения

enter image description here

Мне интересно сравнить это решение, используя численно полученный результат из dde23 на MATLAB, но у меня возникают проблемы с пониманием того, как изменить код по умолчанию, указанный для моей конкретной проблемы. Пока что я изменил код Wiley and Baker Example 23 по умолчанию:

sol = dde23(@ddex1de,[3, 2],@ddex1hist,[0, 3]);
figure;
plot(sol.x,sol.y)
title('MAT 5450 P5');
xlabel('time t');
ylabel('solution y');

function s = ddex1hist(t)
s = ones(3,1);
end

function dydt = ddex1de(t,y,Z)
ylag1 = Z(:,1);
ylag2 = Z(:,2);
dydt = [ ylag1(1)
ylag1(1) + ylag2(2)
y(2)          ];
end

Этот код создает графическое изображение, но я почти абсолютно уверен, что код неправильно адаптирован для моей конкретной проблемы. Я был бы признателен за любую помощь в изменении этого кода для моей проблемы, чтобы я мог сравнить мой ответ, полученный без использования MATLAB, спасибо.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...