穿墙成像的BP算法

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（0L_widthLL）;%在横轴0-6之间划分256个点
%%=== 发射信号参数 ======
c=3.0e8;
B_width=（1）*10^9;%带宽
fc=（1.5）*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/4W_width/4N_line）;%天线位置，-1.5-1.5，划分31
% yL=-d;
%%=== 目标位置 ======
x_tag = -1;
% y_tag=2*tan（（90-30/2）*pi/180）;
y_tag = 3;
tN=length（x_tag）;   %目标个数

%%%========== 回波信号 ==========
pulse=zeros（N_linenT）;
for kk=1:tN  %目标数目
    for k=1:N_line  %天线阵列
        %全遍历
       l1=sqrt（（x_tag（kk）-x）.^2+（y_tag（kk））.^2）;
       l2=sqrt（（nL（k）-x）.^2+d.^2）;
       r_x=sqrt（e_content）*l2+l1;
       r=min（r_x）;

       t_delay=2*r/c;
       pulse（k:）=pulse（k:）+exp（j*2*pi*fc*（t-t_delay））.*（sin（pi*B_width*（t-t_delay））./（pi*（t-t_delay）））;
ll=lins