将DDE23应用于简单的延迟微分方程

时间:2018-11-10 22:49:43

标签: matlab differential-equations

您好,我正在尝试将dde23应用于延迟微分方程

y'(t)= 3y(t-2),当enter image description here在间隔enter image description here上时,历史h(t)= 1。

我已经使用分步方法求解了该方程,以获得分段解

enter image description here

我很想使用从MATLAB上dde23获得的数字结果来比较该解决方案,但是很难理解如何修改针对我的特定问题给出的默认代码。到目前为止,我已经修改了默认的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时得到的答案,谢谢。

0 个答案:

没有答案