maiyadb.m

说明多普勒频率对脉冲压缩输出的影响

% 分析LFM信号在不同多普勒频率时的脉压性能
% 不考虑噪声因素
% 结论：线性调频信号的主旁瓣比会因多普勒频移的增加而降低，但所受到的影响并不明显，
% 即当存在较大的多普勒频率时，同样能获得大的压缩比，也就是说线性调频信号对多普勒频率不敏感。

clear all;
close all;
clc;
%  参数设置 
tao=20e-6;    % 脉宽,s
bw=2e6;       % chirp带宽
fs=8e6;       % 采样率 sample frequency
fd=20e3;      %多普勒频率

% tt=-tao/2+1/fs:1/fs:tao/2;
tt=0:1/fs:tao-1/fs;         % 时间刻度
num_t=length(tt);                       % 一个LFM信号内的采样点数
k=bw/tao;

%产生发射信号
%假设载频为0
st=exp(j*pi*k*tt.^2);
st_fft=fft(st);
figure
subplot(211)
plot(tt,real(st))
xlabel('时间：s')
ylabel('幅度')
title('发射信号')
subplot(212)
plot((0:fs/num_t:fs/2-fs/num_t),abs(st_fft(1:num_t/2)))
xlabel('频率（单位：HZ）')
title('信号频谱')

delay=80e-6;    %假设时延
num_delay=round(delay*fs);
t1=125e-6;
num_t1=round(t1*fs);
tt_all=0:1/fs:t1-1/fs;
sr=exp(j*pi*k*(tt-delay).^2).*exp(j*2*pi*fd*(tt-delay));    %接收信号
% sr=[zeros(1,num_delay),sr,zeros(1,num_t1-num_delay-num_t)];
% figure
% plot(tt_all,real(sr))
figure
plot(tt,real(sr))
xlabel('时间：s')
ylabel('幅度')
title('接收信号')

% match_filter=fliplr(conj(st));
% match_filter_fft=fft(match_filter,1024);
% sr_fft=fft(sr,1024);
% pc_result_fft=match_filter_fft.*sr_fft;
% pc_result=ifft(pc_result_fft,1024);
% figure
% plot(abs(pc_result))

%      pulse compression
match_filter=fliplr(conj(st));   %匹配滤波器的响应
match_filter_fft=fft(match_filter,2*num_t);
sr_fft=fft(sr,2*num_t);
pc_result_fft=match_filter_fft.*sr_fft;     %频域相乘
pc_result=ifft(pc_result_fft,2*num_t);
figure
plot((0:1/fs:2*tao-1/fs),abs(pc_result))
title('回波脉冲压缩处理结果')

pc_result=abs(pc_result);
pc_result=pc_result/max(pc_result);
pc_result=20*log10(pc_result);
figure
% plot((0:1/fs:2*tao-1/fs),pc_result)
% axis([0 2*tao -50 0])
plot(pc_result)
axis([0 320 -50 0])
title('回波脉冲压缩处理结果(dB)')