function mychirp
h=1.25*10^-8; %采樣間隔時(shí)間1/80us
ptwidth=10^-6;%脈沖寬度1us
pcr=2*ptwidth;%脈沖周期2us
tcr=0:h:pcr;
pbwidth=10^7;%chirp信號(hào)帶寬
y=chirp(tcr,0,ptwidth,pbwidth);%產(chǎn)生chirp信號(hào)
pluse=rectpuls(tcr,pcr).*y;%rectpuls產(chǎn)生矩形脈沖，占空比0.5pcr/(tcr-0.5pcr)
circle=10;%脈沖個(gè)數(shù)10
cirpluse=zeros(1,circle*length(pluse));%存放周期脈沖的初始數(shù)組
for i=1:circle
cirpluse(((i-1)*length(pluse)+1):i*length(pluse))=pluse;%得到周期脈沖
end
figure(1);
subplot(2,1,1)
tt=0:(circle*pcr)/((circle*length(pluse))-1):(circle*pcr);
plot(tt,cirpluse);
ffty=fft(cirpluse,2*10^4 )/100;%對(duì)chirp序列進(jìn)行fft變換，幅度減小100倍,2*10^4是FFT點(diǎn)數(shù)
f=tt/h;%各頻率點(diǎn)
subplot(2,1,2)
%plot(ffty);
plot(f(1:floor(length(f)/2)),10*log10(abs(ffty(1:floor(length(f)/2)))));
%可以修改各參數(shù)值，得到不同的信號(hào)