xx.m.txt

我编写的关于求解达芬方程的程序,各位请指教一下？

[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;%subplot(2,1,1)
plot(x1,x2) 
xlabel('displacement ')
ylabel('velocity')
title('phase diagram')
%subplot(2,1,2)
%plot(x1);
%xlabel('time')
%ylabel('displacement ')
%title('waveform diagram')