| 查看: 2215 | 回复: 2 | |||
[交流]
【求助】用Matlab画常微分方程的分支图
|
|
我画分支图常常先差分,然后运行很长时间取最后的几百步就可以了,例如以下 clear; T=1000;dt=0.05;a1=5.0;a2=0.1;b2=2.0;d1=0.4;d2=0.01;e=1.0;r=0; X=zeros(1,T);Y=zeros(1,T);Z=zeros(1,T); X(1)=0.4;Y(1)=0.6;Z(1)=5; hold on for b1=2:0.02:9 for t=2:T X(t)=X(t-1)+(X(t-1)*(1-X(t-1))-a1*X(t-1)*Y(t-1)/(1+b1*X(t-1)))*dt; Y(t)=Y(t-1)+(a1*X(t-1)*Y(t-1)/(1+b1*X(t-1))-a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d1*Y(t-1))*dt; Z(t)=Z(t-1)+(a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d2*Z(t-1))*dt; end plot(b1,Z(T-500:T),'k.','MarkerSize',4) end 但这样做的结果是周期和混沌无法分辨,有人用极值命令做,好像可以解决这方面的问题。 请教一下高手,帮我把上面程序‘plot(b1,Z(T-500:T),'k.','MarkerSize',4) ’这部分用极大值命令代替一下,万分感谢! |
» 猜你喜欢
面上合作单位盖章
已经有5人回复
有多少人是今天查系统知道结果的?
已经有17人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有8人回复
导师吐槽:我怎么摊上了这么个极品研究生!
已经有7人回复
申请删除本帖
已经有7人回复
基金不中,共勉
已经有10人回复
国自然面上复盘~欢迎讨论
已经有11人回复
基金未中,这种答复是模板吗?
已经有6人回复
看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐,对您同样有帮助:
Matlab数值求解二阶常微分方程
已经有9人回复
如何用 数值解法 求 二阶微分方程的本征值
已经有18人回复
matlab求解常微分方程
已经有6人回复
matlab处理常微分方程作图问题
已经有7人回复
求解两个简单非线性常微分方程~~
已经有4人回复
非线性方程组的解法
已经有5人回复
一阶常微分方程组,matlab进行数值求解
已经有3人回复
微分方程组的稳定性分析,急!
已经有10人回复
高手进来用matlab求解常微分方程,十分感谢!
已经有17人回复
matlab 线性常微分方程求解
已经有5人回复
matlab-常微分方程参数估计
已经有12人回复
matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
数值方法求解偏微分方程的问题
已经有4人回复
常微分方程(ODEs)的MATLAB数值解法
已经有119人回复
【求助】matlab解常微分方程,暂停和继续
已经有5人回复
» 抢金币啦!回帖就可以得到:
三无人员面上一把中散金币
+2/1882
发点金币玩
+5/525
面上得中
+1/462
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/199
大连理工大学王铁胜副教授招收2027年入学博士生
+1/189
中国基,散金币!!!!!!
+5/134
东华大学材料科学与工程学院丁奇老师招收2026级硕士研究生
+1/87
欢迎报考中科院苏州纳米所的氢能、燃料电池方向的博士生!(长期有效)
+1/54
欢迎保送和统考中科院苏州纳米所的氢能、燃料电池方向的研究生!(长期有效)
+1/52
中国科学院苏州纳米技术与纳米仿生研究所 仿生自适应智能系统实验室招收联培生
+1/36
国家蛋白质科学中心(北京),招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/15
国家蛋白质科学中心(北京),招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/14
2027材料申博自荐C9本国科大硕
+1/12
新审查指南堵死了"一案双申"的老路,实用新型和发明只能二选一了?
+1/8
项目助理招聘启事
+1/4
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/4
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/4
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/2
细胞治疗全流程关键检测技术与研究策略
+1/1
2楼2010-11-14 02:09:12
3楼2010-11-21 14:25:57









回复此楼
投票: