采用分布式球形阵列的球谐波域声场分离方法

为从测得的混合声场信号中提取出需要的目标声场,提出一种基于分布式球形传声器阵列的声场分离方法。该方法依据声场的球谐波分解,利用阵列各传声器采集到的声场声压信号,获得目标声场与干扰声场的球谐波展开系数,进而估计目标声场。该方法利用声场以不同中心展开的球谐波系数之间的变换关系,直接建立传声器测量声压信号与整体坐标系下声场展开系数的方程,与传统的分布式球形阵列声场分离方法,即先求解局部坐标系下声场展开系数,再变换为整体展开系数的方法,进行比较。分别通过数值仿真和实验说明了提出方法的有效性。结果表明:该方法能够从混合声场中较准确地估计出目标声场,并且在干扰声场能量增大时,保持了较小的声场估计误差,相比于传统方法误差增加更少。     

瞬态波位移场计算方法在相控阵声场模拟中的实验验证

通过实验研究了声束偏转不同角度和声场半径变化时线性相控阵的声场特性,并与基于瞬态波位移解析解的相控阵声场计算方法得到的理论结果进行了对比。首先研究了基于半圆圆弧离面瞬态位移的声场计算方法,并进行了数值计算可靠性的优化。然后建立了一套相控阵瞬态波位移场的测量系统,以半圆形钢板为试件,纵波传感器为接收传感器,测量了偏转角度为-10°,-30°,-50°,-70°,及声场半径为50 mm,75 mm,100 mm和125 mm时线性超声相控阵的声束指向性;同时以基于瞬态波位移解析解的相控阵声场计算方法,计算了相同条件下该相控阵的声束指向性。最后将实验结果与理论结果进行了对比,取得了较好的一致。说明基于瞬态波位移的相控阵声场计算方法可以较好地计算线阵的声场,对于声束偏转较大角度和声场半径变化时也同样有效。      

频率自适应最优权重阵列干涉条纹处理技术

考虑长基线水平阵列波束形成对声场模态呈现的滤波特征,建立了一种随频率改变滤波通带的阵列权矢量设计方法。基于简正波理论,结合声场波束形成特点,在分析了常规波束形成信号LOFAR谱图干涉条纹的清晰度和条纹结构不同于单水听器输出信号条纹现象基础上,采用线性等式约束的二次优化模型,给出了一种频率自适应的最优权重估计。数值仿真验证和试验数据分析表明,利用该方法设计的最优权重阵列波束处理,可以滤波出SRBR(Surface-Reflected Bottom-Reflected)或N-SRBR(Non-Surface-Reflected Bottom-Reflected)模态组成的波束形成信号,信号谱图干涉条纹斜率与理论分析基本一致。      

匀加速运动声源的声场建模及其数值算法研究

针对匀加速运动点声源的声场特性与其运动状态密切相关这一问题,提出匀加速直线运动状态下点声源的声场计算方法。利用此方法建立了匀加速直线运动时点声源的声压模型,对模型中的关键参数声矢量R进行数值解析,并对声压进行数值分析仿真,得出匀加速直线运动时固定接收点的声压数值计算方法。用此方法对固定接收点位置的匀加速点声源声压进行声场建模,结果表明:在声源接近接收者一定距离以后,声压明显增大;在此距离之外,距离对声压的影响不大。     

单层传声器阵列重建相干声场的近场声全息方法研究

针对传声器阵列两侧存在相干声源的非自由声场重建问题,提出基于球面谐波函数扩展近场声全息理论的相干声场重建方法。该方法在已知测量面两侧声源几何位置时,使用单层传声器阵列获取测量面处的声压分布,通过最小二乘法获得与目标声源和干扰噪声源响应对应的最优球波函数扩展项数和最优系数向量,结合测点位置的空间坐标进行声波分解,并分别重建出各声源在测量面上的声压分布。为了验证方法的有效性,分别给出了相干噪声源为球形声源和非球形声源的仿真验证,并在全消声室内对双扬声器产生的相干声场的重建进行了实验验证。结果表明:该方法对球形声源和非球形声源干扰下的声场重建都具有较好的效果,球形声源干扰下的重建精度更高。      

利用宽带声场频率—掠射角干涉结构的深海直达声区目标深度估计方法

针对典型深海环境中宽带声源的深度分辨问题,通过研究深海声场随频率起伏的干涉结构与垂直线阵频域波束输出图中的干涉结构,给出一种直达声区内可区分多水下目标的宽带声源深度估计方法.该方法以近水面目标的射线声场模型为基础,推导出近海面宽带声源接收声场的波束输出表达式,阐明了频域波束输出图中干涉结构与声源深度的对应关系.然后利用...     

高精度球面阵聚焦声源定位方法研究

提出了一种高精度高空间分辨率球面阵聚焦声源定位方法——虚拟源法。该方法通过球面阵波束扫描获得实际声源的空间聚焦谱,并假定各扫描点为虚拟声源,将实际声源聚焦谱看作是全体虚拟源共同作用的结果,由此得到各虚拟源对声场的贡献量,从而可实现声源精确定位。仿真研究分析了频率,阵列孔径,声场模态阶数,信噪比等参数对声源定位性能的影响,并与常规算法进行对比。结果显示,该方法不受频率和阵列孔径的限制,避免了空间“混淆”,能够进行高精度高分辨率声源定位,并具有良好的背景噪声抑制能力。      

A high resolution sonar for sea-bed imaging

A sonar system for making three-dimensional topographic images of the sea-bed is described. The system applies a sectorial scan within each range cell by electronic beam forming. It is designed to resolve 128 beams each separated by 0.5° to cover a 64° sector. The beam forming is facilitated by a real-time Fourier transform of the signal samples from the multielement receiving transducer array. The Fourier transform is implemented as a chirp transform using surface acoustic wave reflective array devices.The design and measured performance of the transducers, and the beam former are discussed in some detail.     

Matched field processing with data-derived modes

The authors demonstrate MFP using data-derived modes and the sound speed profile, using no a priori bottom information. Mode shapes can be estimated directly from vertical line array data, without a priori knowledge of the environment and without using numerical wave field models. However, it is difficult to make much headway with data-derived modes alone, without wave numbers, since only a few modes at a few frequencies may be captured, and only at depths sampled by the array. Using a measured sound speed profile, the authors derive self-consistent, complete sets of modes, wave numbers, and bottom parameters from data-derived modes. Bottom parameters enable modes to be calculated at all frequencies, not just those at which modes were derived from data. This process is demonstrated on SWellEx-96 experiment data. Modes, wave numbers, and bottom parameters are derived from one track and MFP based on this information is demonstrated on another track.     

Accelerated beam tracing algorithm

Determining early specular reflection paths is essential for room acoustics modeling. Beam tracing algorithms have been used to calculate these paths efficiently, thus allowing modeling of acoustics in real-time with a moving listener in simple, or complex but densely occluded, environments with a stationary sound source. In this paper, it is shown that beam tracing algorithms can still be optimized by utilizing the spatial coherence in path validation with a moving listener. Since the precalculations required for the presented technique are relatively fast, the acoustic reflection paths can be calculated even for a moving source in simple cases. Simulations were performed to show how the accelerated algorithm compares with the basic algorithm with varying scene complexity and occlusion. Up to two-orders of magnitude speed-up was achieved.     

线阵波束形成声强缩放方法估算声源的辐射声功率

为了解决波束形成声源识别过程中声源辐射声功率定量计算的问题,给出了阵型简洁、便于组合的线阵声强缩放模型。通过推导线阵的声强缩放系数,建立起线阵波束输出结果与声源辐射声功率之间的换算关系。无论是线阵还是平面阵的声强缩放方法,对于偏离阵列中心位置较远处的声源进行辐射声功率估算时都存在较为明显的误差。通过理论推导和仿真模拟计算,研究了同一单极子点声源在不同位置处的声功率估算偏差随频率、幅度的变化规律,发现该估算偏差只与声源偏离位置有关,而与声源自身的强度信息无关的结论,据此给出了相应的声功率估算修正方法。半消声室实验结果和声压法测量结果对比表明:修正后的线阵声强缩放方法用于中高频声源的辐射声功率计算时,单频声源的估算误差不超过1.0 dB,宽带声源的估算误差不超过1.8 dB。      

An analytical model for a dodecahedron loudspeaker applied to the design of omni-directional loudspeaker arrays

The dodecahedron (12 sided) loudspeaker is used extensively for acoustic measurements as an approximation to an omni-directional sound source, for example in building acoustics and automotive applications. In many instances its use is favoured above other omni-directional source approximations e.g. open pipe because of its greater power. However, analysis by other authors has shown that the radiation deviates from omni-directional above a certain frequency. In this paper an analytical model of the sound pressure field generated by a dodecahedron loudspeaker is constructed. The sound pressure field generated by the model is then compared with measured data and good agreement is demonstrated. Having established and validated the model, it is adapted to form alternative spherical loudspeaker arrays based on Platonic and Archimedean solid geometry and also a number of spherical packing routines. The aim of the analysis is to determine the optimal geometry for an omni-directional source using pistonic loudspeakers.The analysis shows that constructive interference of the pressure field across the spherical baffle surface and not individual loudspeaker pistonic radiation characteristics is the most significant factor with respect to deviations from omni-directional radiation for a complete spherical array. Consequently, loudspeaker elements which physically occupy the spherical baffle surface to their maximum extent provide the most extended frequency range of omni-directional operation. The most optimal spherical array distribution is shown to be the dodecahedron.     

Near field acoustic holography with microphones on a rigid sphere (L)

Spherical near field acoustic holography (spherical NAH) is a technique that makes it possible to reconstruct the sound field inside and just outside a spherical surface on which the sound pressure is measured with an array of microphones. This is potentially very useful for source identification. The sphere can be acoustically transparent or it can be rigid. A rigid sphere is somewhat more practical than an open sphere. However, spherical NAH based on a rigid sphere is only valid if it can be assumed that the sphere has a negligible influence on the incident sound field, and this is not necessarily a good assumption when the sphere is very close to a radiating surface. This Letter examines the matter through simulations and experiments.     

Optical fiber coupler array for multi-point sound field measurements

A pressure-sensitive device using an optical fiber coupler is proposed, and is integrated into an array system to measure sound field distributions. The optical coupling between two optical fibers is modulated by the applied sound field with the use of a membrane mechanically connected to the coupler. In the array system, a large number of the optical coupler sensors are connected two-dimensionally, and sound pressure at each sensor can be read out by switching the input light sources in turn. A prototype of the basic array system based on this principle was made, and its performance was demonstrated. The signal discrimination between each element was confirmed in a 2x2 array; then, using a four-element linear array, we demonstrated the steering of beam based on the phased array method.     

An ultrasonic transducer array for velocity measurement in underwater vehicles

A correlation velocity log (CVL) is an ultrasonic navigation aid for marine applications, in which velocity is estimated using an acoustic transmitter and a receiver array. CVLs offer advantages over Doppler velocity logs (DVLs) in many autonomous underwater vehicle (AUV) applications, since they can achieve high accuracy at low velocities even during hover manoeuvres. DVLs require narrow beam widths, whilst ideal CVL transmitters have wide beam widths. This gives CVLs the potential to use lower frequencies thus permitting operation in deeper water, reducing power requirements for the same depth, or allowing the use of smaller transducers. Moving patterns in the wavefronts across a 2D receiver array are detected by calculating correlation coefficients between bottom reflections from consecutive transmitted pulses, across all combinations of receiver pairings. The position of the peak correlation value, on a surface representing receiver-pairing separations, is proportional to the vessel's displacement between pulses. A CVL aimed primarily for AUVs has been developed. Its acoustical and signal processing design has been optimised through sea trials and computer modelling of the sound field. This computer model is also used to predict how the distribution of the correlation coefficients varies with distance from the peak position. Current work seeks to increase the resolution of the peak estimate using surface fitting methods. Numerical simulations suggest that peak estimation methods significantly improve system precision when compared with simply identifying the position of the maximum correlation coefficient in the dataset. The peak position may be estimated by fitting a quadratic model to the measured data using least squares or maximum likelihood estimation. Alternatively, radial basis functions and Gaussian processes successfully predict the peak position despite variation between individual correlation datasets. This paper summarises the CVL's main acoustical features and signal processing techniques and includes results of sea trials using the device.     

A miniature class V flextensional cymbal transducer with directional beam patterns: the double-driver

The "double-driver" cymbal, a directional class V flextensional transducer, is described in this paper. Its basic structure is a bilaminar piezoelectric disk with metal caps as mechanical transformers and amplifiers. The directivity was accomplished by exciting the double-driver in a combined flexural and bending mode causing the sound pressure to add in one direction and cancel in the opposite direction. The cardioid beam pattern predicted by finite element modeling agreed well with the experimental measurements. A 3 x 3 double-driver array was constructed to demonstrate that under optimal conditions the array can provide a directional beam pattern with a front-to-back ratio of more than 20 dB.     

Acoustic centering of sources measured by surrounding spherical microphone arrays

The radiation patterns of acoustic sources have great significance in a wide range of applications, such as measuring the directivity of loudspeakers and investigating the radiation of musical instruments for auralization. Recently, surrounding spherical microphone arrays have been studied for sound field analysis, facilitating measurement of the pressure around a sphere and the computation of the spherical harmonics spectrum of the sound source. However, the sound radiation pattern may be affected by the location of the source inside the microphone array, which is an undesirable property when aiming to characterize source radiation in a unique manner. This paper presents a theoretical analysis of the spherical harmonics spectrum of spatially translated sources and defines four measures for the misalignment of the acoustic center of a radiating source. Optimization is used to promote optimal alignment based on the proposed measures and the errors caused by numerical and array-order limitations are investigated. This methodology is examined using both simulated and experimental data in order to investigate the performance and limitations of the different alignment methods.     

球形阵列自适应网格分层定位算法

刚性球形传声器阵列可以在无空间模糊的条件下进行球谐域数据处理,基于球形阵列的可控波束形成器导向响应功率(SRP)算法定位精度高,但是计算量大,计算效率低。通过将球面致密的全局网格搜索替换为分层搜索策略可以有效减小SRP算法的计算量。提出MRE-SRP算法保持球形阵列SRP定位精度的同时降低计算量,首先通过球谐域MUSIC(SH-MUSIC)算法判断入射声源的数量减小搜索区域;其次将相对熵模型引入球谐域SRP(SH-SRP)定位算法中,提取网格分层前后的信息增益,设计自适应网格选择判据,实现分层多分辨率网格的精准再细分,从而降低计算量。实验验证了所提出算法的性能,结果显示在单双声源定位中,该算法可以实现较高的定位精度,精准选择分层网格,计算量减少75%以上。     

Presence of nonlinear interference effects as a source of low frequency excitation force in vibro-acoustography

The vibro-acoustography imaging method consists of forming an image of the deformability of a tissue submitted to a low frequency fLF stress field. This sound field can be created locally by means of a focused annular array emitting two primary beams driven at two close frequencies fa and fb = fa + fLF. In the existing literature, the origin of this stress field has been identified as the low frequency radiation pressure of the two primary beams. However, this work intends to show that another contribution to this internal stress is the low frequency field distributed in the object volume and created by the nonlinear interferences of the two primary beams. This nonlinear field was calculated in the case of multiple ring annular arrays and compared with the qLF beam experimentally measured in a water tank. The agreement between the theoretical and experimental curves provides information on the possibility that this nonlinear effect takes place in vibro-acoustography.     

非均匀组织医学超声发射声场仿真

为了分析医学超声在非均匀组织中的分布特性,建立了超声发射声场的计算模型。采用交错网格有限差分法对耦合超声非线性方程进行数值求解,获得了声速、密度及非线性参数非均匀分布情况下人体腹壁组织内的超声声场分布数据。同均匀介质相比:当声速均匀而密度非均匀时,声束仍聚焦良好,焦点处声能下降了1.8 dB;当密度均匀而声速非均匀时,声束发散严重,焦点处声能下降了3.8 dB,下降程度与非均匀组织接近。组织声速在空间分布的非均匀性是导致聚焦声束能量分布畸变的主要原因。