雷达信号处理流程

%%雷达信号基本处理流程 脉压
%%平台：R2016a

clear all;close all;clc;

%% 基本参数设置
f0 = 10e9;       %载频
Tp = 10e-6;      %脉冲宽度
B  = 10e6;       %信号带宽
fs = 100e6;      %采样率
R0 = 3000;       %目标初始距离
c  = 3e8;        %光速
tr = 2*R0/c;     %目标点延时
k  = B/Tp;       %调频斜率
N  = 4096;       %采样点数
t  = （0:N-1）/fs; %采样时间步进

%% 发射信号
S0 = rectpuls（t-Tp/2Tp） .* exp（1i*pi*k*（t-Tp/2）.^2）;
figure（1）;
subplot（211）;plot（t*c/2real（S0））;title（‘发射信号实部‘）;xlabel（‘距离/m‘）;
subplot（212）;plot（t*c/2imag（S0））;title（‘发射信号虚部‘）;xlabel（‘距离/m‘）;
%subplot（313）;plot（abs（fft（S0）））;

%% 发射信号频谱
S0_fft = fft（S0N）;

%% 回波构造
S1 = rectpuls（t-tr-Tp/2Tp） .* exp（1i*pi*k*（t-tr-Tp/2）.^2） .* exp（-1i*2*pi*f0*tr）;
figure（2）;
subplot（211）;plot（t*c/2real（S1））;title（‘回波函数实部‘）;xlabel（‘距离/m‘）;
subplot（212）;plot（t*c/2imag（S1））;title（‘回波函数虚部‘）;xlabel（‘距离/m‘）;
%subplot（313）;plot（abs（fft（S1）））;

%% 回波频谱
f = fs/N * （-N/2:N/2 - 1）;
S1_fft = fft（S1N）;
figure（3）;
subplot（211）;plot（abs（S1_fft））;title（‘回波频谱‘）;xlabel（‘频率‘）;
subplot（212）;plot（fftshift（abs（S1_fft）））;title（‘回波频谱‘）;xlabel（‘频率‘）;

%% 脉冲压缩
S_F = ifft（ S1_fft .* conj（S0_fft））;
figure（4）;
subplot（211）;plot（t*c/2abs（S_F））;title（‘脉压‘）;xlabel（‘距离/m‘）;
subplot（212）;plot（t*c/2db（abs（S_F）/max（S_F）））;title（‘归一化‘）;xlabel（‘距离/m‘）;

%% 频域加窗
win = hamming（410）;
window = [zeros（11843）