%匹配濾波器實(shí)現(xiàn)
clear;
%%%%%%%%%%%%%%%%%%%%%%%%1
t=1:0.01:100
T=50;
s1=rectpuls(t,T);
%%%%%%%%%%%%%%%%%%%%%%%%2
f0=0.5;
w0=2*pi*f0;
s2=s1.*cos(w0*t);
%%%%%%%%%%%%%%%%%%%%%%%%3
k=1;
s3=s1.*cos(w0*t+0.5*k*t.^2);
%%%%%%%%%%%%%%%%%%%%%%%%4
N=100;
tao=T/N;
s4=0;
for index=1:N
    rect=rectpuls(index*tao,T);
    c=round(rand)*2-1;
    s4=s4+c*rect.*cos(w0*t);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%

for s_num=1:4
figure(s_num)
switch s_num
    case 1
        signal=s1;
    case 2
        signal=s2;
    case 3
        signal=s3;
    case 4
        signal=s4;
end

[m,n]=size(signal);
noise = wgn(m,n,10); % white Gaussian noise
signal_noise=signal+noise;
%%%%%%%%%%%%%%%
subplot(221)
plot(t,signal)
title('signal')
xlabel('t');
ylabel('s(t)');
%%%%%%%%%%%%%%%
subplot(222)
plot(t,signal_noise)
title('add Gaussian noise')
xlabel('t');
ylabel('n(t)')
%%%%%%%%%%%%%%%%%%%%%%%%%%
f1=fft(signal_noise);
f2=fft(fliplr(signal));
sig_mat=fliplr(signal);
output=ifft(f1.*f2);
subplot(223)
plot(t,sig_mat)
xlabel('t')
ylabel('h(t)')
title('matched filter')
%%%%%%%%%%%%%%%%%%%

subplot(224)
plot(t,output)
xlabel('t')
ylabel('y(t)')
title('output signal')
end