匹配滤波器仿真matlab代码

clear all; close all; clc; f0 = 30; fs = 8*f0; k=1; t0=2; t=[0:0.01:2-0.01]; t1=[0:0.01:4-0.02]; signal0=3*cos(2*pi*f0*t+pi*k*t.^2)+randn(1,length(t)); h=3*cos(2*pi*f0*(t0-t)+pi*k*(t0-t).^2)+randn(1,length(t)); s0=conv2(signal0,h); figure subplot(211) plot(t,signal0) axis([0,2,-10,10]) xlabel('Time/s'); ylabel('Amplitude/v'); title('Original Signal with Noise'); subplot(212) plot(t1,s0) xlabel('Time/s'); ylabel('Amplitude/v'); title('Signal Passed Matched Filter'); figure subplot(211) signal_spectrum00 = fftshift(fft(signal0,fs)); index_spectrum00 = [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum00,abs(signal_spectrum00)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Original Signal with Noise'); subplot(212) signal_spectrum01 = fftshift(fft(s0,fs)); index_spectrum01= [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum01,abs(signal_spectrum01)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Signal Passed Matched Filter'); signal1=5*ones(1,200)+randn(1,200); s1=conv2(signal1,signal1); figure subplot(211) plot(t,signal1) xlabel('Time/s'); ylabel('Amplitude/v'); title('Original Signal with Noise'); axis([0,2,0,10]) subplot(212) plot(t1,s1) xlabel('Time/s'); ylabel('Amplitude/v'); title('Signal Passed Matched Filter'); figure subplot(211) signal_spectrum10 = fftshift(fft(signal1,fs)); index_spectrum10 = [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum10,abs(signal_spectrum10)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Original Signal with Noise'); subplot(212) signal_spectrum11 = fftshift(fft(s1,fs)); index_spectrum11= [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum11,abs(signal_spectrum11)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Signal Passed Matched Filter'); signal2=3*cos(2*pi*f0*t)+randn(1,length(t)); h=3*cos(2*pi*f0*(t0-t))+randn(1,length(t)); s2=conv2(signal2,h); figure subplot(211) plot(t,signal2) axis([0,2,-10,10]) xlabel('Time/s'); ylabel('Amplitude/v'); title('Original Signal with Noise'); subplot(212) plot(t1,s2) xlabel('Time/s'); ylabel('Amplitude/v'); title('Signal Passed Matched Filter'); figure subplot(211) signal_spectrum20 = fftshift(fft(signal2,fs)); index_spectrum20 = [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum20,abs(signal_spectrum20)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Original Signal with Noise'); subplot(212) signal_spectrum21 = fftshift(fft(s2,fs)); index_spectrum21= [-fs/4:0.5:(fs/2-1)/2]; plot(index_spectrum21,abs(signal_spectrum21)); xlabel('Frequency/Hz'); ylabel('Spectrum'); title('Signal Passed Matched Filter');