24小时热门版块排行榜

返回列表

欧阳球飞

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.6
散金: 25
帖子: 62
在线: 42.6小时
虫号: 3799827
注册: 2015-04-10
专业: 机械工程

[求助] 分岔图问题已有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

回复此楼

» 猜你喜欢

仅有一篇一作相关文章，其余的都是通讯，好申请面上么？已经有6人回复
water research重复率要求已经有3人回复
金属材料论文润色/翻译怎么收费? 已经有118人回复
使用高温共聚焦显微镜原位观察熔盐-陶瓷的高温腐蚀反应过程已经有0人回复
使用高温激光共聚焦显微镜对GCr15轴承钢中碳化物高温行为的原位观察与分析已经有1人回复
求助TiAl合金的B2相pdf卡片已经有0人回复
紧急招博士生报考已经有20人回复
引领原子制造时代，成就未来材料之基--徐州国基原子制造有限公司已经有3人回复
工程热物理方向-考博-自荐已经有0人回复

1楼 2016-04-07 08:53:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

霜小妹二

木虫 (正式写手)

哈哈

应助: 154 (高中生)
金币: 1806.7
红花: 25
帖子: 531
在线: 66小时
虫号: 3237838
注册: 2014-05-27
专业: 机械动力学

【答案】应助回帖

★ ★ ★ ★ ★
欧阳球飞: 金币+5, 感谢回复 2016-04-22 18:38:11

我还原了你的公式。其实就是一个单自由度弹簧阻尼系统受到“非线性的力激励”，然后求其稳态响应（幅频响应曲线），不知道我说的对不对。如果真是这样，我建议还原方程，使用迭代求解方法。迭代方法有很多，你可以查一下。我最近在做非线性动力学问题，我的是“阻尼是高度非线性（幅值的六次多项式函数）”的，采用迭代方法，非常好求解。