%x‘-a*x+b*x^3=Asin（2*pi*f*t）;
%畫出系統輸入輸出的時頻圖和相圖;
%ode4法;
global?a?b?A?f0;
a=1;b=1;A=0.406;f0=0.01;
fs=100;
t0=0;tf=400;
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）*fsr/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;
xout=ode4（@sr11tx0）;
xout=xout‘;
subplot（423）plot（txout）title（‘輸出信號時域圖‘）;
y=fft（xoutN）;
f=（0:N/2-1）*fsr/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（1000:end）xdot（1000:end））title（‘相圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????725??2010-11-16?04:12??bianchidu.m
?????文件????????2053??2010-11-07?23:29??ode4.m
?????文件?????????731??2010-12-14?05:32??odesr.m
?????文件?????????111??2010-07-02?04:39??shuchuxinzaobi.m
?????文件??????????76??2010-12-14?05:29??sr11.m
?????文件?????????693??2010-11-23?05:23??temp.m