%BPA?多點目標?單站SAR
%參數來自表6.1
%date:2010/10/14
clc;
clear?all;
close?all;

%（1）參數定義
%================================================
c=3e8;
j=sqrt（-1）;
pi=3.1416;
fc=5.3e9;
lamda=c/fc;
D=4;??%方位向天線尺寸
Va=150;
Kr=20e12;
Tr=2.5e-6;
sq_ang=3.5/180*pi;

Br=Kr*Tr;
Frfactor=1.2;
Fr=Br*Frfactor;
Ba=0.886*2*Va*cos（sq_ang）/D;
Fafactor=1.2;
Fa=Ba*Fafactor;

R_near=2e4;???%近距離
R_far=R_near+1000;???%遠距離
%R=（R_near+R_far）/2;%assumed?a?fixed?R?for?implement
%Y_min=-50;%both?R_near?and?R_far?are?closest?approach?range
%Y_max=Y_min+100;

La_near=0.886*R_near*lamda/cos（sq_ang）^2/D;%近程點的合成孔徑長度
La_far=0.886*R_far*lamda/cos（sq_ang）^2/D;%遠距離點的合成孔徑長度
Tc_near=-R_near*tan（sq_ang）/Va;??%近距離點的波束中心穿越時刻
Tc_far=-R_far*tan（sq_ang）/Va;??%遠距離點的波束中心穿越時刻
fdc=2*Va*sin（sq_ang）/lamda;??%多普勒中心
Y_min=Va*Tc_far;
Y_max=Y_min+100;


Rmin=sqrt（R_near^2+（Tc_near*Va+La_near/2）^2）;%雷達與目標區域之間的最短距離
Rmax=sqrt（R_far^2+（Tc_far*Va-La_far/2）^2）;%雷達與目標區域之間的最長距離

disp（‘parameters:‘）;
disp（‘minimalslantrange:‘）;disp（Rmin）;
disp（‘maximalslantrange:‘）;disp（Rmax）;
disp（‘rangeresolution:‘）;disp（0.886*（c/2/Br））;
disp（‘azimuthresolution:‘）;disp（0.886*Va/Ba）;
disp（‘dopplercentroidfrequency:‘）;disp（fdc）;
%================================================

%（2）回波模型
%================================================
Nr=（2*Rmax/c-2*Rmin/c+Tr）*Fr;
Nr=2^nextpow2（Nr）;
tr=linspace（-Tr/2+2*Rmin/cTr/2+2*Rmax/cNr）;
Fr=（Nr-1）/（Tr/2+2*Rmax/c-（-Tr/2+2*Rmin/c））;
Na=（（Tc_near+La_near/2/Va）-（Tc_far-La_far/2/Va））*Fa;
Na=2^nextpow2（Na）;
ta=linspace（Tc_far-La_far/2/VaTc_near+La_near/2/VaNa）;
Fa=（Na-1）/（Tc_near+La_near/2/Va-（Tc_far-La_far/2/Va））;

Rpt=[R_near?R_near+500?R_near+1000];%點目標位置
Ypt=[0?0?0];
La=0.886*Rpt*lamda/（cos（sq_ang）^2）/D;%每個目標的合成孔徑長度
Tc=-Rpt*tan（sq_ang）/Va;%每個目標的波束中心穿越時刻
Npt=length（Rpt）;

Y_high=max（Ypt）+50;%%確定成像網格的范圍，以使其包括目標點;
Y_low=min（Ypt）-50;%%成像區域是（R_right-R_left）*（Y_high-Y_low）這一片區域
R_left=R_near-50;
R_right=R_far+50;

disp（‘numberofpointtargets:‘）;disp（Npt）;
disp（‘rangesamplenumber:‘）;disp（Nr）;
disp（‘azimuthsamplenumber:‘）;disp（Na）;
disp（‘rangesamplerate:‘）;disp（Fr）;