24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3231  |  回复: 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


分岔图问题-1
想出的激励频率分岔图.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

霜小妹二

木虫 (正式写手)

哈哈

【答案】应助回帖

★ ★ ★ ★ ★
欧阳球飞: 金币+5, 感谢回复 2016-04-22 18:38:11
我还原了你的公式。其实就是一个单自由度弹簧阻尼系统受到“非线性的力激励”,然后求其稳态响应(幅频响应曲线),不知道我说的对不对。如果真是这样,我建议还原方程,使用迭代求解方法。迭代方法有很多,你可以查一下。我最近在做非线性动力学问题,我的是“阻尼是高度非线性(幅值的六次多项式函数)”的,采用迭代方法,非常好求解。
没事儿就进行交流~
2楼2016-04-13 19:22:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

178202Z081

新虫 (初入文坛)

有人了解知道怎么出带随机扰动的分叉图吗??
3楼2019-03-08 21:50:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 欧阳球飞 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工求调剂! +20 RichLi_ 2026-04-06 20/1000 2026-04-07 09:20 by dick_runner
[考研] 085600材料与化工301分求调剂院校 +18 刺痛jk 2026-04-06 19/950 2026-04-07 00:02 by longlong0822
[考研] 材料求调剂 +15 一样YWY 2026-04-05 15/750 2026-04-06 22:42 by chenzhimin
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 材料调剂 +12 一样YWY 2026-04-05 13/650 2026-04-06 15:38 by lin-da
[考研] 336材料与化工085600求调剂 +9 水星记infp 2026-04-05 12/600 2026-04-06 11:46 by 水星记infp
[考研] 281求调剂 +8 椰子蘑菇 2026-04-06 8/400 2026-04-06 10:38 by 蓝云思雨
[考研] 求调剂 +5 wos666 2026-04-03 5/250 2026-04-06 10:13 by 蓝云思雨
[考研] 材料334求调剂 +19 Eecho# 2026-04-03 19/950 2026-04-06 08:37 by 小小树2024
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 080200学硕,277分,数一104,求带走! +7 瓶子PZ 2026-03-31 7/350 2026-04-05 17:49 by liucky
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 0854求调剂 +4 assdll 2026-04-03 4/200 2026-04-04 22:17 by hemengdong
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 求调剂,一志愿南京师范大学计算机专硕,初试373,六级通过, +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
信息提示
请填处理意见