xiyinpen.m

来自「Poincare映射图像有不同的做法」· M 代码 · 共 27 行

    % Duffing方程仿真图==相轨迹图 
% 方程在duffing.m文件中定义 
[t,x]=ode45(@duffing,[0,2800],[0,1.5]); 
x1=x(:,1); 
x2=x(:,2); 
x1=mod(x1,2*pi); 
x1(x1>pi)=x1(x1>pi)-2*pi; 
figure 
h=plot(x1,x2) 

%%%%%%%%%% 以下部分定义Duffing方程，此文件另存为duffing.m 
function dx=duffing(t,x) 
% 定义达芬Duffing方程： 
%  diff(x,2)+mu*diff(x,1)-x+x^3=F*cos(w*t) 
mu=0.1; 
F=2.09; 
w=0.5; 
dx=[x(2);F*cos(w*t)-mu*x(2)+x(1)-x(1).^3];

    
%function dy = vdp1000(t,y)
%dy = zeros(2,1);    % a column vector
%dy(1) = y(2);
%dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1);
%[T,Y] = ode15s(@vdp1000,[0 3000],[2 0]);
%plot(T,Y(:,1),'-o')