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