24小时热门版块排行榜

返回列表

落月随风

铜虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有188人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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