RDA_imaging.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%            小斜视角（3.5°）
%                RDA
%              点目标仿真
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 截止到2014.10.10. 17:19 p.m.

%%
clear;
close all;
clc;
% --------------------------------------------------------------------
% 定义参数
% --------------------------------------------------------------------
R_nc = 20e3;                % 景中心斜距
Vr = 150;                   % 雷达有效速度
Tr = 2.5e-6;                % 发射脉冲时宽
Kr = 20e12;                 % 距离调频率
f0 = 5.3e9;                 % 雷达工作频率
BW_dop = 80;                % 多普勒带宽
Fr = 60e6;                  % 距离采样率
Fa = 200;                   % 方位采样率
Naz = 1024;                 % 距离线数（即数据矩阵，行数）——这里修改为1024。
Nrg = 320;                  % 距离线采样点数（即数据矩阵，列数）
sita_r_c = （3.5*pi）/180;	% 波束斜视角，3.5 度，这里转换为弧度
c = 3e8;                    % 光速

R0 = R_nc*cos（sita_r_c）;	% 与R_nc相对应的最近斜距，记为R0
Nr = Tr*Fr;                 % 线性调频信号采样点数
BW_range = Kr*Tr;           % 距离向带宽
lamda = c/f0;               % 波长
fnc = 2*Vr*sin（sita_r_c）/lamda;     % 多普勒中心频率，根据公式（4.33）计算。
La_real = 0.886*2*Vr*cos（sita_r_c）/BW_dop; % 方位向天线长度，根据公式（4.36）
beta_bw = 0.886*lamda/La_real;              % 雷达3dB波束     
La = 0.886*R_nc*lamda/La_real;              % 合成孔径长度
a_sr = Fr / BW_range;       % 距离向过采样因子
a_sa = Fa / BW_dop;         % 方位向过采样因子

Mamb = round（fnc/Fa）;       % 多普勒模糊

NFFT_r = Nrg;               % 距离向FFT长度
NFFT_a = Naz;               % 方位向FFT长度

% --------------------------------------------------------------------
% 设定仿真点目标的位置
% 以距离向作为x轴正方向
% 以方位向作为y轴正方向
% -------------------------------------------------------------------- 
delta_R0 = 0;       % 将目标1的波束中心穿越时刻，定义为方位向时间零点。
delta_R1 = 120; 	% 目标1和目标2的方位向距离差，120m
delta_R2 = 50;      % 目标2和目标3的距离向距离差，50m

% 目标1
x1 = R0;            % 目标1的距离向距离
y1 = delta_R0 + x1*tan（sita_r_c）;	% 目标1的方位向距离

% 目标2
x2 = x1;            % 目标2和目标1的距离向距离相同
y2 = y1 + delta_R1; % 目标2的方位向距离
% 目标3
x3 = x2 + delta_R2;                 % 目标3和目标2有距离向的距离差，为delta_R2
y3 = y2 + delta_R2*tan（sita_r_c）;  	% 目标3的方位向距离
% 定义以下数组，便于处理
x_range = [x1x2x3];
y_azimuth = [y1y2y3];

% 计算三个目标各自的波束中心穿越时刻
nc_1 = （y1-x1*tan（sita_r_c））/Vr;    % 目标1的波束中心穿越时刻。
nc_2 = （y2-x2*tan（sita_r_c））/Vr;    % 目标2的波束中心穿越时刻。
nc_3 = （y3-x3*tan（sita_r_c））/Vr;    % 目标3的波束中心穿越时刻。
nc_target = [nc_1nc_2nc_3];       % 定义该数组，便于处理。

%% 
% --------------------------------------------------------------------
% 距离（方位）向时间，频率相关定义
% --------------------------------------------------------------------
% 距离
tr = 2*x1/c + （ -Nrg/2 : （Nrg/2-1） ）/Fr;                % 距离时间轴
fr = （ -NFFT_r/2 : NFFT_r/2-1 ）*（ Fr/NFFT_r ）;          % 距离频率轴
% 方位
ta = （ -Naz/2: Naz/2-1 ）/Fa;                            % 方位时间轴
fa = fnc + （ -NFFT_a/2 : NFFT_a/2-1 ）*（ Fa/NFFT_a ）;	% 方位频率轴

% 生成距离（方位）时间（频率）矩阵
tr_mtx = ones（Naz1）*tr;    % 距离时间轴矩阵，大小：Naz*Nrg
ta_mtx = ta.‘*ones（1Nrg）;  % 方位时间轴矩阵，大小：Naz*Nrg

%% 
% -------------------------------------------------------------------