时滞微分方程求解之三ddesd--变时滞
Posted studyer_domi
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了时滞微分方程求解之三ddesd--变时滞相关的知识,希望对你有一定的参考价值。
参考:
绘制含有时滞的微分方程组的曲线图像的matlab程序 – MATLAB中文论坛
ddesd说明:具有状态依赖时滞的 DDE- MATLAB & Simulink- MathWorks 中国
考虑下面的变时滞微分方程
试绘制上述含有时变时滞微分方程的图像
分析:该方程应该在某个t0时间之后成立,初始值必须是定义在t0之前的一个关于t的单值向量函数phi(t)。
我假设t0=0吧,phi(t)=[1;-1];
matlab程序:
function ddeex
t0 = 0;
tfinal = 5;
tspan = [t0, tfinal];
sol = ddesd(@ddefun,@delay,@hist,tspan);
plot(sol.x,sol.y,'.-')
xlabel('time')
ylabel('value')
legend('x','y')
function d=delay(t,y)
d=[t-sin(t);t-cos(t);t-2;t-1];
function h=hist(t)
h=[1;-1];
function dydt = ddefun(t,y,Z)
dydt = [1-0.05*Z(1,1)-0.04*Z(2,2); 2-0.02*Z(1,3)-0.7*Z(2,4)];
图像如下
以上是关于时滞微分方程求解之三ddesd--变时滞的主要内容,如果未能解决你的问题,请参考以下文章