bcdSR 变尺度随机共振代码

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

%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