SAR点目标RDA实现

RD算法实现，包括低斜视角和大斜视角（SRC方式2、SRC方式3）两种情况，系统固有参数见《Digital Processing of SAR Data》一书。

%%
close all;
clear all;
clc;

Num = 1;                                                                   %目标点数
%% 参数设定
%——————————————机载参数设定——————————————%
j = sqrt（-1）;                     
Theta_rc = 3.5*pi/180;                                                     %波束斜视角（低斜视角）
Vr = 150;                                                                  %雷达有效速度
Rc = 20e3;                                                                 %景中心斜距
R0 = Rc*cos（Theta_rc）;                                                     %最近斜距
c = 3e8;                                                                   %光速
f0 = 5.3e9;                                                                %载频
lambda = c/f0;                                                             %载波波长
Fr = 60e6;                                                                 %距离向采样率
Tr = 2.5e-6;                                                               %发射脉冲时宽
Kr = 20e12;                                                                %距离向调频率
PRF = 100;                                                                 %脉冲重复频率
Fa = PRF;                                                                  %方位向采样率
Nr = 320;                                                                  %距离向采样点数
Na = 256;                                                                  %方位向采样点数
Delta_Fdop = 80;                                                           %多普勒带宽
f_eta_c = 320;                                                             %多普勒中心频率
eta_c = -8.1;                                                              %波束中心偏移时间
La = 0.886*2*Vr*cos（Theta_rc）/Delta_Fdop;                                  %天线孔径
Ka = 2*Vr.^2*cos（Theta_rc）.^2/（lambda*Rc）;                                 %方位向调频率
Ta = 0.886*lambda*Rc/（La*Vr*cos（Theta_rc））;                                %目标照射时间
P_Target = [Vr*eta_c Rc];                                                  %目标位置（斜距平面）
%% 采样时间
%——————————————采样时间———————————————%
tr = 2*Rc/c+（-Nr/2:1:Nr/2-1）/Fr;                                           %距离向采样时间
ta = eta_c+（-Na/2:1:Na/2-1）/PRF;                                           %方位向采样时间
%% 瞬时斜距
R = sqrt（R0.^2+（Vr*ta）.^2）;                                                %瞬时斜距
R = R.‘*ones（1Nr）;                                                        %R是一个Na*Nr（256*320）的矩阵
                                                                           %注意！：“.‘”是非共轭转置，“‘”是共轭转置                                                                                                                                      
%% sinc窗（方位向包络加权）
%——————————sinc窗（方位向包络加权）———————————%
Wa = sinc（La*atan（Vr*（ta-eta_c）/R0）/lambda）.^2;                            %92
Wa = Wa.‘*ones（1Nr）;
%% 矩形窗（距离向包络加权）
%——————————矩形窗（距离向包络加权）———————————%
Wr = zeros（Na Nr）;                                                        %86
for i = 1:Na
    for k = 1:Nr
        if abs（tr（k）-2*R（i）/c）<=（1/2*Tr）                             

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1229509  2015-09-17 22:26  RDA_point\Fig\FIG\RD_point\AC.fig

     文件    1256802  2015-09-17 22:25  RDA_point\Fig\FIG\RD_point\FA.fig

     文件      21052  2015-09-17 22:27  RDA_point\Fig\FIG\RD_point\Map.fig

     文件     519395  2015-09-17 22:23  RDA_point\Fig\FIG\RD_point\Raw.fig

     文件    1254744  2015-09-17 22:25  RDA_point\Fig\FIG\RD_point\RC.fig

     文件    1218865  2015-09-17 22:26  RDA_point\Fig\FIG\RD_point\RCMC.fig

     文件    1256838  2015-09-17 22:26  RDA_point\Fig\FIG\RD_point\Unwrap.fig

     文件    1451834  2015-09-17 22:27  RDA_point\Fig\FIG\RD_point\Upx16.fig

     文件    1076109  2015-09-17 22:38  RDA_point\Fig\FIG\RD_SRC2_point\AC.fig

     文件    1258528  2015-09-17 22:37  RDA_point\Fig\FIG\RD_SRC2_point\FA+SRC2.fig

     文件    1257982  2015-09-17 22:31  RDA_point\Fig\FIG\RD_SRC2_point\FA.fig

     文件      21011  2015-09-17 22:38  RDA_point\Fig\FIG\RD_SRC2_point\Map.fig

     文件     534585  2015-09-17 22:31  RDA_point\Fig\FIG\RD_SRC2_point\Raw.fig

     文件    1255497  2015-09-17 22:31  RDA_point\Fig\FIG\RD_SRC2_point\RC.fig

     文件     984802  2015-09-17 22:37  RDA_point\Fig\FIG\RD_SRC2_point\RCMC.fig

     文件    1258562  2015-09-17 22:37  RDA_point\Fig\FIG\RD_SRC2_point\Unwrap.fig

     文件    1916806  2015-09-17 22:38  RDA_point\Fig\FIG\RD_SRC2_point\Upx16.fig

     文件    1076193  2015-09-17 22:47  RDA_point\Fig\FIG\RD_SRC3_point\AC.fig

     文件    1258762  2015-09-17 22:46  RDA_point\Fig\FIG\RD_SRC3_point\FA.fig

     文件      21147  2015-09-17 22:48  RDA_point\Fig\FIG\RD_SRC3_point\Map.fig

     文件     534839  2015-09-17 22:45  RDA_point\Fig\FIG\RD_SRC3_point\Raw.fig

     文件    1255955  2015-09-17 22:45  RDA_point\Fig\FIG\RD_SRC3_point\RC+SRC3.fig

     文件     984979  2015-09-17 22:47  RDA_point\Fig\FIG\RD_SRC3_point\RCMC.fig

     文件    1258758  2015-09-17 22:46  RDA_point\Fig\FIG\RD_SRC3_point\Unwrap.fig

     文件    1885721  2015-09-17 22:47  RDA_point\Fig\FIG\RD_SRC3_point\Upx16.fig

     文件      22814  2015-09-17 16:37  RDA_point\Fig\JPEG\RD_point\AC.jpg

     文件      24941  2015-09-17 16:34  RDA_point\Fig\JPEG\RD_point\FA.jpg

     文件      63227  2015-09-17 17:46  RDA_point\Fig\JPEG\RD_point\Map.jpg

     文件      43769  2015-09-17 16:32  RDA_point\Fig\JPEG\RD_point\Raw.jpg

     文件      24119  2015-09-17 16:33  RDA_point\Fig\JPEG\RD_point\RC.jpg

............此处省略41个文件信息