資源簡介
該程序需要在主頁面先下載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
- 上一篇:顯示進度,百分比 和網速 多少kb.rar
- 下一篇:對插補技術的理解.docx
評論
共有 條評論