SAR成像bp算法代码

%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）;