matlab随机共振程序

%用随机共振方法检测微弱信号，观察奇倍频现象
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