24小时热门版块排行榜

返回列表

落月随风

铜虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有105人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有116人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复

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];