%x‘-a*x+b*x^3=Asin（2*pi*f*t）;
%畫出系統輸入輸出的時頻圖和相圖;
%ode4法;
global?a?b?A?f0;
a=1;b=1;A=0.3;f0=0.01;
fs=5;
t0=0;tf=800;
dt=1/fs;N=（tf-t0）/dt+1;
t=t0:dt:tf;
xin=A*sin（2*pi*f0*t）;
y1=fft（xinN）;
f1=（0:N/2-1）*fs/N;
yy1=abs（y1）*2/N;
subplot（421）plot（txin）title（‘輸入信號時域圖‘）;
subplot（422）plot（f1yy1（1:（N-1）/2））title（‘輸入信號頻譜圖‘）axis（[00.2502]）;
%t=[t0tf];
x0=0;%exp（-10）;
xout=ode4（@sr11tx0）;
xout=xout‘;
subplot（423）plot（txout）title（‘輸出信號時域圖‘）;
y=fft（xoutN）;
f=（0:N/2-1）*fs/N;
yy=abs（y）*2/N;
subplot（424）plot（fyy（1:（N-1）/2））title（‘輸出信號頻譜圖‘）axis（[00.2502]）;
xdot=（a*xout-b*xout.^3+xin）;
subplot（212）plot（xout（10000:end）xdot（10000:end））title（‘相圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????????76??2010-12-14?05:29??sr11.m

?????文件????????693??2010-11-23?05:23??temp.m

?????文件????????731??2010-12-14?05:32??odesr.m

?????文件????????111??2010-07-02?04:39??shuchuxinzaobi.m

-----------?---------??----------?-----??----

?????????????????1611????????????????????4