matlab随机共振代码

变尺度随机共振代码，几年前用过的，

%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