%用隨機共振方法檢測微弱信號，觀察奇倍頻現象
clear?all
clc
fs=50;??%采樣頻率fs=50*f1
Ts=1/fs;%采樣時間
h=1/fs;?%時間步長
t=0:Ts:4095*Ts;
D=0.2;?%噪聲強度
A=1;%信號幅度

%雙穩態系統參數
a=1;
b=1;
f1=0.1;
s=A*sin（2*pi*f1*t）;???%信號
noise=sqrt（2*D）*randn（size（t））;??%噪聲
x1=s+noise;

xn=RK（abhx1）;%四階龍格庫塔方法求輸出信號
%輸出信號求傅立葉變換求頻譜
y=f