%%?多普勒雷達信號處理過程
%唐黎明?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）?????????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）;??????????????????????????????????