%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%????????????小斜視角（3.5°）
%????????????????RDA
%??????????????點目標仿真
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?截止到2014.10.10.?17:19?p.m.

%%
clear;
close?all;
clc;
%?--------------------------------------------------------------------
%?定義參數(shù)
%?--------------------------------------------------------------------
R_nc?=?20e3;????????????????%?景中心斜距
Vr?=?150;???????????????????%?雷達有效速度
Tr?=?2.5e-6;????????????????%?發(fā)射脈沖時寬
Kr?=?20e12;?????????????????%?距離調(diào)頻率
f0?=?5.3e9;?????????????????%?雷達工作頻率
BW_dop?=?80;????????????????%?多普勒帶寬
Fr?=?60e6;??????????????????%?距離采樣率
Fa?=?200;???????????????????%?方位采樣率
Naz?=?1024;?????????????????%?距離線數(shù)（即數(shù)據(jù)矩陣，行數(shù)）——這里修改為1024。
Nrg?=?320;??????????????????%?距離線采樣點數(shù)（即數(shù)據(jù)矩陣，列數(shù)）
sita_r_c?=?（3.5*pi）/180;	%?波束斜視角，3.5?度，這里轉(zhuǎn)換為弧度
c?=?3e8;????????????????????%?光速

R0?=?R_nc*cos（sita_r_c）;	%?與R_nc相對應的最近斜距，記為R0
Nr?=?Tr*Fr;?????????????????%?線性調(diào)頻信號采樣點數(shù)
BW_range?=?Kr*Tr;???????????%?距離向帶寬
lamda?=?c/f0;???????????????%?波長
fnc?=?2*Vr*sin（sita_r_c）/lamda;?????%?多普勒中心頻率，根據(jù)公式（4.33）計算。
La_real?=?0.886*2*Vr*cos（sita_r_c）/BW_dop;?%?方位向天線長度，根據(jù)公式（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的方位向距離
%?定義以下數(shù)組，便于處理
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];???????%?定義該數(shù)組，便于處理。

%%?
%?--------------------------------------------------------------------
%?距離（方位）向時間，頻率相關定義
%?--------------------------------------------------------------------
%?距離
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

%%?
%?-------------------------------------------------------------------