基于随机共振Runge-Kutta算法的强噪声信号滤波matlab仿真,含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); rng('default') fs=5; %参样频率 f=0.01; %信号频率 Ts=1/fs;%参样时间 h=1/fs; %时间步长 t=0:Ts:4095*Ts; D=0.04; %噪声强度 %双稳态系统参数 a=1; b=1; % s=0.3*sin(2*pi*0.01*t); %信号 s=0.3*sin(2*pi*0.01*t); %信号 x1=s+sqrt(2*D)*randn(size(t)); %噪声 %输入无噪信号傅立叶变换 y=fft(s,4096); pyy=y.*conj(y)/4096; ff=fs*(0:2048)/4096; figure(1); subplot(2,1,1);plot(t,s); title('输入无噪信号');xlabel('时间t/s');ylim([-0.5,0.5]);ylabel('信号幅度A'); subplot(2,1,2);plot(ff,pyy(1:2049)); xlabel('频率f/Hz');ylabel('频谱幅度');xlim([0,0.05]);title('输入无噪信号的频谱'); %输入加噪信号傅立叶变换 y=fft(x1,4096); pyy=y.*conj(y)/4096; ff=fs*(0:2048)/4096; figure(2) subplot(2,1,1);plot(t,x1); title('输入加噪噪信号');xlabel('时间t/s');ylabel('信号幅度A'); subplot(2,1,2);plot(ff,pyy(1:2049)); xlabel('频率f/Hz');ylabel('频谱幅度');xlim([0,0.05]);ylim([0,1500]);title('输入加噪信号的频谱'); %四阶龙格库塔法对双稳态输出信号求解 x=sr(a,b,h,x1); x=sr(a,b,h,x); x(1:100)=[]; x=x-mean(x); %输出信号求傅立叶变换 y=fft(x,4096); py=y.*conj(y)/4096; ff=fs*(0:2048)/4096; figure(3); subplot(2,1,1);plot(t(101:end),x); title('输出信号');xlabel('时间t/s');ylabel('信号幅度A'); subplot(2,1,2);plot(ff,py(1:2049)); xlabel('频率f/Hz');ylabel('频谱幅度');xlim([0,0.05]);ylim([0,5]);title('输出信号的频谱');