function mychirp2
h=1.25*10^-8; %采樣間隔時間1/80us
t=0:h:2*10^-6;%
y=chirp(t,0,10^-6,10^7);%產生chirp信號
pluse=mypluse.*y;
figure(1);
subplot(2,1,1)
plot(t,cirpluse);
ffty=fft(y)/100;%對chirp信號進行fft變換
f=t/h;%各頻率點
subplot(2,1,2)
plot(f(1:floor(length(f)/2)),10*log10(abs(ffty(1:floor(length(f)/2)))));