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