一种多普勒效应消除方法与流程

本发明涉及物理专业中噪声类领域多普勒效应消除方法，更具体地说是在时域—波数域消除多普勒效应的一种方法。

背景技术：

噪声分析是机械零部件状态监测和故障诊断的一个重要方法。但当机械零部件处于移动状态时，其诊断的有效性将大大降低，主要原因就是机械零部件与传声器之间的相对运动引起的多普勒效应。多普勒效应会对采集到的声信号造成幅值畸变和频率偏移，这会对诊断效果造成严重影响。因此如何准确地消除多普勒效应对移动机械零部件的诊断有重要意义。目前消除多普勒效应的方法主要有频域方法和时域方法。频域方法能够准确地消除多普勒效应，但通常只能对稳态信号进行计算。时域方法直接在时域进行计算，不仅可以计算稳态信号，还可对非稳态信号进行计算。目前消除多普勒效应的时域计算方法主要有重采样方法，非线性时间映射方法以及基于完整morse声学模型的去多普勒效应方法。但这些方法在消除多普勒效应时都有一些缺陷，比如重采样方法需要对信号进行插值拟合计算，计算效率较低，非线性时间映射方法和基于完整morse声学模型的去多普勒效应方法需要声源的初始位置信息，使其在实际应用中受到限制。

技术实现要素：

本发明为弥补重采样方法的计算效率低以及非线性时间映射方法和基于完整morse声学模型的去多普勒效应方法需要声源初始位置信息的不足，提供一种高效、无需声源初始位置的多普勒效应消除方法。

本发明为解决技术问题采用的技术方案是：

本发明多普勒效应消除方法的特点是按如下步骤进行：

步骤1、在空间中建立两个坐标系，一个是静止的固定坐标系of(xf,yf,zf)，另一个是随声源以相同速度运动的移动坐标系om(xm,ym,zm)；两个坐标系原点在t＝0时刻位于同一位置；在空间中放置一声源，声源在xfofyf平面内做直线移动；在移动声源辐射声场中设置固定坐标系下的测量面，在测量面上等间隔均匀分布n个测量点；n个测量点中第n个测量点在固定坐标系下坐标为第n个测量点在移动坐标系下坐标为第n个测量点在两个坐标系下坐标的相对关系为：

其中vx表示声源移动速度在xf轴正方向上的分量，vy表示声源移动速度在yf轴正方向上的分量；

步骤2、将n个传声器一一对应地放置在所述n个测量点处进行测量，获得声压信号，包括：

在接收时刻t，各测量点处在固定坐标系下的声压信号一一对应为：

将所有n个测量点处的声压信号pf表征为：

在接收时刻t，各测量点处在移动坐标系下的声压信号一一对应为：

将所有n个测量点在移动坐标系下声压信号pm表征为：

在接收时刻t，各测量点处不含多普勒效应的声压信号一一对应为：

将所有n个测量点不含多普勒效应的声压信号pd表征为：

步骤3、根据n个测量点在两个坐标系下坐标的相对关系，将n个测量点在移动坐标系下的声压信号pm表征为：

步骤4、根据测量点处声压信号在任一坐标系下不变的先验条件pm＝pf，计算获得各测量点含有多普勒效应的声压信号pf与不含多普勒效应的声压信号pd在时域—波数域的相对关系表征为：f(pf)＝exp(-jkxvxt)×exp(-jkyvyt)×f(pd)，其中，f表示空间傅立叶变换运算，exp表示指数运算，kx表示在xf方向上的波数划分，ky表示在yf方向上的波数划分，j表示虚数单位；

步骤5、根据各测量点测得的含有多普勒效应的声压信号pf，结合数值傅立叶变换运算，计算获得各测量点不含多普勒效应的声压信号表征为：

pd＝f^-1[exp(jkxvxt)×exp(jkyvyt)×f(pf)]，

其中f^-1表示空间逆傅立叶变换运算，从而实现多普勒效应的消除。

本发明多普勒效应消除方法的特点也在于：所述n个测量点是在平面上等间隔均匀分布。

本发明多普勒效应消除方法的特点也在于：声源移动速度为固定值或随时间变化。

本发明多普勒效应消除方法的特点也在于：所述声源为一个或多个，声源信号为稳态信号或非稳态信号。

本发明多普勒效应消除方法的特点也在于：所述空间坐标系为三维坐标系或二维坐标系。

与已有技术相比，本发明有益效果体现在：

1、本发明方法无需插值拟合计算，计算效率较高；

2、本发明方法无需声源初始位置的先验信息，能够准确消除未知初始位置声源的多普勒效应。

附图说明

图1为本发明方法中固定坐标系、移动坐标系、声源s及测量点位置示意图；

图2a为在声源移动时测量点a接收时域信号与利用本发明方法消除多普勒效应后的时域信号的对比，图2a中以点表征的曲线为实际测量得到的时域信号，以实线表征的曲线为采用本发明方法消除多普勒效应后的时域信号；

图2b为在声源移动时测量点a接收信号的频谱与利用本发明方法消除多普勒效应后的信号频谱的对比，图2b中以点表征的曲线为实际测量得到的信号频谱，以实线表征的曲线为采用本发明方法消除多普勒效应后的信号频谱；

图2c为在声源移动时测量点b接收时域信号与利用本发明方法消除多普勒效应后的时域信号的对比，图2c中以点表征的曲线为实际测量得到的时域信号，以实线表征的曲线为采用本发明方法消除多普勒效应后的时域信号；

图2d为在声源移动时测量点b接收信号的频谱与利用本发明方法消除多普勒效应后的信号频谱的对比，图2d中以点表征的曲线为实际测量得到的信号频谱，以实线表征的曲线为采用本发明方法消除多普勒效应后的信号频谱。

具体实施方式

本实施例中多普勒效应消除方法是按如下步骤进行：

步骤1、如图1所示，在空间中建立两个坐标系，一个是静止的固定坐标系of(xf,yf,zf)，以f表示固定坐标系，另一个是随声源以相同速度运动的移动坐标系om(xm,ym,zm)，以m表示移动坐标系；两个坐标系原点在t＝0时刻位于同一位置；在空间中放置一声源，声源在xfofyf平面内做直线移动，具体实施中，声源s以速度v做匀速直线移动；在移动声源辐射声场中设置固定坐标系下的测量面，在测量面上等间隔均匀分布n个测量点，具体实施中，测量点组成一个平面阵列，并与xfofyf平面平行；n个测量点中第n个测量点在固定坐标系下坐标为第n个测量点在移动坐标系下坐标为第n个测量点在两个坐标系下坐标的相对关系为：

其中vx表示声源移动速度在xf轴正方向上的分量，vy表示声源移动速度在yf轴正方向上的分量。

步骤2、将n个传声器一一对应地放置在n个测量点处进行测量，获得声压信号，包括：

在接收时刻t，各测量点处在固定坐标系下的声压信号一一对应为：

将所有n个测量点处的声压信号pf表征为：

在接收时刻t，各测量点处在移动坐标系下的声压信号一一对应为：

将所有n个测量点在移动坐标系下声压信号pm表征为：

在接收时刻t，各测量点处不含多普勒效应的声压信号一一对应为：

将所有n个测量点不含多普勒效应的声压信号pd表征为：以d表示消除多普勒效应后的信号。

步骤3、根据n个测量点在两个坐标系下坐标的相对关系，将n个测量点在移动坐标系下的声压信号pm表征为：

步骤4、根据测量点处声压信号在任一坐标系下不变的先验条件pm＝pf，计算获得各测量点含有多普勒效应的声压信号pf与不含多普勒效应的声压信号pd在时域—波数域的相对关系表征为：f(pf)＝exp(-jkxvxt)×exp(-jkyvyt)×f(pd)，其中，f表示空间傅立叶变换运算，具体实施中：

exp表示指数运算，kx表示在xf方向上的波数划分，ky表示在yf方向上的波数划分，j表示虚数单位。

步骤5、根据各测量点测得的含有多普勒效应的声压信号pf，结合数值傅立叶变换运算，计算获得各测量点不含多普勒效应的声压信号表征为：

pd＝f^-1[exp(jkxvxt)×exp(jkyvyt)×f(pf)]，

其中f^-1表示空间逆傅立叶变换运算，具体实施中：

π为圆周率，从而实现多普勒效应的消除。

为了进行空间傅立叶变换运算，n个测量点是在平面上等间隔均匀分布，声源移动速度为固定值或随时间变化，声源为一个或多个，声源信号为稳态信号或非稳态信号，空间坐标系为三维坐标系或二维坐标系。

具体实施中，声源s采用一个单极子，在固定坐标系of(xf,yf,zf)中位于(0.5m,0.5m,0m)，声源移动速度v在xf轴正方向上的分量为vx＝60m/s，在yf轴正方向上的分量为vy＝60m/s。测量阵列位于zf＝0.2m的平面上，阵列包含121×121个等间隔均匀分布的测量点，相邻两个测量点之间距离为0.05m。声源s辐射三个频率组合的正弦信号，其表达式为：

s(t)＝sin(2πf1t)+sin(2πf2t)+sin(2πf3t)(1)

在式(1)中，频率f1＝960hz，f2＝1000hz，f3＝1040hz。时域信号采样频率为10.24khz，采样点数为512。

为检验本发明方法消除多普勒效应的效果，选取了两个测量点，即测量点a和测量点b，其位置分别为a(0.5m,0.5m,0.2m)、b(1m,1m,0.2m)。图2a和图2c分别表示测量点a和测量点b接收时域信号与利用本发明方法消除多普勒效应后的时域信号的对比，图中以点表征的曲线为实际测量得到的时域信号，以实线表征的曲线为采用本发明方法消除多普勒效应后的时域信号；图2b和图2d分别表示测量点a和测量点b接收信号的频谱与利用本发明方法消除多普勒效应后的信号频谱的对比，图中以点表征的曲线为实际测量得到的信号频谱，以实线表征的曲线为采用本发明方法消除多普勒效应后的信号频谱。比较图中的实线和点线可以看出，采用本发明方法可以很好地修正多普勒效应带来的幅值畸变和频率偏移，消除测量信号中的多普勒效应。