变尺度随机共振的源代码

变尺度随机共振的源代码，包括调用的龙哥库塔方法函数代码，百分百管用

%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