24小时热门版块排行榜    

查看: 3229  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 欧阳球飞 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 考研化学学硕调剂,一志愿985 +3 张vvvv 2026-03-15 5/250 2026-03-16 20:25 by 张vvvv
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 0703化学调剂 ,六级已过,有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见