matlab 解延迟微分方程,matlab – 求解延迟微分方程(DDE)系统,约束给出非负解

你有一个棘手的问题,我不确定是否有一步到位的解决方案.我很乐意为愿意提供替代答案的任何人提供荣誉.

根据延迟的长度,一个选项是多次运行方程,每次迭代都将x的旧值传递给最新更新.

例如,假设你的延迟是一个小时.在第一个小时内,运行标记为NonNegative的ode45.将Value与Time参数一起存储到New矩阵中,然后再次运行算法.这次确保添加两个输入参数：旧的解决方案矩阵和旧的时间矩阵

dx = 1/(1+y*z) - x;

dy = (y*z)^2/(1+(y*z)^2) - y;

tindex = find(told>t,1) -1 % find the upper index which best approximates t

X = xold(tindex) + (xold(tindex+1)-xold(tindex))*(t-told(tindex))/(told(tindex+1)-told(tindex)) % or interpolation method of your choosing

dz = X - z;

现在洗,冲洗,重复.注意,X现在是一个准时间相关项,如ode45中的例3中所示.