均勻間距線列陣波束形成器[運用分享]

1、實驗3 均勻間距線列陣波束形成器姓 名：逯仁杰班 級：20120001（12級陳賡1班）學(xué) 號：20120111121.實驗?zāi)康耐ㄟ^本實驗的學(xué)習(xí)，加深對聲納技術(shù)中波束形成和方位估計的概念理解，理解聲納信號處理的基本過程，為今后聲納信號處理的工作和學(xué)習(xí)建立概念、奠定堅實的基礎(chǔ)。2.實驗原理 波束形成器的本質(zhì)是一個空間濾波器。當(dāng)對基陣各基元接收信號作補(bǔ)償處理，使得各基元對某個特定方向上的信號能夠同相相加，獲得一個最大的響應(yīng)輸出（幅度相加）；相應(yīng)的各基元對其它方向的信號非同相相加，產(chǎn)生一定的相消效果的響應(yīng)輸出（對于各基元噪聲相互獨立的情況時功率相加）。這就是波束形成的工作原理。 常用的波束形成方法主

2、要有時延波束形成法和頻域波束形成法。在此基礎(chǔ)上針對不同的陣形、設(shè)計要求以及背景噪聲特性下還發(fā)展了許多波束形成算法。針對不同的陣形時的波束形成方法是指依賴于陣形的特殊性（如直線陣、圓陣、體積陣等）而得到的波束形成算法：如直線陣波束形成法、圓陣波束形成法，體積陣波束形成法等。 針對不同的設(shè)計要求也衍生出多種新型的波束形成算法。當(dāng)對不同的頻率響應(yīng)要求相同的波束寬度時有恒定束寬波束形成法，當(dāng)對波束的旁瓣級有要求時可采用切比雪夫加權(quán)波束形成法。當(dāng)要求對陣列誤差具有寬容性響應(yīng)時失配條件下的波束形成器6,362-382。如果利用噪聲干擾的統(tǒng)計特性有高分辨最小方差無畸變響應(yīng)（MVDR）波束形成法，線性約束最小

3、方差（LCMV）波束形成法，線性約束最小功率（LCMP）波束形成法，自適應(yīng)波束形成法等。 但不管是何種波束形成方法，其目的均是在干擾背景下獲取某個方向的信號或估計信號的方位。下面僅給出時延波束形成和相移波束形成的基本原理。時延波束形成法（時域）相移波束形成法（頻域）3.實驗內(nèi)容（1）仿真等間距直線陣基元接收信號，對所接收信號進(jìn)行延時波束形成，估計目標(biāo)方位；分析波束形成性能。 參數(shù)：陣元數(shù) 16，中心頻率1500Hz，帶寬500Hz，信號脈寬20ms，信噪比20dB。（2）采用頻域波束形成方法對所接收信號進(jìn)行波束形成，估計目標(biāo)方位，分析波束形成性能。參數(shù)不變。4. 實驗結(jié)果及數(shù)據(jù)分析clc;cl

4、ear all;close all; N=input(請輸入陣元數(shù):);f0=1500;B=500;c=1500;T=0.02;d=0.5*c/f0;fs=100*f0;A=0.1;L=3*(N-1)*d/c;t=0:1/fs:L; theta0=input(請輸入目標(biāo)方位角（角度）:);theta0=theta0/180*pi;x=zeros(N,length(t);for k=0:N-1 if(theta00) for k=0:N-1 tao=fix(k*d*sin(theta/180*pi)/c*fs)+1; y(k+1,1:length(t)+1-tao)=x(k+1,tao:end)

5、; end else for k=0:N-1 tao=fix(N-k-1)*d*sin(theta/180*pi)/c*fs)-1; y(k+1,1:length(t)+1+tao)=x(k+1,-tao:end); end end sumy=sum(y); ps(theta+91,:)=sum(sumy.*sumy)/length(sumy);endfigureplot(-90:90,20*log10(ps/max(ps);set(gca,FontSize,20);title(時域波束形成（無噪聲）);xlabel(theta/度);ylabel(輸出/dB);grid on取基元為16，入

6、射角為30度，也的確在30度輸出最大。n=normrnd(0,1,1,length(t);w=2*f0-B/2,f0+B/2/fs;b=fir1(128,w,bandpass);np=filter(b,1,n);pn=std(np)2;ps=A2/2;ks=sqrt(ps/(100*pn);x=zeros(N,length(t);xx=A*cos(2*pi*f0*t)+ks*np;for k=0:N-1 if(theta00) for k=0:N-1 tao=fix(k*d*sin(theta/180*pi)/c*fs )+1; y(k+1,1:length(t)+1-tao)=x(k+1,t

7、ao:end); end else for k=0:N-1 tao=fix(N-k-1)*d*sin(theta/180*pi)/c*fs)-1; y(k+1,1:length(t)+1+tao)=x(k+1,-tao:end); end end sumy=sum(y); ps(theta+91,:)=sum(sumy.*sumy)/length(sumy);endfigure plot(-90:90,20*log10(ps);set(gca,FontSize,20);title(theta方向掃描結(jié)果（有噪聲）);xlabel(theta/度);ylabel(輸出/dB);grid on f

8、ftx=fft(x,fs);pa=fftx(1500,:);phas=atan(real(pa)./imag(pa);for n=1:N; phas(n)=phas(n)+pi*(n-1);endfigurestem(phas);set(gca,FontSize,20);title(1500Hz處各基元接收信號相位關(guān)系);ylabel(phi);px=(0:N-1).*d;alpha=(-90:90)*pi./180;s=exp(-j*2*pi*f0/c*px*sin(alpha);for k=0:N-1 if(theta00) tao=fix(N-k-1)*d*sin(theta0)/c*f

9、s)-1; x(k+1,-tao:end)=xx(1:length(t)+1+tao); Sf=fft(x(k+1,:),fs); R1(k+1)=Sf(f0+1); else tao=fix(k*d*sin(theta0)/c*fs)+1; x(k+1,tao:end)=xx(1:length(t)+1-tao); Sf=fft(x(k+1,:),fs); R1(k+1)=Sf(f0+1); end endR0=R1*s;yy=abs(real(R0);yy=yy./max(yy);figureplot(-alpha.*180./pi,20.*log10(yy);set(gca,FontSize,20);xlabel(theta/度);ylabel(輸出/dB);title(頻域波束形成)grid on可見各基元間相位基本是線性，頻域波束圖同樣能確定入射波的方位角在30度處。因為頻域的相位補(bǔ)償與時域的時延補(bǔ)償作用類似，都使各基元接收實現(xiàn)了同相相加，獲得最大的幅值?，F(xiàn)將基元數(shù)變更為40，其他參數(shù)不變。通過比較觀察，可以發(fā)現(xiàn)，時域波束形成圖中的尖峰數(shù)即為陣元個數(shù)減一，即N-1個，頻域圖中各個尖峰之間會出現(xiàn)小的旁瓣，可以起到增加通帶增益抑制阻帶，也對確定目標(biāo)方位更有益處。改變N，可以明顯看出主波束寬度有所變化，N越大，主波束寬度越小，指向性越好。6.總結(jié)與心