CN111564151B - A narrow-band active noise reduction optimization system for in-vehicle engine order noise - Google Patents

A narrow-band active noise reduction optimization system for in-vehicle engine order noise Download PDF

Info

Publication number
CN111564151B
CN111564151B CN202010401423.XA CN202010401423A CN111564151B CN 111564151 B CN111564151 B CN 111564151B CN 202010401423 A CN202010401423 A CN 202010401423A CN 111564151 B CN111564151 B CN 111564151B
Authority
CN
China
Prior art keywords
signal
internal reference
noise reduction
reference signal
filter
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.)
Expired - Fee Related
Application number
CN202010401423.XA
Other languages
Chinese (zh)
Other versions
CN111564151A (en
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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN202010401423.XA priority Critical patent/CN111564151B/en
Publication of CN111564151A publication Critical patent/CN111564151A/en
Application granted granted Critical
Publication of CN111564151B publication Critical patent/CN111564151B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1785Methods, e.g. algorithms; Devices
    • G10K11/17853Methods, e.g. algorithms; Devices of the filter
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1781Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions
    • G10K11/17813Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms
    • G10K11/17819Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms between the output signals and the reference signals, e.g. to prevent howling
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1785Methods, e.g. algorithms; Devices
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/10Applications
    • G10K2210/128Vehicles
    • G10K2210/1282Automobiles
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3028Filtering, e.g. Kalman filters or special analogue or digital filters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

The invention discloses a narrow-band active noise reduction optimization system for the order noise of an engine in a vehicle, which is used for smoothing an obtained rotating speed signal through an exponential smoothing formula aiming at the condition that the rotating speed of the engine fluctuates in the actual running process of the vehicle, constructing a more stable internal reference signal by using the smoothed rotating speed signal, obtaining an output signal of the narrow-band active noise reduction system according to the internal reference signal, and then updating the weight coefficient of a filter so as to continuously update the output signal. The system can solve the problem that the noise reduction performance of the system is reduced due to the fact that the internal reference signal presents obvious unsteady state characteristics caused by the fluctuation of the rotating speed of the engine in the actual situation, and has the potential of popularization and application in the field of active noise reduction in the vehicle.

Description

一种车内发动机阶次噪声的窄带主动降噪优化系统A narrow-band active noise reduction optimization system for in-vehicle engine order noise

技术领域technical field

本发明涉及车内主动噪声控制技术领域,更具体的是,本发明涉及一种车内发动机阶次噪声的窄带主动降噪优化系统。The invention relates to the technical field of in-vehicle active noise control, and more particularly, the invention relates to a narrow-band active noise reduction optimization system for in-vehicle engine order noise.

背景技术Background technique

传统上,噪声的被动控制主要通过控制噪声源、切断传播途径、保护接收者三种方法进行降噪,但这些方法只对车内波长较短的中、高频噪声能有效降噪,而对于波长较长的低频噪声,要想取得良好的降噪效果,则需要更加厚重的材料即更加高昂的成本。而噪声主动控制(Active Noise Control,ANC)技术是利用声波干涉的原理,通过在场景中合理布置麦克风及次级声源的方式来达到降噪目的的一种方法,它能够在基本不增加目标对象重量的条件下,有效地抑制其噪声水平,进而改善整体环境的声品质,是一种高效、高性价比的噪声控制方法。Traditionally, the passive control of noise mainly reduces noise by three methods: controlling the noise source, cutting off the propagation path, and protecting the receiver. For low-frequency noise with longer wavelengths, in order to achieve a good noise reduction effect, heavier materials are required, that is, higher costs. Active Noise Control (ANC) technology uses the principle of sound wave interference to achieve noise reduction by rationally arranging microphones and secondary sound sources in the scene. It is an efficient and cost-effective noise control method to effectively suppress the noise level of the object under the condition of the weight of the object, thereby improving the sound quality of the overall environment.

当前针对车内ANC系统的研究主要聚焦于发动机阶次噪声消除(Engine OrderCancellation,EOC)与路噪主动控制(Road Noise Cancellation,RNC)两项技术。其中,发动机阶次噪声通常是通过基于滤波-x最小均方(Filter-x Least Mean square,FxLMS)算法的窄带主动降噪(Narrowband ANC,NANC)系统来实现对阶次噪声的选择性控制,利用转速传感器获得转速信号,然后通过阶次噪声的频率与转速之间的关系为窄带主动降噪系统构造内部参考信号。然而,实际情况中发动机的工作状态并非绝对平稳,即便车辆运行在相对平稳的匀速工况下,发动机的转速仍在一定范围波动,这将导致所构造的内部参考信号具有明显的非稳态特性,使得算法的收敛性能和降噪效果下降。而指数平滑法(Exponential Smoothing,ES)恰好主要应用于平滑或预测时间序列,其原理为任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。若能将该方法应用于降低发动机阶次噪声的窄带主动降噪系统中对发动机转速信号进行平滑处理,将有助于构造更加平稳的内部参考信号,进而提升系统的性能、推进实际应用。The current research on in-vehicle ANC systems mainly focuses on two technologies, Engine Order Cancellation (EOC) and Road Noise Cancellation (RNC). Among them, the engine order noise is usually controlled by a narrowband active noise reduction (Narrowband ANC, NANC) system based on the Filter-x Least Mean Square (FxLMS) algorithm to achieve selective control of the order noise. The rotational speed signal is obtained by the rotational speed sensor, and then the internal reference signal is constructed for the narrow-band active noise reduction system through the relationship between the frequency of the order noise and the rotational speed. However, in practice, the working state of the engine is not absolutely stable. Even if the vehicle runs at a relatively stable and constant speed, the engine speed still fluctuates within a certain range, which will cause the constructed internal reference signal to have obvious non-steady-state characteristics. , which reduces the convergence performance and noise reduction effect of the algorithm. Exponential smoothing (ES) happens to be mainly used for smoothing or forecasting time series. The principle is that the exponential smoothing value of any period is the weighted average of the actual observation value of the current period and the exponential smoothing value of the previous period. If the method can be applied to smooth the engine speed signal in the narrow-band active noise reduction system that reduces engine order noise, it will help to construct a more stable internal reference signal, thereby improving the performance of the system and promoting practical applications.

发明内容SUMMARY OF THE INVENTION

本发明的目的是设计开发了一种车内发动机阶次噪声的窄带主动降噪优化系统,针对发动机转速存在波动的情况,通过指数平滑公式对获取的转速信号进行平滑,利用平滑后的转速信号来构造更加平稳的内部参考信号,根据内部参考信号得到窄带主动降噪系统的输出信号,而后更新滤波器的权值系数以对输出信号进行更新,从而改善窄带主动降噪算法的收敛性能和降噪性能。The purpose of the present invention is to design and develop a narrow-band active noise reduction optimization system for in-vehicle engine order noise. In view of the fluctuation of the engine speed, the obtained speed signal is smoothed by the exponential smoothing formula, and the smoothed speed signal is used to smooth the obtained speed signal. To construct a more stable internal reference signal, obtain the output signal of the narrowband active noise reduction system according to the internal reference signal, and then update the weight coefficient of the filter to update the output signal, thereby improving the convergence performance and reduction of the narrowband active noise reduction algorithm. noise performance.

本发明提供的技术方案为:The technical scheme provided by the present invention is:

一种车内发动机阶次噪声的窄带主动降噪优化系统,包括如下步骤:A narrow-band active noise reduction optimization system for in-vehicle engine order noise, comprising the following steps:

步骤一、获取发动机转速信号;Step 1. Obtain the engine speed signal;

步骤二、将所述发动机转速信号转化为平滑转速信号:Step 2: Convert the engine speed signal into a smooth speed signal:

R(t)=λ·R(t-1)+(1-λ)·r(t);R(t)=λ·R(t-1)+(1-λ)·r(t);

式中,R(t)为平滑转速信号,t为时间指数,且t=1、2、3…n,λ为平滑指数,且λ的取值范围为[0.9,1),r(t)为转速信号;In the formula, R(t) is the smooth speed signal, t is the time index, and t=1, 2, 3...n, λ is the smoothing index, and the value range of λ is [0.9, 1), r(t) is the speed signal;

步骤三、根据平滑转速信号获得发动机阶次噪声的频率、第一内部参考信号和第二内部参考信号,由此得到窄带主动降噪优化系统的扬声器输出信号:Step 3: Obtain the frequency of the engine order noise, the first internal reference signal and the second internal reference signal according to the smoothed rotational speed signal, thereby obtaining the speaker output signal of the narrowband active noise reduction optimization system:

Figure BDA0002489617940000021
Figure BDA0002489617940000021

式中,yN(t)为第t时刻的输出信号,i为目标阶次数,且i=1、2、3…q,q是目标窄带成分的角频率总数,

Figure BDA0002489617940000022
为第t时刻的滤波器第一权值系数,
Figure BDA0002489617940000023
为第t时刻的滤波器第二权值系数,xai(t)为第一内部参考信号,xbi(t)为第二内部参考信号。In the formula, y N (t) is the output signal at time t, i is the target order, and i=1, 2, 3...q, q is the total number of angular frequencies of the target narrowband components,
Figure BDA0002489617940000022
is the first weight coefficient of the filter at time t,
Figure BDA0002489617940000023
is the second weight coefficient of the filter at time t, x ai (t) is the first internal reference signal, and x bi (t) is the second internal reference signal.

优选的是,所述发动机阶次噪声的频率满足:Preferably, the frequency of the engine order noise satisfies:

ωi=2πR(t)ηi/60;ω i =2πR(t)η i /60;

式中,ωi为系统参考噪声信号中第i个目标阶次成分对应的角频率,ηi为第i个阶次成分对应的谐波数。In the formula, ω i is the angular frequency corresponding to the i-th target order component in the system reference noise signal, and η i is the harmonic number corresponding to the i-th order component.

优选的是,所述第一内部参考信号满足:Preferably, the first internal reference signal satisfies:

xai(t)=cos(ωit);x ai (t)=cos(ω i t);

式中,xai(t)为第一内部参考信号,ωi为系统参考噪声信号中第i个目标阶次成分对应的角频率,t为时间指数,且t=1、2、3…n。In the formula, x ai (t) is the first internal reference signal, ω i is the angular frequency corresponding to the i-th target order component in the system reference noise signal, t is the time index, and t=1, 2, 3...n .

优选的是,所述第二内部参考信号满足:Preferably, the second internal reference signal satisfies:

xbi(t)=sin(ωit);x bi (t)=sin(ω i t);

式中,xbi(t)为第二内部参考信号。In the formula, x bi (t) is the second internal reference signal.

优选的是,所述滤波器第一权值系数能够自适应更新为:Preferably, the first weight coefficient of the filter can be adaptively updated as:

Figure BDA0002489617940000031
Figure BDA0002489617940000031

式中,

Figure BDA0002489617940000032
为第t时刻的滤波器第一权值系数,
Figure BDA0002489617940000033
为第t+1时刻的滤波器第一权值系数,μN为窄带主动降噪算法的滤波器的权值更新步长,
Figure BDA0002489617940000034
为第一内部参考信号经估计的次级通路滤波后得到的滤波信号,eN(t)为残余噪声。In the formula,
Figure BDA0002489617940000032
is the first weight coefficient of the filter at time t,
Figure BDA0002489617940000033
is the first weight coefficient of the filter at time t+1, μ N is the weight update step size of the filter of the narrowband active noise reduction algorithm,
Figure BDA0002489617940000034
is the filtered signal obtained after the first internal reference signal is filtered by the estimated secondary path, and e N (t) is the residual noise.

优选的是,所述滤波器第二权值系数能够自适应更新为:Preferably, the second weight coefficient of the filter can be adaptively updated as:

Figure BDA0002489617940000035
Figure BDA0002489617940000035

式中,

Figure BDA0002489617940000036
为第t时刻的滤波器第二权值系数,
Figure BDA0002489617940000037
为第t+1次时刻的滤波器第二权值系数,
Figure BDA0002489617940000038
为第二内部参考信号经估计的次级通路滤波后得到的滤波信号。In the formula,
Figure BDA0002489617940000036
is the second weight coefficient of the filter at time t,
Figure BDA0002489617940000037
is the second weight coefficient of the filter at time t+1,
Figure BDA0002489617940000038
is a filtered signal obtained by filtering the estimated secondary path for the second internal reference signal.

优选的是,所述第一内部参考信号经估计的次级通路滤波后得到的滤波信号满足:Preferably, the filtered signal obtained after the first internal reference signal is filtered by the estimated secondary path satisfies:

Figure BDA0002489617940000039
Figure BDA0002489617940000039

式中,

Figure BDA00024896179400000310
为次级通路传递函数的估计的脉冲响应。In the formula,
Figure BDA00024896179400000310
is the estimated impulse response of the transfer function of the secondary path.

Figure BDA00024896179400000311
Figure BDA00024896179400000311

式中,

Figure BDA00024896179400000312
为次级通路传递函数的估计的脉冲响应。In the formula,
Figure BDA00024896179400000312
is the estimated impulse response of the transfer function of the secondary path.

优选的是,所述残余噪声满足:Preferably, the residual noise satisfies:

eN(t)=dN(t)-y′N(t);e N (t)=d N (t)-y′ N (t);

式中,dN(t)为初级期望信号,y′N(t)为第t时刻的输出信号经过次级通路的次级抵消信号。In the formula, d N (t) is the primary desired signal, and y′ N (t) is the secondary cancellation signal of the output signal at the t-th time passing through the secondary path.

优选的是,还包括:Preferably, it also includes:

指数平滑模块,其用于对采集的转速信号进行指数平滑处理,将得到的结果用来构造第一内部参考信号和第二内部参考信号;an exponential smoothing module, which is used to perform exponential smoothing processing on the collected rotational speed signal, and use the obtained result to construct the first internal reference signal and the second internal reference signal;

二权值更新模块,其用于实时更新滤波器的第一权值系数和第二权值系数;Two weight update modules, which are used to update the first weight coefficient and the second weight coefficient of the filter in real time;

预测滤波器模块,其用于接收所述二权值更新模块的结果和计算输出信号;a prediction filter module, which is used for receiving the result of the two-weight updating module and calculating the output signal;

误差合成模块,其用于对初级期望信号和次级抵消信号的相反数进行求和,并且所述误差合成模块能够将所得结果传输至所述二权值更新模块。An error synthesis module for summing the inverses of the primary desired signal and the secondary cancellation signal, and the error synthesis module can transmit the obtained result to the two-weight updating module.

本发明所述的有益效果:The beneficial effects of the present invention:

本发明设计开发的一种车内发动机阶次噪声的窄带主动降噪优化系统,与传统的用于车内的窄带主动降噪系统相比,其优势在于仅通过一个计算量很小的平滑公式,便可以显著地改善算法在发动机转速存在波动的情况下的收敛性能和降噪性能,使得车内的降噪效果更优越。Compared with the traditional narrow-band active noise reduction system used in the vehicle, the narrow-band active noise reduction optimization system designed and developed by the invention has the advantage of only passing a smooth formula with a small amount of calculation. , it can significantly improve the convergence performance and noise reduction performance of the algorithm when the engine speed fluctuates, making the noise reduction effect in the car more superior.

附图说明Description of drawings

图1为本发明所述的一种针对车内发动机阶次噪声的窄带主动降噪优化系统的原理框图。FIG. 1 is a schematic block diagram of a narrow-band active noise reduction optimization system for in-vehicle engine order noise according to the present invention.

图2为本发明所述的对比试验中车速为30km/h的匀速工况下采集的车内噪声参考信号的时域图。FIG. 2 is a time-domain diagram of an in-vehicle noise reference signal collected under a constant speed operating condition with a vehicle speed of 30 km/h in the comparative test according to the present invention.

图3为本发明所述的对比试验中车速为30km/h的匀速工况下采集的同步转速信号图。FIG. 3 is a graph of a synchronous rotational speed signal collected under a constant speed operating condition with a vehicle speed of 30 km/h in the comparative test according to the present invention.

图4为本发明所述的对比试验中车速为30km/h的匀速工况下采集的车内噪声参考信号的频域图。FIG. 4 is a frequency domain diagram of an in-vehicle noise reference signal collected under a constant speed operating condition with a vehicle speed of 30 km/h in the comparative test according to the present invention.

图5为本发明所述的对比试验中车速为60km/h的匀速工况下采集的车内噪声参考信号的时域图。FIG. 5 is a time-domain diagram of an in-vehicle noise reference signal collected under a constant speed operating condition with a vehicle speed of 60 km/h in the comparative test according to the present invention.

图6为本发明所述的对比试验中车速为60km/h的匀速工况下采集的同步转速信号图。FIG. 6 is a graph of a synchronous rotational speed signal collected under a constant speed operating condition with a vehicle speed of 60 km/h in the comparative test according to the present invention.

图7为本发明所述的对比试验中车速为60km/h的匀速工况下采集的车内噪声参考信号的频域图。FIG. 7 is a frequency domain diagram of an in-vehicle noise reference signal collected under a constant speed operating condition with a vehicle speed of 60 km/h in the comparative test according to the present invention.

图8为本发明所述的对比试验中采用的初级通道与次级通道的幅频响应曲线图。FIG. 8 is a graph showing the amplitude-frequency response of the primary channel and the secondary channel used in the comparative test of the present invention.

图9为本发明所述的对比试验中采用的初级通道与次级通道的相频响应曲线图。FIG. 9 is a phase-frequency response curve diagram of the primary channel and the secondary channel used in the comparative test of the present invention.

图10为本发明所述的对比试验中优化窄带主动降噪系统的指数平滑模块对车速为30km/h的匀速工况下采集的转速信号的平滑效果图。FIG. 10 is a graph showing the smoothing effect of the speed signal collected by the exponential smoothing module of the optimized narrowband active noise reduction system in the comparative test according to the present invention on the speed signal collected under a constant speed operating condition with a vehicle speed of 30km/h.

图11为本发明所述的对比试验中对车速为30km/h的匀速工况下采集的车内噪声降噪的频域降噪效果图。FIG. 11 is a frequency domain noise reduction effect diagram of in-vehicle noise reduction collected under a constant speed operating condition with a vehicle speed of 30 km/h in the comparative test according to the present invention.

图12为本发明所述的对比试验中优化窄带主动降噪系统的指数平滑模块对车速为60km/h的匀速工况下采集的转速信号的平滑效果图。12 is a graph showing the smoothing effect of the speed signal collected by the exponential smoothing module of the optimized narrowband active noise reduction system in the comparative test of the present invention on the speed signal collected under a constant speed operating condition with a vehicle speed of 60 km/h.

图13为本发明所述的对比试验中对车速为60km/h的匀速工况下采集的车内噪声降噪的频域降噪效果图。FIG. 13 is a frequency domain noise reduction effect diagram of in-vehicle noise reduction collected under a constant speed operating condition with a vehicle speed of 60km/h in the comparative test according to the present invention.

图14为本发明所述的对比试验中对车速为30km/h的匀速工况下采集的车内噪声降噪的误差信号的均方值(MSE)曲线图。FIG. 14 is a mean square value (MSE) curve diagram of the error signal of in-vehicle noise reduction collected under a constant speed operating condition with a vehicle speed of 30km/h in the comparative test according to the present invention.

图15为本发明所述的对比试验中对车速为60km/h的匀速工况下采集的车内噪声降噪的误差信号的均方值(MSE)曲线图。FIG. 15 is a mean square value (MSE) curve diagram of the error signal of the noise reduction in the vehicle collected under the constant speed operating condition of the vehicle speed of 60km/h in the comparative test according to the present invention.

具体实施方式Detailed ways

下面结合对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。The present invention will be further described in detail below, so that those skilled in the art can implement it with reference to the description.

如图1所示,本发明提供的一种车内发动机阶次噪声的窄带主动降噪优化系统包括:指数平滑模块、二权值更新模块、预测滤波器模块、误差合成模块。其中,所述指数平滑模块用于对采集的转速信号进行指数平滑处理,将得到的结果用来构造第一内部参考信号和第二内部参考信号;所述二权值更新模块,采用自适应迭代公式来实时更新预测滤波器系数,并将所得结果传输到预测滤波器模块;所述预测滤波器模块用于计算输出信号,作为优选的是,所述预测滤波器采用有限脉冲响应滤波器(FIR滤波器);所述误差合成模块对初级期望信号和次级抵消信号的相反数进行求和,并将所得结果传输到二权值更新模块。As shown in FIG. 1 , a narrow-band active noise reduction optimization system for in-vehicle engine order noise provided by the present invention includes: an exponential smoothing module, a two-weight update module, a prediction filter module, and an error synthesis module. Wherein, the exponential smoothing module is used to perform exponential smoothing processing on the collected rotational speed signal, and the obtained results are used to construct the first internal reference signal and the second internal reference signal; the two weight updating module adopts adaptive iteration formula to update the prediction filter coefficients in real time, and transmit the result to the prediction filter module; the prediction filter module is used to calculate the output signal, and preferably, the prediction filter adopts a finite impulse response filter (FIR filter); the error synthesis module sums the inverses of the primary expected signal and the secondary cancellation signal, and transmits the result to the two-weight updating module.

如图1所示,主动噪声控制系统的基本原理是声波的叠加相消,有主动降噪系统的算法迭代计算,通过扬声器发出一列与目标噪声或初级期望噪声幅值相同相位相反的信号。其中,扬声器与目标降噪区域之间的通道的传递函数为次级通路传递函数,它代表通道对声音幅值和相位的影响;

Figure BDA0002489617940000061
是次级通路传递函数的估计,S(z)为次级通路传递函数,并且设定
Figure BDA0002489617940000062
次级通路传递函数代表现实存在的一种对声音信号的影响,它是客观存在,不需要获取,而次级通路传递函数的估计可以通过次级通路建模辨识来获得;dN(t)是初级期望信号,也叫做初级噪声,即车内目标降噪区域的噪声,yN(t)是降噪系统的输出信号,y′N(t)是输出信号经过次级通路的次级抵消信号,eN(t)是初级噪声与次级抵消信号在目标降噪区域叠加后的残余噪声信号,也叫做误差信号,用于反馈至二权值更新模块进行权值更新。As shown in Figure 1, the basic principle of the active noise control system is the superposition and cancellation of sound waves. There is an algorithm iterative calculation of the active noise reduction system, and a series of signals with the same phase and opposite phase as the target noise or the primary expected noise are emitted through the speaker. Among them, the transfer function of the channel between the speaker and the target noise reduction area is the secondary path transfer function, which represents the influence of the channel on the sound amplitude and phase;
Figure BDA0002489617940000061
is the estimate of the transfer function of the secondary path, S(z) is the transfer function of the secondary path, and set
Figure BDA0002489617940000062
The transfer function of the secondary path represents a real influence on the sound signal, it exists objectively and does not need to be acquired, and the estimation of the transfer function of the secondary path can be obtained through the identification of the secondary path modeling; d N (t) is the primary desired signal, also called primary noise, that is, the noise in the target noise reduction area in the vehicle, y N (t) is the output signal of the noise reduction system, and y′ N (t) is the secondary cancellation of the output signal through the secondary path The signal, e N (t) is the residual noise signal after the primary noise and the secondary cancellation signal are superimposed in the target noise reduction area, also called the error signal, which is used for feedback to the two weight update module for weight update.

在窄带主动降噪子系统中,xai(t)为子系统第一内部参考信号,xbi(t)为子系统第二内部参考信号,

Figure BDA0002489617940000063
为第一内部参考信号经估计的次级通路滤波后得到的滤波信号,
Figure BDA0002489617940000064
为第二内部参考信号经估计的次级通路滤波后得到的滤波信号。In the narrowband active noise reduction subsystem, x ai (t) is the first internal reference signal of the subsystem, x bi (t) is the second internal reference signal of the subsystem,
Figure BDA0002489617940000063
is the filtered signal obtained by filtering the estimated secondary path for the first internal reference signal,
Figure BDA0002489617940000064
is a filtered signal obtained by filtering the estimated secondary path for the second internal reference signal.

在针对车内发动机阶次噪声进行主动降噪时,降噪过程如下:When performing active noise reduction for the engine order noise in the car, the noise reduction process is as follows:

首先,对由转速传感器提供的同步信号,即为转速信号进行指数平滑处理得到平滑转速信号:First, perform exponential smoothing on the synchronization signal provided by the speed sensor, that is, the speed signal, to obtain a smooth speed signal:

R(t)=λ·R(t-1)+(1-λ)·r(t);R(t)=λ·R(t-1)+(1-λ)·r(t);

式中,R(t)为平滑转速信号,t为时间指数,且t=1、2、3…n,λ为平滑指数,它是一个常数,且λ的取值范围为[0.9,1),r(t)为转速信号。In the formula, R(t) is the smoothing speed signal, t is the time index, and t=1, 2, 3...n, λ is the smoothing index, which is a constant, and the value range of λ is [0.9, 1) , r(t) is the speed signal.

由此计算目标窄带成分的频率,即发动机阶次噪声的频率为:From this, the frequency of the target narrowband component, that is, the frequency of the engine order noise, is calculated as:

ωi=2πR(n)i/60;ω i =2πR(n)i/60;

式中,ωi为系统参考噪声信号中第i个目标阶次成分对应的角频率,ηi为第i个阶次成分对应的谐波数,i为目标阶次数,且i=1、2、3…q,q是目标窄带成分的角频率总数。In the formula, ω i is the angular frequency corresponding to the i-th target order component in the system reference noise signal, η i is the harmonic number corresponding to the i-th order component, i is the target order, and i=1, 2 , 3...q, q is the total number of angular frequencies of the target narrowband component.

窄带主动降噪子系统据此合成第一内部参考信号和第二内部参考信号,结果如下:The narrowband active noise reduction subsystem synthesizes the first internal reference signal and the second internal reference signal accordingly, and the results are as follows:

xai(t)=cos(ωit);x ai (t)=cos(ω i t);

xbi(t)=sin(ωit);x bi (t)=sin(ω i t);

式中,xai(t)为第一内部参考信号,xbi(t)为第二内部参考信号。In the formula, x ai (t) is the first internal reference signal, and x bi (t) is the second internal reference signal.

第一内部参考信号和第二内部参考信号分别经估计的次级通路滤波后,得到滤波信号,结果如下:After the first internal reference signal and the second internal reference signal are filtered by the estimated secondary paths, respectively, the filtered signals are obtained, and the results are as follows:

Figure BDA0002489617940000071
Figure BDA0002489617940000071

Figure BDA0002489617940000072
Figure BDA0002489617940000072

式中,

Figure BDA0002489617940000073
为第一内部参考信号经估计的次级通路滤波后得到的滤波信号,
Figure BDA0002489617940000074
为第二内部参考信号经估计的次级通路滤波后得到的滤波信号,
Figure BDA0002489617940000075
为次级通路传递函数的估计的脉冲响应,s(t)为次级通道传递函数的脉冲响应。In the formula,
Figure BDA0002489617940000073
is the filtered signal obtained by filtering the estimated secondary path for the first internal reference signal,
Figure BDA0002489617940000074
is the filtered signal obtained by filtering the estimated secondary path for the second internal reference signal,
Figure BDA0002489617940000075
is the estimated impulse response of the transfer function of the secondary channel, and s(t) is the impulse response of the transfer function of the secondary channel.

预测滤波器模块将所得滤波器第一权值系数和滤波器第二权值系数分别与第一内部参考信号和第二内部参考信号进行卷积并求和作为子系统输出信号,表示如下:The prediction filter module convolves the obtained first weight coefficient of the filter and the second weight coefficient of the filter with the first internal reference signal and the second internal reference signal respectively and sums them up as the subsystem output signal, which is expressed as follows:

Figure BDA0002489617940000076
Figure BDA0002489617940000076

式中,yN(t)为t时刻的输出信号,i为目标阶次数,且i=1、2、3…q,q是目标窄带成分的角频率总数,

Figure BDA0002489617940000077
为t时刻的滤波器第一权值系数,
Figure BDA0002489617940000078
为t时刻的滤波器第二权值系数,xai(t)为第一内部参考信号,xbi(t)为第二内部参考信号。In the formula, y N (t) is the output signal at time t, i is the target order, and i=1, 2, 3...q, q is the total number of angular frequencies of the target narrowband components,
Figure BDA0002489617940000077
is the first weight coefficient of the filter at time t,
Figure BDA0002489617940000078
is the second weight coefficient of the filter at time t, x ai (t) is the first internal reference signal, and x bi (t) is the second internal reference signal.

误差合成模块将初级期望信号和次级抵消信号的相反数进行求和,得到误差信号,结果如下:The error synthesis module sums the inverses of the primary desired signal and the secondary cancellation signal to obtain an error signal. The results are as follows:

eN(t)=dN(t)-y′N(t);e N (t)=d N (t)-y′ N (t);

式中,dN(t)为初级期望信号,y′N(t)为第t次输出信号经过次级通路的次级抵消信号。In the formula, d N (t) is the primary desired signal, and y' N (t) is the secondary cancellation signal of the t-th output signal passing through the secondary path.

将误差信号和滤波信号带入二权值更新模块中的自适应迭代公式,进行权值更新。The error signal and the filtered signal are brought into the adaptive iterative formula in the two-weight updating module to update the weights.

二权值更新模块中的自适应迭代公式基于最速下降原理,表示为:The adaptive iterative formula in the two-weight update module is based on the principle of steepest descent and is expressed as:

Figure BDA0002489617940000079
Figure BDA0002489617940000079

Figure BDA0002489617940000081
Figure BDA0002489617940000081

其中,

Figure BDA0002489617940000082
为第t时刻预测滤波器的第一权值系数,
Figure BDA0002489617940000083
为第t+1时刻预测滤波器的第一权值系数;
Figure BDA0002489617940000084
为第t时刻预测滤波器的第二权值系数,
Figure BDA0002489617940000085
为第t+1时刻预测滤波器的第二权值系数,μN为窄带主动降噪子系统的滤波器的权值更新步长,Ja(t)为自适应迭代公式的第一代价函数,Jb(t)为自适应迭代公式的第二代价函数,
Figure BDA0002489617940000086
为第一代价函数的梯度,
Figure BDA0002489617940000087
为第二代价函数的梯度。in,
Figure BDA0002489617940000082
is the first weight coefficient of the prediction filter at time t,
Figure BDA0002489617940000083
is the first weight coefficient of the prediction filter at time t+1;
Figure BDA0002489617940000084
is the second weight coefficient of the prediction filter at time t,
Figure BDA0002489617940000085
is the second weight coefficient of the prediction filter at time t+1, μ N is the weight update step size of the filter of the narrowband active noise reduction subsystem, and J a (t) is the first cost function of the adaptive iterative formula , J b (t) is the second cost function of the adaptive iteration formula,
Figure BDA0002489617940000086
is the gradient of the first cost function,
Figure BDA0002489617940000087
is the gradient of the second cost function.

并且第一代价函数和第二代价函数满足:And the first cost function and the second cost function satisfy:

Ja(t)=Jb(t)=E[eN 2(t)];J a (t)=J b (t)=E[e N 2 (t)];

第一代价函数的梯度和第二代价函数的梯度表示如下:The gradient of the first cost function and the gradient of the second cost function are expressed as follows:

Figure BDA0002489617940000088
Figure BDA0002489617940000088

Figure BDA0002489617940000089
Figure BDA0002489617940000089

其中,eN(t)为误差信号,即误差合成子模块的输出信号。Among them, e N (t) is the error signal, that is, the output signal of the error synthesis sub-module.

可得出预测滤波器的第一权值系数和第二权值系数的更新结果如下:The update results of the first weight coefficient and the second weight coefficient of the prediction filter can be obtained as follows:

Figure BDA00024896179400000810
Figure BDA00024896179400000810

Figure BDA00024896179400000811
Figure BDA00024896179400000811

不断重复上述过程,即可实现目标场景中发动机阶次噪声的有效控制。By repeating the above process continuously, the effective control of the engine order noise in the target scene can be achieved.

为检验本发明所提系统对车内发动机阶次噪声的降噪性能,将本发明中提供的针对车内发动机阶次噪声的优化窄带主动降噪系统与现有技术中的传统窄带主动降噪系统进行对比试验如下:In order to test the noise reduction performance of the system proposed in the present invention for the order noise of the in-vehicle engine, the optimized narrow-band active noise reduction system for the in-vehicle engine order noise provided in the present invention is compared with the traditional narrow-band active noise reduction in the prior art. The system is compared and tested as follows:

如图2-4所示,车速为30km/h的匀速工况下车内噪声以及同步的转速信号的测试结果,图5-7所示,车速为60km/h的匀速工况下车内噪声以及同步的转速信号的测试结果,其中,图2与图5车速为30km/h和60km/h两工况下采集的车内噪声信号,图3与图6为两工况下采集的对应的同步转速信号,图4与图7为两工况下采集的车内噪声信号的频谱图。As shown in Figure 2-4, the test results of in-vehicle noise and synchronous speed signal under the constant speed condition of the vehicle speed of 30km/h, as shown in Figure 5-7, under the constant speed condition of the vehicle speed of 60km/h And the test results of the synchronous speed signal, among which, Figure 2 and Figure 5 are the in-vehicle noise signals collected under the two operating conditions of 30km/h and 60km/h, and Figure 3 and Figure 6 are the corresponding data collected under the two operating conditions. Synchronous speed signal, Figure 4 and Figure 7 are the spectrum diagrams of the in-vehicle noise signals collected under two operating conditions.

实验中,P(z)为初级通路传递函数,S(z)为次级通路传递函数,两者均以阶数为64阶的FIR滤波器表示,两个通路的频率响应如附图8与图9所示,其中,图8为两个通路的幅频响应曲线,图9为两个通路的相频响应曲线。In the experiment, P(z) is the transfer function of the primary path, S(z) is the transfer function of the secondary path, both of which are represented by an FIR filter with an order of 64. The frequency responses of the two paths are shown in Figures 8 and 8. As shown in FIG. 9 , wherein, FIG. 8 is the amplitude-frequency response curve of the two channels, and FIG. 9 is the phase-frequency response curve of the two channels.

本发明所提供的系统中的二权值更新模块和预测滤波器模块中所涉及的滤波器均为阶数为1的FIR滤波器。其中,二权值更新模块中的步长参数在车速为30km/h的工况下为μN=9.6×10-4,在60km/h的工况下对应的步长为μN=1.6×10-3,系统降噪的目标窄带成分频率对应发动机阶次噪声的一阶和二阶,即η1=1,η2=2。The filters involved in the two-weight update module and the prediction filter module in the system provided by the present invention are all FIR filters with an order of 1. Among them, the step size parameter in the two-weight update module is μ N =9.6×10 −4 under the working condition of the vehicle speed of 30km/h, and the corresponding step size under the working condition of 60km/h is μ N =1.6× 10 -3 , the target narrow-band component frequency of system noise reduction corresponds to the first and second orders of engine order noise, ie η 1 =1, η 2 =2.

由图10-13所示,可以从频域上看出,传统窄带主动降噪系统对于实车中发动机阶次噪声的控制并不是十分理想,在目标窄带成分频率附近的频带甚至呈现出声压级升高的现象,降噪性能不佳。相比之下,本发明所提的系统则对发动机阶次噪声的抑制具有更优越的降噪效果,并且不会引起附近频带的噪声声压级的增大。As shown in Figure 10-13, it can be seen from the frequency domain that the traditional narrowband active noise reduction system is not very ideal for the control of engine order noise in real vehicles, and the frequency band near the target narrowband component frequency even shows sound pressure. The phenomenon of level rise, the noise reduction performance is not good. In contrast, the system proposed in the present invention has a superior noise reduction effect on the suppression of engine order noise, and does not cause an increase in the noise sound pressure level of the nearby frequency band.

根据计算,传统主动窄带主动降噪系统在30km/h和60km/h两种匀速工况下分别能实现0.62dB和0.86dB的总降噪量,而本发明所提系统在两工况下分别能实现1.77dB和1.40dB的总降噪量,具有更突出的降噪性能。According to the calculation, the traditional active narrowband active noise reduction system can achieve a total noise reduction of 0.62dB and 0.86dB respectively under two uniform speed conditions of 30km/h and 60km/h, while the system proposed in the present invention can achieve a total noise reduction of 0.62dB and 0.86dB respectively under the two working conditions. It can achieve a total noise reduction of 1.77dB and 1.40dB, with more prominent noise reduction performance.

此外,如图14、图15所示,为两系统误差信号的均方值曲线,可以看出,两系统在收敛性能中的稳态误差方面的表现与在其在频域上的降噪能力基本一致,并且本发明所提系统呈现出更好的稳定性。In addition, as shown in Figure 14 and Figure 15, which are the mean square value curves of the error signals of the two systems, it can be seen that the performance of the two systems in terms of steady-state error in the convergence performance and their noise reduction ability in the frequency domain Basically the same, and the system proposed in the present invention exhibits better stability.

本发明设计开发的车内发动机阶次噪声的窄带主动降噪优化系统,与传统的用于车内的窄带主动降噪系统相比,其优势在于仅通过一个计算量很小的平滑公式,便可以显著地改善算法在发动机转速存在波动的情况下的收敛性能和降噪性能,具体为,通过指数平滑公式对获取的转速信号进行平滑,利用平滑后的转速信号来构造更加平稳的内部参考信号,根据所述内部参考信号得到窄带主动降噪系统的输出信号,而后更新滤波器的权值系数以对输出信号进行更新,使得车内的降噪性能更优越。Compared with the traditional narrow-band active noise reduction system used in the vehicle, the narrow-band active noise reduction optimization system designed and developed by the present invention has the advantage that only a smooth formula with a small amount of calculation is used, and the It can significantly improve the convergence performance and noise reduction performance of the algorithm when the engine speed fluctuates. Specifically, the obtained speed signal is smoothed by the exponential smoothing formula, and the smoothed speed signal is used to construct a more stable internal reference signal. , obtain the output signal of the narrowband active noise reduction system according to the internal reference signal, and then update the weight coefficient of the filter to update the output signal, so that the noise reduction performance in the vehicle is better.

尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的实施例。Although the embodiment of the present invention has been disclosed as above, it is not limited to the application listed in the description and the embodiment, and it can be applied to various fields suitable for the present invention. For those skilled in the art, it can be easily Therefore, the invention is not limited to the specific details and embodiments shown and described herein without departing from the general concept defined by the appended claims and the scope of equivalents.

Claims (4)

1. A narrow-band active noise reduction optimization system for the order noise of an engine in a vehicle is characterized by comprising the following steps:
step one, obtaining an engine rotating speed signal;
step two, converting the engine rotating speed signal into a smooth rotating speed signal:
R(t)=λ·R(t-1)+(1-λ)·r(t);
wherein, r (t) is a smooth rotation speed signal, t is a time index, t is 1, 2, 3 … n, λ is a smooth index, λ has a value range of [0.9, 1 ], and r (t) is a rotation speed signal;
step three, obtaining the frequency of the engine order noise, the first internal reference signal and the second internal reference signal according to the smooth rotating speed signal:
ω i =2πR(t)η i /60;
x ai (t)=cos(ω i t);
x bi (t)=sin(ω i t);
in the formula, omega i For the angular frequency, eta, corresponding to the ith target order component in the system reference noise signal i Harmonic number, x, corresponding to the ith order component ai (t) is the first internal reference signal, x bi (t) is a second internal reference signal;
the loudspeaker output signal of the narrow-band active noise reduction optimization system is obtained by the following steps:
Figure FDA0003749869780000011
in the formula, y N (t) is the output signal at time t, i is the target order number, and i is 1, 2, 3 … q, q is the total number of angular frequencies of the target narrow-band component,
Figure FDA0003749869780000012
the first weight coefficients of the filter at the t-th time,
Figure FDA0003749869780000013
is the second weight coefficient, x, of the filter at the t-th moment ai (t) is the first internal reference signal, x bi (t) is a second internal reference signal;
wherein, the first weight coefficient of the filter can be adaptively updated as:
Figure FDA0003749869780000014
in the formula,
Figure FDA0003749869780000015
the first weight coefficients of the filter at the t-th time,
Figure FDA0003749869780000016
is the first weight coefficient, mu, of the filter at the t +1 th moment N Is a narrow bandThe weight update step size of the filter of the active noise reduction algorithm,
Figure FDA0003749869780000017
a filtered signal obtained by filtering the estimated secondary path of the first internal reference signal, e N (t) is residual noise;
the second weight coefficient of the filter can be adaptively updated as follows:
Figure FDA0003749869780000021
in the formula,
Figure FDA0003749869780000022
the second weight coefficients of the filter at the t-th time,
Figure FDA0003749869780000023
the second weight coefficients of the filter at the t +1 th time,
Figure FDA0003749869780000024
filtering the second internal reference signal by the estimated secondary path to obtain a filtered signal;
wherein, the narrowband of engine order noise in the car initiative noise reduction optimizing system includes:
an exponential smoothing module, configured to perform exponential smoothing on the collected engine speed signal, and use an obtained result to construct the first internal reference signal and the second internal reference signal;
the two weight updating modules are used for updating the first weight coefficient and the second weight coefficient of the filter in real time;
a prediction filter module for receiving the results of the two weight update modules and calculating output signals;
an error synthesis module for summing the inverses of the primary desired signal and the secondary cancellation signal, and the error synthesis module can transmit the result to the two weight update module.
2. The system for narrowband active noise reduction optimization of an in-vehicle engine order noise of claim 1, wherein the first internal reference signal is filtered by the estimated secondary path to obtain a filtered signal satisfying:
Figure FDA0003749869780000025
in the formula,
Figure FDA0003749869780000026
an estimated impulse response of the secondary path transfer function.
3. The system for narrowband active noise reduction optimization of an in-vehicle engine order noise of claim 2, wherein the second internal reference signal is filtered by the estimated secondary path to obtain a filtered signal satisfying:
Figure FDA0003749869780000027
in the formula,
Figure FDA0003749869780000028
an estimated impulse response of the secondary path transfer function.
4. The in-vehicle engine order noise narrowband active noise reduction optimization system of claim 3, wherein the residual noise satisfies:
e N (t)=d N (t)-y′ N (t);
in the formula (d) N (t) is the primary desired signal, y' N And (t) is a secondary cancellation signal of the output signal at the t-th time point through the secondary path.
CN202010401423.XA 2020-05-13 2020-05-13 A narrow-band active noise reduction optimization system for in-vehicle engine order noise Expired - Fee Related CN111564151B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010401423.XA CN111564151B (en) 2020-05-13 2020-05-13 A narrow-band active noise reduction optimization system for in-vehicle engine order noise

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010401423.XA CN111564151B (en) 2020-05-13 2020-05-13 A narrow-band active noise reduction optimization system for in-vehicle engine order noise

Publications (2)

Publication Number Publication Date
CN111564151A CN111564151A (en) 2020-08-21
CN111564151B true CN111564151B (en) 2022-09-23

Family

ID=72073433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010401423.XA Expired - Fee Related CN111564151B (en) 2020-05-13 2020-05-13 A narrow-band active noise reduction optimization system for in-vehicle engine order noise

Country Status (1)

Country Link
CN (1) CN111564151B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112328949B (en) * 2020-10-26 2024-02-27 中科上声(苏州)电子有限公司 Reference signal generation method and device for active noise reduction system of automobile engine
CN113593516B (en) * 2021-07-22 2024-04-02 中国船舶集团有限公司第七一一研究所 Active vibration and noise control method, system, storage medium and ship
CN114530140B (en) * 2022-01-19 2024-12-24 中科上声(苏州)电子有限公司 A vehicle-mounted noise reduction method, system and storage medium for narrowband noise
CN114677997B (en) * 2022-02-14 2024-08-06 中国第一汽车股份有限公司 Real vehicle active noise reduction method and system based on acceleration working condition

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08177452A (en) * 1994-12-20 1996-07-09 Fujitsu Ten Ltd Noise controller
CN106089361A (en) * 2016-06-30 2016-11-09 重庆长安汽车股份有限公司 A kind of car intrinsic motivation active noise reduction system and method
CN106382143A (en) * 2016-12-01 2017-02-08 吉林大学 Active noise reduction device and active noise reduction method based on engine speed
CN206299429U (en) * 2016-12-01 2017-07-04 吉林大学 A kind of active noise reducing device based on engine speed
CN107642426A (en) * 2017-08-31 2018-01-30 清华大学苏州汽车研究院(相城) A kind of automobile engine noise initiative control method and system
CN110335582A (en) * 2019-07-11 2019-10-15 吉林大学 An active noise reduction method suitable for active control of impulse noise
WO2019241657A1 (en) * 2018-06-14 2019-12-19 Harman International Industries, Incorporated Concurrent fxlms system with common reference and error signals

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08177452A (en) * 1994-12-20 1996-07-09 Fujitsu Ten Ltd Noise controller
CN106089361A (en) * 2016-06-30 2016-11-09 重庆长安汽车股份有限公司 A kind of car intrinsic motivation active noise reduction system and method
CN106382143A (en) * 2016-12-01 2017-02-08 吉林大学 Active noise reduction device and active noise reduction method based on engine speed
CN206299429U (en) * 2016-12-01 2017-07-04 吉林大学 A kind of active noise reducing device based on engine speed
CN107642426A (en) * 2017-08-31 2018-01-30 清华大学苏州汽车研究院(相城) A kind of automobile engine noise initiative control method and system
WO2019241657A1 (en) * 2018-06-14 2019-12-19 Harman International Industries, Incorporated Concurrent fxlms system with common reference and error signals
CN110335582A (en) * 2019-07-11 2019-10-15 吉林大学 An active noise reduction method suitable for active control of impulse noise

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Active noise control: a tutorial review";S. M. Kuo;《IEEE》;19991231;全文 *
"New Equalizing Scheme of Active Noise Equalization System in Automobile Cabin";Liang W;《IEEE International Conference on Multimedia & Expo》;20071231;全文 *
"自适应振动噪声主动控制若干关键问题研究";浦玉学;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20160715;全文 *

Also Published As

Publication number Publication date
CN111564151A (en) 2020-08-21

Similar Documents

Publication Publication Date Title
CN111564151B (en) A narrow-band active noise reduction optimization system for in-vehicle engine order noise
CN111402853B (en) A Wide and Narrow Band Hybrid Active Noise Reduction Algorithm Suitable for Vehicles
US10373600B2 (en) Active noise control system
CN105489225B (en) A Feedforward Narrowband Active Noise Control System with On-Line Identification of Secondary Channels
CN101819766B (en) Multi-channel active noise control method for abating noises
Kuo et al. Frequency-domain delayless active sound quality control algorithm
CN106531145B (en) Recurrence active noise control method based on M estimator
CN101802905A (en) Active silencer and method of controlling active silencer
GB2501325A (en) Non-adaptive controller for an ANC system, using coefficients determined from experimental data
CN113078884B (en) Add adaptive algorithm for nonlinear fitting
Durak et al. Adaptive fractional Fourier domain filtering
CN106935238A (en) A kind of active noise balance control method based on the control of optimal critical band amplitude gain
US20230252967A1 (en) Road noise cancellation shaping filters
CN113327570A (en) Narrow-band feedforward type active noise control system and method
CN104575512A (en) Nonlinear narrowband active noise control method based on Volterra filter
CN111627414A (en) Active denoising method and device and electronic equipment
CN109658947A (en) A kind of active noise controlling method of synchronous modeling and control
CN105913836A (en) Fixed-point real time noise reduction method based on DSP
CN109961773B (en) Active noise reduction method for rotary mechanical order noise
CN106782594A (en) One kind suppresses the unmatched method of frequency in arrowband active noise controlling
CN113299260A (en) Active noise reduction method based on EMFNL filter online modeling secondary channel
US20230169949A1 (en) System and method for providing frequency dependent dynamic leakage for a feed forward active noise cancellation (anc)
CN113096629B (en) Relative path virtual sensing method for single-channel feedback active noise control system
Chang et al. Enhanced offline secondary path modeling method for narrowband active noise control system
CN113851104A (en) Feedback active noise control system and method with on-line identification of secondary channel

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220923

CF01 Termination of patent right due to non-payment of annual fee