求解duffing系统用频闪法画分岔图的Matlab程序

雨露均沾

求解duffing系统用频闪法画分岔图的Matlab程序

william

function AA

global mu w F;
mu=0.3;w=1.2;
cF=;
k=0;
for F=cF
    k=k+1;
=ode45(@duffing,,);
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=;
title(s);
end

function dx=duffing(t,x)
global mu w F;
dx=;
杜芬(Duffing)方程 Poincare 截面图
function aa
global mu w F;
mu=0.05;w=1;
F=7.5;
x0=;
tspan = 0:0.1*2*pi:50000;
=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=;
网上找的程序  你自己看看行不行

敷衍会致命

william 发表于 2016-9-28 13:44
function AA

global mu w F;

这个程序好用吗，楼主找到了吗

雨露均沾

william 发表于 2016-9-28 13:44
function AA

global mu w F;

这个之前我也找过，但是好像不太有用欸

敷衍会致命

雨露均沾 发表于 2016-10-12 11:02
这个之前我也找过，但是好像不太有用欸

那你现在用哪个

雨露均沾

还在探索，感觉论坛里现在有好几篇蛮有用的

Catsayer

雨露均沾 发表于 2016-10-12 20:25
还在探索，感觉论坛里现在有好几篇蛮有用的

闪频法是什么意思？

雨露均沾

因为duffing是非自治系统，所以要用频闪法，差不多就是固定一个周期，然后打点的意思

Catsayer

雨露均沾 发表于 2016-10-14 16:19
因为duffing是非自治系统，所以要用频闪法，差不多就是固定一个周期，然后打点的意思

哦哦 一个周期一打点?