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]）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-04-28?11:20??非線性雙穩系統隨機共振程序\
?????文件?????????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