24小时热门版块排行榜

返回列表

落月随风

铜虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有152人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有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];