合成孔径BP成像.m

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