标题:求助:对带时滞的方程组仿真
只看楼主
zhaoyuan325
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2007-5-27
 问题点数:0 回复次数:0 
求助:对带时滞的方程组仿真

对带时滞的方程组仿真,该程序不报错也不显示图片,是工具有问题还是程序有问题,烦高手指点!

sol=dde23(@twodelay,[0.05 2],[1.4;2.5],[0,300]);
figure (1);
plot(sol.x,sol.y(1,:),sol.x,sol.y(2,:))
figure(2);
plot(sol.y(1,:),sol.y(2,:))

function y=twodelay(t,x,Z)
y=ones(2,1);
xlag1=Z(:,1);
xlag2=Z(:,2);
K=0.02;
k=0.1;
bt=1.77;
gm1=0.1;
gm2=2.4;
A=20;
y(1)=-[K+k*x(2)+bt*x(1)]*x(1)+2*exp(-gm1*0.5)*bt*xlag1(1)*xlag1(1);
y(2)=-gm2*x(2)+A*k*xlag2(2)*xlag2(2);

搜索更多相关主题的帖子: 时滞 方程 仿真 
2007-05-31 09:49



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-143746-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.057910 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved