声纳技术实验报告三汇总

1、1. 实验目的(1) 、熟悉波束形成的物理机理和基本原理，熟悉时延波束形成基本方法(2) 、练习MATLAB用于声纳信号处理的方法。2. 基本原理(1)、波束形成的物理机理波束形成的目的是使多阵元构成的基阵经过适当处理得到在预定方向的指 向性，主要内容是多元阵如何形成指向性波束及如何控制波束指向的技术，主要包括发射系统和接收系统的波束形成。发射系统的波束形成可以达到能量聚焦以增大探测距离及定向的目的，常见的如激光、手电筒、凸透镜聚焦、相控阵多普勒计程仪等设备就是应用这个原理。接收系统的波束形成的目的是获得预定方向的信号以控制其他方向的信号和干 扰、准确定向和在形成多个波束时分辨多个目标，像望远

2、镜等就具有类似功能。发射系统的波束形成针对物理空间，而接收系统针对观测空间。本实验的重 点是发射系统的波束形成。Tn_ JNiN厂N 1LjL.d1T2111 10E2(2)、波束形成的基本原理狭义上，波束形成技术是指将按一定几何形状(直线，圆柱等)排列的多元基阵各阵元的输出经过 处理(例如加权、延时、求和等)形成空间指向性的 方法。在广义上，波束形成技术是将一个多元阵经适 当处理时期对某些空间方向的声波具有所需响应的 方法。在物理上等效于空间滤波器，可以滤去空间某 些方位的信号，只让指定方位的信号通过。对于任意阵，其波束形成分析如下：如右图所示，以阵元0为参考点，信号到达其他 各阵元的提前时

3、间为：* 1， 2，I ,- N J.因此，将各阵元接收到的信号分别延时1,川,-N 丄后相加，即可实现同向相加的输出，在该方向的输出幅度最大。波束形成后基阵的指向性图如下所示：总之，波束形成的基本原理可叙述为：对多元阵阵元接收信号进行时延或相 移补偿，使对预定方向的入射信号形成同相相加。延时求和波束形成是将各基元信号进行适当组合，若波束形成器的时间延时与声波入射方向相匹配时，信号产生等相位相加，相对于噪声得到增强。用数字方法实现延迟时，延迟时间只能是采样间隔的整数倍，即d sin 60丄. 山丄 一I = kA =，日0 = arc sin (k Ac / d ), k = 0, 1, 2

4、,川，也一米样间隔，R（刃0.9正确补偿时延0.80.70.60.50.4未补偿时延0.30.20.1020406080I60-80-60-40-20N 波束形成器输出为：b(n,k)八 WjX(nik), k =0, _1, _2，川i三常规波束形成器示意图如下所示:(3) 、波束形成的计算机仿真方法1) 、确定预成波束的方位；2) 、结合阵元位置确定程差;3) 、在不同的波束补偿时延；4) 、将补偿程差后的信号累加；5) 、求功率；6) 、画出不同波束的波束输出结果；7) 、估计方位。3. 实验内容a、产生带时延的阵元接收信号b、产生一定信噪比下的信号c、对信号进行补偿时延d、波束输出e、

5、估计方位4. 实验结果及数据分析首先对时延波束形成、频域波束形成及 DFT 波束形成进行比较： 实验程序：clc;clear;% 参数 c=1500;f0=1.5e3;fs=150e3;T=50e-3;SNR=20;lamda=c./f0;D=lamda./2;N=21;maxdelay=(N-1).*D./c;Nmaxdelay=round(maxdelay.*fs); alpha=(-90:90)*pi./180;px=(0:N-1).*D; theta0=50*pi./180;N_delay0=D.*sin(theta0)./c.*fs;Wn=2*1000 2000/fs; bb=fir

6、1(128,Wn,bandpass);%带通滤波器系数% 信号模型t=0:1./fs:T; sig00=sqrt(2).*cos(2.*pi.*f0.*t);for mm=1:N;N0=round(N_delay0.*(mm-1); sig0(mm,:)=zeros(1,Nmaxdelay.*2-N0) sig00 zeros(1,Nmaxdelay.*1+N0);end len=length(sig0(1,:);for mm=1:N;% 加噪声%噪声能量归一化%噪声与信噪比的转化NN=normrnd(0,1,1,len); noise=filter(bb,1,NN); noise=noise

7、/std(noise);n oise=noise/10A(SNR/20);sig(mm,:)=sig0(mm,:)+noise;end % 时域波束形成tt=(0:len-1)./fs;tt1=(Nmaxdelay+1:len-1-Nmaxdelay)./fs;for mm=1:length(alpha);for nn=1:N;tc0=(nn-1).*D.*sin(alpha(mm)./c;ss(nn,:)=interp1(tt,sig0(nn,:),tt1-tc0); end sss(mm,:)=sum(ss,1); sout(mm)=std(sss(mm,:);endfigureplot(

8、alpha.*180./pi,20.*log10(sout./max(sout)xlabel(方位角 /度)ylabel(波束输出/分贝)title(时域波束形成)grid on% 频域波束形成 vscan=exp(-j*2*pi*f0/c*px*sin(alpha);for ii=1:NSf=fft(sig(ii,:),fs);R1(ii)=Sf(f0+1);endR0=R1*vscan;yy=abs(real(R0); yy=yy./max(yy);figureplot(alpha.*180./pi,20.*log10(yy)xlabel(方位角 /度)ylabel(波束输出/分贝)tit

9、le(频域波束形成)grid on% DFT波束形成for ii=1:NSf=fft(sig(ii,:),fs);R1(ii)=Sf(fO+1);endnfft=1024;Beam_DFT=fftshift(abs(fft(R1, nfft);Beam_DFT=20*log10(Beam_DFT./max(Beam_DFT); nn仁n fft/2+1: nfft,1: nfft/2;angle = asi n(lamda/D*- nfft/2: nfft/2-1/nfft)*180/pi; figureplot(a ngle,Beam_DFT)xlabel(方位角 /度)ylabel(波束输

10、出/分贝)title(DFT波束形成) grid on实验结果；时域波束形成-100-80-60-40-2002040方位角/度60 80 1005-O 5O-25-2O-3贝什岀输束波频域波束形成贝什岀输束波-1C0-10-20-30-40-50)0-80-60-40-20020方位角/度DFT波束形成406080100T/ i11Pa4沁AA 卩1 4Ii11111h1 fI1I 1|11111111111111贝什岀输束波-100-80-60-40-20020406080100方位角/度下面对时延波束形成中参数进行变化比较：主要参数为方位角 乙、信噪比SNR和阵元数N，上述为0 = 50

11、 SN R =20， N = 2 1E =30, SNR =20, N =21时（此时及以下只需改变程序中的相应参数，其它部分没有变化，因此程序不再重述，只给出相应结果）：贝疔出输束波时域波束形成-5-105O 52O-3-40-100-80-60-40-20020406080100方位角/度% =30、， SN R =0, N =21 时:八1 1 1AI1 ”，1 IIIIi ； ii 11/ /.A A i ; I Iba Aiii n i i h：i ii1,II 1 1 i! M!l II 1JI1,; |1I1:VE ff 1F 1ii1I u0-5-10时域波束形成5 0 5 0

12、 51 2 2 3 3 - - - - - 贝疔出输束波-40-100-80-60-40-2002040方位角/度60 80 100E =3 0； , SNR =0, N =40 时:贝什岀输束波1|.1I111 !1 ；|i.pl-!ft H fi山口HII 111!|. .1 1 111 11r11110-5-10-15-20-25-30时域波束形成-35-40-100-80-60-40-2002040方位角/度60 80 100-455. 结论1、 由以上各图可以看出，波束输出都有最大值且都在预定方位角-0处，即在九 角度上基阵的输出最大，具有最大的自然指向性，这也就是波束形成的目的所在。2、由前三个图形可以得到三种波束形成方法的比较如下：相对来说，时域波束形成的主波束幅值突出较明显，但其主波束宽度较宽，在实际应用时会影响精确 度；频域波束形成的主波束宽度最窄，可以达到较高的角度分辨率，但其主波束 在幅值上与其它波束相差不大，不宜用于幅值分辨率小的场合；DFT波束形成的主波束在幅值上与其它波束相差也不大，而且主波束宽度也较宽。3、改变方位角时（即第一图和第四图），可以很清晰的看到在各个不同角度上的 波束输出，这就是波束形成的目