技术领域
[0001] 本
发明涉及物理专业中噪声类领域多普勒效应消除方法，更具体地说是一种消除旋转声源多普勒效应的方法。
背景技术
[0002] 在实际工程应用中，噪声分析能够作为机械产品运行状态监测的重要手段。但是假如零部件处于旋转运动状态，对其工作状态进行监测的准确度会大大降低，导致该现象的主要原因是被监测的机械零部件和
传声器之间存在相对运动引起的多普勒效应，多普勒效应会对采集到的声
信号造成频散现象，最终得到的零部件故障分析结果也不可靠。因此消除多普勒效应对于
旋转机械部件的状态监测有着重要意义。目前用于消除旋转声源导致的多普勒效应的方法主要有频域和时域两种。频域方法是通过对静止的环形传声器阵列所测得的声压
频谱在周向上进行旋转模态分解，计算得到虚拟旋转传声器处的声压频谱，由于虚拟旋转传声器
和声源处于相对静止状态，所以其声压频谱中不包含多普勒效应，使用该方法能够较为准确地消除多普勒效应，但是当用于宽频噪声处理时，需要对每个
频率进行操作，因此导致计算效率低；时域方法主要是插值法，让传声器虚拟地和声源同步旋转，之后判断每一时刻下虚拟传声器
位置，借助于附近的静止传声器测得的声压，插值得到虚拟传声器的声压，即不含多普勒效应的声压结果，但是该方法的
精度与传声器的数量密切相关，需要在环形阵列上布置密集的传声器，否则得到的结果十分不理想。
发明内容
[0003] 本发明是为避免上述
现有技术所存在的不足，提供一种高效、高精度的消除旋转声源多普勒效应的方法。
[0004] 本发明解决技术问题所采用的技术方案是：
[0005] 本发明消除旋转声源多普勒效应的方法的特点是：在旋转声源
辐射声场中设置测量面，利用传声器测量获得测量面中各测点含有多普勒效应的声压信号；在空间中建立静止柱
坐标系和随声源以相同
角速度旋转的旋转柱坐标系；根据静止柱坐标系和旋转柱坐标系之间的相对运动关系，建立各测点在静止柱坐标系和旋转柱坐标系下坐标之间的相对关系；根据测量声压在静止柱坐标系和旋转柱坐标系下相等的先验条件，并结合坐标的相对关系计算出含有多普勒效应的声压角谱与不含多普勒效应的声压角谱之间的对应关系；根据声压角谱之间的对应关系以及各测点含有多普勒效应的声压，计算获得各测点不含多普勒效应的声压，从而消除多普勒效应。
[0006] 本发明旋转声源多普勒效应消除方法的特点是按如下步骤进行：
[0007] 步骤a、在旋转声源辐射声场中设置由N个均匀分布且静止不动的传声器构成的环形测量阵列，所述环形测量阵列所在平面为测量面，所述测量面与声源旋转平面平行，且测量面中心与声源旋转中心的连线垂直于声源旋转平面，利用传声器测量获得t时刻旋转声源所辐射的声压信号，其为含有多普勒效应的声压信号；
[0008] 步骤b、在空间中分别建立静止柱坐标系 和以旋转角速度Ω与声源同步旋转的旋转柱坐标系 所述静止柱坐标系和旋转柱坐标系的原点均为环形测量阵列的中心点，静止柱坐标系和旋转柱坐标系的z轴方向为环形测量阵列的轴线方向，静止柱坐标系和旋转柱坐标系在测量的初始时刻，即t＝0的时刻，是重合的；
[0009] 步骤c、将环形测量阵列中第n个传声器在静止柱坐标系下的坐标表示为将环形测量阵列中第n个传声器在旋转柱坐标系下的坐标表示为根据静止柱坐标系和旋转柱坐标系之间的相对运动关系，获得所述第n个传声器在静止柱坐标系和旋转柱坐标系下坐标之间的相对关系如式(1)：
[0010]
[0011] 步骤d、按如下方式表征各声压信号：
[0012] 在t时刻，静止柱坐标系下N个传声器测得的声压信号一一对应为：将静止柱坐标系下N个传声器测得的声压信号
表征为式(2)：
[0013]
[0014] 在t时刻，旋转柱坐标系下N个传声器测得的声压信号一一对应为：将旋转柱坐标系下N个传声器测得的声压信号
表征为式(3)：
[0015]
[0016] 根据式(1)和式(3)，将N个传声器在旋转柱坐标系下测得的声压信号表征为式(4)：
[0017]
[0018] 步骤e、根据各传声器测得的声压信号在任一坐标系下不变的先验条件，则有式(5)：
[0019]
[0020] 对式(5)进行 角正变换获得静止柱坐标系下含有多普勒效应的声压角谱Ps(rs,m,zs,t)和旋转柱坐标系下的不含多普勒效应的声压角谱Pr(rs,m,zs,t)之间的对应关系如式(6)：
[0021] Ps(rs,m,zs,t)＝Pr(rs,m,zs,t)e-jmΩt (6)
[0022] 式(6)中，j表示虚数单位；m表示周向模态数；
[0023] 步骤f、式(6)中声压角谱Ps(rs,m,zs,t)为已知量，是通过对各传声器测得的声压信号 进行 角正变换所获得，则将式(6)变换为式(7)，式(7)中声压角谱Pr(rs,m,zs,t)为待确定量：
[0024] Pr(rs,m,zs,t)＝Ps(rs,m,zs,t)ejmΩt (7)
[0025] 通过对式(7)进行 角逆变换，获得各测点不含多普勒效应的声压信号如式(8)：
[0026]
[0027] 实现旋转声源多普勒效应的消除。
[0028] 本发明旋转声源多普勒效应消除方法的特点也在于：各时刻声源旋转角速度为已知的匀速变化，或为已知的非匀速变化。
[0029] 本发明旋转声源多普勒效应消除方法的特点也在于：所述声源为一个或多个，且每个声源旋转平面均与所述测量面平行，每个声源旋转中心与测量面中心的连线垂直于声源旋转平面；声源为多个时，所有声源的旋转角速度变化相一致。
[0030] 本发明旋转声源多普勒效应消除方法的特点也在于：所述声源信号为稳态信号或非稳态信号。
[0031] 与已有技术相比，本发明有益效果体现在：
[0032] 1、本发明方法在处理宽频噪声时因采用时域处理方法，可通过一次运算获得整个频带的声场信息，而无需对每个频率进行操作，因此相较于频域处理方法，本发明方法效率高；
[0033] 2、本发明方法因采用解析的 角变换方法，无需使用数值插值方法，因此相较于时域插值方法可以在相同传声器数量条件下得到更为精确的重建结果。
附图说明
[0034] 图1为本发明中静止柱坐标系、旋转柱坐标系、声源以及测点位置示意图；
[0035] 图2a为本发明方法中测点A处声压时域信号图；
[0036] 图2b为本发明方法测点A处声压频域信号图；
[0037] 图3a为本发明方法测点B处声压时域信号图；
[0038] 图3b为本发明方法测点B处声压频域信号图；
[0039] 图4a为本发明方法在0.002s时所有测点的声压信号图；
[0040] 图4b为本发明方法在0.004s时所有测点的声压信号图；
[0041] 图5为本发明方法与时域插值方法在使用相同传声器数量条件下的声压信号重建误差对比图。
具体实施方式
[0042] 本
实施例中消除旋转声源多普勒效应的方法的特点是：在旋转声源辐射声场中设置测量面，利用传声器测量获得测量面中各测点含有多普勒效应的声压信号；在空间中建立静止柱坐标系和随声源以相同角速度旋转的旋转柱坐标系；根据静止柱坐标系和旋转柱坐标系之间的相对运动关系，建立各测点在静止柱坐标系和旋转柱坐标系下坐标之间的相对关系；根据测量声压在静止柱坐标系和旋转柱坐标系下相等的先验条件，并结合坐标的相对关系计算出含有多普勒效应的声压角谱与不含多普勒效应的声压角谱之间的对应关系；根据声压角谱之间的对应关系以及各测点含有多普勒效应的声压，计算获得各测点不含多普勒效应的声压，从而消除多普勒效应。
[0043] 参见图1，本实施例中消除旋转声源多普勒效应的方法按如下步骤进行：
[0044] 步骤a、在旋转声源辐射声场中设置由N个均匀分布且静止不动的传声器构成的环形测量阵列，环形测量阵列所在平面为测量面，测量面与声源旋转平面平行，且测量面中心与声源旋转中心的连线垂直于声源旋转平面，利用传声器测量获得t时刻旋转声源所辐射的声压信号，其为含有多普勒效应的声压信号；
[0045] 步骤b、在空间中分别建立静止柱坐标系 和以旋转角速度Ω与声源同步旋转的旋转柱坐标系 静止柱坐标系和旋转柱坐标系的原点均为环形测量阵列的中心点，静止柱坐标系和旋转柱坐标系的z轴方向为环形测量阵列的轴线方向，静止柱坐标系和旋转柱坐标系在测量的初始时刻，即t＝0的时刻，是重合的；
[0046] 步骤c、将环形测量阵列中第n个传声器在静止柱坐标系下的坐标表示为将环形测量阵列中第n个传声器在旋转柱坐标系下的坐标表示为根据静止柱坐标系和旋转柱坐标系之间的相对运动关系，获得第n个传声器在静止柱坐标系和旋转柱坐标系下坐标之间的相对关系如式(1)：
[0047]
[0048] 步骤d、按如下方式表征各声压信号：
[0049] 在t时刻，静止柱坐标系下N个传声器测得的声压信号一一对应为：将静止柱坐标系下N个传声器测得的声压信号
表征为式(2)：
[0050]
[0051] 在t时刻，旋转柱坐标系下N个传声器测得的声压信号一一对应为：将旋转柱坐标系下N个传声器测得的声压信号
表征为式(3)：
[0052]
[0053] 根据式(1)和式(3)，将N个传声器在旋转柱坐标系下测得的声压信号表征为式(4)：
[0054]
[0055] 步骤e、根据各传声器测得的声压信号在任一坐标系下不变的先验条件，则有式(5)：
[0056]
[0057] 对式(5)进行 角正变换获得静止柱坐标系下含有多普勒效应的声压角谱Ps(rs,m,zs,t)和旋转柱坐标系下的不含多普勒效应的声压角谱Pr(rs,m,zs,t)之间的对应关系如式(6)：
[0058] Ps(rs,m,zs,t)＝Pr(rs,m,zs,t)e-jmΩt (6)
[0059] 式(6)中，j表示虚数单位；m表示周向模态数；
[0060] 步骤f、式(6)中声压角谱Ps(rs,m,zs,t)为已知量，是通过对各传声器测得的声压信号 进行 角正变换所获得，则将式(6)变换为式(7)，式(7)中声压角谱Pr(rs,m,zs,t)为待确定量：
[0061] Pr(rs,m,zs,t)＝Ps(rs,m,zs,t)ejmΩt (7)
[0062] 通过对式(7)进行 角逆变换，获得各测点不含多普勒效应的声压信号如式(8)：
[0063]
[0064] 实现旋转声源多普勒效应的消除。
[0065] 具体实施中，各时刻声源旋转角速度为已知的匀速变化，或为已知的非匀速变化；声源为一个或多个，且每个声源旋转平面均与测量面平行，每个声源旋转中心与测量面中心的连线垂直于声源旋转平面；声源为多个时，所有声源的旋转角速度变化相一致；声源信号为稳态信号或非稳态信号。
[0066] 针对本发明方法按如下方式进行检验
[0067] 环形测量阵列由48个均匀分布且静止不动的传声器组成，测量阵列位于静止柱坐标系zs＝0m的平面上，测量阵列的半径为0.3m，测量阵列的中心即为静止柱坐标系的原心。旋转声源采用一个单极子，单极子在静止柱坐标系中的初始坐标为(0.2m,0rad,0.3m)，单极子以角速度Ω＝12000πrad/s绕zs轴旋转，单极子声源的旋转半径为0.2m，单极子声源的旋转平面位于zs＝0.3m的平面上。声源辐射信号S(t)为正弦信号，如式(9)所表达：
[0068] S(t)＝sin(2πft) (9)
[0069] 在式(9)中，频率f＝100Hz；时域信号
采样频率为12.8kHz，采样点数设定为128。
[0070] 为了检验本发明方法消除旋转声源多普勒效应的效果，选取了两个测点，即测点A(0.3m,0rad,0m)和测点
[0071] 图2a和图3a分别表示为本发明方法测点A和测点B所在位置处声压时域信号图；图2a中，实线即曲线a1为直接测量的声压时域信号，其包含旋转声源的多普勒效应；点线即曲线b1表示旋转声源所辐射的理论声压时域信号；虚线曲线c1表示采用本发明方法消除多普勒效应后的声压时域信号；图3a中，实线即曲线a3表示直接测量的声压时域信号，包含旋转声源的多普勒效应；点线即曲线b3表示旋转声源所辐射的理论声压时域信号，虚线即曲线c3表示采用本发明方法消除多普勒效应后的声压时域信号。由图2a和图3a可见本发明方法在时域内能很好地修正多普勒效应带来的声压幅值畸变，消除测量信号的多普勒效应。
[0072] 图2b和图3b分别表示为本发明方法测点A和测点B所在位置处声压频谱图；图2b中实线即曲线a2表示直接测量的声压频谱，包含旋转声源的多普勒效应；点线即曲线b2表示旋转声源所辐射的理论声压频谱；虚线即曲线c2表示采用本发明方法消除多普勒效应后的频谱；；
[0073] 图3b中实线即曲线a4表示直接测量的声压频谱，包含旋转声源的多普勒效应；点线即曲线b4表示旋转声源所辐射的理论声压频谱，虚线即曲线c4表示采用本发明方法消除多普勒效应后的频谱；由图2b和图3b可见本发明方法在频域内能很好地修正多普勒效应带来的声压幅值畸变和
频率偏移，消除测量信号的多普勒效应。
[0074] 图4a和图4b分别表示为本发明方法在0.002s和0.004s接收时刻各个测点所测量的声压信号；图4a中实线即曲线a5表示直接测量的声压信号，包含旋转声源的多普勒效应；点线即曲线b5表示旋转声源所辐射的理论声压信号；虚线即曲线c5表示采用本发明方法消除多普勒效应后的声压信号；图4b中实线即曲线a6表示直接测量的声压信号，包含旋转声源的多普勒效应；点线即曲线b6表示旋转声源所辐射的理论声压信号；虚线即曲线c6表示采用本发明方法消除多普勒效应后的声压信号；由图4a和图4b可见采用本发明方法消除多普勒效应后的声压信号与旋转声源所辐射的理论声压信号基本一致，因此本发明方法可以在空间域内消除测量信号的多普勒效应。
[0075] 图5表示为本发明方法与时域插值方法在使用相同传声器数量条件下的声压信号重建误差对比图；图5中实线即曲线a7表示采用本发明方法计算得到的重建误差；点线即曲线b7表示采用时域插值方法计算得到的重建误差；由图5可见本发明方法与时域插值方法相比，可以在使用相同传声器数量条件下得到更为精确的重建结果。
