%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%????????????小斜視角（3.5°）
%????????????????RDA
%??????????????點(diǎn)目標(biāo)仿真
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?截止到2014.10.10.?17:19?p.m.

%%
clear;
close?all;
clc;
%?--------------------------------------------------------------------
%?定義參數(shù)
%?--------------------------------------------------------------------
R_nc?=?20e3;????????????????%?景中心斜距
Vr?=?150;???????????????????%?雷達(dá)有效速度
Tr?=?2.5e-6;????????????????%?發(fā)射脈沖時(shí)寬
Kr?=?20e12;?????????????????%?距離調(diào)頻率
f0?=?5.3e9;?????????????????%?雷達(dá)工作頻率
BW_dop?=?80;????????????????%?多普勒帶寬
Fr?=?60e6;??????????????????%?距離采樣率
Fa?=?200;???????????????????%?方位采樣率
Naz?=?1024;?????????????????%?距離線數(shù)（即數(shù)據(jù)矩陣，行數(shù)）——這里修改為1024。
Nrg?=?320;??????????????????%?距離線采樣點(diǎn)數(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相對(duì)應(yīng)的最近斜距，記為R0
Nr?=?Tr*Fr;?????????????????%?線性調(diào)頻信號(hào)采樣點(diǎn)數(shù)
BW_range?=?Kr*Tr;???????????%?距離向帶寬
lamda?=?c/f0;???????????????%?波長
fnc?=?2*Vr*sin（sita_r_c）/lamda;?????%?多普勒中心頻率，根據(jù)公式（4.33）計(jì)算。
La_real?=?0.886*2*Vr*cos（sita_r_c）/BW_dop;?%?方位向天線長度，根據(jù)公式（4.36）
beta_bw?=?0.886*lamda/La_real;??????????????%?雷達(dá)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長度

%?--------------------------------------------------------------------
%?設(shè)定仿真點(diǎn)目標(biāo)的位置
%?以距離向作為x軸正方向
%?以方位向作為y軸正方向
%?--------------------------------------------------------------------?
delta_R0?=?0;???????%?將目標(biāo)1的波束中心穿越時(shí)刻，定義為方位向時(shí)間零點(diǎn)。
delta_R1?=?120;?	%?目標(biāo)1和目標(biāo)2的方位向距離差，120m
delta_R2?=?50;??????%?目標(biāo)2和目標(biāo)3的距離向距離差，50m

%?目標(biāo)1
x1?=?R0;????????????%?目標(biāo)1的距離向距離
y1?=?delta_R0?+?x1*tan（sita_r_c）;	%?目標(biāo)1的方位向距離

%?目標(biāo)2
x2?=?x1;????????????%?目標(biāo)2和目標(biāo)1的距離向距離相同
y2?=?y1?+?delta_R1;?%?目標(biāo)2的方位向距離
%?目標(biāo)3
x3?=?x2?+?delta_R2;?????????????????%?目標(biāo)3和目標(biāo)2有距離向的距離差，為delta_R2
y3?=?y2?+?delta_R2*tan（sita_r_c）;??	%?目標(biāo)3的方位向距離
%?定義以下數(shù)組，便于處理
x_range?=?[x1x2x3];
y_azimuth?=?[y1y2y3];

%?計(jì)算三個(gè)目標(biāo)各自的波束中心穿越時(shí)刻
nc_1?=?（y1-x1*tan（sita_r_c））/Vr;????%?目標(biāo)1的波束中心穿越時(shí)刻。
nc_2?=?（y2-x2*tan（sita_r_c））/Vr;????%?目標(biāo)2的波束中心穿越時(shí)刻。
nc_3?=?（y3-x3*tan（sita_r_c））/Vr;????%?目標(biāo)3的波束中心穿越時(shí)刻。
nc_target?=?[nc_1nc_2nc_3];???????%?定義該數(shù)組，便于處理。

%%?
%?--------------------------------------------------------------------
%?距離（方位）向時(shí)間，頻率相關(guān)定義
%?--------------------------------------------------------------------
%?距離
tr?=?2*x1/c?+?（?-Nrg/2?:?（Nrg/2-1）?）/Fr;????????????????%?距離時(shí)間軸
fr?=?（?-NFFT_r/2?:?NFFT_r/2-1?）*（?Fr/NFFT_r?）;??????????%?距離頻率軸
%?方位
ta?=?（?-Naz/2:?Naz/2-1?）/Fa;????????????????????????????%?方位時(shí)間軸
fa?=?fnc?+?（?-NFFT_a/2?:?NFFT_a/2-1?）*（?Fa/NFFT_a?）;	%?方位頻率軸

%?生成距離（方位）時(shí)間（頻率）矩陣
tr_mtx?=?ones（Naz1）*tr;????%?距離時(shí)間軸矩陣，大小：Naz*Nrg
ta_mtx?=?ta.‘*ones（1Nrg）;??%?方位時(shí)間軸矩陣，大小：Naz*Nrg

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