声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: lihaitao123

[稳定性与分岔] 软弹簧duffing方程不能用ode45么?

[复制链接]
发表于 2019-8-29 16:06 | 显示全部楼层
你好,楼主,我做了一下你这个软刚度非线性,用的是和你这个一样的数据,带入了你说的初值,但是还是出现警告,我把我的程序复制上,楼主你给帮忙看一下和你的区别。
回复 支持 反对
分享到:

使用道具 举报

发表于 2019-8-29 16:09 | 显示全部楼层
lihaitao123 发表于 2014-4-14 14:41
clear all
%数据的初始化
global  gamma omega f

function df=dafen(t,x)

global f2 w  alpha1 c  alpha3

df=[x(2);
    -alpha1*x(1)-alpha3*x(1)^3+f2*cos(w*t)-c*x(2)];
end


clear all
close all
global f2 w  alpha1 c  alpha3
%%%%%%%%%%%%%%%%%%%%%%%
f2=0.35;
alpha1=1;
alpha3=-1
c=0.4;
w=0.8;
tic %% tic,toc用于计时
  range=0.36;
k=0;
YY2=[];
for f2=range   
    period=2*pi/w; %
    step=2*pi/500/w;  %步长。
   y0=[0,-0.5 ];
    k=k+1

    tspan=[0:step:1000*period];
  [t,Y]=ode45(@dafen,tspan,y0);
  nn =floor(length(Y)*0.9);
  YY(k,:)=Y(nn:500:end,1);
    end
% figure
%  plot(range,YY,'k.','markersize',1);
figure
plot(Y(400000:end,1),Y(400000:end,2))
hold on
plot(Y(400000:500:end,1),Y(400000:500:end,2),'r.')
figure
plot(t(450000:end),Y(450000:end,1))
toc  
楼主,帮忙给我看一下,用的你的初值但是还是出现警告





发表于 2019-8-29 16:13 | 显示全部楼层
clear all
close all
global f2 w  alpha1 c  alpha3
f2=0.35;
alpha1=1;
alpha3=-1
c=0.4;
w=0.8;
tic %% tic,toc用于计时
%%%%%%%%%%%%8%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  range=7080%%%  参数范围
  range=0.35:0.001:0.39;
k=0;
YY2=[];
for f2=range   
    period=2*pi/w; %
    step=2*pi/500/w;  %步长。
   y0=[0 -0.50];
    k=k+1
  % discard the first 60 periodic data;
    %除去前面60个周期的数据,并将最后的结果作为下一次积分的初值
    tspan=[0:step:1000*period];
  [t,Y]=ode45(@dafen,tspan,y0);
  nn =floor(length(Y)*0.9);
  YY(k,:)=Y(nn:500:end,1);
    end
figure
  plot(range,YY,'k.','markersize',1);
figure
plot(Y(400000:end,1),Y(400000:end,2))
hold on
plot(Y(400000:500:end,1),Y(400000:500:end,2),'r.')
figure
plot(t(450000:end),Y(450000:end,1))
toc  


function df=dafen(t,x)

global f2 w  alpha1 c  alpha3

df=[x(2);
    -alpha1*x(1)-alpha3*x(1)^3+f2*cos(w*t)-c*x(2)];
end


楼主,帮忙给我看下和你的区别。是我哪里弄错了那
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-3-29 22:10 , Processed in 0.146986 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表