24小时热门版块排行榜

返回列表

落月随风

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 237.6
红花: 1
帖子: 107
在线: 23.3小时
虫号: 1666450
注册: 2012-03-05
性别: MM
专业: 机械动力学

[求助] 紧急求助duffing算子的幅频分岔曲线及整个方程的分岔图matlab代码

求duffing算子的幅频分岔曲线及整个方程的分岔图matlab代码，如有不甚感激！！

回复此楼

» 猜你喜欢

如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有296人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有42人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有1人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复

1楼 2012-04-11 14:58:34

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

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
落月随风: 金币+10, ★★★很有帮助, 我想要分岔过度到混沌的图，不过还是感谢 2012-04-12 21:59:45

以前帮别人写的代码，不知是否你想要的
拷到editor中保存为.m文件，点击run运行
倍周期分岔图

CODE:

function AA

global mu w F;

mu=0.3;w=1.2;

cF=[0.20 0.27 0.289 0.293 0.32 0.35 0.37 0.4 0.65 0.7];

k=0;

for F=cF

    k=k+1;

[t,x]=ode45(@duffing,[0,2000],[2,0.01]);

x1=x(:,1);

count=find(t>1800);

Y=x1(count,:);

t1=t(count,:);

subplot(5,2,k)

plot(t1,Y)

 xlabel('t');

ylabel('x ');

s=['t=0 x=2 y=0.01 F=',num2str(F)];

title(s);

end

function dx=duffing(t,x)

global mu w F;

dx=[x(2);F*cos(w*t)-mu*x(2)+x(1)-x(1).^3];

杜芬(Duffing)方程 Poincare 截面图

CODE:

function aa

global mu w F;

mu=0.05;w=1;

F=7.5;

x0=[1;0.2];

tspan = 0:0.1*2*pi:50000;

[t,y]=ode45(@duffing,tspan,x0);

m=zeros(7000,1);

n=zeros(7000,1);

for i=1:7000

    m(i,1)=y(10*i,1);

    n(i,1)=y(10*i,2);

end

plot(m,n,'k.','markersize',1);

function dx=duffing(t,x)

global mu w F;

dx=[x(2);F*cos(w*t)-mu*x(2)-x(1).^3];