JP2007219262A - Active vibrating noise controlling apparatus - Google Patents

Active vibrating noise controlling apparatus Download PDF

Info

Publication number
JP2007219262A
JP2007219262A JP2006040881A JP2006040881A JP2007219262A JP 2007219262 A JP2007219262 A JP 2007219262A JP 2006040881 A JP2006040881 A JP 2006040881A JP 2006040881 A JP2006040881 A JP 2006040881A JP 2007219262 A JP2007219262 A JP 2007219262A
Authority
JP
Japan
Prior art keywords
noise
vibration
signal
delay time
model
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.)
Granted
Application number
JP2006040881A
Other languages
Japanese (ja)
Other versions
JP4735319B2 (en
Inventor
Mensler Michel
ミシエル メンスレ
Yoshiro Takamatsu
吉郎 高松
Kazuo Nakamura
一男 中村
Kazuto Sato
一人 佐藤
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.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor 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 Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP2006040881A priority Critical patent/JP4735319B2/en
Publication of JP2007219262A publication Critical patent/JP2007219262A/en
Application granted granted Critical
Publication of JP4735319B2 publication Critical patent/JP4735319B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide an active vibrating noise controlling apparatus, capable of reducing only noise caused by vibration of a casing which is a control object, without reducing in-vehicle music and voice of an attendant which is not a control object. <P>SOLUTION: The active vibrating noise controlling apparatus comprises: a case 5 in which vibration propagates from outside to inside as noise; a sensor 7 for detecting the vibration of the casing 5; an actuator 8 for generating control vibration on the casing 5; a signal arithmetic means 31 for generating a control instructing signal to the actuator 8 so that the noise in the case may be reduced; and a noise estimating means 32 for estimating the noise inside the casing by using the output signal of the sensor 7 and the control instruction signal. <P>COPYRIGHT: (C)2007,JPO&INPIT

Description

本発明は、外部から筐体内に侵入する振動又は騒音に制御振動を重ね合わせることで当該振動又は騒音を低減する能動振動騒音制御装置に関する。   The present invention relates to an active vibration noise control device that reduces vibration or noise by superimposing control vibration on vibration or noise that enters the casing from the outside.

車室内騒音を低減することを目的としたフィードバック制御を行う騒音制御装置が特許文献1に開示されている。ここでは、マイクを使って100Hz以下の周波数帯域のロードノイズを低減することを目的としている。低周波ロードノイズは車室内音場の音響モードが原因となっていると仮定をした上で、この音響モードの腹にマイクロフォンを設置しフィードバック制御を行う方式である。
特開2000−322066号公報
A noise control device that performs feedback control for the purpose of reducing vehicle interior noise is disclosed in Patent Document 1. Here, an object is to reduce road noise in a frequency band of 100 Hz or less using a microphone. The low frequency road noise is assumed to be caused by the acoustic mode of the vehicle interior sound field, and a feedback control is performed by installing a microphone on the belly of the acoustic mode.
JP 2000-322066 A

しかしながら、上記の騒音制御装置では、マイクロフォンを用いて車室内騒音を測定するため、制御する必要のない車室内音楽や乗員の声などもロードノイズと共に参照信号として用いてしまい、フィードバック制御により低減してしまうという問題点があった。   However, in the above-described noise control device, vehicle interior noise is measured using a microphone, and therefore interior music that is not required to be controlled, passenger voices, and the like are also used as reference signals together with road noise, and are reduced by feedback control. There was a problem that it was.

また、マイクロフォンを車室内に多数配置するため、マイクロフォン自体のコストやハーネスのコストや重量が増加するという問題点があった。   Further, since a large number of microphones are arranged in the vehicle interior, there is a problem that the cost of the microphone itself and the cost and weight of the harness increase.

さらに、ハーネスを車室内にはりめぐらせて実現していたため、ハーネス長が長くなりシステムの信頼性が悪化してしまうという問題点があった。   Furthermore, since the harness is realized by running around the vehicle interior, there is a problem that the harness length becomes long and the reliability of the system deteriorates.

本発明の特徴は、振動が外部から内部へ騒音として伝達する筐体と、この筐体の振動を検出するセンサと、筐体に制御振動を生成するアクチュエータと、筐体内騒音を低減するように、アクチュエータへの制御指令信号を生成する信号演算手段と、センサの出力信号及び制御指令信号を用いて筐体内騒音を推定する騒音推定手段とを備える能動振動騒音制御装置であることである。   A feature of the present invention is that a housing that transmits vibration from outside to inside as noise, a sensor that detects vibration of the housing, an actuator that generates control vibration in the housing, and noise in the housing are reduced. The active vibration noise control device includes a signal calculation unit that generates a control command signal to the actuator, and a noise estimation unit that estimates the noise in the housing using the output signal of the sensor and the control command signal.

本発明に係わる能動振動騒音制御装置によれば、センサの出力信号及び制御指令信号を用いて筐体内騒音を推定することにより、制御対象とならない車室内音楽や乗員の声などを低減させること無く、制御対象となる筐体の振動による騒音だけを低減させることができる。   According to the active vibration noise control device according to the present invention, by estimating the noise in the housing using the output signal of the sensor and the control command signal, it is possible to reduce vehicle interior music and occupant voice that are not controlled. Only the noise caused by the vibration of the casing to be controlled can be reduced.

以下図面を参照して、本発明の実施の形態を説明する。図面の記載において同一あるいは類似の部分には同一あるいは類似な符号を付している。   Embodiments of the present invention will be described below with reference to the drawings. In the description of the drawings, the same or similar parts are denoted by the same or similar reference numerals.

(第1の実施の形態)
図1を参照して、本発明の第1の実施の形態に係わる能動振動騒音制御装置の全体構成を説明する。第1の実施の形態に係わる能動振動騒音制御装置(AVC:Active Vibration Controller)は、振動が外部から内部へ騒音として伝達する筐体(車体5)と、外部から伝わる筐体の振動を検出するセンサ(加速度センサ7)と、筐体に制御振動を生成するアクチュエータ8と、筐体の内部に伝達する騒音(筐体内騒音)を低減するように、アクチュエータ8への制御指令信号を生成する信号演算手段(コントローラ31)と、加速度センサ7の出力信号及びアクチュエータ8への制御指令信号を用いて筐体内騒音を推定する騒音推定手段(騒音推定部32)とを備える。
(First embodiment)
With reference to FIG. 1, the overall configuration of an active vibration noise control apparatus according to a first embodiment of the present invention will be described. An active vibration noise control device (AVC: Active Vibration Controller) according to the first embodiment detects vibrations of a casing (vehicle body 5) that transmits vibration from outside to inside as noise and a casing that is transmitted from the outside. A sensor (acceleration sensor 7), an actuator 8 that generates control vibration in the housing, and a signal that generates a control command signal to the actuator 8 so as to reduce noise transmitted to the inside of the housing (in-housing noise). Calculation means (controller 31) and noise estimation means (noise estimation unit 32) for estimating the noise in the housing using the output signal of the acceleration sensor 7 and the control command signal to the actuator 8 are provided.

筐体は、外部から加わる力によって振動が生じ、この振動が内部へ騒音として伝達し得るものであればよい。例えば、飛行機、電車、自動車、エレベータなど人間の移動手段がこれに含まれる。特に、自動車は、凹凸が激しい路面にタイヤを接して走行するため、タイヤから伝わる振動の制御が重要な技術課題である。そこで実施の形態では、本発明に係わる能動振動騒音制御装置を自動車に適用した場合を例に取り説明する。つまり、ここでは筐体の一例として車体を挙げて説明する。   Any housing may be used as long as vibration is generated by a force applied from the outside and the vibration can be transmitted to the inside as noise. This includes, for example, human moving means such as airplanes, trains, cars, and elevators. In particular, since an automobile travels with a tire on a road surface with severe irregularities, control of vibration transmitted from the tire is an important technical issue. Therefore, in the embodiment, a case where the active vibration noise control device according to the present invention is applied to an automobile will be described as an example. That is, here, a vehicle body will be described as an example of a housing.

加速度センサ7は、自己に加わる加速度を測定する機能を備えるセンサであり、加速度センサ7を車体上に固定して加速度を測定することにより、測定される加速度の変化から車体の振動を検出することができる。加速度センサ7以外にも、外部から伝わる筐体の振動を検出するセンサとして、センサが固定された位置(座標)の変化から直接振動を検出するセンサであっても構わないし、また、タイヤから伝わる振動を検出するのであれば、タイヤ又はそれが接続された車軸の振動を検出するセンサであっても構わない。図1の例では、加速度センサ7は、車体5のうちボンネットやルーフパネルに比べてタイヤ4に近いフロアパネル6の下面に固定されている。   The acceleration sensor 7 is a sensor having a function of measuring acceleration applied to itself, and by detecting the acceleration by fixing the acceleration sensor 7 on the vehicle body, the vibration of the vehicle body is detected from the change in the measured acceleration. Can do. In addition to the acceleration sensor 7, the sensor for detecting the vibration of the casing transmitted from the outside may be a sensor that directly detects the vibration from a change in the position (coordinates) at which the sensor is fixed, or transmitted from the tire. As long as vibration is detected, it may be a sensor that detects vibration of a tire or an axle to which the tire is connected. In the example of FIG. 1, the acceleration sensor 7 is fixed to the lower surface of the floor panel 6 that is closer to the tire 4 than the hood or roof panel in the vehicle body 5.

アクチュエータ8は、車体5の中のフロアパネル6下に固定され、電気的な信号(制御指令信号)に受けて自らが歪むことでフロアパネル6に歪みを生成する。これにより、車体5に制御振動を生成して外乱による振動又は騒音を相殺する。アクチュエータ8として、結晶に電界を加えると電界に比例した歪みが生じるピエゾ効果(圧電効果)を利用したピエゾ電気アクチュエータ(piezo-electric actuator)を使用することができる。アクチュエータ8の両面に電界を印加することによりアクチュエータ8は歪む。   The actuator 8 is fixed under the floor panel 6 in the vehicle body 5, and generates distortion in the floor panel 6 by receiving distortion in response to an electrical signal (control command signal). Thereby, the control vibration is generated in the vehicle body 5 to cancel the vibration or noise caused by the disturbance. As the actuator 8, a piezo-electric actuator utilizing a piezo effect (piezoelectric effect) that generates a strain proportional to the electric field when an electric field is applied to the crystal can be used. The actuator 8 is distorted by applying an electric field to both surfaces of the actuator 8.

コントローラ31及び騒音推定部32は、コントローラユニット30に含まれている。コントローラユニット30は、加速度センサ7で得られた信号に基づいて車室内騒音を低減する制御指令値を演算する。コントローラユニット30は、加速度センサ7で得られた信号を増幅する信号増幅用アンプ33aと、加速度センサ7の出力信号とアクチュエータ8への制御指令信号に基づいて車室内騒音を推定する騒音推定部32と、騒音推定部32で推定された推定騒音に基づいて車室内騒音を低減するアクチュエータ8への制御指令値を演算するコントローラ31とを備える。   The controller 31 and the noise estimation unit 32 are included in the controller unit 30. The controller unit 30 calculates a control command value for reducing vehicle interior noise based on the signal obtained by the acceleration sensor 7. The controller unit 30 includes a signal amplification amplifier 33 a that amplifies the signal obtained by the acceleration sensor 7, and a noise estimation unit 32 that estimates vehicle interior noise based on the output signal of the acceleration sensor 7 and the control command signal to the actuator 8. And a controller 31 that calculates a control command value for the actuator 8 that reduces vehicle interior noise based on the estimated noise estimated by the noise estimation unit 32.

車室内騒音の原因は代表的なものとして、エンジン10の振動に起因する騒音、走行時に路面の凹凸の影響がタイヤから進入することに起因するロードノイズと呼ばれる騒音、走行時に空気の気流によって発生する風切音、などがある。また、タイヤから侵入するロードノイズは、タイヤ及びフロアパネルの振動を介して、これが車室内に伝達する。つまり、フロアパネル6の振動が車体5の内部(車室内)に騒音として伝達する。本実施形態において、騒音推定部32はフロアパネル6の振動に起因する車室内騒音を主に推定することを目的とするため、フロアパネル6下に加速度センサ7を配置して車室内騒音を推定する。したがって、上記例の車室内騒音のうち、フロアパネル振動に起因するロードノイズを中心に制御を行う。   Typical causes of vehicle interior noise are noise caused by the vibration of the engine 10, noise called road noise caused by road surface unevenness entering from tires during driving, and airflow during driving Wind noise, etc. In addition, road noise entering from the tire is transmitted to the vehicle interior via vibration of the tire and the floor panel. That is, the vibration of the floor panel 6 is transmitted as noise to the inside of the vehicle body 5 (vehicle interior). In the present embodiment, the noise estimation unit 32 mainly estimates the vehicle interior noise caused by the vibration of the floor panel 6. Therefore, the acceleration sensor 7 is disposed under the floor panel 6 to estimate the vehicle interior noise. To do. Therefore, control is performed mainly on road noise caused by floor panel vibration among the vehicle interior noise in the above example.

図2に示すように、路面の凹凸の影響により車両の4つのタイヤ4が振動し、タイヤ4から車体に進入したロードノイズの主成分となる振動は、まず車軸及びびサスペンションが取り付けられている点からメンバと呼ばれる剛性の高い梁状の部材に進入する。その後、メンバによって囲まれたフロアパネル6と呼ばれる比較的剛性の低い板状の部材に伝達し、フロアパネル6が振動する。さらに、フロアパネル6の膜振動により車室内9の空気振動が引き起こされ、車室内9に共振現象を起こす。このために、車室内9でロードノイズとしての騒音が聞こえる。フロアパネル6の他にルーフパネルや窓ガラスが振動することにより騒音が発生するが、サスペンション取付け部から主に進入するロードノイズの大部分は、フロアパネル6の振動が原因となっていることが分かっている。このようなタイヤ4と車室内9を結ぶロードノイズの侵入経路上に加速度センサ7が配置されている。   As shown in FIG. 2, the four tires 4 of the vehicle vibrate due to the road surface unevenness, and the main component of road noise entering the vehicle body from the tire 4 is first attached to the axle and suspension. From a point, it enters a highly rigid beam-like member called a member. Thereafter, the transmission is transmitted to a plate-like member called a floor panel 6 surrounded by the members and having relatively low rigidity, and the floor panel 6 vibrates. Furthermore, the air vibration in the vehicle interior 9 is caused by the membrane vibration of the floor panel 6, causing a resonance phenomenon in the vehicle interior 9. For this reason, noise as road noise is heard in the passenger compartment 9. Although noise is generated when the roof panel and the window glass vibrate in addition to the floor panel 6, most of the road noise that mainly enters from the suspension mounting portion is caused by the vibration of the floor panel 6. I know it. The acceleration sensor 7 is disposed on the road noise intrusion path connecting the tire 4 and the vehicle interior 9.

一般的なタイヤ4からのロードノイズの進入経路に加え、フロアパネル6に貼付されたアクチュエータ8による制御振動の伝達経路が存在する。アクチュエータ8が発生した制御振動はロードノイズと同様にフロアパネル6を振動させ、そのフロアパネル6の膜振動により車室内9の空気振動が起こり、車室内9への制御騒音となる。一方、アクチュエータ8から発生する振動は加速度センサ7にも伝播する。したがって、加速度センサ7で計測されるフロアパネル6の振動は、タイヤ4から車体に進入する振動のうち加速度センサ7の取付け位置を通る成分と、ピエゾアクチュエータ8が発生する制御振動のうち同位置を通る成分とが重ね合わされた振動となる。   In addition to a road noise approach path from a general tire 4, there is a control vibration transmission path by an actuator 8 attached to the floor panel 6. The control vibration generated by the actuator 8 vibrates the floor panel 6 in the same way as road noise, and air vibration in the vehicle interior 9 occurs due to membrane vibration of the floor panel 6, and becomes control noise to the vehicle interior 9. On the other hand, vibration generated from the actuator 8 also propagates to the acceleration sensor 7. Therefore, the vibration of the floor panel 6 measured by the acceleration sensor 7 has the same position among the vibration that enters the vehicle body from the tire 4 and the component that passes through the mounting position of the acceleration sensor 7 and the control vibration that the piezo actuator 8 generates. It becomes a vibration that is superposed on the passing component.

図2中の範囲9は車室内空間の一例として、乗員耳位置を含む球体状の領域を指す。この空間での騒音は、タイヤ4から車体に入力された騒音と、アクチュエータ8が発生する制御振動とが重ね合わされた車室内騒音になる。コントローラ31は車室内領域9での騒音が低減されるような制御指令信号をアクチュエータ8に供給して制御振動を発生させる。路面の凹凸に起因するフロアパネル6の振動とアクチュエータ8が生成する制御振動の両者が互いに打ち消しあうように重ね合わさって車室内領域9での騒音は低減される。   A range 9 in FIG. 2 indicates a spherical region including an occupant ear position as an example of a vehicle interior space. The noise in this space becomes vehicle interior noise in which noise input from the tire 4 to the vehicle body and control vibration generated by the actuator 8 are superimposed. The controller 31 supplies a control command signal for reducing noise in the vehicle interior region 9 to the actuator 8 to generate control vibration. Noise in the vehicle interior region 9 is reduced by superimposing the vibration of the floor panel 6 caused by the road surface unevenness and the control vibration generated by the actuator 8 so as to cancel each other.

図1のコントローラ31は、騒音推定部32で推定された音圧(車室内騒音)を入力として、車室内空間9での騒音を低減するような制御指令信号を演算する機能を備える固定係数のFIRもしくはIIRフィルタである。例えば、コントローラ31は以下の手順に従って設計をすればよい。   The controller 31 in FIG. 1 has a fixed coefficient having a function of calculating a control command signal that reduces the noise in the vehicle interior space 9 by using the sound pressure (vehicle interior noise) estimated by the noise estimation unit 32 as an input. FIR or IIR filter. For example, the controller 31 may be designed according to the following procedure.

まず、「アクチュエータ8への制御指令信号から車室内空間9での騒音までの伝達特性」のモデル化を行う。モデル化の手順は以下のとおりである。まず、アクチュエータ8とマイクロフォンをそれぞれフロアパネル6と車室内空間9に設置する。アクチュエータ8へあらかじめ生成しておいたM系列信号を入力し、入力信号とマイクロフォン信号を計測する。得られた入出力信号に対し、例えば Lennart Ljung,“System Identification.Second edition.Thory for the user”,Prentice Hall PTR,1999に記載のシステム同定手法を用いることにより、例えば状態空間モデル
x(n+1) = Ax(n) + Bu(n)
y(n) = Cx(n) + Du(n)
の形式でモデル化がなされる。ここで、x(t)はモデルの状態変数、u(t)はピエゾアクチュエータへの入力信号、y(t)は車室内空間9での音圧を表し、マトリクスA、B、C、Dはそれぞれモデルのパラメータである。また、nはサンプル時刻を表す。
First, modeling of “transfer characteristics from a control command signal to the actuator 8 to noise in the vehicle interior space 9” is performed. The modeling procedure is as follows. First, the actuator 8 and the microphone are installed in the floor panel 6 and the vehicle interior space 9, respectively. An M-sequence signal generated in advance is input to the actuator 8, and the input signal and the microphone signal are measured. By using the system identification method described in, for example, Lennart Ljung, “System Identification. Second edition. Thorough for the user”, Parent Hall PTR, 1999, for example, a state space model.
x (n + 1) = Ax (n) + Bu (n)
y (n) = Cx (n) + Du (n)
Is modeled in the form Here, x (t) is a model state variable, u (t) is an input signal to the piezoelectric actuator, y (t) is a sound pressure in the vehicle interior space 9, and matrices A, B, C, and D are Each is a model parameter. N represents the sample time.

次に、図3に示すブロック線図を用いて制御系設計の方法を説明する。ここで、コントローラ34は信号Bが入力され、アクチュエータ8への制御指令信号を出力する。制御指令信号は伝達関数モデル35に入力される。伝達関数モデル35は、上の手順でモデル化した「アクチュエータ8の制御指令信号から車室内空間9での騒音までの伝達特性」を示すモデルである。伝達関数モデル35から出力される信号と信号A(外乱)とが重ねあわされて信号B(車室内騒音)となる。図3の信号Aから信号Bまでの伝達特性が所望の周波数帯域において小さくなるようにコントローラ34を設計すればよい。設計法としては、例えば、PID制御、H2制御、H∞制御などを用いることができる。また、
細江、荒木、「制御系設計、H∞制御とその応用」、朝倉書店、1994
須田、「PID制御」、朝倉書店、1992
に記載されている設計法を用いることで、数1の形式のI次動的コントローラKi,i=1,・・・,Iを設計することが出来る。

Figure 2007219262
Next, a control system design method will be described with reference to the block diagram shown in FIG. Here, the controller 34 receives the signal B and outputs a control command signal to the actuator 8. The control command signal is input to the transfer function model 35. The transfer function model 35 is a model indicating “transfer characteristics from the control command signal of the actuator 8 to the noise in the vehicle interior space 9” modeled in the above procedure. The signal output from the transfer function model 35 and the signal A (disturbance) are overlapped to form a signal B (vehicle interior noise). The controller 34 may be designed so that the transfer characteristic from the signal A to the signal B in FIG. 3 becomes small in a desired frequency band. As a design method, for example, PID control, H2 control, H∞ control, or the like can be used. Also,
Hosoe, Araki, “Control System Design, H∞ Control and its Applications”, Asakura Shoten, 1994
Suda, “PID Control”, Asakura Shoten, 1992
By using the design method described in the above, it is possible to design the I-order dynamic controller Ki, i = 1,.
Figure 2007219262

図4に示すように、騒音推定部32は、加速度センサ7で得られたフロアパネルの振動(出力信号)と、アクチュエータ8への制御指令信号とが入力される。騒音推定部32は、加速度センサ7が取り付けられた位置における振動から車室内空間9の騒音までの伝達特性を表す第1の伝達モデル14と、アクチュエータ8への制御指令信号から車室内空間9の騒音までの伝達特性を表す第2の伝達モデル15と、アクチュエータ8への制御指令信号から加速度センサ7の振動までの伝達特性を表す第3の伝達モデル16とをデータとして記憶している。第1乃至第3の伝達モデル14、15、16はそれぞれ図1の伝達経路1、2、3の数学モデルとなっている。なお、第1乃至第3の伝達モデル14、15、16は上記のモデル化手順に従うことで得ることが出来る。   As shown in FIG. 4, the noise estimation unit 32 receives the floor panel vibration (output signal) obtained by the acceleration sensor 7 and the control command signal to the actuator 8. The noise estimation unit 32 includes a first transmission model 14 representing transmission characteristics from vibration at a position where the acceleration sensor 7 is attached to noise in the vehicle interior space 9, and a control command signal to the actuator 8 in the vehicle interior space 9. The second transfer model 15 representing the transfer characteristic up to noise and the third transfer model 16 representing the transfer characteristic from the control command signal to the actuator 8 to the vibration of the acceleration sensor 7 are stored as data. The first to third transmission models 14, 15, and 16 are mathematical models of the transmission paths 1, 2, and 3 in FIG. The first to third transfer models 14, 15, and 16 can be obtained by following the above modeling procedure.

具体的に説明すると、第1の伝達モデル14は、加速度センサ7が取り付けられた位置におけるフロアパネル振動が車室内空間9まで伝達するときの伝達特性を表す。第2の伝達モデル15は、アクチュエータ8への入力信号が、アクチュエータ8へ入力されてアクチュエータが生成する制御振動が車室内空間9まで伝達するときの伝達特性を表す。第3の伝達モデル16は、アクチュエータ8への入力信号が、アクチュエータ8へ入力されてアクチュエータが生成する制御振動が加速度センサ7が設置された位置まで伝達するときの伝達特性を表す。   More specifically, the first transmission model 14 represents a transmission characteristic when the floor panel vibration at the position where the acceleration sensor 7 is attached is transmitted to the vehicle interior space 9. The second transfer model 15 represents a transfer characteristic when an input signal to the actuator 8 is input to the actuator 8 and a control vibration generated by the actuator is transmitted to the vehicle interior space 9. The third transmission model 16 represents a transmission characteristic when an input signal to the actuator 8 is input to the actuator 8 and a control vibration generated by the actuator is transmitted to a position where the acceleration sensor 7 is installed.

次に、騒音推定部32によって車室内騒音を推定するアルゴリズムの例を説明する。   Next, an example of an algorithm for estimating the vehicle interior noise by the noise estimation unit 32 will be described.

先ず、第1のアルゴリズムにおいて、騒音推定部32は、加速度センサ7の出力信号を「そのまま」第1の伝達モデル14へ入力し、制御指令信号を第2の伝達モデル15へ入力する。そして、両モデル14、15から得られる出力値の和を、車室内騒音(推定音圧:SPL)として出力する。   First, in the first algorithm, the noise estimation unit 32 inputs the output signal of the acceleration sensor 7 “as is” to the first transmission model 14 and inputs the control command signal to the second transmission model 15. Then, the sum of output values obtained from both models 14 and 15 is output as vehicle interior noise (estimated sound pressure: SPL).

第2のアルゴリズムにおいて、騒音推定部32は、加速度センサ7の出力信号のうちアクチュエータ8が発生する振動の成分を減算し、これを第1の伝達モデル14へ入力する。そして、第1の伝達モデル14から得られる出力と、制御指令信号を第2の伝達モデル15へ入力して得られる出力との和を車室内騒音として出力する。このように、加速度センサ7の出力信号のうちアクチュエータ8が発生する振動の成分を減算することにより、加速度センサ7の出力値のうちアクチュエータ8による制御振動を除去して路面の凹凸によるロードノイズだけを抽出することができ、ロードノイズ(外乱)の測定精度が向上する。   In the second algorithm, the noise estimation unit 32 subtracts the vibration component generated by the actuator 8 from the output signal of the acceleration sensor 7 and inputs this to the first transmission model 14. Then, the sum of the output obtained from the first transmission model 14 and the output obtained by inputting the control command signal to the second transmission model 15 is output as vehicle interior noise. In this way, by subtracting the vibration component generated by the actuator 8 from the output signal of the acceleration sensor 7, the control vibration by the actuator 8 is removed from the output value of the acceleration sensor 7, and only road noise due to road surface unevenness is obtained. And the measurement accuracy of road noise (disturbance) is improved.

第2のアルゴリズムを具体的に説明すると、図4に示すように、騒音推定部32は、加速度センサ7の出力信号から第3の伝達モデル16の出力を減算した信号を第1の伝達モデル14へ入力して得られる出力と、制御指令信号を第2の伝達モデル15へ入力して得られる出力との和を車室内騒音として出力する。   The second algorithm will be described specifically. As shown in FIG. 4, the noise estimation unit 32 uses a signal obtained by subtracting the output of the third transfer model 16 from the output signal of the acceleration sensor 7 as the first transfer model 14. The sum of the output obtained by inputting to and the output obtained by inputting the control command signal to the second transfer model 15 is output as vehicle interior noise.

図5のフローチャートを参照して、図4に示す騒音推定部32による騒音推定手順(第2のアルゴリズム)の詳細を説明する。   Details of the noise estimation procedure (second algorithm) by the noise estimation unit 32 shown in FIG. 4 will be described with reference to the flowchart of FIG.

(イ)まず、ステップS101で加速度センサ7の出力信号を騒音推定部32へ入力する。これを信号Cとする。   (A) First, the output signal of the acceleration sensor 7 is input to the noise estimation unit 32 in step S101. This is signal C.

(ロ)ステップS102に進み、アクチュエータ8への制御指令信号を騒音推定部32へ入力する。これを信号Dとする。   (B) Proceeding to step S102, the control command signal for the actuator 8 is input to the noise estimation unit 32. This is signal D.

(ハ)ステップS103に進み、信号Dに第3の伝達モデル16を乗算する。これを信号Eとする。   (C) Proceeding to step S103, the signal D is multiplied by the third transfer model 16. This is signal E.

(ニ)ステップS104に進み、信号Cから信号Eを減算する。これを信号Fとする。ステップS105に進み、信号Fに第1の伝達モデル14を乗算する。これを信号Gとする。   (D) Proceeding to step S104, the signal E is subtracted from the signal C. This is signal F. In step S105, the signal F is multiplied by the first transfer model 14. This is a signal G.

(ホ)ステップS106に進み、信号Dに第2の伝達モデル15を乗算する。これを信号Hとする。   (E) Proceeding to step S106, the signal D is multiplied by the second transfer model 15. This is signal H.

(へ)ステップS107に進み、信号Gと信号Hを加算する。これを信号Jとする。最後に、ステップS108にて信号Jを車室内空間位置9での推定音圧として出力する。   (F) Proceed to step S107, and the signal G and the signal H are added. This is signal J. Finally, in step S108, the signal J is output as the estimated sound pressure at the vehicle interior space position 9.

次に、上記した騒音推定部32を用いて車室内騒音を推定することによる効果を説明する。図6は実際に測定した運転席右耳元位置での音圧と、騒音推定装置を用いて推定した同位置での音圧を比較した図である。なお、運転席右耳元位置は車室内位置9の一例である。ここで、実線42は実際に車室内位置9でマイクロフォンにより測定した騒音(実測音圧)であり、点線43は騒音推定部32による推定音圧である。図6に示すように、推定音圧は実測音圧とほぼ同等の値を示している。よって、マイクロフォンを用いることなく、騒音推定部32によって車室内騒音を精度良く推定できることがわかる。   Next, the effect by estimating vehicle interior noise using the noise estimation part 32 mentioned above is demonstrated. FIG. 6 is a diagram comparing the sound pressure at the driver's right ear position actually measured with the sound pressure at the same position estimated using the noise estimation device. The driver's seat right ear position is an example of the vehicle interior position 9. Here, the solid line 42 is the noise (actually measured sound pressure) actually measured by the microphone at the vehicle interior position 9, and the dotted line 43 is the estimated sound pressure by the noise estimating unit 32. As shown in FIG. 6, the estimated sound pressure has a value substantially equal to the actually measured sound pressure. Therefore, it can be seen that the vehicle interior noise can be accurately estimated by the noise estimation unit 32 without using a microphone.

図7は、図4の第3の伝達モデル16の有無による実測音圧と推定音圧とのコヒーレンスの変化を示す。実線46は、図4において第3の伝達モデル16がある場合(第2のアルゴリズム)、つまり、加速度センサ7の出力信号のうちアクチュエータ8が発生する振動の成分を減算した信号を第1の伝達モデル14へ入力する場合を示す。点線47は、第3の伝達モデル16がない場合(第1のアルゴリズム)、つまり、加速度センサ7の出力信号をそのまま第1の伝達モデルへ14入力する場合を示す。図7は、騒音推定部32により推定される推定音圧とマイクロフォンにより測定した実測音圧との間のコヒーレンス(Coherence)関数を計算した結果を示している。ここで、コヒーレンス関数とは、2つの信号間の相関を表し、推定精度を表す尺度であり、信号x(t)とy(t)の間のコヒーレンス関数は数2により定義される。

Figure 2007219262
FIG. 7 shows a change in coherence between the actually measured sound pressure and the estimated sound pressure due to the presence or absence of the third transfer model 16 of FIG. The solid line 46 shows the case where the third transmission model 16 is present in FIG. 4 (second algorithm), that is, the first transmission signal obtained by subtracting the vibration component generated by the actuator 8 from the output signal of the acceleration sensor 7. The case where it inputs to the model 14 is shown. A dotted line 47 indicates a case where there is no third transmission model 16 (first algorithm), that is, a case where the output signal of the acceleration sensor 7 is input to the first transmission model as it is. FIG. 7 shows a result of calculating a coherence function between the estimated sound pressure estimated by the noise estimation unit 32 and the actually measured sound pressure measured by the microphone. Here, the coherence function is a measure representing the correlation between two signals and representing the estimation accuracy, and the coherence function between the signals x (t) and y (t) is defined by Equation 2.
Figure 2007219262

ここで、Cxyは信号x(t)と信号y(t)との間のコヒーレンス、Sxyは信号x(t)と信号y(t)との間のクロススペクトル密度、Sx、Syはそれぞれ信号x(t)、信号y(t)のスペクトル密度を表す。実測音圧と推定音圧とが完全に一致するときコヒーレンス関数の値は1となり、2つの信号が無相関であるときコヒーレンス関数は0を値としてとる。 Here, C xy is the coherence between the signal x (t) and the signal y (t), S xy is the cross spectral density between the signal x (t) and the signal y (t), S x , S y Represents the spectral densities of the signal x (t) and the signal y (t), respectively. The value of the coherence function is 1 when the actually measured sound pressure and the estimated sound pressure completely match, and the value of the coherence function is 0 when the two signals are uncorrelated.

図7を見ると分かるように、第3の伝達モデル16を入れた場合は第3の伝達モデル16を入れない場合に比べ、コヒーレンス関数が高くなっている。よって、第3の伝達モデル16を入れることによって、より精度の高い音圧推定が実現されることが分かる。   As can be seen from FIG. 7, the coherence function is higher when the third transfer model 16 is inserted than when the third transfer model 16 is not inserted. Therefore, it can be seen that by inserting the third transmission model 16, more accurate sound pressure estimation is realized.

以上説明したように、本発明の第1の実施の形態によれば、以下の作用効果が得られる。   As described above, according to the first embodiment of the present invention, the following operational effects can be obtained.

騒音推定部32を設けることにより、車室内9にマイクロフォンを設置しないので、乗員の声や音楽を低減することなく、車外から進入する振動に起因した車室内騒音のみを低減する制御を行うことが出来る。また、マイクを用いて制御する場合に比べ、センサとアクチュエータ間距離が短く無駄時間が少ないため、制御効果が向上する。   Since the microphone is not installed in the vehicle interior 9 by providing the noise estimation unit 32, it is possible to perform control to reduce only vehicle interior noise caused by vibration entering from outside the vehicle without reducing the voice and music of the occupant. I can do it. In addition, the control effect is improved because the distance between the sensor and the actuator is short and the dead time is short as compared with the case of controlling using a microphone.

また、第1及び第2の伝達モデル14、15を備えることにより、当該第1及び第2の伝達モデル14、15を用いて推定を行うことができ、過渡応答を回避して推定精度が良くなる。この結果、騒音低減制御の効果が高まる。   In addition, by providing the first and second transfer models 14 and 15, estimation can be performed using the first and second transfer models 14 and 15, and a transient response is avoided and the estimation accuracy is good. Become. As a result, the effect of noise reduction control is enhanced.

具体的には、第1及び第2の伝達モデル14、15へ入力して得られる出力値との和を車内騒音とすることにより、過渡応答を回避して推定精度が良くなる。   Specifically, by making the sum of the output values obtained by inputting to the first and second transfer models 14 and 15 into the vehicle interior noise, a transient response is avoided and the estimation accuracy is improved.

第1及び第2の伝達モデル14、15に加えて、第3の伝達モデル16を更に備えることにより、アクチュエータ8の制御指令信号が大きくなった場合にも加速度センサ7の出力信号のうち車外から進入する振動成分(外乱)を正確に計測することが出来るので、車室内音圧の推定精度が高まる。また、加速度センサ7とアクチュエータ8の配置自由度が増す。以上の結果、騒音低減制御の効果が高まる。   In addition to the first and second transmission models 14 and 15, the third transmission model 16 is further provided, so that even when the control command signal of the actuator 8 becomes large, the output signal of the acceleration sensor 7 from outside the vehicle. Since the entering vibration component (disturbance) can be accurately measured, the estimation accuracy of the vehicle interior sound pressure is increased. Further, the degree of freedom of arrangement of the acceleration sensor 7 and the actuator 8 is increased. As a result, the effect of noise reduction control is enhanced.

アクチュエータ8から振動までの第3の伝達モデル16のみを用いるため、より簡便でより精度の良い振動推定部32を形成することが出来る。この結果、騒音低減制御の効果が高まる。   Since only the third transmission model 16 from the actuator 8 to the vibration is used, the vibration estimation unit 32 that is simpler and more accurate can be formed. As a result, the effect of noise reduction control is enhanced.

筐体として車体を例に取り、加速度センサとアクチュエータを車体のフロアパネル上に配置することにより、車室内騒音を除去することができる。   By taking a vehicle body as an example of the case and disposing the acceleration sensor and the actuator on the floor panel of the vehicle body, it is possible to remove vehicle interior noise.

(第2の実施の形態)
第2の実施の形態では、本発明をフィードフォワード型能動振動制御装置に適用した場合の例を説明する。
(Second Embodiment)
In the second embodiment, an example in which the present invention is applied to a feedforward active vibration control device will be described.

先ず、図8〜図10を参照して、一般的に知られている適応フィルタを用いたフィードフォワード型能動振動制御装置について説明し、その後、図11及び図12を参照して本実施形態に係わるフィードフォワード型能動振動制御装置について説明する。   First, a feedforward active vibration control apparatus using a generally known adaptive filter will be described with reference to FIGS. 8 to 10, and then the present embodiment will be described with reference to FIGS. 11 and 12. The feedforward active vibration control apparatus will be described.

図8に示すように、従来から知られている適応フィルタを用いたフィードフォワード型能動振動制御装置は、車体に入力される路面の凹凸の影響を測定する加速度センサ39と、車室内空間9の騒音を検出するマイクロフォン40と、フロアパネル上に配置されたアクチュエータ8と、車室内所定空間9での騒音を低減するためにアクチュエータ8への制御指令信号を計算して、アクチュエータ8の動作を制御するコントローラユニット36とを備える。   As shown in FIG. 8, a conventionally known feedforward active vibration control device using an adaptive filter includes an acceleration sensor 39 that measures the influence of road surface unevenness input to the vehicle body, and a vehicle interior space 9. Control the operation of the actuator 8 by calculating a control command signal to the actuator 8 in order to reduce the noise in the predetermined space 9 in the vehicle interior 40 and the microphone 8 that detects noise, the actuator 8 disposed on the floor panel. And a controller unit 36.

コントローラユニット36は、マイクロフォン54からの信号を増幅するアンプ33aと、マイクロフォン54からの信号に基づいて適応フィルタ38を制御する適応則演算部37と、加速度センサ39の出力信号を増幅するアンプ33cと、適応フィルタ38の出力を増幅するアンプ33bとを備える。   The controller unit 36 includes an amplifier 33a that amplifies the signal from the microphone 54, an adaptive law calculation unit 37 that controls the adaptive filter 38 based on the signal from the microphone 54, and an amplifier 33c that amplifies the output signal of the acceleration sensor 39. And an amplifier 33b for amplifying the output of the adaptive filter 38.

図9は図8のシステム全体のブロック線図である。コントローラユニット36は、通常のコンピュータシステムもしくは電子回路で実現されている。第1の振動伝播経路41は加速度センサ39の出力信号から車室内空間9での騒音までの伝達関数を表し、第2の振動伝播経路42はアクチュエータ8への制御指令信号から車室内空間9での騒音までの伝達関数を表す。   FIG. 9 is a block diagram of the entire system of FIG. The controller unit 36 is realized by a normal computer system or an electronic circuit. The first vibration propagation path 41 represents a transfer function from the output signal of the acceleration sensor 39 to the noise in the vehicle interior space 9, and the second vibration propagation path 42 is transmitted from the control command signal to the actuator 8 in the vehicle interior space 9. Represents the transfer function up to noise.

以下、図9に従って図8のシステムにおけるアクチュエータ8への制御指令信号の演算法を説明する。加速度センサ39で得られた信号をx(n)、路面から車体に入力される外乱が車室内空間9に作る騒音をy(n)、ピエゾアクチュエータ8への制御指令信号をu(n)、マイクロフォン54で得られた誤差信号をe(n)とおく。ここで、nはサンプル時刻を表す。また、適応フィルタ38、第1の振動伝播経路41及び第2の振動伝播経路42のモデルをそれぞれ数3、数4、数5とする。

Figure 2007219262
Figure 2007219262
Figure 2007219262
Hereinafter, a method of calculating a control command signal to the actuator 8 in the system of FIG. 8 will be described with reference to FIG. X (n) is a signal obtained by the acceleration sensor 39, y (n) is a noise generated in the vehicle interior space 9 by a disturbance input from the road surface to the vehicle body, u (n) is a control command signal to the piezoelectric actuator 8. The error signal obtained by the microphone 54 is set as e (n). Here, n represents the sample time. Further, the models of the adaptive filter 38, the first vibration propagation path 41, and the second vibration propagation path 42 are represented by Equation 3, Equation 4, and Equation 5, respectively.
Figure 2007219262
Figure 2007219262
Figure 2007219262

ここで、z-iはZ変換のオペレータを表す。なお、wiがサンプル時刻nに依存しているのは、可変係数フィルタを使用しているからである。適応フィルタ38の出力がアクチュエータ8への入力信号として計算される。このとき、信号y(n)とu(n)はそれぞれ数6及び数7として計算される。

Figure 2007219262
Figure 2007219262
Here, z -i represents an operator of Z conversion. The reason why w i depends on the sampling time n is that a variable coefficient filter is used. The output of the adaptive filter 38 is calculated as an input signal to the actuator 8. At this time, the signals y (n) and u (n) are calculated as Equation 6 and Equation 7, respectively.
Figure 2007219262
Figure 2007219262

以上の式に基づいて計算されたu(n)がアクチュエータ8へ入力される。   U (n) calculated based on the above formula is input to the actuator 8.

一方、適応則演算部37ではフィルタ係数wi(n)の更新が以下の演算によって行われる。まず、評価関数Jを数8と定義する。

Figure 2007219262
On the other hand, the adaptive law calculation unit 37 updates the filter coefficient w i (n) by the following calculation. First, the evaluation function J is defined as Equation 8.
Figure 2007219262

一般的に知られたLMS(Least Mean Square)アルゴリズムでは、この評価関数Jがwi(n)に対して最小になるように適用フィルタを更新していく。wi(n)の更新則は数9で与えられる。ここで、KeとKyは設計パラメータである。

Figure 2007219262
In a generally known LMS (Least Mean Square) algorithm, the applied filter is updated so that the evaluation function J is minimized with respect to w i (n). The update rule for w i (n) is given by Equation 9. Here, Ke and Ky are design parameters.
Figure 2007219262

また、r(n-i)は数10で定義する。

Figure 2007219262
Also, r (ni) is defined by Equation 10.
Figure 2007219262

以上の更新演算が適応則演算部37で行われる。   The above update calculation is performed by the adaptive law calculation unit 37.

次に、図10のフローチャートを参照して、適応フィルタ38で行われる演算手順を説明する。   Next, the calculation procedure performed by the adaptive filter 38 will be described with reference to the flowchart of FIG.

(い)先ずステップS201において、加速度センサ39及びマイクロフォン40の出力信号をサンプル時間n毎に取得する。ステップS202に進み、加速度センサ39の出力信号x(n)に対して、数7の計算式にて演算することにより、アクチュエータへの指令信号u(n)を求める。   (Ii) First, in step S201, output signals of the acceleration sensor 39 and the microphone 40 are acquired every sampling time n. Proceeding to step S202, a command signal u (n) to the actuator is obtained by calculating the output signal x (n) of the acceleration sensor 39 using the formula (7).

(ろ)ステップS203に進み、マイクロフォンの出力信号e(n)と、加速度センサ39の出力信号x(n)を用いて、数9及び数10に従って、更新されたフィルタ係数wi(n)を計算する。ステップS204に進み、ステップS203で計算されたフィルタ係数wi(n)を更新する。 (B) Proceeding to step S203, using the output signal e (n) of the microphone and the output signal x (n) of the acceleration sensor 39, the updated filter coefficient w i (n) is calculated according to Equations 9 and 10. calculate. Proceeding to step S204, the filter coefficient w i (n) calculated at step S203 is updated.

(は)ステップS205に進み、ステップS202で計算されたアクチュエータへの制御指令信号u(n)を出力する。ステップS201〜S205をサンプル時間n毎に行うことでフィードフォワード制御が行われ、車室内空間9での騒音を低減することができる。   In step S205, the control command signal u (n) to the actuator calculated in step S202 is output. Feed forward control is performed by performing steps S201 to S205 every sample time n, and noise in the vehicle interior space 9 can be reduced.

次に、本実施形態に係わるフィードフォワード型能動振動騒音制御装置について説明する。   Next, the feedforward active vibration noise control apparatus according to this embodiment will be described.

図11に示すように、第2の実施の形態に係わるフィードフォワード型能動振動制御装置では、図8に比べて、車室内空間9に配置されるマイクロフォンの替わりに第2の加速度センサ7がフロアパネルの所定の箇所に貼り付けられている。また、アンプ33aと適応則演算部37との間に騒音推定部32が配置され、適応フィルタ38からの出力信号が騒音推定部32に入力される。つまり、騒音推定部32には、第2の加速度センサ7から得られた信号にアンプ33aを介して得られる信号と、アクチュエータ8への制御指令信号とが入力される。ここで、騒音推定部32は図3に示した例と同様の装置を用いればよい。騒音推定部32の出力信号は車室内空間9での推定音圧である。よって、適応則演算部37は、図8のマイクロフォンによる実測音圧の代りに、騒音推定部32による推定音圧を入力とする。適応フィルタ38では、加速度センサ39から得られた外乱信号に基づいて車室内空間9での騒音を低減する制御指令信号が演算される。   As shown in FIG. 11, in the feedforward type active vibration control apparatus according to the second embodiment, the second acceleration sensor 7 is a floor instead of the microphone arranged in the vehicle interior space 9 as compared with FIG. Affixed to a predetermined location on the panel. In addition, a noise estimation unit 32 is disposed between the amplifier 33 a and the adaptive law calculation unit 37, and an output signal from the adaptive filter 38 is input to the noise estimation unit 32. That is, the signal obtained from the second acceleration sensor 7 through the amplifier 33 a and the control command signal to the actuator 8 are input to the noise estimation unit 32. Here, the noise estimation unit 32 may use a device similar to the example shown in FIG. The output signal of the noise estimation unit 32 is an estimated sound pressure in the vehicle interior space 9. Therefore, the adaptive law calculating unit 37 receives the estimated sound pressure by the noise estimating unit 32 instead of the actually measured sound pressure by the microphone of FIG. The adaptive filter 38 calculates a control command signal for reducing noise in the vehicle interior space 9 based on the disturbance signal obtained from the acceleration sensor 39.

図12を参照して、図11のコントローラユニット36の演算手順を説明する。   With reference to FIG. 12, the calculation procedure of the controller unit 36 of FIG. 11 will be described.

(A)先ずステップS301において、車外からタイヤを介して車体へ伝達する振動を計測するために用いる第1の加速度センサ39の出力信号と、車室内騒音を推定するために用いる第2の加速度センサ7の出力信号とをそれぞれ取得する。   (A) First, in step S301, the output signal of the first acceleration sensor 39 used for measuring vibration transmitted from the outside of the vehicle via the tire to the vehicle body, and the second acceleration sensor used for estimating vehicle interior noise. 7 output signals are obtained.

(B)ステップS302に進み、第1の加速度センサ39の出力信号x(n)に対して、数7の演算を行うことで、アクチュエータ8への制御指令信号u(n)を求める。ステップS303に進み、第2の加速度センサ7の出力信号とアクチュエータ8への制御指令信号u(n)を用いて、図5に示すアルゴリズムによる騒音推定演算を行うことで、車室内空間9での推定音圧(推定騒音)を計算する。   (B) Proceeding to step S302, the control command signal u (n) to the actuator 8 is obtained by performing the calculation of Expression 7 on the output signal x (n) of the first acceleration sensor 39. Proceeding to step S303, using the output signal of the second acceleration sensor 7 and the control command signal u (n) to the actuator 8 to perform noise estimation calculation by the algorithm shown in FIG. Calculate the estimated sound pressure (estimated noise).

(C)ステップS304に進み、ステップS303で計算された推定騒音と、第1の加速度センサ39の出力信号x(n)を用いて、数9及び数10の演算式に従って、更新されたフィルタ係数wi(n)を計算する。ステップS305に進み、ステップS304で計算されたフィルタ係数wi(n)を更新する。 (C) Proceeding to step S304, using the estimated noise calculated in step S303 and the output signal x (n) of the first acceleration sensor 39, the updated filter coefficients according to the mathematical expressions of Equations 9 and 10 Calculate w i (n). Proceeding to step S305, the filter coefficient w i (n) calculated at step S304 is updated.

(D)ステップS306に進み、ステップS202で計算されたアクチュエータ8への制御指令信号u(n)を出力する。ステップS301〜S306をサンプル時間n毎に行うことでフィードフォワード制御が行われ、車室内空間9での騒音が低減される。   (D) Proceeding to step S306, the control command signal u (n) to the actuator 8 calculated at step S202 is output. Feed forward control is performed by performing steps S301 to S306 every sample time n, and noise in the vehicle interior space 9 is reduced.

このように、騒音推定部32により車室内空間9における騒音を推定することにより、図8に示した従来型と比較してマイクロフォンを装備する必要が無く、フィードフォワード型能動振動制御装置を実現することができる。従って、第1の実施の形態と同様に、車室内空間9にマイクロフォンを設置しないので、乗員の声や音楽を低減することなく、車外から進入する振動に起因した車室内騒音のみを低減する制御を行うことが出来る。また、マイクを用いて制御する場合に比べ、加速度センサ−アクチュエータ間距離が短く無駄時間が少ないため、制御効果が向上する。   Thus, by estimating the noise in the vehicle interior space 9 by the noise estimation unit 32, it is not necessary to equip a microphone as compared with the conventional type shown in FIG. 8, and a feedforward type active vibration control device is realized. be able to. Therefore, as in the first embodiment, since no microphone is installed in the vehicle interior space 9, control is performed to reduce only vehicle interior noise caused by vibrations entering from outside the vehicle without reducing the voice and music of the occupant. Can be done. In addition, the control effect is improved because the distance between the acceleration sensor and the actuator is short and the dead time is short as compared with the case where control is performed using a microphone.

なお、説明を簡便にするため、加速度センサの数を第1の実施の形態では1つ、第2の実施の形態では2つとし、アクチュエータの数を1つとしたシステムを例に取り説明したが、センサやアクチュエータを複数備えたシステムに対しても同様の議論が可能であり、本発明を適用することができる。すなわち、加速度センサ及びアクチュエータの数は本発明の特許請求の範囲を狭めるものではない。   In order to simplify the description, the number of acceleration sensors is one in the first embodiment, two in the second embodiment, and a system with one actuator is described as an example. The same discussion can be made for a system including a plurality of sensors and actuators, and the present invention can be applied. That is, the number of acceleration sensors and actuators does not narrow the scope of the claims of the present invention.

実際に、r個の加速度センサおよびs個のアクチュエータを用いてt個の車室内空間における騒音低減制御を行う場合、図1及び図11中の騒音推定部32は図13に示す構成に変更すればよい。図13において、信号線横のr、s、tの記号はその信号線のチャンネル数を表している。第1の伝達モデル14は加速度センサの出力信号から車室内空間9の騒音までのt行r列の伝達関数マトリクスであり、第2の伝達モデル15はアクチュエータへの制御指令信号から車室内空間9の騒音までのt行s列の伝達関数マトリクスであり、第3の伝達モデル16はアクチュエータへの制御指令信号から加速度センサの取付け位置の振動までのr行s列の伝達関数マトリクスである。このように、加速度センサおよびアクチュエータが複数になることで、出力信号や制御指令信号が伝達する信号線を複数チャンネルにして、騒音推定部32内にあらかじめ保持しておく第1乃至第3の伝達モデル14〜16をマトリクス形式にする。このことにより、複数センサ及びアクチュエータの場合にも、第1乃至第2の実施の形態を容易に拡張できる。   Actually, when performing noise reduction control in t vehicle interior spaces using r acceleration sensors and s actuators, the noise estimation unit 32 in FIGS. 1 and 11 is changed to the configuration shown in FIG. That's fine. In FIG. 13, symbols r, s, and t beside a signal line indicate the number of channels of the signal line. The first transfer model 14 is a transfer function matrix of t rows and r columns from the output signal of the acceleration sensor to the noise in the vehicle interior space 9, and the second transfer model 15 is a vehicle interior space 9 from the control command signal to the actuator. The third transfer model 16 is an r row and s column transfer function matrix from a control command signal to the actuator to a vibration at the mounting position of the acceleration sensor. As described above, since there are a plurality of acceleration sensors and actuators, the signal lines for transmitting the output signal and the control command signal are made into a plurality of channels, and the first to third transmissions are held in advance in the noise estimation unit 32. The models 14 to 16 are made into a matrix form. Thus, the first and second embodiments can be easily expanded even in the case of a plurality of sensors and actuators.

(第3の実施の形態)
第3の実施の形態に係わる能動振動騒音制御装置を説明する前に、先ず、図14(a)、図14(b)及び図15を参照して、当該能動振動騒音制御装置が搭載される車両の構造について説明する。
(Third embodiment)
Before explaining the active vibration noise control apparatus according to the third embodiment, first, the active vibration noise control apparatus is mounted with reference to FIG. 14 (a), FIG. 14 (b) and FIG. The structure of the vehicle will be described.

図14(a)及び14(b)に示すように、アクチュエータが固定されている車体5のフロアパネル6は、3つの主要な部分に分けることができる。具体的には、フロアパネル6は、乗員が座るシートの下に位置する平板状のキャビンフロア17と、燃料タンクが固定されるタンクパネル18と、スペアタイヤ及びトランクの下が位置するスペアタイヤパネル19とを備える。キャビンフロア17は、エンジンルームとタンクパネル18の間に位置している。スペアタイヤパネル19は、タンクパネル18と車体5の後端の間に位置している。   As shown in FIGS. 14A and 14B, the floor panel 6 of the vehicle body 5 to which the actuator is fixed can be divided into three main parts. Specifically, the floor panel 6 includes a flat cabin floor 17 positioned under a seat on which an occupant sits, a tank panel 18 to which a fuel tank is fixed, and a spare tire panel positioned under a spare tire and a trunk. 19. The cabin floor 17 is located between the engine room and the tank panel 18. The spare tire panel 19 is located between the tank panel 18 and the rear end of the vehicle body 5.

図15に示すように、図14(b)のフロアパネル6の上には、車両の進行方向に向かって縦方向に伸びる縦メンバ20a〜20dと、横方向に伸びる横メンバ21a〜21dとが配置されている。フロアパネル6は車体を軽量化するべく薄い平板状に形成されているため、フロアパネル6だけでは十分な剛性が得られない。フロアパネル6よりも剛性が高い棒状の縦メンバ20a〜20d及び横メンバ21a〜21dをフロアパネル6上に固定して、フロアパネル6の剛性を高めている。縦メンバ20a〜20d及び横メンバ21a〜21dの厚みはフロアパネル6よりも厚く形成されている。縦メンバ20a〜20d及び横メンバ21a〜21dは、フロアパネル6の一方の面(表面又は裏面)側のみに固定されていても、両方の面(表裏面)に固定されていても構わない
以下、本発明の第3の実施の形態に係わる能動振動騒音制御装置について説明する。一般的に、路面からタイヤを介して車体に伝播する振動及び騒音は、波としての性質を有するため、フロアパネルの振動が車室内空間の騒音まで伝達するのに、ある一定の遅延時間を要する。つまり、振動及び騒音は、フロアパネルから乗員の耳まで瞬間的に伝達するものではなく、一定の遅延時間をかけて伝達する。この遅延時間によって騒音推定部32の騒音推定性能は大きく左右される。騒音推定部32は、第1の伝達モデルを用いて、加速度センサの取り付け位置における振動から車室内空間の騒音を推定し、第2の伝達モデルを用いて、アクチュエータへの制御指令信号から車室内空間の騒音を推定している。よって、フロアパネルから車室内空間までの遅延時間はこの騒音推定部32に含まれることが望ましい。第3の実施の形態では、フロアパネルから車室内空間まで波として伝達する振動及び騒音の遅延時間を考慮した車室内騒音を推定する装置及び方法について説明する。
As shown in FIG. 15, on the floor panel 6 in FIG. 14B, there are vertical members 20 a to 20 d extending in the vertical direction toward the traveling direction of the vehicle, and horizontal members 21 a to 21 d extending in the horizontal direction. Has been placed. Since the floor panel 6 is formed in a thin flat plate shape to reduce the weight of the vehicle body, sufficient rigidity cannot be obtained with the floor panel 6 alone. The bar-shaped vertical members 20a to 20d and the horizontal members 21a to 21d having higher rigidity than the floor panel 6 are fixed on the floor panel 6 to increase the rigidity of the floor panel 6. The vertical members 20 a to 20 d and the horizontal members 21 a to 21 d are formed thicker than the floor panel 6. The vertical members 20a to 20d and the horizontal members 21a to 21d may be fixed to only one surface (front surface or back surface) side of the floor panel 6 or may be fixed to both surfaces (front and back surfaces). Next, an active vibration noise control apparatus according to the third embodiment of the present invention will be described. In general, vibration and noise propagating from the road surface to the vehicle body via the tire have properties as waves. Therefore, it takes a certain delay time for the vibration of the floor panel to be transmitted to the noise in the passenger compartment. . That is, vibration and noise are not transmitted instantaneously from the floor panel to the passenger's ears, but are transmitted with a certain delay time. The noise estimation performance of the noise estimation unit 32 is greatly influenced by this delay time. The noise estimation unit 32 estimates the noise in the vehicle interior space from the vibration at the mounting position of the acceleration sensor using the first transmission model, and uses the second transmission model to calculate the vehicle interior from the control command signal to the actuator. Estimates the noise in the space. Therefore, it is desirable that the delay time from the floor panel to the vehicle interior space is included in the noise estimation unit 32. In the third embodiment, an apparatus and method for estimating vehicle interior noise in consideration of vibration and noise delay times transmitted as waves from the floor panel to the vehicle interior space will be described.

図16に示すように、車室内空間にマイクロフォン40を配置して車室内騒音を実測してフロアパネル上のアクチュエータへの制御指令信号をコントロールするフィードバック型の能動振動騒音制御装置(NRC:Noise Reduction Control)が一般的に知られている。この装置は、車室内における1又は2以上の測定ポイントにおける騒音(音圧レベルSPL:Sound Pressure Level [dB])が低減されるように、アクチュエータ8を用いて車体のフロアパネルに制御振動を生成することを目的としている。1又は2以上のアクチュエータ8をフロアパネルの所定の位置に設置して制御振動をフロアパネル上に生成する。そして、1又は2以上のマイクロフォン40を車室内の所定の位置に設置して車室内騒音を測定する。マイクロフォン40により測定された音圧レベルSPLは、フィードバック型のコントローラ43に入力される。コントローラ43として、例えば、一般的なH∞制御を用いることができる。外乱によるフロアパネルの振動及びアクチュエータ8による制御振動は、それぞれ所定の遅延時間をかけて車室内騒音として伝達され、マイクロフォン40により音圧レベルとして測定される。よって、図16に示す制御システムには、物理的に存在する総ての遅延時間を含んだ構成を形成していることになる。ここで、総ての遅延時間には、外乱(ここでは、路面の凹凸によりタイヤを介して車体に侵入する振動及び騒音を示す)から車室内空間までの第1の遅延時間と、アクチュエータから車室内空間までの第2の遅延時間とを含まれる。   As shown in FIG. 16, a feedback type active vibration noise control device (NRC: Noise Reduction) that arranges a microphone 40 in a vehicle interior space, measures vehicle interior noise, and controls a control command signal to an actuator on the floor panel. Control) is generally known. This device uses the actuator 8 to generate control vibrations on the floor panel of the vehicle body so that noise (sound pressure level SPL: Sound Pressure Level [dB]) at one or more measurement points in the passenger compartment is reduced. The purpose is to do. One or more actuators 8 are installed at predetermined positions on the floor panel to generate a control vibration on the floor panel. Then, one or two or more microphones 40 are installed at predetermined positions in the vehicle interior, and vehicle interior noise is measured. The sound pressure level SPL measured by the microphone 40 is input to the feedback type controller 43. As the controller 43, for example, general H∞ control can be used. The vibration of the floor panel due to the disturbance and the control vibration by the actuator 8 are transmitted as vehicle interior noise over a predetermined delay time, and are measured by the microphone 40 as the sound pressure level. Therefore, the control system shown in FIG. 16 has a configuration including all physically existing delay times. Here, all delay times include the first delay time from disturbance (in this case, vibration and noise entering the vehicle body via the tire due to road surface irregularities) to the vehicle interior space, and from the actuator to the vehicle. And a second delay time to the indoor space.

これに対して、図17は、室内騒音をマイクロフォンで直接測定することなく、車体の振動から車室内騒音までの伝達モデルを用いて車室内騒音を推定する能動振動騒音制御装置を示す。図17の装置では、騒音推定部32を仮想的なマイクロフォンとして備える。騒音推定部32が備える伝達関数と実際の伝達システム間に狂いが無ければ、図16と図17のシステムは正確に一致し、マイクロフォン40を騒音推定部32で置き換えることが可能である。   On the other hand, FIG. 17 shows an active vibration noise control apparatus that estimates vehicle interior noise using a transfer model from vehicle body vibration to vehicle interior noise without directly measuring room noise with a microphone. In the apparatus of FIG. 17, the noise estimation unit 32 is provided as a virtual microphone. If there is no error between the transfer function provided in the noise estimation unit 32 and the actual transfer system, the systems in FIG. 16 and FIG. 17 match exactly, and the microphone 40 can be replaced with the noise estimation unit 32.

コントローラ31として、例えば、一般的なH∞制御を用いることができる。アクチュエータの制御指令信号Vと車室内音圧SPL間の開ループにおける伝達モデルは、一般的な伝達関数の形式で求めることができる。また、複数の制御指令信号V及び複数の車室内音圧SPLを含むシステムにおいても容易に図16及び図17の装置を適用可能であることは前述したとおりである。開ループからなるシステムは、数11で表される。ここで、sはラプラス変換における変数である。

Figure 2007219262
As the controller 31, for example, general H∞ control can be used. An open-loop transfer model between the actuator control command signal V and the vehicle interior sound pressure SPL can be obtained in the form of a general transfer function. Further, as described above, the apparatus shown in FIGS. 16 and 17 can be easily applied to a system including a plurality of control command signals V and a plurality of vehicle interior sound pressures SPL. A system consisting of an open loop is expressed by Equation 11. Here, s is a variable in Laplace transform.
Figure 2007219262

コントローラ31は、ループ整形による設計法(the loop shaping design procedure)を用いて計算することができる。この計算手法は、D.McFarlane and K.Golver,IEEE Transactions on Automatic Control.vol.37,no.6,June 1992,pp.759−769,“A Loop Shaping Design Procedure Using H∞ Synthesis” に開示されている。ここでは、数12に示す計算式によりCを求めている。

Figure 2007219262
The controller 31 can perform the calculation using a loop shaping design procedure. This calculation method is described in D.A. McFarlane and K.M. Golver, IEEE Transactions on Automatic Control. vol. 37, no. 6, June 1992, pp. 759-769, “A Loop Shaping Design Procedure Using H∞ Synthesis”. Here, C∞ is obtained by the calculation formula shown in Equation 12.
Figure 2007219262

ここで、Gsは、W1とWで重み付けされたシステムの伝達関数であって、整形関数(the shaping function)と呼ぶ。

Figure 2007219262
Here, Gs is a transfer function of the system weighted by W 1 and W 2 and is called a shaping function.
Figure 2007219262

このような条件式を満たす図16に示すコントローラCは、数14で表される。

Figure 2007219262
The controller C shown in FIG. 16 that satisfies such a conditional expression is expressed by Equation 14.
Figure 2007219262

ここで、W1とWは、通常、制御対象とする振動及び騒音の周波数帯を選択することができるフィルターの伝達関数である。[Ms]とNsは数15に示すように、標準化されたGsのコプライム因子である。なお、パラメータεは、安定許容範囲と呼ばれ、通常0.2〜0.3が適当である。

Figure 2007219262
Here, W 1 and W 2 are usually transfer functions of a filter that can select a frequency band of vibration and noise to be controlled. [Ms] and Ns are standardized Gs coprime factors as shown in Equation 15. The parameter ε is called a stable allowable range, and 0.2 to 0.3 is usually appropriate.
Figure 2007219262

次に、図17の騒音推定部32について説明する。騒音推定部Eは、車室内騒音SPLの伝達におけるマイクロフォン40に置き換えられる。騒音推定部Eは、実際には、車室以内騒音[SPL]の推定を行うため、騒音推定部Eが使用する伝達モデルはできる限り正確であることが望まれる。以後、推定された車室以内騒音(音圧)を[SPL]と表記し、マイクロフォンで実際に測定された車室内騒音をSPLと表記する。   Next, the noise estimation unit 32 in FIG. 17 will be described. The noise estimation unit E is replaced with the microphone 40 in the transmission of the vehicle interior noise SPL. Since the noise estimation unit E actually estimates the in-cabin noise [SPL], it is desirable that the transmission model used by the noise estimation unit E is as accurate as possible. Hereinafter, the estimated in-cabin noise (sound pressure) is denoted as [SPL], and the cabin noise actually measured by the microphone is denoted as SPL.

図18に示すように、第3の実施の形態に係わるコントローラユニット30は、騒音推定部32とコントローラ31とを備える。プラント5は、外乱d及びアクチュエータの制御振動から車室内騒音SPLiまで振動及び騒音が伝達する筐体の一例としての車体を示す。騒音推定部32は、外乱dとして加速度センサから7の出力信号とコントローラ31からプラント5内のアクチュエータ8へ送信される制御指令信号とを受信して車室内騒音[SPL]iを推定する。コントローラ31は、車室内騒音[SPL]iに基づいて、プラントから出力される車室内騒音SPLiが軽減されるようにアクチュエータ8への制御指令信号を制御する。   As shown in FIG. 18, the controller unit 30 according to the third embodiment includes a noise estimation unit 32 and a controller 31. The plant 5 shows a vehicle body as an example of a housing that transmits vibration and noise from disturbance d and actuator control vibration to vehicle interior noise SPLi. The noise estimation unit 32 receives the output signal 7 from the acceleration sensor and the control command signal transmitted from the controller 31 to the actuator 8 in the plant 5 as the disturbance d, and estimates the vehicle interior noise [SPL] i. Based on the vehicle interior noise [SPL] i, the controller 31 controls the control command signal to the actuator 8 so that the vehicle interior noise SPLi output from the plant is reduced.

図19に示すように、図18の騒音推定部32は、図17の加速度センサ7が取り付けられた位置から車室内騒音SPLを推定するポイントまでの伝達特性を表す第1の伝達モデルGdと、アクチュエータ8への制御指令信号から車室内騒音SPLまでの伝達特性を示す第2の伝達モデルGaとをデータとして備える。騒音推定部Eは、加速度センサ7の出力信号を第1の伝達モデルGdへ入力して得られる出力値と、制御指令信号を第2の伝達モデルGaへ入力して得られる出力値との和を車室内騒音[SPL]iとして出力する。   As shown in FIG. 19, the noise estimation unit 32 in FIG. 18 includes a first transfer model Gd representing transfer characteristics from a position where the acceleration sensor 7 in FIG. 17 is attached to a point for estimating the vehicle interior noise SPL, A second transfer model Ga indicating transfer characteristics from the control command signal to the actuator 8 to the vehicle interior noise SPL is provided as data. The noise estimation unit E adds the output value obtained by inputting the output signal of the acceleration sensor 7 to the first transfer model Gd and the output value obtained by inputting the control command signal to the second transfer model Ga. Is output as vehicle interior noise [SPL] i.

ここで、第2の伝達モデルGaは、アクチュエータ8が配置されたフロアパネル位置から車室内騒音を推定するポイントまでの伝達特性を示すものであっても構わない。   Here, the second transmission model Ga may indicate a transmission characteristic from the floor panel position where the actuator 8 is disposed to a point where the vehicle interior noise is estimated.

第1の伝達モデルGd及び第2の伝達モデルGaは、振動が伝達するプラント5を特定することにより求められ、実際のプラント5の物理的な総ての特徴を含有していることが望ましい。よって、第1及び第2の伝達モデルには、それぞれフロアパネルから車室内騒音まで騒音及び振動が波として伝達するために必要な時間(遅延時間)が含まれていることが望ましい。   The first transmission model Gd and the second transmission model Ga are obtained by specifying the plant 5 to which the vibration is transmitted, and preferably include all the physical characteristics of the actual plant 5. Therefore, it is desirable that the first and second transmission models include time (delay time) necessary for transmitting noise and vibration as waves from the floor panel to the vehicle interior noise.

騒音推定部32aは、図19に示すように、アクチュエータへの制御指令信号Vと、外乱を示す信号αが入力される。制御指令信号Vは、結果的には、アクチュエータ8が設置されている位置におけるフロアパネル6の振動に相当する。また、外乱を示す信号αは、加速度センサ7により測定されたフロアパネルの振動に相当する。   As shown in FIG. 19, the noise estimation unit 32a receives a control command signal V to the actuator and a signal α indicating disturbance. As a result, the control command signal V corresponds to the vibration of the floor panel 6 at the position where the actuator 8 is installed. The signal α indicating the disturbance corresponds to the vibration of the floor panel measured by the acceleration sensor 7.

騒音推定部32aが十分な精度によって車室内騒音[SPL]を推定する場合、図17に示す能動振動制御装置は、図16に示すマイクロフォンで車室内騒音SPLを実測するタイプの能動振動制御装置に正確に一致する。図17のシステムの利点は、車室内に配置するマイクロフォンが不要となる点であり、ロードノイズと共に、ラジオの音、乗員の声、風の音なども制御振動によりかき消されてしまう危険性が無くなる。   When the noise estimation unit 32a estimates the vehicle interior noise [SPL] with sufficient accuracy, the active vibration control device shown in FIG. 17 is an active vibration control device of the type that actually measures the vehicle interior noise SPL with the microphone shown in FIG. Match exactly. The advantage of the system of FIG. 17 is that a microphone disposed in the passenger compartment is not required, and there is no risk that the radio noise, the occupant's voice, the wind noise, etc. will be drowned out by the control vibration as well as the road noise.

しかしながら、騒音推定部32aは実際のプラント5が有する遅延時間を含んでいるため、コントローラの高い制御効率は期待できない。そこで、第3の実施の形態では、より良い制御効率、つまり騒音除去効率を実現するため、新たな騒音推定部Eの構造を提案する。即ち、遅延時間を補償してプラント5に実際に生じる遅延時間よりも短い遅延時間によって車室内騒音SPLを推定することにより、より良い制御効率を実現する。   However, since the noise estimation unit 32a includes the delay time that the actual plant 5 has, high control efficiency of the controller cannot be expected. Therefore, in the third embodiment, a new structure of the noise estimation unit E is proposed in order to realize better control efficiency, that is, noise removal efficiency. That is, better control efficiency is realized by compensating the delay time and estimating the vehicle interior noise SPL based on a delay time shorter than the delay time actually generated in the plant 5.

時間遅延を含まないプラントは、遅延時間を含むプラントにくらべてより単純な構造を有しており、制御効率が高く、得られる制御ゲインも大きい。これに対して、遅延時間を含むプラントの場合、騒音制御が難しく、制御ゲインを高くすることができないため、制御効率は通常低くなってしまう。そこで、プラントの異なる遅延時間を計算及び評価し、これを騒音推定部Eから取り除く。こうして得られた新たな騒音推定部[E]は、単に車室内騒音SPLを推定するのではなく、推定される車室内騒音SPLを予測することができる。   A plant that does not include a time delay has a simpler structure than a plant that includes a delay time, and has high control efficiency and a large control gain. On the other hand, in the case of a plant including a delay time, noise control is difficult and the control gain cannot be increased, so that the control efficiency is usually low. Therefore, a different delay time of the plant is calculated and evaluated, and this is removed from the noise estimation unit E. The new noise estimation unit [E] obtained in this way can predict the estimated vehicle interior noise SPL rather than simply estimating the vehicle interior noise SPL.

図20は、遅延時間を考慮した騒音推定部32bの構成を示す。外乱に関する情報αは、先ず第1の伝達モデル14bに入力され、その後、遅延時間モデル44に入力される。一方、アクチュエータへの制御指令信号Vは第2の伝達モデル15bに入力され、その後、遅延時間モデル45へ入力される。   FIG. 20 shows a configuration of the noise estimation unit 32b in consideration of the delay time. The information α regarding the disturbance is first input to the first transmission model 14 b and then input to the delay time model 44. On the other hand, the control command signal V to the actuator is input to the second transmission model 15 b and then input to the delay time model 45.

遅延時間補償τを用いて遅延時間Ta及びTdを短縮(補償)する方法を説明する。騒音推定部[E]は、数16に示すようなGa(s)及びGd(s)とを備える。

Figure 2007219262
A method for shortening (compensating) the delay times Ta and Td using the delay time compensation τ will be described. The noise estimation unit [E] includes Ga (s) and Gd (s) as shown in Equation 16.
Figure 2007219262

数16において、第1の伝達モデル[G]a(s)及び第2の伝達モデル[G]d(s)は、それぞれGa(s)及びGd(s)の遅延時間の項を除去したものに相当する。図21では、1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLとを含むシステムを示している。   In Equation 16, the first transfer model [G] a (s) and the second transfer model [G] d (s) are obtained by removing the delay time terms of Ga (s) and Gd (s), respectively. It corresponds to. FIG. 21 shows a system including one actuator 8, one disturbance (acceleration sensor 7), and one vehicle interior noise SPL.

騒音推定部[E]は数16に示したGa(s)及びGd(s)を使用するため、何らハードウェアの変更なしに、遅延時間補償τを用いて容易に遅延時間Ta及びTdを短縮することができる。この特徴が、マイクロフォンの代りに騒音推定部が出力する信号を用いる利点である。遅延時間Tdが遅延時間Ta以上である場合、遅延時間補償τは遅延時間Taに等しくなる。この場合、アクチュエータは、時間遅延が存在しないかのような振る舞いをする。ここで[E]は、遅延時間を補償した騒音推定部を示す。   Since the noise estimation unit [E] uses Ga (s) and Gd (s) shown in Equation 16, the delay times Ta and Td can be easily shortened using the delay time compensation τ without any hardware change. can do. This feature is an advantage of using a signal output from the noise estimation unit instead of the microphone. When the delay time Td is equal to or longer than the delay time Ta, the delay time compensation τ is equal to the delay time Ta. In this case, the actuator behaves as if there is no time delay. Here, [E] indicates a noise estimation unit that compensates for the delay time.

騒音推定部[E]を用いてコントローラユニットCを演算する。コントローラユニットCの演算は、前述したとおりであるが、数11の代りに、数17を用いる。

Figure 2007219262
The controller unit C is calculated using the noise estimation unit [E]. The calculation of the controller unit C is as described above, but the equation 17 is used instead of the equation 11.
Figure 2007219262

図22は、遅延時間を補償した騒音推定部32cの構成を示す。外乱に関する信号は先ず第1の伝達モデル14bに入力され、その後、τにより遅延時間補償された遅延時間モデル44bに入力される。一方、アクチュエータへの制御指令信号Vは第2の伝達モデル15bに入力され、その後、τにより遅延時間補償された遅延時間モデル45bへ入力される。遅延時間モデル44b、45bの出力信号を重ね合わせることにより車室内騒音[SPL]を推定することができる。   FIG. 22 shows a configuration of the noise estimation unit 32c that compensates for the delay time. A signal related to the disturbance is first input to the first transfer model 14b, and then input to the delay time model 44b compensated for the delay time by τ. On the other hand, the control command signal V to the actuator is input to the second transmission model 15b and then input to the delay time model 45b compensated for the delay time by τ. Vehicle interior noise [SPL] can be estimated by superimposing the output signals of the delay time models 44b and 45b.

図23は、1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLからなるシステムを示し、伝達経路上には、遅延時間補償された遅延時間モデルだけを示している。   FIG. 23 shows a system composed of one actuator 8, one disturbance (acceleration sensor 7) and one vehicle interior noise SPL, and only a delay time model compensated for the delay time is shown on the transmission path.

遅延時間補償された遅延時間Ta及びTdの計算又は推定には3つの方法がある。   There are three methods for calculating or estimating the delay times Ta and Td compensated for the delay time.

第1に、遅延時間は、伝達する距離の計算により推定することができる。時間遅延は本質的には空気中の音速に依存している。例えば、アクチュエータ8が車室内騒音SPLを測定する箇所から距離da離れた場所に配置されている場合、空気中の音速Vsを用いて数18により表すことができる。

Figure 2007219262
First, the delay time can be estimated by calculating the distance traveled. The time delay is essentially dependent on the speed of sound in the air. For example, when the actuator 8 is arranged at a location da away from the location where the vehicle interior noise SPL is measured, it can be expressed by Equation 18 using the sound velocity Vs in the air.
Figure 2007219262

第2に、遅延時間は信号処理の速度から推定することができる。入力信号と出力信号を測定することにより、システムの特定が可能となる。したがって、入力信号と出力信号の時間差を計算することができる。   Second, the delay time can be estimated from the speed of signal processing. The system can be specified by measuring the input signal and the output signal. Therefore, the time difference between the input signal and the output signal can be calculated.

第3に、システムを特定して伝達モデルG(s)を取得した後、伝達モデルGの不安定なゼロとともに遅延時間を推定できる。   Third, after acquiring the transfer model G (s) by specifying the system, the delay time can be estimated together with the unstable zero of the transfer model G.

遅延時間Tを含む一般的な伝達関数G(s)は、数19により得られる。

Figure 2007219262
A general transfer function G (s) including the delay time T is obtained by Equation 19.
Figure 2007219262

伝達関数G(s)は遅延時間Tを含むため、理想的な伝達関数N(s)/D(s)は、N(s)の中に不安定なゼロを含み、結果的に安定したゼロも含む。N(s)は数20に示すように分解される。

Figure 2007219262
Since the transfer function G (s) includes a delay time T, the ideal transfer function N (s) / D (s) includes an unstable zero in N (s), resulting in a stable zero. Including. N (s) is decomposed as shown in Equation 20.
Figure 2007219262

[N]は安定零点を有するN(s)の一部を示し、N'は不安定零点を含むN(s)の他の一部を示す。よって、数19は、数21になる。

Figure 2007219262
[N] represents a part of N (s) having a stable zero, and N ′ represents another part of N (s) including an unstable zero. Therefore, Equation 19 becomes Equation 21.
Figure 2007219262

数21は数22に変形することができる。

Figure 2007219262
Equation 21 can be transformed into Equation 22.
Figure 2007219262

[N](s)は安定化された[N](-s)であり、例えば、[N](-s)の根は、[N](s)の根に相対する。

Figure 2007219262
[N] (s) is a stabilized [N] (-s), for example, the root of [N] (-s) is relative to the root of [N] (s).
Figure 2007219262

[G]は、分子が安定な伝達関数G(s)であり、全域通過伝達関数[N](-s)/[N](s)は、遅延時間項の近似で用いられる周知のパデ近似(Pade approximation)を表す。遅延時間項を数24に示す。

Figure 2007219262
[G] is a transfer function G (s) in which the numerator is stable, and the all-pass transfer function [N] (− s) / [N] (s) is a well-known Padé approximation used in the approximation of the delay term. (Pade approximation). The delay time term is shown in Equation 24.
Figure 2007219262

G'(s)とG(s)の極は同じである。不安定なゼロが伝達関数の中に無いため、G'(s)は「最小位相(minimum phase)」であり、非最小位相(non-minimum phase)伝達関数と呼ばれるG(s)とは異なる。   The poles of G ′ (s) and G (s) are the same. Since there is no unstable zero in the transfer function, G ′ (s) is a “minimum phase” and is different from G (s) called a non-minimum phase transfer function. .

次に、遅延時間補償された騒音推定部Eを有するコントローラの制御効率を解析する。最初に、実際にマイクロフォンを用いて制御効率を解析する。このとき、騒音推定部Eは使用しない。次に、遅延時間補償が成されていない騒音推定部Eを使用したコントローラの制御効率を解析する。最後に、遅延時間補償が成された騒音推定部Eを使用したコントローラの制御効率を解析する。なお、1つのアクチュエータと、1つの外乱(加速度センサ)と、1つの車室内騒音SPLからなるシステムについて上記の分析を実施する。   Next, the control efficiency of the controller having the noise estimation unit E compensated for the delay time is analyzed. First, the control efficiency is actually analyzed using a microphone. At this time, the noise estimation unit E is not used. Next, the control efficiency of the controller using the noise estimation unit E that is not compensated for the delay time is analyzed. Finally, the control efficiency of the controller using the noise estimation unit E for which the delay time compensation has been performed is analyzed. In addition, said analysis is implemented about the system which consists of one actuator, one disturbance (acceleration sensor), and one vehicle interior noise SPL.

先ず、マイクロフォンを用いたシステムについて制御効率を解析する。図24は、車室内騒音SPLを実際に測定するマイクロフォンを備える能動振動制御装置の構成を示す。車室内騒音SPLは数25により得られる。

Figure 2007219262
First, the control efficiency of a system using a microphone is analyzed. FIG. 24 shows a configuration of an active vibration control apparatus including a microphone that actually measures the vehicle interior noise SPL. The vehicle interior noise SPL is obtained by Equation 25.
Figure 2007219262

コントローラCがアクチュエータを駆動するために出力する制御指令信号Vは、数26により得られる。

Figure 2007219262
A control command signal V output for the controller C to drive the actuator is obtained by Equation 26.
Figure 2007219262

数25と数26の計算式を組み合わせると、SPLは数27になる。

Figure 2007219262
Combining Equations 25 and 26, SPL becomes Equation 27.
Figure 2007219262

更に変形すると、数28となる。

Figure 2007219262
Further deformation yields Equation 28.
Figure 2007219262

上記のSPL(s)は伝達関数H(s)の中には表れない。しかし、分母及び分子は、それぞれ遅延時間Ta及びTdを含む。この典型的な形式は制御が難しいことで知られており、具体的には、コントローラCのゲインを大きくすることができず、発散する危険性が高い。   The above SPL (s) does not appear in the transfer function H (s). However, the denominator and numerator include delay times Ta and Td, respectively. This typical form is known to be difficult to control. Specifically, the gain of the controller C cannot be increased, and there is a high risk of divergence.

次に、遅延時間補償のない仮想的マイクロフォンを備えるシステムについて制御効率を解析する。   Next, the control efficiency is analyzed for a system including a virtual microphone without delay time compensation.

図25は、仮想的なマイクロフォン(騒音推定部)であって遅延時間補償のないものを備える能動振動制御装置を示す。騒音推定部32a内の各伝達モデル14、15は、例えばプラント5の特定により推定される。数27に示す車室内騒音SPLは同じである。しかし、コントローラCが出力するアクチュエータ8を駆動するための制御指令信号Vは数29となる。

Figure 2007219262
FIG. 25 shows an active vibration control apparatus including a virtual microphone (noise estimation unit) without delay time compensation. The transmission models 14 and 15 in the noise estimation unit 32a are estimated by specifying the plant 5, for example. The vehicle interior noise SPL shown in Equation 27 is the same. However, the control command signal V for driving the actuator 8 output from the controller C is expressed by Equation 29.
Figure 2007219262

ここでSPLは数30で表される。

Figure 2007219262
Here, SPL is expressed by Equation 30.
Figure 2007219262

よって、制御指令信号Vは数31となる。

Figure 2007219262
Therefore, the control command signal V is expressed by Equation 31.
Figure 2007219262

そして、数15のSPL(s)は数32になる。

Figure 2007219262
Then, SPL (s) in Expression 15 becomes Expression 32.
Figure 2007219262

伝達関数Hの中には時間遅延は表れていないが、分母と分子にそれぞれ遅延時間が表れている。また、上記したような実際のマイクロフォンを使用したシステムにおける問題点がここでも存在するため、コントローラCの制御効率は低く留まってしまう。よって、次に示すように、遅延時間補償された騒音推定部Eが望まれる。   Although no time delay appears in the transfer function H, delay times appear in the denominator and numerator, respectively. Further, since there is a problem in the system using the actual microphone as described above, the control efficiency of the controller C remains low. Therefore, a noise estimation unit E that is compensated for delay time is desired as shown below.

最後に、遅延時間補償のある仮想的マイクロフォンを備えるシステムについて制御効率を解析する。   Finally, control efficiency is analyzed for a system with a virtual microphone with delay time compensation.

図26は、仮想的なマイクロフォン(騒音推定部32c)であって遅延時間補償があるものを備える能動振動制御装置を示す。ここでの制御指令信号Vは、数33により表される。

Figure 2007219262
FIG. 26 shows an active vibration control apparatus including a virtual microphone (noise estimation unit 32c) having a delay time compensation. The control command signal V here is expressed by Equation 33.
Figure 2007219262

推定される車室内騒音[SPL]は、数34により表される。

Figure 2007219262
The estimated vehicle interior noise [SPL] is expressed by Equation 34.
Figure 2007219262

数25のSPL(s)は、数35となる。

Figure 2007219262
SPL (s) in Expression 25 is expressed in Expression 35.
Figure 2007219262

[H](s)の分母は、H(s)及びH(s)より小さい遅延時間を有する。遅延時間補償τと遅延時間Taが等しい場合、時間遅延は完全に削除され、SPL(s)は数36となる。

Figure 2007219262
The denominator of [H] (s) has a delay time smaller than H (s) and H (s). When the delay time compensation τ and the delay time Ta are equal, the time delay is completely deleted, and SPL (s) is expressed by Equation 36.
Figure 2007219262

<第3の実施の形態の変形例>
第3の実施の形態では、1つのアクチュエータ8と、1つの外乱と、1つの車室内騒音SPLとを備えるシステムにおける騒音推定部Eについて説明したが、図22及び図23に示した能動振動制御装置は、2以上のアクチュエータ、外乱、車室内騒音SPLを備えるシステムに対して適用可能である。つまり、2以上のアクチュエータ及び加速度センサをフロアパネル上に配置し、車室内の複数箇所における騒音について騒音推定を実施するシステムについても容易に拡張することができる。
<Modification of Third Embodiment>
In the third embodiment, the noise estimation unit E in the system including one actuator 8, one disturbance, and one vehicle interior noise SPL has been described. However, the active vibration control illustrated in FIGS. The apparatus is applicable to a system including two or more actuators, a disturbance, and a vehicle interior noise SPL. That is, it is possible to easily expand a system in which two or more actuators and acceleration sensors are arranged on the floor panel and noise estimation is performed for noise at a plurality of locations in the vehicle interior.

具体的には、Na個のアクチュエータとNd個の外乱(加速度センサ)を用いた場合の、遅延時間補償されたj番目の車室内騒音[SPL]jは、数37により与えられる。

Figure 2007219262
Specifically, when using Na actuators and Nd disturbances (acceleration sensors), the delay time-compensated j-th vehicle interior noise [SPL] j is given by Equation 37.
Figure 2007219262

ここで、Viはi番目のアクチュエータへの制御指令信号を示し、[G]a,ijはi番目のアクチュエータからj番目の車室内騒音SPLまでの遅延時間を含まない伝達関数を示し、e-s(Ta,ij-τij)はi番目のアクチュエータからj番目の車室内騒音SPLまでの遅延時間であって、遅延時間補償τijにより補償されたものを示す。[G]d,kjはk番目の外乱からj番目の車室内騒音SPLまでの遅延時間を含まない伝達関数を示し、e-s(Td,kj-τij)はk番目の外乱からj番目の車室内騒音SPLまでの遅延時間であって、遅延時間補償τijにより補償されたものを示す。Dk(s)は、k番目の外乱信号、つまり加速度センサから出力される信号を示す。なお、Nd個の外乱を含むシステムにおける車室内騒音SPLijを求める騒音推定部Eの構成を図29に示す。 Here, Vi represents a control command signal to the i-th actuator, [G] a, ij represents a transfer function not including a delay time from the i-th actuator to the j-th vehicle interior noise SPL, and e −. s (Ta, ij-τij) represents a delay time from the i-th actuator to the j-th vehicle interior noise SPL, which is compensated by the delay time compensation τij. [G] d, kj represents a transfer function that does not include a delay time from the kth disturbance to the jth vehicle interior noise SPL, and e −s (Td, kj−τij) represents the jth disturbance from the kth disturbance. This is the delay time to the vehicle interior noise SPL and compensated by the delay time compensation τij. Dk (s) represents a kth disturbance signal, that is, a signal output from the acceleration sensor. FIG. 29 shows the configuration of the noise estimation unit E for obtaining the vehicle interior noise SPLij in a system including Nd disturbances.

<遅延時間補償の効率の実験例>
遅延時間補償された騒音推定部Eとこれを含むコントローラユニットCについて説明する。
<Experimental example of delay compensation efficiency>
The noise estimation unit E compensated for the delay time and the controller unit C including the same will be described.

車両の運転者の右耳の位置を騒音制御対象とする。1つのピエゾアクチュエータ8をフロアパネル6のタンクパネル19上に載置し、外乱を制御対象位置にリンクした箇所を振動するシェーカーにより形成する。外乱信号はガウシアン分布を形成している。コントローラCは、110Hzを中心周波数とする周波数帯の騒音を除去することを目的としている。実験結果を図27及び図28に示す。なお、図27及び図28に示す実験結果は、同じシステム(プラント5)について得られたものであるが、騒音推定部E及びコントローラCは互いに異なる。   The position of the right ear of the driver of the vehicle is the noise control target. One piezo actuator 8 is placed on the tank panel 19 of the floor panel 6, and is formed by a shaker that vibrates a location where the disturbance is linked to the control target position. The disturbance signal forms a Gaussian distribution. The controller C is intended to remove noise in a frequency band having a center frequency of 110 Hz. The experimental results are shown in FIGS. 27 and 28 are obtained for the same system (plant 5), but the noise estimation unit E and the controller C are different from each other.

先ず、図27の実験において、騒音推定部Eは総ての時間遅延であって遅延時間補償されていないものを有する。騒音除去効率は最大で15dBに到達し、騒音除去が有効に機能している周波数帯は、およそ107〜120Hzである。   First, in the experiment of FIG. 27, the noise estimation unit E has all the time delays and no delay time compensation. Noise removal efficiency reaches a maximum of 15 dB, and the frequency band in which noise removal is functioning effectively is approximately 107 to 120 Hz.

次に、図28に示す実験において、騒音推定部[E]は、遅延時間Taが遅延時間補償τに等しい場合の遅延時間を有する。騒音除去の最大値は18dBであり、騒音除去が有効に機能している周波数帯は、およそ80〜130Hzである。   Next, in the experiment shown in FIG. 28, the noise estimation unit [E] has a delay time when the delay time Ta is equal to the delay time compensation τ. The maximum value of noise removal is 18 dB, and the frequency band in which noise removal is functioning effectively is approximately 80 to 130 Hz.

このような実験例により、遅延時間補償の利点が示された。つまり、遅延時間補償を加えることにより、騒音除去が有効に機能する周波数帯を広げることができ、且つ、騒音除去効率を高めることができる。   Such an experimental example showed the advantage of delay time compensation. In other words, by adding the delay time compensation, it is possible to widen the frequency band in which noise removal functions effectively, and to improve noise removal efficiency.

以上説明したように、本発明の第3の実施の形態によれば、以下の作用効果が得られる。   As described above, according to the third embodiment of the present invention, the following operational effects can be obtained.

騒音推定部32は、センサが取り付けられた位置における振動から車室内騒音までの第1の遅延時間Tdと、制御指令信号から車室内騒音までの第2の遅延時間Taとを考慮して、第1の遅延時間後及び第2の遅延時間後の車室内騒音を推定する。これにより、より現実に近い車室内騒音の推定を行うことが出来る。   The noise estimation unit 32 takes into account the first delay time Td from the vibration at the position where the sensor is attached to the vehicle interior noise, and the second delay time Ta from the control command signal to the vehicle interior noise. The vehicle interior noise after the first delay time and after the second delay time is estimated. As a result, the vehicle interior noise can be estimated more realistically.

第2の遅延時間Taは、制御指令信号から車室内騒音までの遅延時間のうち最も短い最短遅延時間である。これにより、遅延時間を補償して車体に実際に生じる遅延時間よりも短い遅延時間によって車室内騒音SPLを推定することにより、より良い制御効率を実現することができる。   The second delay time Ta is the shortest delay time among the delay times from the control command signal to the vehicle interior noise. Thus, better control efficiency can be realized by compensating the delay time and estimating the vehicle interior noise SPL based on a delay time shorter than the delay time actually generated in the vehicle body.

第1の遅延時間Tdは、第2の遅延時間Taを補正することにより求められる。遅延時間の計算負担が軽減され、容易なアルゴリズムにより遅延時間を計算することができる。   The first delay time Td is obtained by correcting the second delay time Ta. The calculation load of the delay time is reduced, and the delay time can be calculated by an easy algorithm.

アクチュエータが複数存在する場合、第2の遅延時間Taは各アクチュエータについて個別に設定される。これにより、複数の振動伝達経路のそれぞれについての遅延時間を設定することができる。   When there are a plurality of actuators, the second delay time Ta is individually set for each actuator. Thereby, the delay time about each of a plurality of vibration transmission paths can be set.

(第4の実施の形態)
上記第1乃至第3の実施の形態に係わる能動振動騒音制御装置によれば、ロードノイズをその他の音(ラジオ、乗員の声、風、雨など)から分離することができるようになるが、騒音推定部Eにより推定された騒音[SPL]と実際の騒音SPLとの間での不一致が問題となる場合がある。
(Fourth embodiment)
According to the active vibration noise control apparatus according to the first to third embodiments, road noise can be separated from other sounds (radio, passenger voice, wind, rain, etc.). Inconsistency between the noise [SPL] estimated by the noise estimation unit E and the actual noise SPL may be a problem.

そこで、第4の実施の形態では、実際にマイクロフォンを用いて騒音SPLを測定して騒音推定部Eに含まれる伝達モデルを更新する技術について説明する。ここで問題となるのは、どのタイミングで伝達モデルの更新を行うかである。   Therefore, in the fourth embodiment, a technique for actually measuring the noise SPL using a microphone and updating the transmission model included in the noise estimation unit E will be described. The problem here is at what timing the transmission model is updated.

騒音推定部Eが備える伝達モデルが完全に実際の伝達モデルに一致している場合、推定された騒音[SPL]は、実際の騒音SPLと正確に一致する。しかし、以下に示す2つの理由から、両者は必ずしも一致しない。   When the transmission model included in the noise estimation unit E completely matches the actual transmission model, the estimated noise [SPL] exactly matches the actual noise SPL. However, for the following two reasons, they do not necessarily match.

先ず、騒音推定部が有する2つのタイプのエラーがある。路面の振動には直接アクセスすることができず、フロアパネルの振動を用いて近似されている。更に、フロアパネルの振動測定点の数は限られているため、測定されるデータに路面の振動についての総ての情報が含まれていない。以上が1点目のエラーである。そして、2点目のエラーは、騒音推定部Eで使用している伝達モデルは、実際のシステムの対応物に過ぎず、伝達モードの数、伝達モデルの正確性には限界がある点である。   First, there are two types of errors that the noise estimator has. Road surface vibrations are not directly accessible and are approximated using floor panel vibrations. Furthermore, since the number of vibration measurement points on the floor panel is limited, the measured data does not include all information about road vibration. This is the first error. The second error is that the transmission model used in the noise estimation unit E is only a counterpart of the actual system, and there are limits to the number of transmission modes and the accuracy of the transmission model. .

そして、実際の騒音SPLは、路面の振動となんら関係のないその他の音(ラジオ、乗員の声、風、雨など)と混合されているため、推定される騒音と実際の騒音との間の関連性は低い。   Since the actual noise SPL is mixed with other sounds (radio, occupant's voice, wind, rain, etc.) that have nothing to do with road vibration, there is a difference between the estimated noise and the actual noise. Relevance is low.

そこで、騒音推定部E内の伝達モデルを更新する必要が生じる。ここでは、コヒーレンスを用いた更新を説明する。この更新は2つの段階に分類できる。第1の段階は、コヒーレンス関数を用いた伝達モデルの更新を行うことを決定する段階である。第2の段階は、騒音推定部Eの出力と実際のマイクロフォンの出力とを用いて、実際に伝達モデルの更新を行う段階である。   Therefore, it is necessary to update the transmission model in the noise estimation unit E. Here, updating using coherence will be described. This update can be classified into two stages. The first stage is a stage in which it is determined to update the transfer model using the coherence function. The second stage is a stage where the transmission model is actually updated using the output of the noise estimation unit E and the output of the actual microphone.

更新を決定する第1の段階に関して、従来では、更新時を外的な出来事に基づいて定めていた。例えば、特開2000−120768号公報では、ドアが閉まる前を更新時としていた。そして、予め定めた外的な事項が生じたときに限り、伝達モデルの更新を実施する。   With respect to the first stage of determining the update, conventionally, the update time is determined based on an external event. For example, in Japanese Patent Application Laid-Open No. 2000-120768, the time before the door is closed is set as the update time. The transmission model is updated only when a predetermined external matter occurs.

これに対して本実施の形態では、システムが自動的に更新を行う時及び更新を行わない時を自動的に決定する。具体的には、騒音推定部が複数のポイントで推定した車室内騒音[SPL]と実際にマイクロフォンを用いて測定した車室内騒音SPLとを比較して更新時を決定する。   On the other hand, in this embodiment, the system automatically determines when to update and when not to update. Specifically, the vehicle interior noise [SPL] estimated by the noise estimation unit at a plurality of points is compared with the vehicle interior noise SPL actually measured using a microphone to determine the update time.

制御対象となる周波数帯における推定音圧[SPL]と実測音圧SPLとの間のコヒーレンス関数Cs(f)を計算する。コヒーレンス関数が大きい場合、例えば0.9以上、理想的には1.0近くである場合、推定音圧[SPL]と実測音圧SPLとが線形的に関連して、ロードノイズに混合されるその他の音(風、乗員の声等)がないことを意味している。よって、この場合、騒音推定部の伝達モデルの更新を有効に実施することができる。   A coherence function Cs (f) between the estimated sound pressure [SPL] and the measured sound pressure SPL in the frequency band to be controlled is calculated. When the coherence function is large, for example, 0.9 or more, ideally close to 1.0, the estimated sound pressure [SPL] and the actually measured sound pressure SPL are linearly related and mixed with road noise. It means that there is no other sound (wind, occupant's voice, etc.). Therefore, in this case, the transmission model of the noise estimation unit can be effectively updated.

騒音推定部の伝達モデルGと実際のシステムとの不一致は、コヒーレンス関数Cs(f)に影響を与えず、実測音圧SPLに含まれ得るロードノイズ以外の音だけがコヒーレンス関数Cs(f)に影響を与える。例えば、コヒーレンス関数Cs(f)が1.0よりも低くなり、ほとんど0に近くなってしまう。推定音圧[SPL]と実測音圧SPLが正確に一致していないことは、コヒーレンス関数Cs(f)は影響を受けない。   The discrepancy between the transfer model G of the noise estimation unit and the actual system does not affect the coherence function Cs (f), and only sounds other than road noise that can be included in the measured sound pressure SPL are included in the coherence function Cs (f). Influence. For example, the coherence function Cs (f) becomes lower than 1.0 and is almost close to 0. The fact that the estimated sound pressure [SPL] and the measured sound pressure SPL do not exactly match does not affect the coherence function Cs (f).

推定音圧[SPL]と実測音圧SPL間のコヒーレンス関数Cs(f)を計算する利点の1つは、フロアパネル上の総てのポイントにおける振動αiと実測音圧SPL間のコヒーレンス関数を測定する場合に比べて、計算に対する負担及び記憶容量に対する負担が極端に軽減されることである。   One advantage of calculating the coherence function Cs (f) between the estimated sound pressure [SPL] and the measured sound pressure SPL is to measure the coherence function between the vibration αi and the measured sound pressure SPL at all points on the floor panel. Compared with the case where it does, the burden with respect to a calculation and the burden with respect to a storage capacity are reduced extremely.

図30は、n個の第1の伝達モデル14−1〜14−nと、m個の第2の伝達モデル15−1〜15−mとを備える騒音推定部32を示す。n個の加速度センサがα1〜αnの振動を測定し、m個のアクチュエータ(への制御指令信号V1〜Vm)が能動騒音制御に使用される。   FIG. 30 illustrates a noise estimation unit 32 including n first transmission models 14-1 to 14-n and m second transmission models 15-1 to 15-m. n acceleration sensors measure vibrations α1 to αn, and m actuators (control command signals V1 to Vm) are used for active noise control.

推定音圧[SPL]が使用されるとき1つのコヒーレンス関数Cs(f)を使用できる。

Figure 2007219262
One coherence function Cs (f) can be used when the estimated sound pressure [SPL] is used.
Figure 2007219262

ここで、xは入力信号(推定音圧[SPL])を示し、yは出力信号(実測音圧SPL)を示す。Pxx(f)は信号x(t)の自己スペクトル、Pyy(f)は信号y(t)の自己スペクトル、Pxy(f)は周波数fにおける信号x(t)と信号y(t)の交差スペクトルを示す。これら総てのスペクトルは例えばWelch法により計算できる。信号は、例えば所定のサンプリング周波数にて1秒間測定される。   Here, x represents an input signal (estimated sound pressure [SPL]), and y represents an output signal (measured sound pressure SPL). Pxx (f) is the self spectrum of the signal x (t), Pyy (f) is the self spectrum of the signal y (t), and Pxy (f) is the cross spectrum of the signal x (t) and the signal y (t) at the frequency f. Indicates. All these spectra can be calculated, for example, by the Welch method. The signal is measured for 1 second at a predetermined sampling frequency, for example.

複数の入力信号が使用されるとき、数39に示すように、複数のコヒーレンス関数Cmが必要である。計算に対する負荷及び必要とされる記憶容量は、1つのコヒーレンス関数の場合に比べて非常に大きくなる。

Figure 2007219262
When a plurality of input signals are used, a plurality of coherence functions Cm are required as shown in Equation 39. The computational burden and the required storage capacity are very large compared to a single coherence function.
Figure 2007219262

ここで、Sxy(f)とSxx(f)は、数40に示すとおりである。

Figure 2007219262
Here, Sxy (f) and Sxx (f) are as shown in Equation 40.
Figure 2007219262

kは入力信号の数を示し、制御が行われずn個の加速度センサが使用されているだけである時はnであり、制御が行われn個の加速度センサとm個のアクチュエータが使用されている時はn+mである。複数のコヒーレンス関数Cmを計算しなくても、1つのコヒーレンス関数Csについて計算すれば十分であることは明らかである。よって、実際にマイクロフォンで測定された実測音圧SPLと仮想的マイクロフォン(騒音推定部)により推定された推定音圧[SPL]とから1つのコヒーレンス関数Cs(f)を計算する。   k indicates the number of input signals, and is n when control is not performed and only n acceleration sensors are used, and control is performed and n acceleration sensors and m actuators are used. N + m. Obviously, it is sufficient to calculate one coherence function Cs without calculating a plurality of coherence functions Cm. Therefore, one coherence function Cs (f) is calculated from the actually measured sound pressure SPL actually measured by the microphone and the estimated sound pressure [SPL] estimated by the virtual microphone (noise estimation unit).

一旦、更新の決定が成されると、更新そのものは極めてありふれた方法で行われる。しかしながら、制御が行われている時であっても更新は行うことができる。なぜなら、推定音圧[SPL]の一部はアクチュエータによる制御振動に相当し(これをSPLaとする)、図30に示すようにこれは仮想的なマイクロフォンの中に含まれている。推定音圧[SPL]は数41で表される。

Figure 2007219262
Once the update decision is made, the update itself is done in a very common way. However, updating can be performed even when control is being performed. This is because part of the estimated sound pressure [SPL] corresponds to the control vibration by the actuator (this is referred to as SPLa), and this is included in the virtual microphone as shown in FIG. The estimated sound pressure [SPL] is expressed by Equation 41.
Figure 2007219262

制御を行わない時、アクチュエータへの制御指令信号V1〜Vmはすべてゼロである。よって、SPLaの項は無くなる。制御が行われていない時であっても更新は行うことができる。   When the control is not performed, the control command signals V1 to Vm to the actuator are all zero. Therefore, the SPLa term disappears. Updates can be made even when control is not being performed.

更新モデルK(s)について説明する。更新モデルK(s)は、推定騒音[SPL]のうちフロアパネルの振動に相当するSPLαに対して適用される。これは以下の2つの理由からである。   The update model K (s) will be described. The update model K (s) is applied to SPLα corresponding to floor panel vibration in the estimated noise [SPL]. This is for the following two reasons.

第1に、制御が行われていない場合、信号SPLaはゼロにほぼ等しい。よって、信号SPLaに相当する部分の更新は適切ではない。   First, when no control is taking place, signal SPLa is approximately equal to zero. Therefore, the update of the portion corresponding to the signal SPLa is not appropriate.

第2に、推定騒音[SPL]のうち信頼できない結果を表しているのはSPLαである。この部分はフロアパネルの幾つかのポイントにおける振動αiを使用し、車室内でロードノイズとなる路面の振動をあらわしている。更に、第1の伝達モデルGαiは、システムの特定による実験データから求められる。よって、第1の伝達モデルGαiと実際のシステムとの間には避けることができない違いが存在する。   Second, SPLα represents an unreliable result among the estimated noise [SPL]. This part uses vibration αi at several points on the floor panel, and represents road vibration that becomes road noise in the passenger compartment. Further, the first transfer model Gαi is obtained from experimental data obtained by specifying the system. Thus, there is an unavoidable difference between the first transfer model Gαi and the actual system.

この2つの理由から、更新モデルK(s)は信号SPLαの部分に適用される。よって、推定音圧[SPL]は数42になる。

Figure 2007219262
For these two reasons, the update model K (s) is applied to the part of the signal SPLα. Therefore, the estimated sound pressure [SPL] is expressed by Equation 42.
Figure 2007219262

図31は、更新モデル64を取り入れたモデル更新部61の構成を示すブロック図である。   FIG. 31 is a block diagram illustrating a configuration of the model update unit 61 that incorporates the update model 64.

次に、更新モデルKを計算する演算部63について説明する。判定部62が更新を決定した時、数43に示すように、実測音圧SPLを用いて更新モデルKを計算する。

Figure 2007219262
Next, the calculation unit 63 that calculates the update model K will be described. When the determination unit 62 determines to update, the update model K is calculated using the measured sound pressure SPL as shown in Equation 43.
Figure 2007219262

これを変形すると、数44となる。

Figure 2007219262
When this is transformed, the following equation 44 is obtained.
Figure 2007219262

更新モデルKは、実測音圧SPLと騒音推定部の路面振動に関する部分SPLαとの間の伝達関数としてみることができる。更新モデルK(s)は、単入力−単出力の周知の特定化技術により計算可能である。   The update model K can be viewed as a transfer function between the measured sound pressure SPL and the portion SPLα related to the road surface vibration of the noise estimation unit. The update model K (s) can be calculated by a well-known specification technique of single input-single output.

図31に示すように、第3の実施の形態では、モデル更新部61が、第1の伝達モデルGαの後段に付加されている。モデル更新部61は、判定部62と、演算部63とを備える。判定部62は、推定音圧[SPL]と実測音圧SPL間のコヒーレンス関数を計算して更新モデルKの更新を行うか否かを判断する。演算部63は、判定部62が更新を決定した場合、更新モデルKを演算及び更新する。   As shown in FIG. 31, in the third embodiment, a model update unit 61 is added to the subsequent stage of the first transfer model Gα. The model update unit 61 includes a determination unit 62 and a calculation unit 63. The determination unit 62 determines whether or not to update the update model K by calculating a coherence function between the estimated sound pressure [SPL] and the actually measured sound pressure SPL. The calculation unit 63 calculates and updates the update model K when the determination unit 62 determines to update.

更新のアルゴリズムは以下のとおりである。(1)先ず、サンプリング時間n毎にマイクロフォンを用いて実測音圧SPLを測定する。(2)T秒毎に実測音圧SPLと推定音圧[SPL]間のコヒーレンス関数Cを計算する。(3)所定の閾値Ctに対してコヒーレンス関数Cを比較して更新の是非を決定する。CがCtよりも大きければ、更新を決定して次の段階に進み、CがCt以下であれば更新は行わない。(4)制御が行われている場合、実測音圧SPL、SPLα及びSPLaとを用いて更新モデルK(s)を更新する。なお、更新モデルK(s)はゲイン1へ初期化される。   The update algorithm is as follows. (1) First, the measured sound pressure SPL is measured using a microphone at every sampling time n. (2) The coherence function C between the measured sound pressure SPL and the estimated sound pressure [SPL] is calculated every T seconds. (3) The coherence function C is compared with a predetermined threshold Ct to determine whether to update. If C is larger than Ct, update is determined and the process proceeds to the next stage. If C is equal to or less than Ct, no update is performed. (4) When control is performed, the update model K (s) is updated using the measured sound pressures SPL, SPLα, and SPLa. Note that the update model K (s) is initialized to gain 1.

上記アルゴリズムの具体例を示す。ここでは、制御は行われていない。つまり、SPLaはゼロである。加速度センサの数は1であり、SPLの測定ポイントは1つである。よって、システムは、単入力−単出力(Single-Input Single-Output)モデルで表される。   A specific example of the above algorithm will be shown. Here, no control is performed. That is, SPLa is zero. The number of acceleration sensors is 1, and there is one SPL measurement point. Therefore, the system is represented by a single-input single-output model.

1組のデータを用いて、第1の伝達モデルGαが計算され、騒音推定部に装備される。第1の伝達モデルGαへ加速度センサからの出力信号αが入力されて推定音圧[SPL]を出力する。   A first transfer model Gα is calculated using one set of data and is installed in the noise estimation unit. The output signal α from the acceleration sensor is input to the first transmission model Gα, and the estimated sound pressure [SPL] is output.

他の1組のデータを用いて、図32に示すように、推定音圧[SPL]と実測音圧SPLを比較する。推定音圧[SPL]と実測音圧SPLが非常に近くてもゲインは正確に一致しない。よって、騒音推定つまり制御効率を改善するために更新を実施する。   As shown in FIG. 32, the estimated sound pressure [SPL] is compared with the actually measured sound pressure SPL using another set of data. Even if the estimated sound pressure [SPL] and the actually measured sound pressure SPL are very close, the gains do not exactly match. Therefore, updating is performed to improve noise estimation, that is, control efficiency.

前述したように、制御対象の周波数帯においてコヒーレンス関数Cが所定の閾値Ctよりも大きい場合、更新モデルKの更新を行う。なお、理論的な閾値Ctは1にほぼ等しい。   As described above, when the coherence function C is larger than the predetermined threshold Ct in the frequency band to be controlled, the update model K is updated. The theoretical threshold Ct is approximately equal to 1.

図33は、コヒーレンス関数Cの計算結果の一例を示す。ここで制御対象の周波数帯は、80〜230Hzである。この制御対象周波数帯においてコヒーレンス関数Cは1に非常に近いため、更新モデルKの更新を行うことができることが分かる。なお、0〜80Hzと250〜1000Hzでのコヒーレンス関数はとても低い。よって、Kの計算は数44により行うことができる。   FIG. 33 shows an example of the calculation result of the coherence function C. Here, the frequency band to be controlled is 80 to 230 Hz. Since the coherence function C is very close to 1 in this control target frequency band, it can be seen that the update model K can be updated. Note that the coherence function at 0 to 80 Hz and 250 to 1000 Hz is very low. Therefore, the calculation of K can be performed by Equation 44.

以上説明したように、本発明の第4の実施の形態によれば、以下の作用効果が得られる。   As described above, according to the fourth embodiment of the present invention, the following operational effects can be obtained.

筐体内騒音を測定するマイクロフォンを配置して、マイクロフォンにより測定された車室内騒音と騒音推定部32により推定された車室内騒音とを用いて第1の伝達モデルを更新する。これにより、騒音推定部Eにより推定された騒音[SPL]と実際の騒音SPLとが一致しないという課題を解決することができる。   A microphone for measuring the noise in the housing is arranged, and the first transmission model is updated using the vehicle interior noise measured by the microphone and the vehicle interior noise estimated by the noise estimation unit 32. Thereby, the problem that the noise [SPL] estimated by the noise estimation unit E does not match the actual noise SPL can be solved.

マイクロフォンにより測定された車室内騒音と騒音推定部32により推定された車室内騒音との間のコヒーレンス関数Cの値が所定値以上であるときに限り、第1の伝達モデルを更新する。更新の必要がある場合に限り更新を行うため、効率的に更新を行うことができる。   The first transfer model is updated only when the value of the coherence function C between the vehicle interior noise measured by the microphone and the vehicle interior noise estimated by the noise estimation unit 32 is greater than or equal to a predetermined value. Since the update is performed only when the update is necessary, the update can be performed efficiently.

予め定めた周波数帯において、マイクロフォンにより測定された車室内騒音SPLと騒音推定部32により推定された車室内騒音[SPL]との間のコヒーレンス関数Cの値が所定値以上であるときに限り、第1の伝達モデルを更新する。対象となるコヒーレンス関数Cの周波数帯を限定することにより、より効率的な更新を行うことができる。   Only when the value of the coherence function C between the vehicle interior noise SPL measured by the microphone and the vehicle interior noise [SPL] estimated by the noise estimation unit 32 is a predetermined value or more in a predetermined frequency band. Update the first transmission model. By limiting the frequency band of the target coherence function C, more efficient updating can be performed.

上記のように、本発明は、4つの実施の形態によって記載したが、この開示の一部をなす論述及び図面はこの発明を限定するものであると理解すべきではない。この開示から当業者には様々な代替実施の形態、実施例及び運用技術が明らかとなろう。即ち、本発明はここでは記載していない様々な実施の形態等を包含するということを理解すべきである。したがって、本発明はこの開示から妥当な特許請求の範囲に係る発明特定事項によってのみ限定されるものである。   As described above, the present invention has been described with reference to the four embodiments. However, it should not be understood that the description and drawings constituting a part of this disclosure limit the present invention. From this disclosure, various alternative embodiments, examples and operational techniques will be apparent to those skilled in the art. That is, it should be understood that the present invention includes various embodiments not described herein. Therefore, the present invention is limited only by the invention specifying matters according to the scope of claims reasonable from this disclosure.

本発明の第1の実施の形態に係わる能動振動騒音制御装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an active vibration noise control apparatus according to a first embodiment of the present invention. 路面の凹凸の影響によるタイヤ及び車体の振動および車室内での騒音の主な伝播経路を示す斜視図である。It is a perspective view which shows the main propagation path | route of the vibration of the tire and vehicle body by the influence of the unevenness | corrugation of a road surface, and the noise in a vehicle interior. 制御系設計のためのブロック線図である。It is a block diagram for control system design. 図1の騒音推定部32の構成を示すブロック線図である。It is a block diagram which shows the structure of the noise estimation part 32 of FIG. 図4の騒音推定部32による騒音推定手順を示すフローチャートである。It is a flowchart which shows the noise estimation procedure by the noise estimation part 32 of FIG. 実際に測定した車室内騒音と騒音推定部32により推定した車室内騒音を比較するグラフである。6 is a graph comparing vehicle interior noise actually measured with vehicle interior noise estimated by the noise estimation unit 32. 図4の騒音推定部32における第3の伝達モデル16の有無による実測音圧と推定音圧とのコヒーレンスの変化を示すグラフである。It is a graph which shows the change of the coherence of the measurement sound pressure by the presence or absence of the 3rd transmission model 16 in the noise estimation part 32 of FIG. 一般的な適応フィルタを用いたフィードフォワード型能動振動制御装置の全体構成を示すブロック図である。It is a block diagram which shows the whole structure of the feedforward type | mold active vibration control apparatus using a general adaptive filter. 図8のシステム全体の信号伝達系を示すブロック図である。It is a block diagram which shows the signal transmission system of the whole system of FIG. 図8の適応フィルタ38で行われる演算手順を示すフローチャートである。It is a flowchart which shows the calculation procedure performed with the adaptive filter 38 of FIG. 第2の実施の形態に係わるフィードフォワード型能動振動制御装置の全体構成を示すブロック図である。It is a block diagram which shows the whole structure of the feedforward type | mold active vibration control apparatus concerning 2nd Embodiment. 図11のコントローラユニット36の演算手順を示すフローチャートである。12 is a flowchart showing a calculation procedure of the controller unit 36 of FIG. 11. 複数個の加速度センサ及びアクチュエータを使用した場合の騒音推定部32の構成を示すブロック図である。It is a block diagram which shows the structure of the noise estimation part 32 at the time of using a some acceleration sensor and actuator. 図14(a)は第3の実施の形態に係わる能動振動騒音制御装置が搭載される車両を示す断面図であり、図14(b)は車両のフロアパネル部分を示す平面図である。FIG. 14A is a cross-sectional view showing a vehicle on which the active vibration and noise control apparatus according to the third embodiment is mounted, and FIG. 14B is a plan view showing a floor panel portion of the vehicle. 図14(b)のフロアパネル上に配置されたメンバを示す平面図である。It is a top view which shows the member arrange | positioned on the floor panel of FIG.14 (b). 車室内にマイクロフォンを配置して車室内騒音を実測するタイプの能動振動騒音制御装置を示すブロック図である。It is a block diagram which shows the active vibration noise control apparatus of the type which arrange | positions a microphone in a vehicle interior and measures vehicle interior noise. フロアパネルの振動から車室内騒音を推測するタイプの能動振動騒音制御装置を示すブロック図である。It is a block diagram which shows the active vibration noise control apparatus of the type which estimates vehicle interior noise from the vibration of a floor panel. 第3の実施の形態に係わる、騒音推定部32を用いたコントローラユニット30の構成を示すブロック図である。It is a block diagram which shows the structure of the controller unit 30 using the noise estimation part 32 concerning 3rd Embodiment. 図18の騒音推定部32の内部構成を示すブロック図である。It is a block diagram which shows the internal structure of the noise estimation part 32 of FIG. 時間遅延を考慮した騒音推定部32bの内部構成を示すブロック図である。It is a block diagram which shows the internal structure of the noise estimation part 32b which considered the time delay. 1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLからなるシステムを示すブロック図である。It is a block diagram which shows the system which consists of one actuator 8, one disturbance (acceleration sensor 7), and one vehicle interior noise SPL. 時間遅延を補償した騒音推定部32cの内部構成を示すブロック図である。It is a block diagram which shows the internal structure of the noise estimation part 32c which compensated the time delay. 1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLからなるシステムを示すブロック図である。It is a block diagram which shows the system which consists of one actuator 8, one disturbance (acceleration sensor 7), and one vehicle interior noise SPL. 実際のマイクロフォンを備える能動振動制御装置を示すブロック図である。It is a block diagram which shows an active vibration control apparatus provided with an actual microphone. 仮想的なマイクロフォン(騒音推定部)であって遅延時間補償のないものを備える能動振動制御装置を示すブロック図である。It is a block diagram which shows an active-vibration control apparatus provided with a virtual microphone (noise estimation part) without a delay time compensation. 仮想的なマイクロフォン(騒音推定部)であって遅延時間補償があるものを備える能動振動制御装置を示すブロック図である。It is a block diagram which shows an active vibration control apparatus provided with what is a virtual microphone (noise estimation part) and has delay time compensation. 遅延時間補償されていない場合の周波数ごとの音圧レベルを示すグラフであり、実線が能動制御を行っている場合、点線が能動制御を行っていない場合を示す。It is a graph which shows the sound pressure level for every frequency when delay time compensation is not carried out, and shows the case where a solid line is performing active control, and the case where a dotted line is not performing active control. 遅延時間補償されている場合の周波数ごとの音圧レベルを示すグラフであり、実線が能動制御を行っている場合、点線が能動制御を行っていない場合を示す。It is a graph which shows the sound pressure level for every frequency when delay time compensation is carried out, and the case where a solid line is performing active control, and the case where a dotted line is not performing active control are shown. Nd個の外乱を含むシステムにおける車室内騒音SPLijを求める騒音推定部Eの構成を示す。The structure of the noise estimation part E which calculates | requires the vehicle interior noise SPLij in the system containing Nd disturbance is shown. n個の第1の伝達モデル14−1〜14−nと、m個の第2の伝達モデル15−1〜15−mとを備える騒音推定部32を示すブロック図である。It is a block diagram which shows the noise estimation part 32 provided with n 1st transmission models 14-1 to 14-n and m 2nd transmission models 15-1 to 15-m. 更新モデル64を取り入れたモデル更新部61の構成を示すブロック図である。It is a block diagram which shows the structure of the model update part 61 which took in the update model 64. FIG. 推定音圧[SPL]と実測音圧SPLを比較するグラフである。It is a graph which compares estimated sound pressure [SPL] and measured sound pressure SPL. コヒーレンス関数Cの計算結果の一例を示すグラフである。6 is a graph illustrating an example of a calculation result of a coherence function C.

符号の説明Explanation of symbols

4…タイヤ
5…車体(筐体)
6…フロアパネル
7、39…加速度センサ(センサ)
8…アクチュエータ
9…車室内空間
10…エンジン
14、14b…第1の伝達モデル
15、15b…第2の伝達モデル
16…第3の伝達モデル
17…キャビンフロア
18…タンクパネル
19…スペアタイヤパネル
20a〜20d…縦メンバ
21a〜21d…横メンバ
30、36…コントローラユニット
31、34、43…コントローラ(信号演算手段)
32、32、32a、32b、32c…騒音推定部(騒音推定手段)
33a、33b、33c…アンプ
35…伝達関数モデル
37…適応則演算部
38…適応フィルタ
40、54…マイクロフォン
41…第1の振動伝播経路
42…第2の振動伝播経路
44、44b…第1の遅延時間モデル
45、45b…第2の遅延時間モデル
61…モデル更新部
62…判定部
63…演算部
64…更新モデル
4 ... Tire 5 ... Car body (housing)
6 ... Floor panel 7, 39 ... Acceleration sensor (sensor)
DESCRIPTION OF SYMBOLS 8 ... Actuator 9 ... Interior space 10 ... Engine 14, 14b ... 1st transmission model 15, 15b ... 2nd transmission model 16 ... 3rd transmission model 17 ... Cabin floor 18 ... Tank panel 19 ... Spare tire panel 20a ˜20d... Vertical member 21a to 21d .. Horizontal member 30, 36... Controller unit 31, 34, 43... Controller (signal calculation means)
32, 32, 32a, 32b, 32c ... noise estimation unit (noise estimation means)
33a, 33b, 33c ... Amplifier 35 ... Transfer function model 37 ... Adaptive law calculation unit 38 ... Adaptive filter 40,54 ... Microphone 41 ... First vibration propagation path 42 ... Second vibration propagation path 44,44b ... First Delay time model 45, 45b ... Second delay time model 61 ... Model update unit 62 ... Determination unit 63 ... Calculation unit 64 ... Update model

Claims (13)

振動が外部から内部へ騒音として伝達する筐体と、
前記筐体の振動を検出するセンサと、
前記筐体に制御振動を生成するアクチュエータと、
筐体内騒音を低減するように、前記アクチュエータへの制御指令信号を生成する信号演算手段と、
前記センサの出力信号及び前記制御指令信号を用いて前記筐体内騒音を推定する騒音推定手段と
を備えることを特徴とする能動振動騒音制御装置。
A housing that transmits vibration from outside to inside as noise,
A sensor for detecting vibration of the housing;
An actuator that generates control vibration in the housing;
Signal calculation means for generating a control command signal to the actuator so as to reduce noise in the housing;
An active vibration noise control device comprising: noise estimation means for estimating the noise in the housing using the output signal of the sensor and the control command signal.
前記騒音推定手段は、前記センサが取り付けられた位置における振動から前記筐体内騒音までの伝達特性を表す第1の伝達モデルと、前記制御指令信号から前記筐体内騒音までの伝達特性を示す第2の伝達モデルとを記憶することを特徴とする請求項1記載の能動振動騒音制御装置。   The noise estimation means includes a first transfer model representing a transfer characteristic from vibration to a noise in the casing at a position where the sensor is attached, and a second transfer characteristic from the control command signal to the noise in the casing. The active vibration and noise control apparatus according to claim 1, wherein a transmission model of the vibration is stored. 前記騒音推定手段は、前記センサの出力信号を前記第1の伝達モデルへ入力して得られる出力値と、前記制御指令信号を前記第2の伝達モデルへ入力して得られる出力値との和を前記筐体内騒音として出力することを特徴とする請求項2記載の能動振動騒音制御装置。   The noise estimation means is a sum of an output value obtained by inputting the output signal of the sensor to the first transfer model and an output value obtained by inputting the control command signal to the second transfer model. The active vibration noise control apparatus according to claim 2, wherein: 前記騒音推定手段は、前記センサの出力信号のうち前記アクチュエータが発生する振動の成分を減算した信号を前記第1の伝達モデルへ入力して得られる出力と、前記制御指令信号を前記第2の伝達モデルへ入力して得られる出力との和を前記筐体内騒音として出力することを特徴とする請求項2記載の能動振動騒音制御装置。   The noise estimation means outputs an output obtained by inputting a signal obtained by subtracting a vibration component generated by the actuator from the output signal of the sensor to the first transmission model, and the control command signal as the second control signal. 3. The active vibration noise control apparatus according to claim 2, wherein a sum of outputs obtained by inputting to the transmission model is output as the in-casing noise. 前記騒音推定手段は、前記制御指令信号から前記センサが取り付けられた位置における振動までの伝達特性を示す第3の伝達モデルを更に備え、
前記センサの出力信号から前記第3の伝達モデルの出力を減算した信号を前記第1の伝達モデルへ入力して得られる出力と、前記制御指令信号を前記第2の伝達モデルへ入力して得られる出力との和を前記筐体内騒音として出力することを特徴とする請求項2記載の能動振動騒音制御装置。
The noise estimation means further includes a third transfer model indicating transfer characteristics from the control command signal to vibration at a position where the sensor is attached,
An output obtained by inputting a signal obtained by subtracting the output of the third transfer model from the output signal of the sensor to the first transfer model, and an input obtained by inputting the control command signal to the second transfer model. The active vibration noise control apparatus according to claim 2, wherein the sum of the output and the output is output as the noise in the housing.
前記筐体は車体であり、前記センサと前記アクチュエータは前記車体のフロアパネル上に配置することを特徴とする請求項1乃至5のいずれか1項に記載の能動振動騒音制御装置。   The active vibration noise control device according to claim 1, wherein the housing is a vehicle body, and the sensor and the actuator are arranged on a floor panel of the vehicle body. 前記騒音推定手段は、前記センサが取り付けられた位置における振動から前記筐体内騒音までの第1の遅延時間と、前記制御指令信号から前記筐体内騒音までの第2の遅延時間とを考慮して、前記第1の遅延時間後及び前記第2の遅延時間後の筐体内騒音を推定することを特徴とする請求項1乃至6いずれか1項に記載の能動振動騒音制御装置。   The noise estimation means takes into account a first delay time from vibration at the position where the sensor is mounted to noise in the housing, and a second delay time from the control command signal to noise in the housing. The active vibration noise control apparatus according to any one of claims 1 to 6, wherein the noise in the housing after the first delay time and after the second delay time is estimated. 前記第2の遅延時間は、前記制御指令信号から前記筐体内騒音までの遅延時間のうち最も短い最短遅延時間であることを特徴とする請求項7記載の能動振動騒音制御装置。   8. The active vibration noise control device according to claim 7, wherein the second delay time is a shortest shortest delay time among delay times from the control command signal to the in-casing noise. 前記第1の遅延時間は、前記第2の遅延時間を補正することにより求められることを特徴とする請求項7又は8記載の能動振動騒音制御装置。   The active vibration noise control apparatus according to claim 7 or 8, wherein the first delay time is obtained by correcting the second delay time. 前記アクチュエータが複数存在する場合、前記第2の遅延時間は各アクチュエータについて個別に設定されることを特徴とする請求項7乃至9いずれか1項に記載の能動振動騒音制御装置。   10. The active vibration noise control device according to claim 7, wherein when there are a plurality of actuators, the second delay time is individually set for each actuator. 11. 前記筐体内騒音を測定するマイクロフォンと、
前記マイクロフォンにより測定された前記筐体内騒音と前記騒音推定手段により推定された前記筐体内騒音とを用いて前記第1の伝達モデルを更新するモデル更新部と
を更に備えることを特徴とする請求項1乃至10いずれか1項に記載の能動振動騒音制御装置。
A microphone for measuring the noise in the housing;
The apparatus further comprises: a model updating unit that updates the first transmission model using the in-casing noise measured by the microphone and the in-casing noise estimated by the noise estimation unit. The active vibration noise control apparatus according to any one of 1 to 10.
前記モデル更新部は、前記マイクロフォンにより測定された前記筐体内騒音と前記騒音推定手段により推定された前記筐体内騒音との間のコヒーレンス関数の値が所定値以上であるときに限り、前記第1の伝達モデルを更新することを特徴とする請求項11記載の能動振動騒音制御装置。   The model update unit is configured so that the first update is performed only when a value of a coherence function between the noise in the casing measured by the microphone and the noise in the casing estimated by the noise estimation unit is equal to or greater than a predetermined value. The active vibration and noise control apparatus according to claim 11, wherein the transmission model is updated. 前記モデル更新部は、予め定めた周波数帯において、前記マイクロフォンにより測定された前記筐体内騒音と前記騒音推定手段により推定された前記筐体内騒音との間のコヒーレンス関数の値が所定値以上であるときに限り、前記第1の伝達モデルを更新することを特徴とする請求項11記載の能動振動騒音制御装置。   The model updating unit has a coherence function value between the in-casing noise measured by the microphone and the in-casing noise estimated by the noise estimating means in a predetermined frequency band equal to or greater than a predetermined value. 12. The active vibration noise control device according to claim 11, wherein the first transmission model is updated only occasionally.
JP2006040881A 2006-02-17 2006-02-17 Active vibration noise control device Expired - Fee Related JP4735319B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006040881A JP4735319B2 (en) 2006-02-17 2006-02-17 Active vibration noise control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006040881A JP4735319B2 (en) 2006-02-17 2006-02-17 Active vibration noise control device

Publications (2)

Publication Number Publication Date
JP2007219262A true JP2007219262A (en) 2007-08-30
JP4735319B2 JP4735319B2 (en) 2011-07-27

Family

ID=38496630

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006040881A Expired - Fee Related JP4735319B2 (en) 2006-02-17 2006-02-17 Active vibration noise control device

Country Status (1)

Country Link
JP (1) JP4735319B2 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101018783B1 (en) 2009-07-24 2011-03-03 한국과학기술원 Apparatus and method for noise control
KR101109035B1 (en) 2009-07-24 2012-01-31 한국과학기술원 Apparatus and method for noise control
KR20150105305A (en) * 2013-01-10 2015-09-16 르노 에스.아.에스. System and corresponding method for controlling the rotation speed of an electric motor of a motor vehicle
WO2015196434A1 (en) * 2014-06-26 2015-12-30 华为技术有限公司 Noise reduction method and apparatus, and mobile terminal
WO2018099890A1 (en) * 2016-11-29 2018-06-07 Brose Fahrzeugteile Gmbh & Co. Kommanditgesellschaft, Bamberg Method for operating an apparatus for suppressing structure-borne noise of an auxiliary assembly of a motor vehicle
CN110024025A (en) * 2016-11-23 2019-07-16 哈曼国际工业有限公司 Dynamic stability control system based on coherence
EP3144928B1 (en) * 2015-09-15 2021-03-24 Harman Becker Automotive Systems GmbH Active road noise control
CN113124052A (en) * 2021-04-16 2021-07-16 中国航空发动机研究院 Method for controlling unbalance vibration of electromagnetic bearing-rotor system and electronic equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08272376A (en) * 1995-04-03 1996-10-18 Nippondenso Co Ltd Method and device for reducing car noise
JPH0972375A (en) * 1995-07-04 1997-03-18 Nippon Soken Inc Active damping device
JP2001295622A (en) * 2000-04-18 2001-10-26 Fuji Xerox Co Ltd Active type muffler
JP2005225348A (en) * 2004-02-12 2005-08-25 Fuji Heavy Ind Ltd In-cabin noise reducing device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08272376A (en) * 1995-04-03 1996-10-18 Nippondenso Co Ltd Method and device for reducing car noise
JPH0972375A (en) * 1995-07-04 1997-03-18 Nippon Soken Inc Active damping device
JP2001295622A (en) * 2000-04-18 2001-10-26 Fuji Xerox Co Ltd Active type muffler
JP2005225348A (en) * 2004-02-12 2005-08-25 Fuji Heavy Ind Ltd In-cabin noise reducing device

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101018783B1 (en) 2009-07-24 2011-03-03 한국과학기술원 Apparatus and method for noise control
KR101109035B1 (en) 2009-07-24 2012-01-31 한국과학기술원 Apparatus and method for noise control
KR20150105305A (en) * 2013-01-10 2015-09-16 르노 에스.아.에스. System and corresponding method for controlling the rotation speed of an electric motor of a motor vehicle
KR102279041B1 (en) * 2013-01-10 2021-07-19 르노 에스.아.에스. System and corresponding method for controlling the rotation speed of an electric motor of a motor vehicle
JP2016512676A (en) * 2013-01-10 2016-04-28 ルノー エス.ア.エス. System and corresponding method for controlling the rotational speed of an electric motor of a motor vehicle
US10089972B2 (en) 2014-06-26 2018-10-02 Huawei Technologies Co., Ltd. Noise reduction method and apparatus, and mobile terminal
WO2015196434A1 (en) * 2014-06-26 2015-12-30 华为技术有限公司 Noise reduction method and apparatus, and mobile terminal
EP3144928B1 (en) * 2015-09-15 2021-03-24 Harman Becker Automotive Systems GmbH Active road noise control
CN110024025A (en) * 2016-11-23 2019-07-16 哈曼国际工业有限公司 Dynamic stability control system based on coherence
JP2020501178A (en) * 2016-11-23 2020-01-16 ハーマン インターナショナル インダストリーズ インコーポレイテッド Coherence-based dynamic stability control system
JP7008701B2 (en) 2016-11-23 2022-01-25 ハーマン インターナショナル インダストリーズ インコーポレイテッド Coherence-based dynamic stability control system
CN110024025B (en) * 2016-11-23 2023-05-23 哈曼国际工业有限公司 Dynamic stability control system based on coherence
WO2018099890A1 (en) * 2016-11-29 2018-06-07 Brose Fahrzeugteile Gmbh & Co. Kommanditgesellschaft, Bamberg Method for operating an apparatus for suppressing structure-borne noise of an auxiliary assembly of a motor vehicle
CN113124052A (en) * 2021-04-16 2021-07-16 中国航空发动机研究院 Method for controlling unbalance vibration of electromagnetic bearing-rotor system and electronic equipment
CN113124052B (en) * 2021-04-16 2022-08-23 中国航空发动机研究院 Method for controlling unbalance vibration of electromagnetic bearing-rotor system and electronic equipment

Also Published As

Publication number Publication date
JP4735319B2 (en) 2011-07-27

Similar Documents

Publication Publication Date Title
JP4735319B2 (en) Active vibration noise control device
US7904212B2 (en) Noise estimating device and noise estimating method
JP4857907B2 (en) Noise control device and noise control method
JP5040163B2 (en) Noise reduction apparatus and method
JP4857928B2 (en) Noise control device and noise control method
JP4765410B2 (en) Active vibration noise control device
JP2006213297A (en) Active noise and vibration control device and method
JP4940792B2 (en) Noise control device and noise control method
JP2007296886A (en) Noise reducing device and method
Oh et al. Development of mass producible ANC system for broad-band road noise
JP2010202073A (en) Apparatus and method for controlling noise
JP2010188871A (en) Active vibration noise control device and transfer characteristic measuring method
JP3028977B2 (en) Active noise control device
JP2009073417A (en) Apparatus and method for controlling noise
US11741926B2 (en) Echo cancelation
JP2010054962A (en) Noise control device and noise control method
JP2010202162A (en) Active vibration and noise control device
Zafeiropoulos et al. State-of-the-art digital road noise cancellation by Harman
JP2010202136A (en) Active vibration and noise control device
Oh et al. Development of an active road noise control system
JP2010202081A (en) Active vibration and noise control device
JP4910419B2 (en) Actuator optimal placement device
JP2009107475A (en) Noise reduction device and its method
JP2827603B2 (en) Active noise control device
JP2009220731A (en) Noise controller for vehicle

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20081222

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110127

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110329

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110411

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140513

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees