您好,我正在尝试将dde23
应用于延迟微分方程
y'(t)= 3y(t-2),当在间隔上时,历史h(t)= 1。
我已经使用分步方法求解了该方程,以获得分段解
我很想使用从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时得到的答案,谢谢。