%计算Duffing方程主程序
clc
clear all
tic
%定义系统参数
r=2; %阻尼比
F0=58.5; %激励力
omega=1.2; %频率
k=1.8;
Ts=0.01; %取样时间
T0=500; %
t0=[0:Ts:T0];
x0=[0.001 0.001 0.001];
options=odeset('RelTol',1e-6,'AbsTol',[1e-8 1e-8 1e-8]);
%**************计算**********************************
[t,x]=ode45('Duffing_equation',t0,x0,options,r,F0,omega,k);
pp=10000; %画图点数
%***************吸引子的相图***********************************
figure(1);
plot(x(end-pp:end,1),x(end-pp:end,2));
xlabel('x1','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
ylabel('x2','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
hold on
figure(2);
plot(x(end-pp:end,1),x(end-pp:end,3));
xlabel('x1','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
ylabel('x3','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
hold on
figure(3);
plot(x(end-pp:end,2),x(end-pp:end,3));
xlabel('x2','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
ylabel('x3','fontsize',20,'fontname','times?new?roman','FontAngle','italic');
toc
xlswrite('g:\test.xls',x,'sheet1'); %保存数据
评论20