随机共振的所有实例

关于随机共振的所有实例，有势阱含义的动画显示，有随机共振系统检测微弱信号，还有输入输出的信噪比计算，非常好的源程序

clear all
clc
fs=5;   %参样频率
f=0.1; %信号频率
Ts=1/fs;%参样时间
h=1/fs; %时间步长
t=0:Ts:4095*Ts;
D=0.02; %噪声强度，方差
A=0.01;
%双稳态系统参数
% a=0.4;
% b=0.3;
a=1;
b=1;
s=A*sin（2*pi*f*t）;         %信号
x1=s+sqrt（2*D）*randn（size（t））;  %噪声
%输入无噪信号傅立叶变换
y=fft（s4096）; 
pyy=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;
figure（1）;
subplot（211）;plot（ts）;
title（‘输入无噪信号‘）;xlabel（‘时间t/s‘）;ylim（[-0.50.5]）;ylabel（‘信号幅度A‘）;
subplot（212）;plot（ffpyy（1:2049））;
xlabel（‘频率f/Hz‘）;ylabel（‘频谱幅度‘）;xlim（[00.05]）;title（‘输入无噪信号的频谱‘）;
%输入加噪信号傅立叶变换
y=fft（x14096）; 
pyy=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;
figure（2）
subplot（211）;plot（tx1）;
title（‘输入加噪噪信号‘）;xlabel（‘时间t/s‘）;ylabel（‘信号幅度A‘）;
subplot（212）;plot（ffpyy（1:2049））;
xlabel（‘频率f/Hz‘）;ylabel（‘频谱幅度‘）;xlim（[00.05]）;ylim（[01500]）;title（‘输入加噪信号的频谱‘）;
%四阶龙格库塔法对双稳态输出信号求解
x=sr（abhx1）;
%输出信号求傅立叶变换
y=fft（x4096）; 
py=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;
figure（3）;
subplot（211）;plot（tx）;
title（‘输出信号‘）;xlabel（‘时间t/s‘）;ylabel（‘信号幅度A‘）;
subplot（212）;plot（ffpy（1:2049））;
xlabel（‘频率f/Hz‘）;ylabel（‘频谱幅度‘）;xlim（[00.05]）;ylim（[01500]）;title（‘输出信号的频谱‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-01-22 14:57  all_files_sr\
     文件        1280  2014-01-21 21:39  all_files_sr\detection.m
     文件          88  2014-01-21 15:52  all_files_sr\f.m
     文件         267  2014-01-21 15:51  all_files_sr\oula.m
     文件         386  2014-01-21 16:07  all_files_sr\powerf.m
     文件        2665  2014-01-21 20:56  all_files_sr\shijing.m
     文件        1076  2008-11-26 16:28  all_files_sr\snr.m
     文件         325  2014-01-21 16:03  all_files_sr\sr.m
     文件          72  2014-01-21 20:21  all_files_sr\test.m
     文件         156  2014-01-21 15:53  all_files_sr\xinhao1.m