【信号检测】基于matlab随机共振弱信号检测【含Matlab源码 3588期】.zip

%% 随机共振stochastic resonace进行弱信号的检测的程序 clear all close all; clc; %% % % *******************************set parameter********************** f=9.251; % freqency of the signal Ts=1/f; h=1/f; D=0.31; a=1; b=1; fai=0.512*pi; % phase of the signal snr=-20; % SNR in dB(logarithmic) t=0:0.001:32; % time vector s=cos(2*pi*f*t+fai); % signal with the Amplify of 1 % % *******************************nosie generation********************** % tr=sqrt(2*D)*randn(size(t)); Amplifytr=max(abs(tr)); ampsig=Amplifytr*10^(snr/20); s=ampsig*s; % % *******************************signal with nosie generation*********** % x1=s+tr; % %% y=fft(x1,4096); pyy=y.*conj(y)/4096; ff=f*(0:2048)/4096; figure(1) subplot(2,2,1) plot(t,x1,'k') title('输入信号sn') xlabel('时间') ylabel('幅度') subplot(2,2,2) plot(ff,pyy(1:2049),'k') xlabel('频率 ') ylabel('幅度') title('输入信号sn的频谱') axis([-0.1 0.1 0 1000]) x=sr(a,b,h,x1); y=fft(x,4096); py=y.*conj(y)/4096; ff=f*(0:2048)/4096; subplot(2,2,3) plot(t,x,'k') title('输出信号x') xlabel('时间') ylabel('幅度') subplot(2,2,4) plot(ff,py(1:2049),'k') xlabel('频率y') ylabel('幅度') title('输出信号x的频谱') axis([-0.1 0.1 0 1000])