TW=50;
p=10;
s=dchirp(TW,p);
h=conj(s(end:-1:1));
noise=sqrt(10)*randn(size(s));
s=s+noise;
s=delay(s,-150)+delay(s,200)+s;
T=16e-6;
N=p*TW;
t=[-T/2:T/N:T/2];t(end)=[];
subplot(311)
plot(t*1e6,real(s)),grid on
title('混有噪声的Chirp信号')
xlabel('时间(\mus)'),ylabel('幅度')
ts=[-N/2:N/2];ts(end)=[];th=ts;
[y,ty]=conv_m(s,ts,h,th);
ty=ty/(p*TW/T);
subplot(313)
plot(ty*1e6,abs(y)),grid on
title('噪声下Chirp信号的匹配滤波器输出')
xlabel('时间(\mus)'),ylabel('幅度')
subplot(312)
plot(t*1e6,noise),grid on
title('噪声')
xlabel('时间(\mus)'),ylabel('幅度')