mychirp.m

线性调频信号及频谱在matlab中的应用

function mychirp
h=1.25*10^-8; %采样间隔时间1/80us
ptwidth=10^-6;%脉冲宽度1us
pcr=2*ptwidth;%脉冲周期2us
tcr=0:h:pcr;
pbwidth=10^7;%chirp信号带宽
y=chirp(tcr,0,ptwidth,pbwidth);%产生chirp信号
pluse=rectpuls(tcr,pcr).*y;%rectpuls产生矩形脉冲，占空比0.5pcr/(tcr-0.5pcr)
circle=10;%脉冲个数10
cirpluse=zeros(1,circle*length(pluse));%存放周期脉冲的初始数组
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;%对chirp序列进行fft变换，幅度减小100倍,2*10^4是FFT点数
f=tt/h;%各频率点
subplot(2,1,2)
%plot(ffty);
plot(f(1:floor(length(f)/2)),10*log10(abs(ffty(1:floor(length(f)/2)))));
%可以修改各参数值，得到不同的信号