PD.m

%% 多普勒雷达信号处理过程
%唐黎明 2017-11-27
%%
clc;
close all;
clear;
%% 参数设置

fs = 80e6;                                                                 %载波频率（Hz）
f0 = 100e6;                                                                %中频（Hz）
B = 5e6;                                                                   %脉冲带宽（Hz）
ts =5e-6;                                                                  %脉冲时宽（s）
prt = 100e-6;                                                              %脉冲重复时间间隔（s） prf = 10kHz
num = 32;
N = round（fs*ts）;                                                          %1个脉冲的采样点数
n = -N/2:1:N/2-1;
NN = round（fs*prt）;                                                        %1个脉冲重复间隔内的采样点数
NN_total = NN*num;                                                         %32个脉冲重复间隔内的采样点数
n_total = -0.5*NN_total:1:0.5*NN_total-1;
t_total = ts/N*n_total;
t = ts/N*n;                                                                %1个脉冲对应的时间范围                                                            
fd = 2.5e3;                                                                %动目标多普勒频率（Hz）
k = B/ts;                                                                  %线性调频波斜率
%% 线性调频回波

s_back0 = exp（1i*（2*pi*1/2*k*t.^2））;                                       %基带
s_back2 = [zeros（10.5*（NN-N））s_back0zeros（10.5*（NN-N））];               %补零
s_back3 = real（repmat（s_back21num）.*exp（1i*（2*pi*（fd+f0）*t_total）））;     %32个脉冲
s_back = s_back3;
% s_back = awgn（s_back3-15‘measured‘）;                                     %加入高斯白噪声SNR=-15dB
% s_back = round（s_back/max（abs（s_back））*（2^15-1））;                          %量化    ---?

%量化导出s_back
s_back_16 = round（（s_back/max（abs（s_back）））*（2^15-1））;
for i = 1:256000
    if s_back_16（i） < 0
        s_back_16（i） = s_back_16（i） + 2^16;
    end
end
s_back_eprt = fopen（‘s_back.txt‘‘w‘）;
fprintf（s_back_eprt‘%x\r\n‘s_back_16）;
fclose（s_back_eprt）;

figure（1）;
subplot（211）;
plot（real（s_back3））;
title（‘线性调频波回波时域图像（不含噪声）‘）;
grid on;
subplot（212）;
plot（real（fft（s_back3）））;
title（‘线性调频波回波频域（不含噪声）‘）;
grid on;

figure（2）
subplot（211）;
plot（s_back）;   
title（‘线性调频波回波时域图像（含噪声）‘）;
subplot（212）;
plot（abs（fft（s_back）））;
title（‘线性调频波回波频域图像（含噪声）‘）;
%% 下变频
s_dds = exp（-1i*2*pi*f0*t_total）;