标题:怎么画时滞微分方程的logistic映射的分叉图
取消只看楼主
冰点水lc
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2015-10-29
 问题点数:0 回复次数:0 
怎么画时滞微分方程的logistic映射的分叉图
请教各位牛人,小女子初学matlab,现遇到如下问题,还请各位多多指教,具体问题如下:
有一个时滞微分方程如图所示:[img=0,1]file:///C:/Users/lc-pc/AppData/Roaming/Tencent/Users/1289882567/QQ/WinTemp/RichOle/7L[_IV0BBB(%7D12LJL@S8@_9.png[/img]
    现在需要时间延迟tau不停变化,从0到20,从而画出x抽为tau,y轴分别为S,L,B。。。等的分叉图。我写得微分方程函数如下,后面的就不知道怎么写了。急求各位大神告知啊!!
function dy=ddefun1(t,y,z)
mu1=0.001;mu2=0.002;mu=0.003;beta1=0.3;beta2=0.62;beta3=0.35;gamma1=0.1;deta1=0.01;deta2=0.05;deta3=0.1;gamma2=0.65;gamma3=0.08;alpha=0.5;eta=0.15;
S=y(1);
L=y(2);
B=y(3);
Q=y(4);
R=y(5);
tauL=z(:,1);
dy=zeros(5,1);
  dy(1)=mu1  + gamma1*L + gamma2*B + gamma3*R - (beta1 *L + beta2*B + beta3)*S - deta1*S - eta* S - mu* S;
  dy(2)=(beta1*L +beta2*B + beta3 )*S - (gamma1    + mu  + deta2  + eta )*L- alpha*tauL(2);
  dy(3)=alpha*tauL(2) - (gamma2  + mu  + deta3  + eta )*B;
  dy(4)=deta1*S + deta2*L + deta3* B - (eta  + mu )*Q;
  dy(5)=mu2  + eta *(S + L + B + Q) - (mu  + gamma3 )*R;
end


[此贴子已经被作者于2015-10-29 22:43编辑过]

搜索更多相关主题的帖子: function matlab 小女子 如图所示 
2015-10-29 22:23



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




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

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