clear;close?all;clc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%?信號基本參數
C?=?3e8;?????????????????????????????????????????????????%?光速
Fc?=?35.75e9;??????????????????????????????????????????%?載頻
lambda?=?C/Fc;??????????????????????????????????????%?波長
H?=?1000;???????????????????????????????????????????????%?高度
angleD?=?72;?????????????????????????????????????????????%?下視角
Rs?=?H/cosd（angleD）;?????????????????????????????????%?參考距離
Nr?=?1024;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%?方位維度
Fa?=?2000;??????????????????????????????????????????????%?方位頻率
PRT?=?1/Fa?;???????????????????????????????????????????%?方位PRT
angleA?=?0.95;????????????????????????????????????????%?方位向波束寬度
Lsar?=（angleA）*pi/180*Rs;??????????????????????%?合成孔徑長度
V?=?46;????????????????????????????????????????????????%?平臺運動速度
fdc?=?0;????????????????????????????????????????????????%（正側視方位中心=0）?
Ka?=?-2*V.^2/lambda/Rs;?????????????????????%?方位調頻斜率?
Tsar?=?Lsar/V;
Ba?=?abs（Ka*Tsar）;
fdM?=?2*V/lambda;%*sind（angleA）;
da?=?V*PRT?%?方位間隔
Xmin?=?-20;
Xmax?=?20;
%?方位全孔徑
Na?=?ceil（（Xmax?-?Xmin?+?Lsar）/da）;
Na?=?2^nextpow2（Na）;
PRT?=?（Xmax?-?Xmin?+?Lsar）/V/Na;
Fa?=?1/PRT;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%?距離維度
Fr?=?900e6;?%?距離采樣頻率??
Br?=?700e6;?%?距離帶寬
Tr?=?10e-6;??%?距離時寬?
Kr?=?Br/Tr;?%?調頻斜率
angleR?=?9.5;?%?距離向波束寬度?
W?=?（Nr-1）*C/Fr/2;?%?測繪帶寬度?341.3333
Rnear?=?Rs?-?W/2;%（Nr/2）*C/Fr/2/2;?%?近點距離?3150.7
Rfar?=?Rs?+?W/2;%（Nr/2-1）*C/Fr/2/2;?%?遠點距離?3321.4
Rmin?=??Rnear?+?Tr/2*C/2;
Rmax?=??Rfar?-?Tr/2*C/2;
dr?=?C/2/Fr?%?距離間隔
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?測繪帶信息
Ymin?=?sqrt（Rmin^2?-?H^2）;?%?2384.7?%測繪帶寬度受到Nr控制，不是原本的視角范圍
Ymax?=?sqrt（Rmax^2?-?H^2）;?%?4246.8
Yc?=?（Ymin?+?Ymax）/2;?????????????????????????%?3077.7
Xc?=?（Xmin?+?Xmax）/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?時頻域ff?tt
tt?=?linspace（?2*Rnear/C?2*Rfar/C??Nr）;
td?=?linspace（（Xmin?-?Lsar/2）/V（Xmax??+?Lsar/2）/VNa）;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?參考坐標
positionReal?=[zeros（1Na）;td*V;H*ones（1Na）];%參考值?（全孔徑）%?Y?X?Z
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?目標
%?Target?=?[Xc?Yc?0;];%單點目標
Target?=?[Xc?Yc?-?50*dr?0;
??????????Xc?Yc?+?200*dr?0];
??????????????????????
[Ntarget~]=size（Target）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?回波信號
????Srnm?=?0;
????R?=?zeros（NtargetNa）;
for?i?=?1:Ntarget
????positionTarget?=?Target（i:）;?%?當前目標位置信息
????R（i:）?=?sqrt（（positionReal（2:）-positionTarget（1））.^2?+?（positionReal（1:）-?positionTarget（2））.^2?+?（positionReal（3:）-positionTarget（3））.^2）;
????tau=2*R（i:）/C;?%?時延
????Tfast?=?tt‘?-?tau;????
????iphase1?=?pi*Kr*Tfast.^2;??%?距離維chirp
????iphase2?=-（4*pi/lambda）*（ones（Nr1）*R（i:））;?%?距離調制相位
????iphase?=?iphase1?+?iphase2;???
????Srnm=Srnm?+?exp（1i*iphase）.*（-Lsar/2<（td.*V-?positionTarget（1））&（td.*V?-?positionTarget（1））end
???