雷达脉冲压缩仿真程序

tp=40e-6;%脉压信号时宽
B=5e+6;%脉压信号带宽

D=tp*B; %压缩比

fs=10e+6;%采样频率
pulse_width=1/fs; %脉压后宽度–(1us,150m)—0.2us-30m–
PC_N=round(tp*fs);%chirp信号采样点数

%%%%%%% LFM signals %%%%%%%%
t=[-tp/2:tp/PC_N:tp/2]; %采样时间点

TranS=exp(jpiB/tp.*t.^2); % 1 x PC_N,LFM 信号

figure(5)
plot(real(TranS));
title(‘LFM 信号实部’);

figure(8)
plot(imag(TranS));
title(‘LFM 信号虚部’);

Hs=TranS;

Hs_length=length(Hs);
for no=1:Hs_length %LFM信号的时间倒置实现
Hs(no)=TranS(Hs_length-no+1);
end

%%%%%%% echol signals %%%%%%%%
RecS=[zeros(1,50),TranS,zeros(1,50)];% 构建回波测试信号
% RecS=TranS;

figure(10)
plot(abs(RecS));
title(‘回波信号幅度值’);
%%%%%%% 时域脉冲压缩 %%%%%%%%
PulComS=conv(RecS,conj(Hs));

figure(20)
% plot(20log10(abs(PulComS)));
plot(20log10(abs(PulComS)/max(abs(PulComS))));

axis([3800 4200 -50 0])
% ylim([-50,0]);
ylabel(‘dB’);
title(‘时域实现：未加权脉冲压缩结果’);

%%%%%%% 频域脉冲压缩 %%%%%%%%
[PulComS2]=pulse_com(B,fs,Hs,RecS);

figure(30)
% plot(20log10(abs(PulComS)));
plot(20log10(abs(PulComS2)/max(abs(PulComS2))));
% axis([3800 4200 -50 0])
ylabel(‘dB’);
ylim([-50,0]);
title(‘频域实现：加权抑制副瓣的脉冲压缩结果’);