LFM脉冲压缩MATLAB程序

function LFM（TB）

mu = B/T;
num = B*T;
points = linspace（-T/2T/2num）;
 
 real = cos（pi*mu.*points.^2）;
 imagepart = sin（pi*mu*points.^2）;
 st = exp（sqrt（-1）*pi*mu*points.^2）;
 %无加权
 ht = exp（-sqrt（-1）*pi*mu*points.^2）;
 yt = conv（stht）;
 yt（num）
 ytchu_db = 20*log10（abs（yt）/abs（yt（num）））;
 fft0 = fftshift（abs（fft（st）））;
 sampling_interval = T /num;
freqlimit = 0.5 / sampling_interval;
freq = linspace（-freqlimitfreqlimitnum）;

%海明加权
for j=1:1:num
 hthm（j） = exp（-sqrt（-1）*pi*mu*points（j）.^2）*（ 0.08 + 0.92*（cos（pi*points（j）/T））^2）;
end
 ythm = conv（sththm）;
 ythmchu_db = 20*log10（abs（ythm）/abs（ythm（num）））;



figure（1）
plot（pointsreal）
%axis（[-T/2 T/2 -1 1]）
grid
xlabel（‘T = 100us B = 1MHz 时间t/s‘）;
ylabel（‘实部‘）;

figure（2）
plot（pointsimagepart）
%axis（[-T/2 T/2 -1 1]）
grid
xlabel（‘T = 100us B = 1MHz  时间t/s‘）;
ylabel（‘虚部‘）;

figure（3）
plot（freqfft0）
%axis（[-T/2 T/2 -1 1]）
gr