CN109194306B - 一种量化汽车噪声调制问题的方法及装置 - Google Patents

一种量化汽车噪声调制问题的方法及装置 Download PDF

Info

Publication number
CN109194306B
CN109194306B CN201810989074.0A CN201810989074A CN109194306B CN 109194306 B CN109194306 B CN 109194306B CN 201810989074 A CN201810989074 A CN 201810989074A CN 109194306 B CN109194306 B CN 109194306B
Authority
CN
China
Prior art keywords
noise
envelope
band
frequency band
determining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810989074.0A
Other languages
English (en)
Other versions
CN109194306A (zh
Inventor
廖祥凝
庞剑
杨亮
靳红英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN201810989074.0A priority Critical patent/CN109194306B/zh
Publication of CN109194306A publication Critical patent/CN109194306A/zh
Application granted granted Critical
Publication of CN109194306B publication Critical patent/CN109194306B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
    • H03H17/0213—Frequency domain filters using Fourier transforms
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/06—Non-recursive filters
    • H03H17/0614—Non-recursive filters using Delta-modulation

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种量化汽车噪声调制问题的方法及装置，以解决现有技术中对有调制的噪声的问题频段需要通过滤波回放手段进行确定以及对噪声调制问题的优化程度需要依靠工程师进行主观判断，导致耗费的时间较多的问题。本发明的量化汽车噪声调制问题的方法，包括：对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；对多个第一噪声分别对应的包络线进行傅立叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度。

Description

一种量化汽车噪声调制问题的方法及装置
技术领域
本发明涉及汽车NVH性能领域，具体是一种量化汽车噪声调制问题的方法及装置。
背景技术
汽车NVH性能是提升汽车品质感的重要因素，近年来受到了广大汽车厂商的重视。整车开发过程中，有时会出现一些噪声的调制问题。采用拟声词对这些调制现象进行描述，如加速“咳咳”声，燃油脉动“嗒嗒”声，排气尾管“框框”声等。这些问题噪声的共同点在于声音间歇存在，并且间歇时间具有规律性。以某工装车恒转速2000rpm爬坡时出现“咳咳”声为例，通过滤波回放手段发现“咳咳”声与300-400Hz的频谱特征有关，如图1所示。由图1可得，300-400Hz频带内出现间隔约16.7Hz(发动机半阶次频率)的特征带，同时由300-400Hz频带内的时域声信号可看出，声信号幅值呈周期性变化，形成了幅值调制现象，使得人耳主观感受上产生了“咳咳”声。并且可推断出“咳咳声”的时间间隔为0.06s，频率间隔约16.7Hz，与300-400Hz频谱特征带的频率间隔一致。
目前整车开发过程中，主要采用滤波回放手段来锁定汽车噪声发生调制的问题频段，这种方式需要耗费较多的时间。另外，在不同解决方案对有调制的噪声进行优化后，需要工程师进行主观评估不同解决方案对噪声调制问题的优化程度，由于不同的人具有不同的主观感受，若期望多个工程师的主观评估具有有效性，需要多个工程师共同进行主观评估，耗费了较多的人力资源和时间。
发明内容
本发明的目的在于提供一种量化汽车噪声调制问题的方法及装置，以解决现有技术中对有调制的噪声的问题频段需要通过滤波回放手段进行确定以及对噪声调制问题的优化程度需要依靠工程师进行主观判断，导致耗费时间较多的问题。
本发明的技术方案为：
本发明提供了一种量化汽车噪声调制问题的方法，包括：
对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；
对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；
对多个第一噪声分别对应的包络线进行傅立叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；
根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度。
优选地，对多个第一噪声分别进行希尔伯特变换，获得所述第一噪声对应的包络线的步骤包括：
通过第一公式：
Figure BDA0001780358960000021
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，Aij coS(2πfijt+φij)为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，φij为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
优选地，对多个第一噪声分别对应的包络线进行傅里叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度的步骤包括：
通过第二公式：
Figure BDA0001780358960000022
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure BDA0001780358960000025
其中，FFT[ej(tk：tk+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure BDA0001780358960000023
为函数向下取整；
通过第三公式：
Figure BDA0001780358960000024
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure BDA0001780358960000031
通过第四公式：
Figure BDA0001780358960000032
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
优选地，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段的步骤包括：
确定主要幅值调制度最高的一个目标频段，将所述目标频段确定为原始有调制的噪声的问题频段。
优选地，所述方法还包括：
根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
优选地，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率的步骤包括：
根据第五公式：
fmj(tk)=f，当
Figure BDA0001780358960000033
确定所述主要幅值调制度对应的调制频率fmj(tk)。
根据本发明的另一方面，本发明还提供了一种量化汽车噪声问题的装置，包括：
第一获取模块，用于对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；
第二获取模块，用于对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；
第三获取模块，用于对多个第一噪声分别对应的包络线进行傅里叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；
第一确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度。
优选地，第二获取模块具体用于：
通过第一公式：
Figure BDA0001780358960000041
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，Aij cos(2πfijt+φij)为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，φij为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
优选地，第三获取模块具体用于：
通过第二公式：
Figure BDA0001780358960000042
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure BDA0001780358960000043
其中，FFT[ej(tk：tk+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure BDA0001780358960000044
为函数向下取整；
通过第三公式：
Figure BDA0001780358960000045
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure BDA0001780358960000046
通过第四公式：
Figure BDA0001780358960000047
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
优选地，第一确定模块包括：
第一确定单元，用于确定主要幅值调制度最高的一个目标频段，将所述目标频段确定为原始有调制的噪声的问题频段。
优选地，所述装置还包括：
第二确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
优选地，根第二确定模块具体用于：
根据第五公式：
fmj(tk)=f，当
Figure BDA0001780358960000051
确定所述主要幅值调制度对应的调制频率fmj(tk)。
本发明的有益效果为：
本发明提供的方法，主要包括了多频段的滤波处理、滤波信号的包络线求取和包络线主要幅值调制度及调制频率的计算这三种技术手段，能快速锁定噪声产生调制的问题频段及其调制程度，为不同解决方案对调制问题的优化力度提供了准确有效的对比手段。问题频段和调制频率的确定有助于快速梳理噪声调制问题的机理，提供解决思路，规避了滤波回放手段的繁冗性。不同解决方案对调制问题的优化力度的有效对比，避免了问题频带三维频谱图对比的片面性和多个工程师主观判断的耗时性问题，提高了整车研发的效率。
附图说明
图1为某工装车以恒转速2000rpm爬坡时产生的噪声的频谱特征图；
图2为本发明的方法的流程示意图；
图3为对有调制的噪声进行多频段带通滤波的示意图；
图4为频段位于300-400Hz带通滤波后的第一噪声对应的包络线的示意图；
图5为有调制的噪声的主要幅值调制度对比图；
图6为300-400Hz频段的调制频率图；
图7为有调制的噪声整改后和整改前的主要幅值调制度对比图。
具体实施方式
参照图2，本发明提供了一种量化汽车噪声调制问题的方法，包括：
步骤1，对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声。
其中，对于有调制的噪声来说，其是在某种特定工况下产生的，通过布置到驾驶舱内的声压采集装置采集，例如，车辆以恒转速2000rpm爬坡时，在副驾驶员右耳处采集到的“咳咳”声。
首先，按照设计的频段生成对应的带通滤波器，例如，如表格1所划分的频段生成对应的带通滤波器。
序号 中心频率 上止频率 下止频率 序号 中心频率 上止频率 下止频率
1 50 0 100 25 1850 1710 1990
2 100 50 150 26 1990 1850 2130
3 150 100 200 27 2150 1990 2310
4 200 150 250 28 2310 2150 2470
5 250 200 300 29 2500 2310 2690
6 300 250 350 30 2690 2500 2880
7 350 300 400 31 2900 2675 3125
8 400 350 450 32 3125 2900 3350
9 450 395 505 33 3400 3125 3675
10 505 450 560 34 3675 3400 3950
11 570 510 630 35 4000 3650 4350
12 630 570 690 36 4350 4000 4700
13 700 630 770 37 4800 4350 5250
14 770 700 840 38 5250 4800 5700
15 840 765 915 39 5800 5250 6350
16 915 840 990 40 6350 5800 6900
17 1000 920 1080 41 7000 6350 7650
18 1080 1000 1160 42 7650 7000 8300
19 1170 1075 1265 43 8500 7600 9400
20 1265 1170 1360 44 9400 8500 10300
21 1370 1265 1475 45 10500 9250 11750
22 1475 1370 1580 46 11750 10500 13000
23 1600 1480 1720 47 13500 11750 15250
24 1720 1600 1840 48 15250 13500 17000
表1
在生成带通滤波器后，基于卷积计算对有原始调制的噪声进行多频段的带通滤波处理，并获得每一频段带通滤波处理后的第一噪声，在图3中的x(t)为原始有调制的噪声，xj(t)为第j个频段带通滤波后的第一噪声信号。
步骤2，对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线。
由图1可得，原始有调制的噪声在问题频段的包络线发生了周期性波动，是一种典型的幅值调制现象，其包络线可表达为一常数项与多个余弦信号的叠加形式。具体的，对于每一第一噪声分别对应的包络线是利用希尔伯特变换进行求解获得的，具体的，该步骤2包括：通过第一公式：
Figure BDA0001780358960000061
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，
Figure BDA0001780358960000071
为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，
Figure BDA0001780358960000072
为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
如图4所示，图4给出了位于300-400Hz频段带通滤波后的第一噪声对应的包络线。
步骤3，对多个第一噪声分别对应的包络线进行傅里叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度。
由第一公式可以看出，相比于常数项，某余弦信号的幅值越大，包络线波动越大，噪声调制问题越明显。因此为表征噪声调制问题的调制程度，计算包络线的幅值调制度，其步骤如下。
将包络线进行快速傅里叶变换(Fast Fourier Transform，FFT)，如第二公式所示，通过第二公式：
Figure BDA0001780358960000073
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure BDA0001780358960000074
其中，FFT[ej(tk：tk+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure BDA0001780358960000075
为函数向下取整；
通过第三公式：
Figure BDA0001780358960000076
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure BDA0001780358960000077
通过第四公式：
Figure BDA0001780358960000078
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
其中，在步骤3中，每一第一噪声对应的包络线的时间间隔设置为0.5s，即Δt＝0.5，然后，通过上述第二公式至第四公式依次进行求解，确定每一频段带通滤波后的包络线在任意时刻的主要幅值调制度。
步骤4，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度。
在通过步骤3计算后，可以获得每一频段带通滤波后的第一噪声对应的包络线的主要幅值调制度，将位于不同区间内的主要幅值调制度分别赋予不同颜色，这样，在图5中，可以通过图中显示的颜色确定不同频段对应的主要幅值调制度，并且，将主要幅值调制度最高的一个目标频段确定为本发明所需的有调制的噪声的问题频段。用户在确定问题频段时，可以通过图片上显示的颜色的不同来确定区间，例如，在图5中，在频段位于300-400Hz的这一频段对应的主要幅值调制度颜色显示为白色时，确定该频段对应的主要幅值调制度为最高，因而将300-400Hz这一频段确定为有调制的噪声的问题频段。
在确定出原始有调制的噪声机理后，提出解决方案，为了确定解决方案的整改程度，需要对整改前和整改后的主要幅度调制度进行对比分析，由于主要幅值调制度是通过不同的颜色进行显示的，在确定解决方案整改程度时，只需要根据整改前和整改后的主要幅值调制度的颜色变化即可快速确定噪声调制问题的优化程度。
优选地，所述方法还包括：
步骤5，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
优选地，步骤5包括：
根据第五公式：
fmj(tk)=f，当
Figure BDA0001780358960000081
确定所述主要幅值调制度对应的调制频率fmj(tk)。
由图5和图6可得某工装车恒转速2000rpm爬坡时，300-400Hz出现了调制频率为发动机半阶次、幅值调制度大于40％的调制现象，产生了明显的“咳咳”声。根据出现调制的问题频段及其调制频率，结合发动机各零部件的模态信息，推断出“咳咳”声的机理为：在发动机一缸激励(由调制频率确定)下，发动机轮系端引起共振(由出现调制的问题频段确定)。基于“咳咳”声的机理，提出更换皮带的方案。在对皮带进行更换后，工程师主观评价“咳咳”声明显减弱，采用本发明方法评估解决方案对调制问题的优化力度，结果如图7所示。由图7可看出，更换皮带的方案使得问题频段的主要幅值调制度明显降低，与工程师主观评价感受一致。
本发明提供的方法，主要包括了多频段的滤波处理、滤波信号的包络线求取和包络线主要幅值调制度及调制频率的计算这三种技术手段，能快速锁定噪声产生调制的问题频段及其调制程度，为不同解决方案对调制问题的优化力度提供了准确有效的对比手段。问题频段和调制频率的确定有助于快速梳理噪声调制问题的机理，提供解决思路，规避了滤波回放手段的繁冗性。不同解决方案对调制问题的优化力度的有效对比，避免了问题频带三维频谱图对比的片面性和多个工程师主观判断的耗时性问题，提高了整车研发的效率。
根据本发明的另一方面，本发明还提供了一种量化汽车噪声问题的装置，包括：
第一获取模块，用于对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；
第二获取模块，用于对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；
第三获取模块，用于对多个第一噪声分别对应的包络线进行傅里叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；
第一确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度。
优选地，第二获取模块具体用于：
通过第一公式：
Figure BDA0001780358960000091
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，
Figure BDA0001780358960000092
为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，
Figure BDA0001780358960000093
为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
优选地，第三获取模块具体用于：
通过第二公式：
Figure BDA0001780358960000101
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure BDA0001780358960000102
其中，FFT[ej(tk：tk+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure BDA0001780358960000103
为函数向下取整；
通过第三公式：
Figure BDA0001780358960000104
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure BDA0001780358960000105
通过第四公式：
Figure BDA0001780358960000106
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
优选地，第一确定模块包括：
第一确定单元，用于确定主要幅值调制度最高的一个目标频段，将所述目标频段确定为原始有调制的噪声的问题频段。
优选地，所述装置还包括：
第二确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
优选地，根第二确定模块具体用于：
根据第五公式：
fmj(tk)=f，当
Figure BDA0001780358960000107
确定所述主要幅值调制度对应的调制频率fmj(tk)。
本发明提供的装置，主要包括了多频段的滤波处理、滤波信号的包络线求取和包络线主要幅值调制度及调制频率的计算这三种技术手段，能快速锁定噪声产生调制的问题频段及其调制程度，为不同解决方案对调制问题的优化力度提供了准确有效的对比手段。问题频段和调制频率的确定有助于快速梳理噪声调制问题的机理，提供解决思路，规避了滤波回放手段的繁冗性。不同解决方案对调制问题的优化力度的有效对比，避免了问题频带三维频谱图对比的片面性和多个工程师主观判断的耗时性问题，提高了整车研发的效率。

Claims (10)

1.一种量化汽车噪声调制问题的方法，其特征在于，包括：
对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；
对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；
对多个第一噪声分别对应的包络线进行傅里叶 变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；
根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度；
对多个第一噪声分别对应的包络线进行傅里叶变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度的步骤包括：
通过第二公式：
Figure FDA0003384982840000011
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure FDA0003384982840000012
其中，FFT[ej(tk：tkt+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure FDA0003384982840000013
为函数向下取整；
通过第三公式：
Figure FDA0003384982840000014
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure FDA0003384982840000015
通过第四公式：
Figure FDA0003384982840000016
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
2.根据权利要求1所述的方法，其特征在于，对多个第一噪声分别进行希尔伯特变换，获得所述第一噪声对应的包络线的步骤包括：
通过第一公式：
Figure FDA0003384982840000021
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，
Figure FDA0003384982840000022
为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，
Figure FDA0003384982840000023
为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
3.根据权利要求1所述的方法，其特征在于，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段的步骤包括：
确定主要幅值调制度最高的一个目标频段，将所述目标频段确定为原始有调制的噪声的问题频段。
4.根据权利要求1所述的方法，其特征在于，所述方法还包括：
根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
5.根据权利要求4所述的方法，其特征在于，根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率的步骤包括：
根据第五公式：
fmj(tk)＝f，当
Figure FDA0003384982840000024
确定所述主要幅值调制度对应的调制频率fmj(tk)。
6.一种量化汽车噪声问题的装置，其特征在于，包括：
第一获取模块，用于对原始有调制的噪声进行多频段带通滤波处理，获得带通滤波后的多个第一噪声；
第二获取模块，用于对多个第一噪声分别进行希尔伯特变换，获得多个第一噪声分别对应的包络线；
第三获取模块，用于对多个第一噪声分别对应的包络线进行傅里叶 变换，获得每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度；
第一确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定原始有调制的噪声的问题频段及其调制程度；
第三获取模块具体用于：
通过第二公式：
Figure FDA0003384982840000031
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻的快速傅里叶变换结果
Figure FDA0003384982840000032
其中，FFT[ej(tk：tk+Δt)]为快速傅里叶变换，Δt为第一噪声对应的包络线的时间间隔，t0为第一噪声对应的包络线的起始时刻，tend为第一噪声对应的包络线的终止时刻，Np为所述第一噪声对应的包络线中的时间间隔的数目，
Figure FDA0003384982840000033
为函数向下取整；
通过第三公式：
Figure FDA0003384982840000034
获得第j个频段带通滤波后的第一噪声对应的包络线在tk时刻且频率为f的幅值调制度
Figure FDA0003384982840000035
通过第四公式：
Figure FDA0003384982840000036
获得第j个频段带通滤波后的包络线在tk时刻的主要幅值调制度。
7.根据权利要求6所述的装置，其特征在于，第二获取模块具体用于：
通过第一公式：
Figure FDA0003384982840000037
获得第j个频段带通滤波后的第一噪声对应的包络线ej(t)，其中，xj(t)为第j个频段带通滤波后的第一噪声，Hilbert(xj(t))为对第j个频段带通滤波后的第一噪声进行希尔伯特变换，A0j为所述第一噪声对应包络线的常数项，
Figure FDA0003384982840000038
为所述第一噪声对应的包络线中第i个余弦信号，Aij为所述余弦信号的幅值，fij为所述余弦信号的频率，
Figure FDA0003384982840000039
为所述余弦信号的相位，N为ej(t)中余弦信号的数目。
8.根据权利要求6所述的装置，其特征在于，第一确定模块包括：
第一确定单元，用于确定主要幅值调制度最高的一个目标频段，将所述目标频段确定为原始有调制的噪声的问题频段。
9.根据权利要求6所述的装置，其特征在于，所述装置还包括：
第二确定模块，用于根据每一所述第一噪声对应的包络线在不同时刻的主要幅值调制度，确定所述主要幅值调制度对应的调制频率，并根据所述调制频率及所述问题频段确定原始有调制的噪声的机理。
10.根据权利要求9所述的装置，其特征在于，根第二确定模块具体用于：
根据第五公式：
fmj(tk)＝f，当
Figure FDA0003384982840000041
确定所述主要幅值调制度对应的调制频率fmj(tk)。
CN201810989074.0A 2018-08-28 2018-08-28 一种量化汽车噪声调制问题的方法及装置 Active CN109194306B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810989074.0A CN109194306B (zh) 2018-08-28 2018-08-28 一种量化汽车噪声调制问题的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810989074.0A CN109194306B (zh) 2018-08-28 2018-08-28 一种量化汽车噪声调制问题的方法及装置

Publications (2)

Publication Number Publication Date
CN109194306A CN109194306A (zh) 2019-01-11
CN109194306B true CN109194306B (zh) 2022-04-08

Family

ID=64916418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810989074.0A Active CN109194306B (zh) 2018-08-28 2018-08-28 一种量化汽车噪声调制问题的方法及装置

Country Status (1)

Country Link
CN (1) CN109194306B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112304632A (zh) * 2020-07-08 2021-02-02 重庆长安汽车股份有限公司 一种描述人耳感知的瞬态调制评价方法
CN113595528B (zh) * 2021-07-30 2023-08-04 重庆长安汽车股份有限公司 一种基于幂函数的自适应变步长lms滤波器及其实现方法
CN113758713B (zh) * 2021-08-09 2023-06-23 重庆长安汽车股份有限公司 粗糙声频段自适应识别方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4784343B2 (ja) * 2006-03-03 2011-10-05 トヨタ自動車株式会社 車両騒音音質制御装置
US8428957B2 (en) * 2007-08-24 2013-04-23 Qualcomm Incorporated Spectral noise shaping in audio coding based on spectral dynamics in frequency sub-bands
DE102008021360A1 (de) * 2008-04-29 2009-11-05 Siemens Aktiengesellschaft Verfahren und Vorrichtung zum Erkennen eines Lagerschadens
CN102608419B (zh) * 2012-03-31 2014-07-16 上海无线电设备研究所 具有噪声抑制性能的自适应瞬时频率测量方法
CN104568444B (zh) * 2015-01-28 2017-02-22 北京邮电大学 变转速火车滚动轴承故障特征频率提取方法

Also Published As

Publication number Publication date
CN109194306A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109194306B (zh) 一种量化汽车噪声调制问题的方法及装置
CN106596123B (zh) 设备故障诊断的方法、装置及系统
Cerrato Automotive sound quality–powertrain, road and wind noise
CN106840378B (zh) 汽车啸叫噪声的评价方法及系统
DE102016208048A1 (de) Automobile, diagnosesysteme und verfahren zur erzeugung von diagnosedaten für automobile
CN103362598B (zh) 一种汽车噪声源分离方法
CN102475554B (zh) 一种利用声品质指导车内声学包装的方法
CN106052849A (zh) 一种汽车车内非平稳异常噪声源识别方法
CN102288285A (zh) 一种单通道振动信号的盲源分离方法
CN105181126A (zh) 加速工况车内轰鸣音传递路径的检测方法
CN202522318U (zh) 一种汽车振动噪声检测系统
CN101163867A (zh) 诊断排气后处理元件的存在的方法和该方法用于车辆车载诊断的应用
CN103443426A (zh) 用于诊断内燃机的增压系统的方法
CN111595590A (zh) 怠速噪声优化方法、装置、设备及存储介质
CN111933105B (zh) 一种汽车噪声控制装置及其控制方法
CN115077690A (zh) 内燃机周期性脉动噪声评价方法
CN111189646A (zh) 车辆nvh自诊断方法、装置、车辆以及控制器和介质
Ajovalasit et al. Analysis of variations in diesel engine idle vibration
Zheng et al. Research on mechanism diagnosis of an idling abnormal noise of automobile engine
CN110021288A (zh) 车辆主动降噪控制方法及其系统
CN111551371A (zh) 车内噪声传递路径测试方法和车内噪声传递路径测试系统
CN112326267B (zh) 确定加速粗糙声影响结果方法及系统
CN112304632A (zh) 一种描述人耳感知的瞬态调制评价方法
CN114689167A (zh) 一种确定车内间歇性低频声问题所在位置的方法及装置
CN113008572B (zh) 用于评价n类汽车车内噪声的响度面积图生成系统和方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant