Здравствуйте, я пытаюсь применить dde23
к дифференциальному уравнению с задержкой
y '(t) = 3y (t - 2) с историей h (t) = 1, когда на интервале .
Я уже решил это уравнение, используя метод шагов для получения кусочного решения
Мне интересно сравнить это решение, используя численно полученный результат из 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, спасибо.