一维距离像（HRRP）及FFT成像对比

基于CS的一维距离像（HRRP）及FFT成像对比

%% ISAR回波模拟+基于压缩感知的一位距离像
% Author：潘小义
clear all;
clc;
close all;
format long;
echo off;
% % 目标几何中心是原点，且是刚体运动的转动中心
Omiga = 0.05*2*pi;%刚体旋转速度
% d = [3 3;-3 3;-3 -3;3 -3];%刚体点相对于旋转中心的坐标 此处目标是定义在目标坐标系中
% d = [-10 2;-9 2;-8 2;-7 2;-6 2;-5 2;-4 2;-3 2;-2 2;-2 3;-2 4;-2 5;-2 6;-1 6;0 3;0 4;0 5;0 6;0 2;1 2;2 2;2 3;2 4;2 5;2 6;2 7;2 8;2 9;2 10;3 10;4 3;4 4;4 5;4 6;4 7;4 8;4 9;4 10;4 2;5 2;6 2;7 2;8 2;9 2;10 2;...
%     -9 1;-8 0;-7 -1;10 -1;10 0;10 1;...
%     -6 -2;-5 -2;-4 -2;-3 -2;-2 -2;-1 -2;0 -2;1 -2;2 -2;3 -2;4 -2;5 -2;6 -2;7 -2;8 -2;9 -2;10 -2];
d = [-10 2;-8 2;-6 2;-4 2;-2 2;-2 4;-2 6;-1 6;0 4;0 6;0 2;2 2;2 4;2 6;2 8;2 10;3 10;4 4;4 6;4 8;4 10;4 2;6 2;8 2;10 2;...
    -9 1;-8 0;-7 -1;10 0;...
    -6 -2;-4 -2;-2 -2;0 -2;2 -2;4 -2;6 -2;8 -2;10 -2];
figure
plot（d（:1）d（:2）‘*‘）;xlabel（‘cross-range/m‘‘Fontsize‘20）;ylabel（‘down-range/m‘‘Fontsize‘20）;axis（[-15 15 -15 15 ]）
set（gca‘Fontsize‘20）;
[RowgColg] = size（d）;
theta_o = zeros（1Rowg）;%刚体点和目标几何中心相对初始角度

for jj=1:Rowg
    if d（jj1）>0
        if d（jj2）>0
            theta_o（jj） = atan（d（jj2）/d（jj1））;
        elseif d（jj2）<0
            theta_o（jj） = 2*pi+atan（d（jj2）/d（jj1））;
        elseif d（jj2）==0
            theta_o（jj） = 0;
        end
    elseif d（jj1）<0
        if d（jj2）>0
            theta_o（jj） = pi+atan（d（jj2）/d（jj1））;
        elseif d（jj2）<0
            theta_o（jj） = pi+atan（d（jj2）/d（jj1））;
        elseif d（jj2）==0
            theta_o（jj） = pi;
        end
    elseif d（jj1）==0
        if d（jj2）>0
            theta_o（jj） = pi/2;
        elseif d（jj2）<0
            theta_o（jj） = 3*pi/2;
        elseif d（jj2）==0
            theta_o（jj） = 0;
        end
    end
end
% 以ISAR视线构建坐标系
f0 = 10e9;%载频
c = 3e8;%
lamda = c/f0;
PW = 100e-6;%单个脉冲持续时间
PRF = 1000;
PRI = 1/PRF;%重复周期
M = 128;%方位向采样数
Tp = PRI*M;%总时长
BW = 500e6;%脉内调制带宽
k = BW/PW;%调制斜率
N = 32;%每一个脉冲采样点数（距离单元数）
deltat = PW/（N-1）;%采样间隔
t_all = 0:deltat:Tp;%全时间
t = 0:deltat:PW;%快时间

R0R = 40000;%单站：雷达距离目标几何中心位置
R0J = 40000;% 雷达发射机和接收机的距离
alpha = （0）/180*pi;%双基地角
RRJ = sqrt（R0R^2+R0J^2-2*R0R*R0J*cos（alpha））;%接收站距离雷达发射机位置

det_r = c/2/BW;%距离分辨率
det_a = lamda/2/Omiga/Tp/cos（alpha/2）;%单站ISAR方位向分辨率
for ii=1:M%刚体
     tm = （ii-1）*PRI;%慢时间
    for jj = 1:Rowg
        R（jjii） = R0R+sqrt（d（jj1）^2+d（jj2）^2）*sin（theta_o（jj）+Omiga*tm）+R0J-sqrt（d（jj1）^2+d（jj2）^2）*cos（3*pi/2-theta_o（jj）-Omiga*tm-alpha）;
    end  
end
figure
plot（R‘）;
%% 回波
sr = zeros（MN）;%全回波
A = 1;%刚体散射点反射强度
for jj=1:M
    %% 刚体点回波
    for ii=1:Rowg
        tdelay（jj:） = t-R（iijj）/c;
        sr（jj:） = sr（jj:）+A*exp（sqrt（-1）*2*pi*（f0*tdelay（jj:）+0.5*k.*tdelay（jj:）.*tdelay（jj:）））;
    end    
end
figure;imagesc（abs（sr））
sr = sr/max（max（abs（sr）））;%归一化
%% Dechirping处理获取一维距离像
%参考信号模拟
HRRP = zeros（MN）;
tref = t-（R0R+R0J）/c;%参考距离对应延时
sref  = exp（sqrt（-1）*2*pi*（f0*tref+0.5*k.*tref.*tref））;%参考信号
for ii=1:M
    sif（ii:） = sr（ii:）.*conj（sref）;
end

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

     文件       4783  2013-01-28 17:10  ISARCSSim_dr.m

----------- ---------  ---------- -----  ----

                 4783                    1