基于后投影算法穿墙雷达成像

该程序需要在主页面先下载wall_delay程序，该程序在墙后建立了长为2米，宽为0.6米的矩形，利用后投影算法以及快速时延估计进行成像。并且计算最佳带宽，最佳中心频率，阵列天线数量，陈列天线间隔。使得成像质量很好。缺点在于：计算量偏大。

clear all;
close all;
clc;

%%====== 参数设置 =====
%%=== 成像区域 ========
L_width=6;   
W_width=6;
LL=128*2;
WW=128*2;
xx=linspace（-W_width/2W_width/2WW）;%纵轴在-3到3之间划分128*2个点
x=linspace（-W_width/2W_width/2WW*2）;%墙离散成点，用于遍历
yy=linspace（1L_width+1LL）;%在横轴0-6之间划分256个点
%%=== 发射信号参数 ======
c=3.0e8;
B_width=（1）*10^9;%带宽
fc=（0.75）*10^9;  %载频
tRange=10e-8; %脉冲持续时间
Rmax=tRange*c/2; %脉冲最远传播距离
nT=1024*1;%采样点数
fs=nT/tRange;%采样频率
t=linspace（0tRangenT）; %时间划分
%%=== 墙壁设置====
d=0.24;
e_content=8.6;
%%=== 阵列参数 =======单发单收
N_line=31;
nL=linspace（-W_width/2W_width/2N_line）;%天线位置，-1.5-1.5，划分31
% yL=-d;
%%=== 目标位置 ======
x_tag = linspace（-1113）;       %设置目标长度要使得以最小分辨率增长
% y_tag=2*tan（（90-30/2）*pi/180）;
tN=length（x_tag）;   %目标个数


immage=zeros（WWLL）;    %初始化成像矩阵

for i = 0:c/（2*B_width）:0.6      %c/（2*B_width）为间隔最小分辨距离 i为每次y增加长度
    
         y_tag = linspace（3+i3+i13）;   %设置目标的宽度
         
        %%%========== 回波信号 ==========
        pulse=zeros（N_linenT）;  %初始化天线采样数据

    for kk=1:tN  %目标数目
        for k=1:N_line  %天线阵列
        %全遍历
       l11=sqrt（（x_tag（kk）-x）.^2+（y_tag（kk））.^2）;
       l22=sqrt（（nL（k）-x）.^2+d.^2）;
       r_x=sqrt（e_content）*l22+l11;
       r=min（r_x）;

       %t_delay=2*r/c;
       t_delay1= wall_delay（x_tag（kk）y_tag（kk）nL（k）de_content）;
       pulse（k:）=pulse（k:）+exp（j*2*pi*fc*（t-t_delay1））.*（sin（pi*B_width*（t-t_delay1））./（pi*（t-t_delay1）））;
ll=linspace（0RmaxnT）;
figure（1）
stairs（llpulse（1:））; %函数用于绘制阶梯状图
xlabel（‘距离/m‘）
ylabel（‘回波幅度‘）
    end
end

figure（2）
bianll=linspace（0RmaxnT）;
imagesc（ll1:N_lineabs（pulse））;
xlabel（‘距离y/m‘）
ylabel（‘阵列编号‘）
colorbar
%最低点大概为3.8的原因，3为目标点，Y负半轴位墙体0.24米，0.24*sqrt（Er）=0.8;3+0.24*sqrt（Er）=3.8;

%%%========== BP成像 =============
    h = hamming（WW）;%设置汉明窗加权系数
    pxiel1=zeros（WWLL）;%256*256
    pxiel2=zeros（WWLL）;
    pxiel3=zeros（WWLL）;
    for nk=1:length（nL） %天线位置，-1.5-1.5，划分31
      for k=1:WW      %1-256
        for kk=1:LL
                l1=sqrt（（xx（k）-x）.^2+yy（kk）.^2）;%xx（k）在-3到3之间划分128*2个点，yy（kk）在横轴0-6之间划分256个点
                l2=sqrt（（nL（nk）-x）.^2+d.^2）;%nL（nk）天线位置，-1.5-1.5，划分31
                r_pxiel_x=sqrt（e_content）*l2+l1;
            
                r_pxiel=min（r_pxiel_x）;
                %t_delay=2*r_pxiel/c;
                t_delay= wall_delay（xx（k）yy（kk）nL（nk）de_content）;
                n=round（t_delay*nT/（tRange））+1; %nT采样点数1024   %tRange脉冲持续时间
                pxiel1（kkk）=pxiel1（kkk）+pulse（nkn）;  %无加权处理       
                pxiel2（kkk）=pxiel2（kkk）+abs（pulse（nkn））;%采用等权值进行加权处理
                pxiel3（kkk）=pxiel3（kkk）+abs（pulse（nkn））.*h（kk）;%采用汉明加权进行处理
        end
      end
    end
    immage = immage + pxiel3;
end
%成像
figure（3）;
imagesc（yyxxabs（immage））;
ylabel（‘cross range/m‘）;
xlabel（‘range/m‘）;
title（‘汉明加权输出样本成像‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3169  2018-05-10 16:51  bp_wall_imagining.m