%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?20:12??變尺度隨機共振\bianchidu.m

?????文件???????2053??2010-11-08?15:29??變尺度隨機共振\ode4.m

?????文件????????731??2010-12-14?21:32??變尺度隨機共振\odesr.m

?????文件????????111??2010-07-02?19:39??變尺度隨機共振\shuchuxinzaobi.m

?????文件?????????76??2010-12-14?21:29??變尺度隨機共振\sr11.m

?????文件????????693??2010-11-23?21:23??變尺度隨機共振\temp.m

?????目錄??????????0??2011-03-12?21:15??變尺度隨機共振

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

?????????????????4389????????????????????7