传声器阵列特征值滤波去噪方法

作为二阶统计量的互谱矩阵(CSM)是声学成像算法的核心输入量。为增强传声器阵列的去噪表现,研究了互谱矩阵特征值滤波的机理,并提出了两种新型的特征值滤波方法的设计准则:(1)声源互谱矩阵的Stein无偏风险估计(SURE收缩),即基于SURE准则的特征值软阈值收缩;(2)进一步提高声源互谱矩阵EYM (Eckart-Young-Mirsky)估计误差的最优收缩(Opt-Shrinkage)方法,即对声源特征值所含有的噪声进一步去除。然后,在3000个快照数,0 dB信噪比(SNR),100个传声器的环境下进行仿真。通过仿真,与经典的MUSIC方法比较原声源信号与去噪后信号的互谱矩阵对角线误差,其中MUSIC方法的互谱矩阵对角线误差为74.15%,SURE收缩为41.97%,最优收缩为20.62%。并通过改变快照数、声源数与SNR,对比了不同特征值滤波去噪方法的效果。从仿真结果上看,在少声源(少于40个)工况下SURE收缩方法的去噪效果优于MUSIC方法与最优收缩方法;在声源数超过42个后,最优收缩方法取得更显著的效果。最后,在3个声源、60个传声器数、-5 dB信噪比的声源定位实验中,SURE收缩与最优收缩均获得了相比MUSIC方法更好的去噪效果。研究表明在现有传声器阵列特征值滤波方法的基础上,对特征值的进一步处理可以得到更好的去噪效果。     

椭圆截面侵彻弹体爆炸特性试验研究

为研究椭圆截面侵彻弹体的爆炸特性,设计并开展了静爆威力外场试验。将质量为255 kg的弹体竖立于木质托弹架上,质心距地面高度为2 m,采用试验引信起爆弹体装药。通过航拍无人机实时拍摄整个爆炸过程,在长轴和短轴方向布置扇形效应钢板以获取破片数量及穿甲率,采用超压传感器测量距弹轴7、10和12 m处的冲击波超压,并对弹体爆炸后的宏观景象以及火球、破片和冲击波超压特性进行了详细分析。结果表明,火球演化形貌与破片散布区域关于弹体长轴和短轴呈对称分布;火球演化过程分为快速成长阶段、高温稳定阶段以及自由扩散阶段,火球尺寸在爆炸后41.7 ms达到最大,短轴和长轴方向的最大尺寸分别为21.86、19.29 m,且火球在长轴方向发生了明显的二次膨胀;短轴方向的破片尺寸小、数量多、穿甲能力强,而长轴方向的破片特性恰好相反;冲击波超压峰值、冲量及速度均随传播距离增大而不断减小。综合试验结果对比分析,认为椭圆截面侵彻弹体的非轴对称结构和非均匀壁厚对爆炸特性影响较大,是造成火球形貌及破片非轴对称分布的根本原因。     

微波消解-ICP-AES测定食品塑料包装中钛、铅、铬和镉

采用微波消解技术,电感耦合等离子体原子发射光谱法(ICP-AES)同时测定食品塑料包装中钛、铅、铬和镉.研究了消解程序和消解试剂,建立了塑料样品的微波消解方法.测定了三种不同的食品包装用塑料,方法回收率为73.3%~101.7%,相对标准偏差为0.3%~10.7%,钛、铅、铬和镉在塑料中的检出限分别为1.50、7.38...     

一种预测扩散声场在非规则形状刚性壁面附近干涉图样的方法

扩散声场会在反射边界附近形成干涉图样,研究方法包括平面波模型、简正模态分析、渐进模态分析等,但仅适用于尺度远大于声波波长的矩形声腔。提出一种预测扩散声场在非规则刚性壁面结构附近形成的干涉图样的数值方法,表明结构附近“受挡”声压的互谱矩阵取决于:(1)假定该结构在自由空间中振动辐射声音时其表面法向振速到表面及场点声压的边界元系数矩阵;(2)假定结构置于自由空间中且表面刚性时,点声源辐射声波入射到结构表面上产生的散射声场的边界元系数矩阵;(3)扩散声场均方声压。仿真表明,该途径预测的干涉图样与理论值完全吻合。该预测方法还可用于混响环境下声源附近直达声压均方值的空间分布估计,为混响环境下设备的声源定位提供帮助。     

HELS法在循环平稳声场全息重建中的理论与实验研究

Helmholtz 方程最小二乘法利用一组球面波基函数拟合声源产生的声场,根据重建和实际声压的误差最小原则,利用最小二乘法确定基函数展开的项数以及对应的权重系数,该方法具有计算效率高和需要测点少的优点,在实际工程中有很大的实用性.Helmholtz 方程最小二乘法和其他近场声全息方法一样都是针对平稳声场,对非平稳声场的分析很少.对于实际工程中经常遇到的一类特殊非平稳声场——循环平稳声场,现有的技术多以单通道信号分析为主,其高阶统计量在故障诊断领域应用较广.分析了循环平稳声场中Helmholtz方程最小二乘 关键词： 声全息 循环平稳 Helmholtz 方程 球面波     

n阶常系数线性微分方程的通解公式

ｎ阶常系数线性微分方程的通解公式曹跃颖，吴海军，李彦庆（合肥工业大学机械系９５级，合肥２３００Ｏ９）１．引言众所周知，对ｎ阶常系数线性非齐次微分方程：ｙ（ｎ）＋ｐ１ｙ（ｎ－１）＋ｐ２ｙ（ｎ－２）＋…＋ｐｎｙ＝ｆ（ｘ）（１）其中ｐ１，ｐ２，…，ｐｎ为实...     

激发圆形管道内多阶声模态的扬声器阵列控制方法

声模态发生器是通过控制扬声器阵列在管道内激发声模态波的一种装置。为了解决在管道内同时激发多个声模态的问题,研究了激发圆形管道内多阶声模态的扬声器阵列控制方法。采用轴向多圈布置的声源阵列,并调节各个声源的幅值和相位,实现同时激发包括径向声模态在内的多个声模态。同时考虑声源的周向位置和轴向位置信息,建立各个声源与多个目标模态系数之间的线性关系,运用最小二乘法求解得到激发目标多模态所需各个声源的复强度（包括幅值和相位）,所研制的高阶模态发生器以计算的声源复强度为输入量,采用数字信号系统控制扬声器输出的幅值和相位,用于实现管道内声源激发,该模态激发过程无需针对特定模态优化声源的位置。实验结果表明,所研制的模态发生器可精确激发单个或多个声模态,且目标模态系数信噪比几乎都大于10 dB。     

基于Burton-Miller边界积分方程的二维声学波动问题对角形式快速多极子边界元及其应用

论述了二维声学问题的快速多极子边界元（FMBEM）方程及实现步骤．概述了核函数展开理论,并对FMBEM的4个重要组成部分:源点矩计算、源点矩转移、源点矩至本地展开转移、本地展开转移进行了详细的描述．提出了一种有利于四叉树建立的数据结构．推导了一种比直接数值计算更精确、稳定和高效的解析源点矩计算公式．数值算例验证了FMBEM的正确性和高效性．最后,使用FMBEM对轨道二维声学辐射模型进行了模拟计算．     

利用O型环抑制圆柱绕流涡脱落的数值研究

圆柱绕流涡脱落诱发较大的振动和声,如何有效地抑制值得关注。利用大涡模拟技术求解了Navier—Stokes方程,得到了涡脱落频率,升力脉动幅值及平均阻力系数。计算表明二维模拟不能体现流动基本特征．三维计算与实验吻合较好。为了抑制涡脱落,在直径为D的圆柱表面装入间距为1D．直径为0．0167D的O型环。通过升力、速度谱分析以及柱向横截面流场分析可知．在光滑圆柱外表面加入O型环能诱发流体边界层分离,有效地抑制涡脱落现象．升力脉动和观测点速度脉动幅值几乎完全消失,阻力系数也略微降低,适合在实际工程中采用。     

水下圆柱壳结构表面振速插值方法及其辐射声场预测

提出一种圆柱壳表面速度的全局插值算法,表面速度函数周向采用Fourier级数、轴向采用分段Lagrange函数表示,利用圆柱壳结构表面有限位置处的测量数据,计算得到速度函数的级数展开式系数,然后对表面其它位置处的振速进行插值,避免了传统的局部插值和离散产生的误差。数值算例验证了该插值方法能通过有限点位置处的数据,获取连续光滑的表面速度分布,并能满足圆柱壳辐射声场预测的需求。进而,在实验条件下,使用该插值方法对双层圆柱壳水下辐射声场进行了预测,实验结果表明通过该插值方法得到壳体表面速度分布,利用快速边界元方法计算得到的预测声场与实测声场有很好的一致性,验证了该插值方法在圆柱壳辐射声场预测中的正确性和工程适用性。