| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 2050 人关注过本帖
标题:怎么画时滞微分方程的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
快速回复:怎么画时滞微分方程的logistic映射的分叉图
数据加载中...
 
   



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

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