均勻間距線列陣波束形成器[運(yùn)用分享]

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

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

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

4、ear all;close all; N=input(請(qǐng)輸入陣元數(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(請(qǐng)輸入目標(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(時(shí)域波束形成（無噪聲）);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處各基元接收信號(hào)相位關(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度處。因?yàn)轭l域的相位補(bǔ)償與時(shí)域的時(shí)延補(bǔ)償作用類似，都使各基元接收實(shí)現(xiàn)了同相相加，獲得最大的幅值?，F(xiàn)將基元數(shù)變更為40，其他參數(shù)不變。通過比較觀察，可以發(fā)現(xiàn)，時(shí)域波束形成圖中的尖峰數(shù)即為陣元個(gè)數(shù)減一，即N-1個(gè)，頻域圖中各個(gè)尖峰之間會(huì)出現(xiàn)小的旁瓣，可以起到增加通帶增益抑制阻帶，也對(duì)確定目標(biāo)方位更有益處。改變N，可以明顯看出主波束寬度有所變化，N越大，主波束寬度越小，指向性越好。6.總結(jié)與心