求解耦合时滞微分方程
求解耦合时滞微分方程
提问于 2016-07-25 21:42:38
我正在尝试在势能中模拟粒子的运动，并结合反馈。我遇到的问题是给反馈增加了延迟。因此，在Matlab中工作时，我有一个基本的脚本，如下所示：
clear all; close all;
qx = .3;
qy = .4;
qz = .5;
om_D = 1.2e8;
ep = 1e-10;
t = 0:1e-11:1e-5;
[t,rss]=ode45(@ion_test, t, [1e-7 .45 2e-7 -.45 0.5e-7 .45],[], qx, qy, qz, om_D, ep);其中，ion_test.m是一个函数文件，其内容如下：
function dr = ion_test(t,r, qx, qy, qz, om_D, ep)
dr = zeros(6,1);
dr(1) = r(2);
dr(2) = (2*qx*cos(om_D.*t) + ep.*2*r(3)*r(5));
dr(3) = r(4);
dr(4) = (2*qy*cos(om_D.*t) + ep.*2*r(1)*r(5));
dr(5) = r(6);
dr(6) = (2*qz*cos(om_D.*t) + ep.*2*r(1)*r(3));
end如你所见，r的每个坐标的微分方程是耦合的。我需要补充的是，每个方程都能够调用以前时间的其他坐标值，而不仅仅是它们的当前值。这样做的目的是构建延迟反馈机制。非常感谢您的帮助！
Stack Overflow用户
发布于 2016-07-30 08:10:48
我想我的问题的答案基本上是我需要使用dde23而不是ode45。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/38569393
复制相关文章
相似问题
腾讯云开发者
