随机共振用于弱信号检测的matlab程序

随机共振描述了过阻尼布朗粒子受周期性信号和随机噪声的共同作用下,在非线性双稳态系统中所发生的跃迁现象. 随机共振可用于弱信号的检测，编写的matlab程序，a_b_f是主程序

clear all
clc
fs=5;
Ts=1/fs;
h=1/fs;   %时间步长
t=0:Ts:4095*Ts;
D=0.31;
a=1;
b=1;
s=0.3*sin（2*pi*0.01*t）;
x1=s+sqrt（2*D）*randn（size（t））;


y=fft（x14096）; 
pyy=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;


figure（1）
subplot（221）
plot（tx1‘k‘）
title（‘Input signal sn‘）
xlabel（‘Time‘）
ylabel（‘Amplitude‘）
subplot（222）
plot（ffpyy（1:2049）‘k‘）
xlabel（‘Frequency ‘）
ylabel（‘Amplitude‘）
title（‘Spectrum of input signal sn‘）
axis（[-0.1 0.1 0 1000]）

x=sr（abhx1）;

y=fft（x4096）; 
py=y.*conj（y）/4096;
ff=fs*（0:2048）/4096;


subplot（223）
plot（tx‘k‘）
title（‘Output signal x‘）
xlabel（‘Time‘）
ylabel（‘Amplitude‘）
subplot（224）
plot（ffpy（1:2049）‘k‘）
xlabel（‘Frequency‘）
ylabel（‘Amplitude‘）
title（‘Spectrum of output signal x‘）
axis（[-0.1 0.1 0 1000]）

% [tfrtf]=tfrsp（x‘）;
% [tfrrtfrhat]=tfrrmsc（x‘）;
% % colorbar;
% % figure（1）
% % contour（tabs（f）tfr）;
% figure（5）
% contour（tabs（f）*10000tfr）;
% axis（[0 1024 0 250]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        982  2008-06-11 10:18  随机共振弱信号检测\a_b_f.m

     文件      24576  2008-06-11 16:28  随机共振弱信号检测\sr.doc

     文件        647  2007-10-26 16:35  随机共振弱信号检测\sr.m

     目录          0  2011-04-28 11:20  随机共振弱信号检测

----------- ---------  ---------- -----  ----

                26205                    4