TW=64;p=3;
s=dchirp(TW,p);
h=conj(s( end : -1: 1));
noise=sqrt(10)*randn(size(s));
s=s+noise;
figure(1)
T=16e-6;
N=p*TW;
t=[-T/2:T/N:T/2]; t(end)=[];
plot(t*1e6, real(s)), grid on
title('The Chirp signal in noisy')
xlabel('t(\mus)'),ylabel('The magnitude of s(n)')
figure(2)
ts=[-N/2:N/2]; ts(end)=[]; th=ts;
[y,ty]=conv_ma(s,ts,h,th);
y=conv(s,h);
ty=ty/(p*TW/T);
plot(ty*1e6, abs(y)), grid on
title('The output from the matched to Chirp signal')
xlabel('t(\mus)'),ylabel('The magnitude of y(n)')