| 查看: 3198 | 回复: 2 | ||
[求助]
分岔图问题 已有1人参与
|
|
想得到系统随激励频率wh变化的分岔图，一直出不对，不知道程序哪里出了问题，请精通的大神指点一二。 感激不尽。程序、方程见附件。 主要问题：1.分岔图问题，图出不对，请帮忙看一下程序；2.程序中(1-e*cos(wn*t))*fx，如何绘制图它随时间变化的曲线。 调用函数： function dx = sucheng0328(t,x,b,Fm,Fa,wn,s,e) if x(1)>b fx=x(1)-b;%b 为无量纲后的比值 else if x(1)<-b fx=x(1)+b; else fx=0; end end dx = [x(2); Fm + Fa*cos(wn*t)-2*s*x(2)-(1-e*cos(wn*t))*fx; ]; 求分岔图的程序 tic clear all; clc b = 1; Fm = 0.05; Fa = 0.1; wn = 0.01:0.01:2; s = 0.07; e = 0.1; options=odeset('RelTol',1e-7); tt=2*pi./wn; hold on for j=1:length(wn) [t,y]=ode45(@sucheng0328,0:tt/100:80*tt,[0;0],options,b,Fm,Fa,wn(j),s,e); i = 5000:100:6000; plot(wn(j),y(i,1),'k-') end toc  方程及已知条件.png  想出的激励频率分岔图.png |
回复此楼» 猜你喜欢
【26年入学】计算机视觉/电子信息/诚招博士生 光电教育部重点实验室
已经有2人回复
-
MOCVD可逆气相反应求助
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有180人回复
-
【26年入学】红外半导体/电子信息/嵌入式视觉/诚招博士生 光电教育部重点实验室
已经有0人回复
-
华北电力大学王鹏课题组诚招博士/硕士
已经有3人回复
-
CSC难上岸，25年东京大学，早稻田，京都大学多Lab大量招募日本全奖博士生！
已经有6人回复
-
【26年入学】红外半导体/电子信息/嵌入式视觉/诚招博士生 光电教育部重点实验室
已经有0人回复
-
读博入学申请
已经有1人回复
-
考核制博士自荐
已经有30人回复
-
2026 博士自荐-机器人机构学方向
已经有12人回复
-
湖南大学JQ团队段辉高教授课题组招光学、算法博士啦！2026级，2个名额，等你加入！
已经有11人回复
霜小妹二
木虫 (正式写手)
哈哈
- 应助: 154 (高中生)
- 金币: 1806.7
- 红花: 25
- 帖子: 531
- 在线: 66小时
- 虫号: 3237838
- 注册: 2014-05-27
- 专业: 机械动力学
2楼2016-04-13 19:22:50
178202Z081
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 60.8
- 散金: 5
- 红花: 1
- 帖子: 47
- 在线: 37.3小时
- 虫号: 10056727
- 注册: 2018-08-28
- 专业: 动力学与控制
3楼2019-03-08 21:50:55