JP2007219262A - Active vibrating noise controlling apparatus - Google Patents
Active vibrating noise controlling apparatus Download PDFInfo
- 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
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
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以下の周波数帯域のロードノイズを低減することを目的としている。低周波ロードノイズは車室内音場の音響モードが原因となっていると仮定をした上で、この音響モードの腹にマイクロフォンを設置しフィードバック制御を行う方式である。
しかしながら、上記の騒音制御装置では、マイクロフォンを用いて車室内騒音を測定するため、制御する必要のない車室内音楽や乗員の声などもロードノイズと共に参照信号として用いてしまい、フィードバック制御により低減してしまうという問題点があった。 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
筐体は、外部から加わる力によって振動が生じ、この振動が内部へ騒音として伝達し得るものであればよい。例えば、飛行機、電車、自動車、エレベータなど人間の移動手段がこれに含まれる。特に、自動車は、凹凸が激しい路面にタイヤを接して走行するため、タイヤから伝わる振動の制御が重要な技術課題である。そこで実施の形態では、本発明に係わる能動振動騒音制御装置を自動車に適用した場合を例に取り説明する。つまり、ここでは筐体の一例として車体を挙げて説明する。 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
アクチュエータ8は、車体5の中のフロアパネル6下に固定され、電気的な信号(制御指令信号)に受けて自らが歪むことでフロアパネル6に歪みを生成する。これにより、車体5に制御振動を生成して外乱による振動又は騒音を相殺する。アクチュエータ8として、結晶に電界を加えると電界に比例した歪みが生じるピエゾ効果(圧電効果)を利用したピエゾ電気アクチュエータ(piezo-electric actuator)を使用することができる。アクチュエータ8の両面に電界を印加することによりアクチュエータ8は歪む。
The
コントローラ31及び騒音推定部32は、コントローラユニット30に含まれている。コントローラユニット30は、加速度センサ7で得られた信号に基づいて車室内騒音を低減する制御指令値を演算する。コントローラユニット30は、加速度センサ7で得られた信号を増幅する信号増幅用アンプ33aと、加速度センサ7の出力信号とアクチュエータ8への制御指令信号に基づいて車室内騒音を推定する騒音推定部32と、騒音推定部32で推定された推定騒音に基づいて車室内騒音を低減するアクチュエータ8への制御指令値を演算するコントローラ31とを備える。
The
車室内騒音の原因は代表的なものとして、エンジン10の振動に起因する騒音、走行時に路面の凹凸の影響がタイヤから進入することに起因するロードノイズと呼ばれる騒音、走行時に空気の気流によって発生する風切音、などがある。また、タイヤから侵入するロードノイズは、タイヤ及びフロアパネルの振動を介して、これが車室内に伝達する。つまり、フロアパネル6の振動が車体5の内部(車室内)に騒音として伝達する。本実施形態において、騒音推定部32はフロアパネル6の振動に起因する車室内騒音を主に推定することを目的とするため、フロアパネル6下に加速度センサ7を配置して車室内騒音を推定する。したがって、上記例の車室内騒音のうち、フロアパネル振動に起因するロードノイズを中心に制御を行う。
Typical causes of vehicle interior noise are noise caused by the vibration of the
図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
一般的なタイヤ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
図2中の範囲9は車室内空間の一例として、乗員耳位置を含む球体状の領域を指す。この空間での騒音は、タイヤ4から車体に入力された騒音と、アクチュエータ8が発生する制御振動とが重ね合わされた車室内騒音になる。コントローラ31は車室内領域9での騒音が低減されるような制御指令信号をアクチュエータ8に供給して制御振動を発生させる。路面の凹凸に起因するフロアパネル6の振動とアクチュエータ8が生成する制御振動の両者が互いに打ち消しあうように重ね合わさって車室内領域9での騒音は低減される。
A
図1のコントローラ31は、騒音推定部32で推定された音圧(車室内騒音)を入力として、車室内空間9での騒音を低減するような制御指令信号を演算する機能を備える固定係数のFIRもしくはIIRフィルタである。例えば、コントローラ31は以下の手順に従って設計をすればよい。
The
まず、「アクチュエータ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
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
次に、図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を設計することが出来る。
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,.
図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
具体的に説明すると、第1の伝達モデル14は、加速度センサ7が取り付けられた位置におけるフロアパネル振動が車室内空間9まで伝達するときの伝達特性を表す。第2の伝達モデル15は、アクチュエータ8への入力信号が、アクチュエータ8へ入力されてアクチュエータが生成する制御振動が車室内空間9まで伝達するときの伝達特性を表す。第3の伝達モデル16は、アクチュエータ8への入力信号が、アクチュエータ8へ入力されてアクチュエータが生成する制御振動が加速度センサ7が設置された位置まで伝達するときの伝達特性を表す。
More specifically, the
次に、騒音推定部32によって車室内騒音を推定するアルゴリズムの例を説明する。
Next, an example of an algorithm for estimating the vehicle interior noise by the
先ず、第1のアルゴリズムにおいて、騒音推定部32は、加速度センサ7の出力信号を「そのまま」第1の伝達モデル14へ入力し、制御指令信号を第2の伝達モデル15へ入力する。そして、両モデル14、15から得られる出力値の和を、車室内騒音(推定音圧:SPL)として出力する。
First, in the first algorithm, the
第2のアルゴリズムにおいて、騒音推定部32は、加速度センサ7の出力信号のうちアクチュエータ8が発生する振動の成分を減算し、これを第1の伝達モデル14へ入力する。そして、第1の伝達モデル14から得られる出力と、制御指令信号を第2の伝達モデル15へ入力して得られる出力との和を車室内騒音として出力する。このように、加速度センサ7の出力信号のうちアクチュエータ8が発生する振動の成分を減算することにより、加速度センサ7の出力値のうちアクチュエータ8による制御振動を除去して路面の凹凸によるロードノイズだけを抽出することができ、ロードノイズ(外乱)の測定精度が向上する。
In the second algorithm, the
第2のアルゴリズムを具体的に説明すると、図4に示すように、騒音推定部32は、加速度センサ7の出力信号から第3の伝達モデル16の出力を減算した信号を第1の伝達モデル14へ入力して得られる出力と、制御指令信号を第2の伝達モデル15へ入力して得られる出力との和を車室内騒音として出力する。
The second algorithm will be described specifically. As shown in FIG. 4, the
図5のフローチャートを参照して、図4に示す騒音推定部32による騒音推定手順(第2のアルゴリズム)の詳細を説明する。
Details of the noise estimation procedure (second algorithm) by the
(イ)まず、ステップS101で加速度センサ7の出力信号を騒音推定部32へ入力する。これを信号Cとする。
(A) First, the output signal of the acceleration sensor 7 is input to the
(ロ)ステップS102に進み、アクチュエータ8への制御指令信号を騒音推定部32へ入力する。これを信号Dとする。
(B) Proceeding to step S102, the control command signal for the
(ハ)ステップS103に進み、信号Dに第3の伝達モデル16を乗算する。これを信号Eとする。
(C) Proceeding to step S103, the signal D is multiplied by the
(ニ)ステップ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
(ホ)ステップS106に進み、信号Dに第2の伝達モデル15を乗算する。これを信号Hとする。
(E) Proceeding to step S106, the signal D is multiplied by the
(へ)ステップ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
次に、上記した騒音推定部32を用いて車室内騒音を推定することによる効果を説明する。図6は実際に測定した運転席右耳元位置での音圧と、騒音推定装置を用いて推定した同位置での音圧を比較した図である。なお、運転席右耳元位置は車室内位置9の一例である。ここで、実線42は実際に車室内位置9でマイクロフォンにより測定した騒音(実測音圧)であり、点線43は騒音推定部32による推定音圧である。図6に示すように、推定音圧は実測音圧とほぼ同等の値を示している。よって、マイクロフォンを用いることなく、騒音推定部32によって車室内騒音を精度良く推定できることがわかる。
Next, the effect by estimating vehicle interior noise using the
図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により定義される。
ここで、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
以上説明したように、本発明の第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
また、第1及び第2の伝達モデル14、15を備えることにより、当該第1及び第2の伝達モデル14、15を用いて推定を行うことができ、過渡応答を回避して推定精度が良くなる。この結果、騒音低減制御の効果が高まる。
In addition, by providing the first and
具体的には、第1及び第2の伝達モデル14、15へ入力して得られる出力値との和を車内騒音とすることにより、過渡応答を回避して推定精度が良くなる。
Specifically, by making the sum of the output values obtained by inputting to the first and
第1及び第2の伝達モデル14、15に加えて、第3の伝達モデル16を更に備えることにより、アクチュエータ8の制御指令信号が大きくなった場合にも加速度センサ7の出力信号のうち車外から進入する振動成分(外乱)を正確に計測することが出来るので、車室内音圧の推定精度が高まる。また、加速度センサ7とアクチュエータ8の配置自由度が増す。以上の結果、騒音低減制御の効果が高まる。
In addition to the first and
アクチュエータ8から振動までの第3の伝達モデル16のみを用いるため、より簡便でより精度の良い振動推定部32を形成することが出来る。この結果、騒音低減制御の効果が高まる。
Since only the
筐体として車体を例に取り、加速度センサとアクチュエータを車体のフロアパネル上に配置することにより、車室内騒音を除去することができる。 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
コントローラユニット36は、マイクロフォン54からの信号を増幅するアンプ33aと、マイクロフォン54からの信号に基づいて適応フィルタ38を制御する適応則演算部37と、加速度センサ39の出力信号を増幅するアンプ33cと、適応フィルタ38の出力を増幅するアンプ33bとを備える。
The
図9は図8のシステム全体のブロック線図である。コントローラユニット36は、通常のコンピュータシステムもしくは電子回路で実現されている。第1の振動伝播経路41は加速度センサ39の出力信号から車室内空間9での騒音までの伝達関数を表し、第2の振動伝播経路42はアクチュエータ8への制御指令信号から車室内空間9での騒音までの伝達関数を表す。
FIG. 9 is a block diagram of the entire system of FIG. The
以下、図9に従って図8のシステムにおけるアクチュエータ8への制御指令信号の演算法を説明する。加速度センサ39で得られた信号をx(n)、路面から車体に入力される外乱が車室内空間9に作る騒音をy(n)、ピエゾアクチュエータ8への制御指令信号をu(n)、マイクロフォン54で得られた誤差信号をe(n)とおく。ここで、nはサンプル時刻を表す。また、適応フィルタ38、第1の振動伝播経路41及び第2の振動伝播経路42のモデルをそれぞれ数3、数4、数5とする。
ここで、z-iはZ変換のオペレータを表す。なお、wiがサンプル時刻nに依存しているのは、可変係数フィルタを使用しているからである。適応フィルタ38の出力がアクチュエータ8への入力信号として計算される。このとき、信号y(n)とu(n)はそれぞれ数6及び数7として計算される。
以上の式に基づいて計算されたu(n)がアクチュエータ8へ入力される。
U (n) calculated based on the above formula is input to the
一方、適応則演算部37ではフィルタ係数wi(n)の更新が以下の演算によって行われる。まず、評価関数Jを数8と定義する。
一般的に知られたLMS(Least Mean Square)アルゴリズムでは、この評価関数Jがwi(n)に対して最小になるように適用フィルタを更新していく。wi(n)の更新則は数9で与えられる。ここで、KeとKyは設計パラメータである。
また、r(n-i)は数10で定義する。
以上の更新演算が適応則演算部37で行われる。
The above update calculation is performed by the adaptive
次に、図10のフローチャートを参照して、適応フィルタ38で行われる演算手順を説明する。
Next, the calculation procedure performed by the
(い)先ずステップS201において、加速度センサ39及びマイクロフォン40の出力信号をサンプル時間n毎に取得する。ステップS202に進み、加速度センサ39の出力信号x(n)に対して、数7の計算式にて演算することにより、アクチュエータへの指令信号u(n)を求める。
(Ii) First, in step S201, output signals of the
(ろ)ステップ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
(は)ステップ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
次に、本実施形態に係わるフィードフォワード型能動振動騒音制御装置について説明する。 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
図12を参照して、図11のコントローラユニット36の演算手順を説明する。
With reference to FIG. 12, the calculation procedure of the
(A)先ずステップS301において、車外からタイヤを介して車体へ伝達する振動を計測するために用いる第1の加速度センサ39の出力信号と、車室内騒音を推定するために用いる第2の加速度センサ7の出力信号とをそれぞれ取得する。
(A) First, in step S301, the output signal of the
(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
(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
(D)ステップS306に進み、ステップS202で計算されたアクチュエータ8への制御指令信号u(n)を出力する。ステップS301〜S306をサンプル時間n毎に行うことでフィードフォワード制御が行われ、車室内空間9での騒音が低減される。
(D) Proceeding to step S306, the control command signal u (n) to the
このように、騒音推定部32により車室内空間9における騒音を推定することにより、図8に示した従来型と比較してマイクロフォンを装備する必要が無く、フィードフォワード型能動振動制御装置を実現することができる。従って、第1の実施の形態と同様に、車室内空間9にマイクロフォンを設置しないので、乗員の声や音楽を低減することなく、車外から進入する振動に起因した車室内騒音のみを低減する制御を行うことが出来る。また、マイクを用いて制御する場合に比べ、加速度センサ−アクチュエータ間距離が短く無駄時間が少ないため、制御効果が向上する。
Thus, by estimating the noise in the vehicle
なお、説明を簡便にするため、加速度センサの数を第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
(第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
図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
図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
これに対して、図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
コントローラ31として、例えば、一般的なH∞制御を用いることができる。アクチュエータの制御指令信号Vと車室内音圧SPL間の開ループにおける伝達モデルは、一般的な伝達関数の形式で求めることができる。また、複数の制御指令信号V及び複数の車室内音圧SPLを含むシステムにおいても容易に図16及び図17の装置を適用可能であることは前述したとおりである。開ループからなるシステムは、数11で表される。ここで、sはラプラス変換における変数である。
コントローラ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∞を求めている。
ここで、Gsは、W1とW2で重み付けされたシステムの伝達関数であって、整形関数(the shaping function)と呼ぶ。
このような条件式を満たす図16に示すコントローラCは、数14で表される。
ここで、W1とW2は、通常、制御対象とする振動及び騒音の周波数帯を選択することができるフィルターの伝達関数である。[Ms]とNsは数15に示すように、標準化されたGsのコプライム因子である。なお、パラメータεは、安定許容範囲と呼ばれ、通常0.2〜0.3が適当である。
次に、図17の騒音推定部32について説明する。騒音推定部Eは、車室内騒音SPLの伝達におけるマイクロフォン40に置き換えられる。騒音推定部Eは、実際には、車室以内騒音[SPL]の推定を行うため、騒音推定部Eが使用する伝達モデルはできる限り正確であることが望まれる。以後、推定された車室以内騒音(音圧)を[SPL]と表記し、マイクロフォンで実際に測定された車室内騒音をSPLと表記する。
Next, the
図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
図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
ここで、第2の伝達モデルGaは、アクチュエータ8が配置されたフロアパネル位置から車室内騒音を推定するポイントまでの伝達特性を示すものであっても構わない。
Here, the second transmission model Ga may indicate a transmission characteristic from the floor panel position where the
第1の伝達モデルGd及び第2の伝達モデルGaは、振動が伝達するプラント5を特定することにより求められ、実際のプラント5の物理的な総ての特徴を含有していることが望ましい。よって、第1及び第2の伝達モデルには、それぞれフロアパネルから車室内騒音まで騒音及び振動が波として伝達するために必要な時間(遅延時間)が含まれていることが望ましい。
The first transmission model Gd and the second transmission model Ga are obtained by specifying the
騒音推定部32aは、図19に示すように、アクチュエータへの制御指令信号Vと、外乱を示す信号αが入力される。制御指令信号Vは、結果的には、アクチュエータ8が設置されている位置におけるフロアパネル6の振動に相当する。また、外乱を示す信号αは、加速度センサ7により測定されたフロアパネルの振動に相当する。
As shown in FIG. 19, the
騒音推定部32aが十分な精度によって車室内騒音[SPL]を推定する場合、図17に示す能動振動制御装置は、図16に示すマイクロフォンで車室内騒音SPLを実測するタイプの能動振動制御装置に正確に一致する。図17のシステムの利点は、車室内に配置するマイクロフォンが不要となる点であり、ロードノイズと共に、ラジオの音、乗員の声、風の音なども制御振動によりかき消されてしまう危険性が無くなる。
When the
しかしながら、騒音推定部32aは実際のプラント5が有する遅延時間を含んでいるため、コントローラの高い制御効率は期待できない。そこで、第3の実施の形態では、より良い制御効率、つまり騒音除去効率を実現するため、新たな騒音推定部Eの構造を提案する。即ち、遅延時間を補償してプラント5に実際に生じる遅延時間よりも短い遅延時間によって車室内騒音SPLを推定することにより、より良い制御効率を実現する。
However, since the
時間遅延を含まないプラントは、遅延時間を含むプラントにくらべてより単純な構造を有しており、制御効率が高く、得られる制御ゲインも大きい。これに対して、遅延時間を含むプラントの場合、騒音制御が難しく、制御ゲインを高くすることができないため、制御効率は通常低くなってしまう。そこで、プラントの異なる遅延時間を計算及び評価し、これを騒音推定部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
遅延時間補償τを用いて遅延時間Ta及びTdを短縮(補償)する方法を説明する。騒音推定部[E]は、数16に示すようなGa(s)及びGd(s)とを備える。
数16において、第1の伝達モデル[G]a(s)及び第2の伝達モデル[G]d(s)は、それぞれGa(s)及びGd(s)の遅延時間の項を除去したものに相当する。図21では、1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLとを含むシステムを示している。
In
騒音推定部[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
騒音推定部[E]を用いてコントローラユニットCを演算する。コントローラユニットCの演算は、前述したとおりであるが、数11の代りに、数17を用いる。
図22は、遅延時間を補償した騒音推定部32cの構成を示す。外乱に関する信号は先ず第1の伝達モデル14bに入力され、その後、τにより遅延時間補償された遅延時間モデル44bに入力される。一方、アクチュエータへの制御指令信号Vは第2の伝達モデル15bに入力され、その後、τにより遅延時間補償された遅延時間モデル45bへ入力される。遅延時間モデル44b、45bの出力信号を重ね合わせることにより車室内騒音[SPL]を推定することができる。
FIG. 22 shows a configuration of the
図23は、1つのアクチュエータ8と1つの外乱(加速度センサ7)と1つの車室内騒音SPLからなるシステムを示し、伝達経路上には、遅延時間補償された遅延時間モデルだけを示している。
FIG. 23 shows a system composed of one
遅延時間補償された遅延時間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により表すことができる。
第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により得られる。
伝達関数G(s)は遅延時間Tを含むため、理想的な伝達関数N(s)/D(s)は、N(s)の中に不安定なゼロを含み、結果的に安定したゼロも含む。N(s)は数20に示すように分解される。
[N]は安定零点を有するN(s)の一部を示し、N'は不安定零点を含むN(s)の他の一部を示す。よって、数19は、数21になる。
数21は数22に変形することができる。
[N](s)は安定化された[N](-s)であり、例えば、[N](-s)の根は、[N](s)の根に相対する。
[G]は、分子が安定な伝達関数G(s)であり、全域通過伝達関数[N](-s)/[N](s)は、遅延時間項の近似で用いられる周知のパデ近似(Pade approximation)を表す。遅延時間項を数24に示す。
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により得られる。
コントローラCがアクチュエータを駆動するために出力する制御指令信号Vは、数26により得られる。
数25と数26の計算式を組み合わせると、SPLは数27になる。
更に変形すると、数28となる。
上記の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となる。
ここでSPLは数30で表される。
よって、制御指令信号Vは数31となる。
そして、数15のSPL(s)は数32になる。
伝達関数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により表される。
推定される車室内騒音[SPL]は、数34により表される。
数25のSPL(s)は、数35となる。
[H](s)の分母は、H(s)及びH(s)より小さい遅延時間を有する。遅延時間補償τと遅延時間Taが等しい場合、時間遅延は完全に削除され、SPL(s)は数36となる。
<第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
具体的には、Na個のアクチュエータとNd個の外乱(加速度センサ)を用いた場合の、遅延時間補償されたj番目の車室内騒音[SPL]jは、数37により与えられる。
ここで、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
先ず、図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
第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
推定音圧[SPL]が使用されるとき1つのコヒーレンス関数Cs(f)を使用できる。
ここで、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つのコヒーレンス関数の場合に比べて非常に大きくなる。
ここで、Sxy(f)とSxx(f)は、数40に示すとおりである。
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で表される。
制御を行わない時、アクチュエータへの制御指令信号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になる。
図31は、更新モデル64を取り入れたモデル更新部61の構成を示すブロック図である。
FIG. 31 is a block diagram illustrating a configuration of the
次に、更新モデルKを計算する演算部63について説明する。判定部62が更新を決定した時、数43に示すように、実測音圧SPLを用いて更新モデルKを計算する。
これを変形すると、数44となる。
更新モデル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
更新のアルゴリズムは以下のとおりである。(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
以上説明したように、本発明の第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
マイクロフォンにより測定された車室内騒音と騒音推定部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
予め定めた周波数帯において、マイクロフォンにより測定された車室内騒音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
上記のように、本発明は、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.
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 ...
6 ...
DESCRIPTION OF
32, 32, 32a, 32b, 32c ... noise estimation unit (noise estimation means)
33a, 33b, 33c ...
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.
前記センサの出力信号から前記第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の伝達モデルを更新するモデル更新部と
を更に備えることを特徴とする請求項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.
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)
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)
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 |
-
2006
- 2006-02-17 JP JP2006040881A patent/JP4735319B2/en not_active Expired - Fee Related
Patent Citations (4)
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)
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 |