JP4002498B2

JP4002498B2 - Car interior sound simulation device and order component removal filter

Info

Publication number: JP4002498B2
Application number: JP2002302999A
Authority: JP
Inventors: 正男小松
Original assignee: Ono Sokki Co Ltd
Current assignee: Ono Sokki Co Ltd
Priority date: 2002-10-17
Filing date: 2002-10-17
Publication date: 2007-10-31
Anticipated expiration: 2022-10-17
Also published as: JP2004138788A

Description

【０００１】
【発明の属する技術分野】
本発明は、自動車室内の音をシミュレートする車室内音シミュレート装置およびその車室内音シミュレート装置に好適な次数成分除去フィルタに関する。
【０００２】
【従来の技術】
自動車の車室内音の低減化が問題となっており、その車室内音を能動的に低減化することが試みられている（例えば、特許文献１参照）。
【０００３】
自動車の車室内音が低減化されたかどうかを評価するためには、自動車の車室内の音を走行状態に応じて再現又はシミュレートする必要がある。このような要請から、自動車の車室内音の評価方法の一つとして、室内の音を車の走行状態に応じて再現又はシミュレートする方法が、重要な手法となってきている。
【０００４】
この時の音の再現方法の一つとして、音の成分を、エンジンの回転等の回転速度に比例した次数成分とそれ以外のランダム音とに分離し、車のおかれた、それぞれの状況に応じ、各々の成分音を合成して、実際の走行状況に近い音を再現する方法がある。次数成分の音が支配的であるが、ランダム音も音をリアルに再現するために必要な成分である。この方法に於いては、次数成分とランダム成分とを分離し、音を個々に分析加工した後、合成、再現する方法が一つの重要な手法となっている。
【０００５】
分離の方法として、従来は、音にＦＦＴ等の演算を施し、周波数軸上で次数成分を除去加工し、逆ＦＦＴ演算により、残った次数成分以外の音をランダム成分の音として取り出す手法が主たる方法である。
【０００６】
しかし、信号のデータを固定周波数でサンプルした時、基本次数が低い周波数から高い周波数へと（１０倍程度）変わる場合は、ＦＦＴの周波数分解能を低い周波数に合わせる必要がある。この時、分解能を低い周波数に合わせると、分析の信号の長さを長く取る必要が出てくる。又高い周波数に合わせると低い方の周波数の除去が分解能の不足によって、正しく出来ないという問題がある。
【０００７】
一方、上記問題を回避する方法として、信号を所定の回転角度毎にサンプルし、ＦＦＴで演算した場合は、ＦＦＴの分解能は必要最低限を選択できる。しかし、このように回転角度を基準としてサンプルすると、扱う周波数の上限が回転速度によって限定される。たとえば自動車のエンジン回転は、低回転１０００ｒｐｍから高回転７０００ｒｐｍと適応回転範囲が広い。ＦＦＴ演算の分解能とサンプル上限周波数が災いして、高い回転では、２０ｋＨｚ程度の音までサンプル対応できるが、低い回転の場合は上限周波数範囲が下がって、音の高い成分がカットオフされてしまうという問題がある。
【０００８】
上記問題を回避する方法、つまり次数成分を回転の変化に追従して除去する方法としては、各次数成分に対応してノッチフィルタを用意し、次数成分毎に成分除去を行うか、フィルタを必要数準備し、同時に複数の次数成分を除去する方法がある。しかし同時に行うとコスト的に、もしくは、演算回路の規模から考えても不利であり、更に効率の良い方法が求められている。
【０００９】
【特許文献１】
特開平５−３３３８８３号公報
【００１０】
【発明が解決しようとする課題】
本発明は、上記事情に鑑み、小さい回路規模かつ少ない演算数で、ランダム成分と次数成分を十分な分解能で分解する機能を備えた車室内音シミュレート装置、およびその車室内音シミュレート装置に好適な次数成分除去フィルタを提供することを目的とする。
【００１１】
【課題を解決するための手段】
上記目的を達成する本発明の車室内音シミュレート装置は、自動車室内の音のピックアップして音信号を得るとともに、その音をピックアップした時点における、その音の発生に寄与する回転体の回転速度を含むその自動車の走行状態をピックアップし、その音をピックアップして得た音信号をランダム成分と上記回転体の回転速度に応じた次数成分とに分離し、音の再生にあたり、上記次数成分に基づいて上記回転体の所望の回転速度に応じた次数成分を生成し、その生成された次数成分と、上記ランダム成分に基づいて得られるランダム成分とを合成して出力する車室内音シミュレート装置において、
上記音信号を制御信号に応じた遅延量だけ遅延する可変遅延部と、その可変遅延部により遅延される前の音信号からその可変遅延部で遅延された後の音信号を減算する信号減算部と、上記回転体の回転速度に基づいて、可変遅延部の遅延量を、信号減算部で上記回転体の回転に起因する次数成分が除去されたランダム成分が得られる遅延量に制御する制御信号を生成して可変遅延部に伝達する遅延制御部とを有する次数成分除去フィルタを備えたことを特徴とする。
【００１２】
ここで、「上記次数成分（回転体の回転速度に応じた次数成分）に基づいて、上記回転体の所望の回転速度に応じた次数成分を生成し」は、音をピックアップして音信号を得た時点における回転体の回転速度に応じた次数成分をそのままその回転体の所望の回転速度に応じた次数成分に変換してもよく、あるいは所望の回転速度に応じた次数成分への変換前あるいは変換後において次数成分に加工を加えてもよいことを意味している。
【００１３】
また、「上記ランダム成分（自動車室内の音をピックアップして得た音信号中のランダム成分）に基づいて得られるランダム成分」は、その音信号から抽出したままのランダム成分であってもよく、その抽出したランダム成分に加工を加えた後のランダム成分であってもよいことを意味している。
【００１４】
本発明の車室内音シミュレート装置は、上記回転体の、ある回転速度での車室内の音をピックアップした音信号に基づいて、その音をピックアップしたときの回転体の回転速度と同一あるいは異なる所望の回転速度での音を作り出すことができ、例えば複数の回転速度それぞれでピックアップして得た音信号に基づいてそれら複数の回転速度の間をつなぐことにより、例えばエンジンの回転数が徐々に上昇しているときの音を作り出すことができる。また、次数成分あるいはランダム成分を任意に加工して加工後の音をシミュレートすることもできる。
【００１５】
ここで、上記本発明の車室内音シミュレート装置は、上記の次数成分除去フィルタを備えたものであるため、従来と比べ小さい回路規模あるいは速い演算速度で、音信号を高精度にランダム成分と次数成分に分離することができる。
【００１６】
また、上記目的を達成する本発明の次数成分除去フィルタは、所定の回転体の回転速度に応じた次数成分とその回転体の回転速度とは独立のランダム成分とを含む信号を入力して、次数成分を除去することによりランダム成分を抽出する次数成分除去フィルタにおいて、
上記信号を、制御信号に応じた遅延量だけ遅延する可変遅延部と、その可変遅延部により遅延される前の信号からその可変遅延部で遅延された後の信号を減算する信号減算部と、上記回転体の回転速度に基づいて、可変遅延部の遅延量を、信号減算部で上記回転体の回転に起因する次数成分が除去された信号が得られる遅延量に制御する制御信号を生成して可変遅延部に伝達する遅延制御部とを備えたことを特徴とする。
【００１７】
本発明の次数成分除去フィルタは、いわゆる櫛形のフィルタの一種である。櫛形フィルタは整数倍の高次成分を同時に除去することがが可能である。しかし従来の使用法による櫛形フィルタはサンプル周期の整数倍であったため、除去可能な周波数成分は、定まった基本波と整数倍の周波数の除去になってしまっている。本発明は、回転速度から求めた遅延量だけ信号を遅延し、遅延前の信号から遅延後の信号を減算することでフィルタを実現し、少ない演算で基本波及びその整数倍の高次成分を一括して除去し、ランダム成分を取り出すことができる。
【００１８】
【発明の実施の形態】
以下、本発明の実施形態について説明する。
【００１９】
図１は、本発明の車室内音シミュレート装置の一実施形態を示すブロック図である。
【００２０】
自動車の室内の音がセンサ１０でピックアップされてアナログの音信号Ｓが生成され、Ａ／Ｄ変換器２０でディジタルの音信号に変換される。
【００２１】
また、センサ１１，１２により、この自動車のエンジンやそのエンジンからの動力を伝達する部材等の回転体の回転速度Ｒと、この自動車の走行速度Ｔがそれぞれピックアップされる。
【００２２】
Ａ／Ｄ変換器２０から出力されたディジタルの音信号は、信号分離部３０により、センサ１１により回転速度がピックアップされた回転体の、その回転速度に応じた次数成分と、その回転体の回転速度とは無関係のランダム成分とに分離される。信号分離部３０で分離された音信号の次数成分やランダム成分は、センサ１０でその音信号をピックアップした時の、センサ１１，１２でピックアップされた回転速度Ｒおよび走行速度Ｔに対応づけられて信号記憶部４０に記憶される。
【００２３】
音信号を収録する段階においては、上記の音収録が回転体の複数の回転速度と自動車の複数の走行速度との各組合せごとに行なわれる。
【００２４】
信号記憶部４０に記憶された音信号（次数成分とランダム成分）は、所望のシミュレートを行なうために必要に応じて加工される。
【００２５】
音の再生に当たっては、例えば時間的に変化する走行速度信号Ｔ’と回転速度信号Ｒ’が入力され、信号記憶部４０からは、そこに記憶されている複数の音信号（次数成分とランダム成分）の中から現在入力されている走行速度信号Ｔ’と回転速度信号Ｒ’との組合せに近い１つもしくは複数の音信号（次数成分とランダム成分）が読み出され、信号変換部５０に入力される。信号変換部５０では、信号記憶部４０から読み出された音信号中の次数成分が、現在入力されている回転速度信号Ｒ’に応じた回転速度の次数成分に変換されてランダム成分と合成される。信号記憶部４０から複数の音信号が読み出されるときは、上記の合成前あるいは合成後において線形補間演算等によりそれら複数の音信号の中間的な音信号（次数成分は現在の回転速度信号Ｓ’に応じた次数成分に変換された後の音信号）が生成される。
【００２６】
信号変換部５０で得られた走行速度信号Ｔ’および回転速度信号Ｒ’に応じた音信号は、Ｄ／Ａ変換器６０でアナログの音信号に変換されて、スピーカ７０から音として空間に放出される。
【００２７】
図２は、図１の信号分離部３０を構成する次数成分除去フィルタを示すブロック図である。
【００２８】
この次数成分除去フィルタ３００は、可変遅延部３１と、信号減算部３２と、遅延制御部３３とを備えている。
【００２９】
可変遅延部３１は、Ａ／Ｄ変換器２０（図１参照）から出力されたディジタルの音信号を入力し、その音信号を、遅延制御部３３で生成された制御信号Ｃに応じた遅延量だけ遅延して出力するものである。
【００３０】
また、信号減算部３２は、可変遅延部３１で遅延される前の音信号からその可変遅延部３１で遅延された後の音信号を減算するものである。
【００３１】
また、遅延制御部３３には、センサ１１でピックアップされた回転速度Ｒが入力され、遅延制御部３３では、その回転速度Ｒに基づいて、可変遅延部３１の遅延量を、信号減算部３２で、その回転速度Ｒで回転する回転体の回転に起因する次数成分が除去されたランダム成分が得られる遅延量に制御するための制御信号Ｃが生成されて可変遅延部３１に入力される。
【００３２】
こうすることにより、この図２に示す次数成分除去フィルタ３００に入力される入力信号が（Ａ）のような周波数特性を有するものである場合、この次数成分除去フィルタ３００では、（Ｂ）の出力信号の周波数特性に見られるように、直流成分、および回転体の回転速度Ｒに起因する基本周波数およびその整数倍の高周波が一括に除去される。
【００３３】
尚、次数成分は、図１の信号分離部３０において、Ａ／Ｄ変換器２０から出力された音信号から、図２の次数成分除去フィルタ３００の信号減算部３２から出力されたランダム成分を減算することにより、抽出される。
【００３４】
次に、図２に示す次数成分除去フィルタの特性を理論的に説明する。
【００３５】
図３は、図２の次数成分除去フィルタの理論上の周波数特性を示す図である。
【００３６】
ここでは、入力をｘ（ｔ）、出力をｙ（ｔ）とし、遅延をＺ^-nで表す。このフィルタの特性はＺ変換を用いると
Ｙ（ｚ）＝（１−Ｚ^-n）Ｘ（ｚ）
Ｈ（ｚ）は、
Ｈ（ｚ）＝１−Ｚ^-n
となり、
Ｚ^-1＝ｅｘｐ（−ｊｗｎＴ）
とすると、Ｔがサンプル周期の時、ｗは２π／ｎＴの周期関数となる。
【００３７】
ｗがπ／Ｔの時、
ｈ（ｗ）＝１−ｃｏｓ（ｎｗＴ）十ｊ・ｓｉｎ（ｎｗＴ）
ゲイン特性Ｇ（ｗ）は

となり、ｎに逆比例して、任意の基本周波数とその高調波を除去できる。
【００３８】
Ｇ（ｗ）はｎｗＴが２πのときサンプル周波数の特性になり、π以上は周波数が折り返しになる。サンプル周波数／２に対しＤＣ＋ｎ／２本の次数成分を除去できる。
【００３９】
従って基本周波数をｆ０としたとき、ｎは、
ｎ＝（サンプル周波数）／ｆ０
サンプル周波数＝１／Ｔ
で求められる。
【００４０】
以上の理論的な考察をグラフに表わすと、図３のようになる。
【００４１】
図４は、図２に示す次数成分除去フィルタの実際の入出力信号波形の一例を示す図である。
【００４２】
図４（Ａ）は入力信号波形を示したものであり、ランダム成分と、回転体の回転速度に応じた次数成分が含まれている。
【００４３】
図４（Ｂ）は、図２の次数成分除去フィルタの出力信号を示したものであり、次数成分が除去された信号が得られている。
【００４４】
次に、従来技術との比較結果について説明する。ここでは、次数成分除去の方式として本発明方式の他に下記の（１）〜（３）の３通りを考え、比較した。
【００４５】
（１）サンプル周波数固定のＦＦＴ方式の場合
（２）サンプル周波数可変のＦＦＴ方式の場合
（３）ノッチフィルタの組み合せた場合
（４）本発明方式を用いた場合
表１は、その比較結果を示す。
【００４６】
【表１】

【００４７】
ここでは、自動車の車室内音に適用した時の場合を考え、演算負荷、メモリ容量、及び応答性について比較を行う。車室内騒音にもっとも大きく影響するエンジン回転速度範囲は一般的に６００ｒｐｍから６０００ｒｐｍである。よって次数成分の除去は、基本周波数範囲が１０Ｈｚ〜１００Ｈｚ間で基本波とその高調波を除く事である。この条件で比較を行なう。
【００４８】
（１）信号を固定周波数で行いＦＦＴ演算する場合
一般的にＦＦＴで周波数分析を行う場合は、ハニングウインドウが使用される。この場合、除去したい周波数がＦＦＴ周波数分析ライン上の場合でも、両サイドのラインにセンターのラインの１／２の大きさのスペクトル信号が表れる。従って次数スペクトル間は最低でも周波数分解能の５倍以上が必要である。つまり除去すべき周波数の基本周波数の周期の５倍以上のデータの長さが分離、除去するために必要になる。
【００４９】
下限の周波数は１０Ｈｚなので、必要なＦＦＴの分解能は１０Ｈｚ／５＝２Ｈｚ（データ長は５００ｍＳ）となる。ＦＦＴ分析結果は５００ｍｓの平均のスペクトルとみなせるので応答は２Ｈｚと考える。
【００５０】
又周波数の上限は、取扱い信号を可聴周波数とし、１０ｋＨｚまでとした時、エリアジング防止を考慮したＦＦＴアナライザ等では一般的にサンプル周波数は１０ｋＨｚ×２．５６＝２５．６ｋＨｚとなる。必要な周波数分解能２Ｈｚとサンプル周波数２５．６ｋＨｚとにより必要なＦＦＴ演算のデータ点数ｎは、
２Ｈｚ×２ⁿ＞２５．６ｋＨｚから
ｎ≧１４となり、
ＦＦＴの演算点数は２¹⁴＝１６３８４ポイント以上が必要である。
【００５１】
ＦＦＴ演算における四則演算数は、１６３８４点の場合、時間軸に戻し、音を発生させるので次数除去の為にＦＦＴを実施し、周波数軸上で次数成分を除去し、その後、周波数軸のデータを時間軸に戻すため、逆ＦＦＴ演算を行う。従って１６３８４／２×（１４＋１４）×２回のバタフライ演算が必要となる。上式の最後の×２は、ハニングウインドウをデータに乗じた時点でデータのはじめと終わりにデータが０になる事に対するハニングウインドウの影響を除去するためである。
【００５２】
ここで、バタフライ演算の四則演算数は、１回のバタフライ演算につき乗算４回、和差４回とすると、１点あたり２８×４乗算、２８×４和差となる。又必要なメモリ数は１６３８４リアル、イマジナリのデータメモリが必要なので３２７６８メモリとなる。
【００５３】
（２）サンプル周波数可変のＦＦＴ方式の場合
回転に合わせサンプルを可変にして（いわゆるトラッキングと言われる方法）でデータサンプルして１０２４点のＦＦＴを実施した場合、演算負荷に関して、ＦＦＴの演算の基本であるバタフライ演算は、
ＦＦＴ１０２４／２×１０バタフライ
ＩＮＶＦＦＴ１０２４／２×１０バタフライ
ハニング時５０％のオーバーラップが必要なので ×２
１０２４点に対して１０２４／２×（（１０十１０）×２）
１点につき２０バタフライ演算が必要
１バタフライにつき掛算４回、和差４回が必要なので
（８０の乗算と８０の和差）／データ
出力される信号の計算のサイクルは２Ｈｚ〜２０Ｈｚとなる。
【００５４】
必要なメモリ数は、１０２４（リアルパート）＋１０２４（イマジナリパート用）＝２０４８メモリとなる。
【００５５】
但しサンプルを可変にするため、実施にあたり、上記条件の場合では１回転に１０２４／５サンプルのためのパルスの発生回路と、折り返し防止のために、回転に応じてカットオフ特性が変わるフィルタが必要となる。
【００５６】
（３）ノッチフィルタの組み合せの場合
デジタルのノッチフィルタ１構成につき、メモリ３、掛算３、和２となる。基本周波数が１００Ｈｚの時、１０ｋＨｚまで除去する場合１００の次数除去になり、メモリ３００、掛算３００、和２００となる。基本周波数が低い時は、更に増加する。基本周波数の応答１０Ｈｚ〜１００Ｈｚ、他のノッチフィルタはその除去周波数と考えられる。
【００５７】
（４）ＦＩＲフィルタを利用して補間を実施した場合
本発明方式では、基本周波数の周期で信号の除去を逐次行なっていくので、ＦＦＴ方式に比べ、上記条件で実施した場合は、最低でも５倍良好な応答性が得られる。
【００５８】
又、サンプルに含まれる全ての高調波を除去可能となる。１点につき２０ステップのＦＩＲ処理と引き算、１点につき２０の乗算と２０の和差及び１回の差を行なう。応答性は、基本周波数の周期だけ信号を遅延させた信号間の差の演算なので、回転速度に比例し１０Ｈｚ〜１００Ｈｚとなる。
【００５９】
また、データメモリは遅延分２５．６ｋＨｚサンプルで１０Ｈｚの周期の遅延を実現する点数が必要なので、２５６００÷１０＝２５６０ポイント分のデータメモリが必要となる。
【００６０】
本発明の方式は、データメモリ関しては、他の方式と比べ優位性はないが、表１に有るように演算数が少ないのが大きな特徴である。２つのＦＦＴの方式に比べて応答も性もよい。ノッチフィルタの組み合せも、応答性がよく、メモリが少ないが、本方式と比べると、１データあたりの演算点数が多くなり、ノッチフィルタ毎に係数を管理するため、操作が煩雑になり、補間にＦＩＲのフィルタを利用する方法を用いると、本方式の方が非常に容易に構成しやすい。
【００６１】
【発明の効果】
以上、説明したように、本発明によれば小さい回路規模かつ少ない演算数で、信号をランダム成分と次数成分に十分に分解能で分離することができる。
【図面の簡単な説明】
【図１】本発明の車室内音シミュレート装置の一実施形態を示すブロック図である。
【図２】図１の信号分離部を構成する次数成分除去フィルタを示すブロック図である。
【図３】図２の次数成分除去フィルタの理論上の周波数特性を示す図である。
【図４】図２に示す次数成分除去フィルタの実際の入出力信号波形の一例を示す図である。
【符号の説明】
１０，１１，１２センサ
２０Ａ／Ｄ変換器
３０信号分離部
３１可変遅延部
３２信号減算部
３３遅延制御部
４０信号記憶部
５０信号変換部
６０Ｄ／Ａ変換器
７０スピーカ
３００次数成分除去フィルタ[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a vehicle interior sound simulating device that simulates sound in a vehicle interior and an order component removal filter suitable for the vehicle interior sound simulating device.
[0002]
[Prior art]
Reduction of vehicle interior noise has become a problem, and attempts have been made to actively reduce vehicle interior noise (see, for example, Patent Document 1).
[0003]
In order to evaluate whether or not the vehicle interior sound has been reduced, it is necessary to reproduce or simulate the sound in the vehicle interior in accordance with the running state. From such a demand, a method of reproducing or simulating the sound in the room according to the running state of the car has become an important technique as one of the evaluation methods of the sound in the car interior of the automobile.
[0004]
As one of the sound reproduction methods at this time, the sound component is separated into the order component proportional to the rotational speed such as the engine speed and the other random sound, and the situation where the car is placed Accordingly, there is a method of reproducing the sound close to the actual driving situation by synthesizing each component sound. The order component sound is dominant, but the random sound is also a component necessary for realistic reproduction of the sound. In this method, an important method is to separate the order component and the random component, individually analyze and process the sound, and then synthesize and reproduce it.
[0005]
Conventionally, as a separation method, a method is mainly used in which sound such as FFT is performed on the sound, the order component is removed on the frequency axis, and the sound other than the remaining order component is extracted as a random component sound by inverse FFT operation. Is the method.
[0006]
However, when the signal data is sampled at a fixed frequency, if the basic order changes from a low frequency to a high frequency (about 10 times), it is necessary to match the FFT frequency resolution to a low frequency. At this time, if the resolution is adjusted to a low frequency, it is necessary to increase the length of the analysis signal. Further, when the frequency is adjusted to a high frequency, there is a problem that the lower frequency cannot be correctly removed due to lack of resolution.
[0007]
On the other hand, as a method of avoiding the above problem, when a signal is sampled at every predetermined rotation angle and calculated by FFT, the minimum resolution of FFT can be selected. However, when sampling is performed based on the rotation angle as described above, the upper limit of the frequency to be handled is limited by the rotation speed. For example, the engine rotation of an automobile has a wide adaptive rotation range from a low rotation of 1000 rpm to a high rotation of 7000 rpm. The resolution of the FFT operation and the sample upper limit frequency are damaged, and at high rotations, it is possible to sample up to a sound of about 20 kHz. However, at low rotations, the upper limit frequency range is lowered and high sound components are cut off. There's a problem.
[0008]
As a method of avoiding the above problem, that is, a method of removing the order component following the change in rotation, a notch filter is prepared corresponding to each order component and the component is removed for each order component or a filter is required. There are methods of preparing several and simultaneously removing a plurality of order components. However, if it is performed simultaneously, it is disadvantageous in terms of cost or from the scale of the arithmetic circuit, and a more efficient method is required.
[0009]
[Patent Document 1]
Japanese Patent Laid-Open No. 5-333833
[Problems to be solved by the invention]
In view of the above circumstances, the present invention provides a vehicle interior sound simulation apparatus having a function of decomposing a random component and an order component with sufficient resolution with a small circuit scale and a small number of operations, and the vehicle interior sound simulation apparatus. It is an object of the present invention to provide a suitable order component removal filter.
[0011]
[Means for Solving the Problems]
The vehicle interior sound simulating apparatus of the present invention that achieves the above object obtains a sound signal by picking up sound in an automobile interior, and the rotational speed of a rotating body that contributes to the generation of the sound at the time of picking up the sound. The sound signal obtained by picking up the sound of the vehicle including the sound is separated into a random component and an order component corresponding to the rotational speed of the rotating body, and the sound is reproduced with the order component. Based on the desired rotational speed of the rotator, and the vehicle interior sound simulator for synthesizing and outputting the generated order component and a random component obtained based on the random component In
A variable delay unit that delays the sound signal by a delay amount corresponding to the control signal, and a signal subtracting unit that subtracts the sound signal delayed by the variable delay unit from the sound signal before being delayed by the variable delay unit. And a control signal that controls the delay amount of the variable delay unit based on the rotation speed of the rotating body to a delay amount that obtains a random component from which the order component due to the rotation of the rotating body is removed by the signal subtracting unit. And an order component removal filter having a delay control unit that generates and transmits it to the variable delay unit.
[0012]
Here, “the order component corresponding to the desired rotational speed of the rotating body is generated based on the order component (the order component corresponding to the rotating speed of the rotating body)” is used to pick up the sound and The order component corresponding to the rotational speed of the rotating body at the time of acquisition may be converted directly to the order component corresponding to the desired rotational speed of the rotating body, or before conversion to the order component corresponding to the desired rotational speed. Alternatively, it means that the order component may be processed after the conversion.
[0013]
In addition, the “random component obtained based on the random component (random component in the sound signal obtained by picking up the sound in the automobile interior)” may be a random component that has been extracted from the sound signal, This means that the extracted random component may be a random component after being processed.
[0014]
The vehicle interior sound simulation apparatus according to the present invention is based on a sound signal obtained by picking up the sound in the vehicle interior at a certain rotational speed of the rotating body, and is the same as or different from the rotational speed of the rotating body when the sound is picked up. A sound can be produced at a desired rotational speed. For example, by connecting between the plurality of rotational speeds based on a sound signal obtained by picking up at each of the plurality of rotational speeds, for example, the engine speed is gradually increased. Can produce sounds when rising. Further, it is possible to simulate the processed sound by arbitrarily processing the order component or the random component.
[0015]
Here, the vehicle interior sound simulation apparatus of the present invention includes the above-described order component removal filter, so that a sound signal can be converted into a random component with high accuracy at a smaller circuit scale or faster operation speed than conventional ones. It can be separated into order components.
[0016]
In addition, the order component removal filter of the present invention that achieves the above object inputs a signal including an order component corresponding to the rotation speed of a predetermined rotating body and a random component independent of the rotation speed of the rotating body, In an order component removal filter that extracts random components by removing order components,
A variable delay unit that delays the signal by a delay amount according to the control signal, a signal subtracting unit that subtracts a signal delayed by the variable delay unit from a signal before being delayed by the variable delay unit, Based on the rotation speed of the rotator, a control signal is generated to control the delay amount of the variable delay unit to a delay amount from which the signal subtraction unit can obtain a signal from which the order component due to the rotation of the rotator is removed. And a delay control unit for transmitting to the variable delay unit.
[0017]
The order component removal filter of the present invention is a kind of so-called comb filter. The comb filter can simultaneously remove integer multiples of higher order components. However, since the comb filter according to the conventional method is an integral multiple of the sample period, the removable frequency component is the removal of the fixed fundamental wave and the integral multiple of the frequency. The present invention realizes a filter by delaying a signal by a delay amount obtained from the rotation speed, and subtracting the delayed signal from the signal before the delay, and the fundamental wave and its higher-order component that is an integral multiple of the fundamental wave with fewer operations. Random components can be extracted by removing all at once.
[0018]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described.
[0019]
FIG. 1 is a block diagram showing an embodiment of a vehicle interior sound simulation apparatus according to the present invention.
[0020]
The sound in the vehicle interior is picked up by the sensor 10 to generate an analog sound signal S, which is converted into a digital sound signal by the A / D converter 20.
[0021]
The sensors 11 and 12 pick up the rotational speed R of the rotating body such as the engine of the automobile and a member for transmitting power from the engine, and the running speed T of the automobile, respectively.
[0022]
The digital sound signal output from the A / D converter 20 is converted into the order component corresponding to the rotational speed of the rotating body picked up by the sensor 11 by the signal separation unit 30 and the rotation of the rotating body. It is separated into random components that are independent of speed. The order components and random components of the sound signal separated by the signal separation unit 30 are associated with the rotational speed R and the traveling speed T picked up by the sensors 11 and 12 when the sound signal is picked up by the sensor 10. It is stored in the signal storage unit 40.
[0023]
In the stage of recording the sound signal, the above sound recording is performed for each combination of a plurality of rotational speeds of the rotating body and a plurality of traveling speeds of the automobile.
[0024]
The sound signal (order component and random component) stored in the signal storage unit 40 is processed as necessary in order to perform a desired simulation.
[0025]
For sound reproduction, for example, a traveling speed signal T ′ and a rotational speed signal R ′ that change with time are input, and a plurality of sound signals (order components and random components) stored therein are input from the signal storage unit 40. ), One or a plurality of sound signals (order component and random component) close to the combination of the currently input traveling speed signal T ′ and the rotational speed signal R ′ are read out and input to the signal conversion unit 50. Is done. In the signal conversion unit 50, the order component in the sound signal read from the signal storage unit 40 is converted into the order component of the rotational speed corresponding to the currently input rotational speed signal R ′ and synthesized with the random component. The When a plurality of sound signals are read from the signal storage unit 40, an intermediate sound signal (order component is the current rotational speed signal S ′) by linear interpolation calculation or the like before or after the above synthesis. Sound signal after being converted into an order component in accordance with.
[0026]
The sound signal corresponding to the traveling speed signal T ′ and the rotational speed signal R ′ obtained by the signal conversion unit 50 is converted into an analog sound signal by the D / A converter 60 and emitted from the speaker 70 into the space as sound. Is done.
[0027]
FIG. 2 is a block diagram showing an order component removal filter constituting the signal separation unit 30 of FIG.
[0028]
The order component removal filter 300 includes a variable delay unit 31, a signal subtraction unit 32, and a delay control unit 33.
[0029]
The variable delay unit 31 receives the digital sound signal output from the A / D converter 20 (see FIG. 1), and delays the sound signal according to the control signal C generated by the delay control unit 33. Output with a delay of only.
[0030]
The signal subtracting unit 32 subtracts the sound signal after being delayed by the variable delay unit 31 from the sound signal before being delayed by the variable delay unit 31.
[0031]
Further, the rotational speed R picked up by the sensor 11 is input to the delay control section 33, and the delay control section 33 determines the delay amount of the variable delay section 31 based on the rotational speed R by the signal subtraction section 32. A control signal C for controlling the amount of delay to obtain a random component from which the order component resulting from the rotation of the rotating body rotating at the rotation speed R is removed is generated and input to the variable delay unit 31.
[0032]
Thus, when the input signal input to the order component removal filter 300 shown in FIG. 2 has a frequency characteristic as shown in (A), the order component removal filter 300 outputs (B). As can be seen from the frequency characteristics of the signal, the DC component, the fundamental frequency resulting from the rotational speed R of the rotating body, and the integer multiple of the fundamental frequency are removed at once.
[0033]
The order component is obtained by subtracting the random component output from the signal subtracting unit 32 of the order component removing filter 300 in FIG. 2 from the sound signal output from the A / D converter 20 in the signal separation unit 30 in FIG. To be extracted.
[0034]
Next, the characteristics of the order component removal filter shown in FIG. 2 will be theoretically described.
[0035]
FIG. 3 is a diagram illustrating theoretical frequency characteristics of the order component removal filter of FIG.
[0036]
Here, the input is x (t), the output is y (t), and the delay is represented by Z ⁻ⁿ . The characteristics of this filter are Y (z) = (1−Z ⁻ⁿ ) X (z) when Z conversion is used.
H (z) is
H (z) = 1−Z ⁻ⁿ
And
Z ⁻¹ = exp (−jwnT)
Then, when T is a sample period, w is a periodic function of 2π / nT.
[0037]
When w is π / T,
h (w) = 1−cos (nwT) + j · sin (nwT)
The gain characteristic G (w) is

Thus, an arbitrary fundamental frequency and its harmonics can be removed in inverse proportion to n.
[0038]
G (w) has the characteristics of the sample frequency when nwT is 2π, and the frequency is turned over when π or more. DC + n / 2 order components can be removed with respect to the sample frequency / 2.
[0039]
Therefore, when the fundamental frequency is f0, n is
n = (sample frequency) / f0
Sample frequency = 1 / T
Is required.
[0040]
The above theoretical consideration is represented in a graph as shown in FIG.
[0041]
FIG. 4 is a diagram showing an example of actual input / output signal waveforms of the order component removal filter shown in FIG.
[0042]
FIG. 4A shows an input signal waveform, which includes a random component and an order component corresponding to the rotational speed of the rotating body.
[0043]
FIG. 4B shows an output signal of the order component removal filter of FIG. 2, and a signal from which the order component has been removed is obtained.
[0044]
Next, a comparison result with the prior art will be described. Here, the following three types (1) to (3) were considered and compared in addition to the method of the present invention as the method of removing the order component.
[0045]
(1) In case of FFT system with fixed sample frequency (2) In case of FFT system with variable sample frequency (3) In case of combination of notch filters (4) In case of using the system of the present invention Table 1 shows the comparison result .
[0046]
[Table 1]

[0047]
Here, considering the case where it is applied to the interior sound of an automobile, the calculation load, the memory capacity, and the response are compared. The engine speed range that has the greatest influence on vehicle interior noise is generally 600 rpm to 6000 rpm. Therefore, the removal of the order component is to remove the fundamental wave and its harmonics when the fundamental frequency range is between 10 Hz and 100 Hz. Comparison is performed under these conditions.
[0048]
(1) When performing FFT calculation with a signal at a fixed frequency In general, when performing frequency analysis with FFT, a Hanning window is used. In this case, even when the frequency to be removed is on the FFT frequency analysis line, a spectrum signal having a size half that of the center line appears on both side lines. Therefore, at least 5 times the frequency resolution is required between the order spectra. That is, a data length of 5 times or more of the period of the fundamental frequency of the frequency to be removed is necessary for separation and removal.
[0049]
Since the lower limit frequency is 10 Hz, the necessary FFT resolution is 10 Hz / 5 = 2 Hz (data length is 500 mS). Since the FFT analysis result can be regarded as an average spectrum of 500 ms, the response is considered to be 2 Hz.
[0050]
When the upper limit of the frequency is an audible frequency up to 10 kHz, the sample frequency is generally 10 kHz × 2.56 = 25.6 kHz in an FFT analyzer or the like considering aliasing prevention. The required number n of FFT calculation data points with a required frequency resolution of 2 Hz and a sample frequency of 25.6 kHz is:
From 2 Hz × 2 ⁿ > 25.6 kHz, n ≧ 14,
The number of FFT calculation points needs to be 2 ¹⁴ = 16384 points or more.
[0051]
In the case of 16384 points, the number of four arithmetic operations in the FFT operation is returned to the time axis and generates a sound, so the FFT is performed to remove the order, the order component is removed on the frequency axis, and then the frequency axis data is In order to return to the time axis, an inverse FFT operation is performed. Therefore, 16384/2 × (14 + 14) × 2 butterfly operations are required. The last x2 in the above equation is for removing the influence of the Hanning window on the data becoming zero at the beginning and end of the data when the data is multiplied by the Hanning window.
[0052]
Here, the number of four arithmetic operations for the butterfly operation is 28 × 4 multiplications and 28 × 4 sum differences per point, where 4 multiplications and 4 sum differences are made per butterfly operation. The number of necessary memories is 16384 real, and since imaginary data memory is required, the number is 32768 memories.
[0053]
(2) In the case of the FFT method with variable sample frequency When the sample is made variable in accordance with the rotation (so-called tracking method) and data sampling is performed and 1024-point FFT is performed, the basics of the FFT calculation with respect to the calculation load The butterfly operation is
FFT 1024/2 × 10 butterfly INVFFT 1024/2 × 10 50% overlap is required for butterfly hanning × 2
1024/2 x 1024 points ((10 + 10) x 2)
20 butterfly operations are required for each point. Since 4 multiplications and 4 sums / differences are required for each butterfly (80 multiplications and 80 sums / differences), the cycle of calculation of the data output signal is 2 Hz to 20 Hz.
[0054]
The number of necessary memories is 1024 (real part) +1024 (for imaginary part) = 2048 memories.
[0055]
However, in order to make the samples variable, in the case of the above conditions, a pulse generation circuit for 1024/5 samples per rotation and a filter whose cut-off characteristics change according to the rotation are necessary to prevent aliasing. It becomes.
[0056]
(3) In the case of a combination of notch filters, the memory 3, the multiplication 3, and the sum 2 are obtained for each digital notch filter 1 configuration. When the fundamental frequency is 100 Hz, when removing up to 10 kHz, the order of 100 is removed, and the memory 300, the multiplication 300, and the sum 200 are obtained. It increases further when the fundamental frequency is low. The response of the fundamental frequency is 10 Hz to 100 Hz, and other notch filters are considered as the removal frequency.
[0057]
(4) When interpolation is performed using an FIR filter In the method of the present invention, signal removal is sequentially performed at the period of the fundamental frequency. Two times better responsiveness can be obtained.
[0058]
Moreover, it becomes possible to remove all the harmonics contained in the sample. FIR processing and subtraction of 20 steps per point, 20 multiplications per point, 20 sum differences and 1 difference are performed. Responsiveness is a calculation of a difference between signals obtained by delaying the signal by the period of the fundamental frequency, and is proportional to the rotation speed and becomes 10 Hz to 100 Hz.
[0059]
Further, since the data memory needs a point that realizes a delay of 10 Hz with a delay of 25.6 kHz samples, a data memory of 25600 ÷ 10 = 2560 points is required.
[0060]
The method of the present invention is not superior to other methods in terms of data memory, but has a major feature that the number of operations is small as shown in Table 1. Compared with the two FFT methods, the response and characteristics are good. The combination of notch filters also has good responsiveness and little memory, but compared to this method, the number of calculation points per data increases and the coefficients are managed for each notch filter, which complicates operations and makes interpolation difficult. If a method using an FIR filter is used, this method is much easier to configure.
[0061]
【The invention's effect】
As described above, according to the present invention, a signal can be separated into a random component and an order component with sufficient resolution with a small circuit scale and a small number of operations.
[Brief description of the drawings]
FIG. 1 is a block diagram showing an embodiment of a vehicle interior sound simulation apparatus according to the present invention.
2 is a block diagram showing an order component removal filter that constitutes the signal separation unit of FIG. 1; FIG.
FIG. 3 is a diagram illustrating theoretical frequency characteristics of the order component removal filter of FIG. 2;
4 is a diagram showing an example of actual input / output signal waveforms of the order component removal filter shown in FIG. 2. FIG.
[Explanation of symbols]
10, 11, 12 Sensor 20 A / D converter 30 Signal separation unit 31 Variable delay unit 32 Signal subtraction unit 33 Delay control unit 40 Signal storage unit 50 Signal conversion unit 60 D / A converter 70 Speaker 300 Order component removal filter

Claims

A sound signal is obtained by picking up the sound in the automobile room, and at the time of picking up the sound, the traveling state of the car including the rotational speed of the rotating body contributing to the generation of the sound is picked up, and the sound is picked up. The sound signal obtained in this manner is separated into a random component and an order component corresponding to the rotational speed of the rotating body, and when reproducing sound, an order component corresponding to a desired rotational speed of the rotating body is obtained based on the order component. In the vehicle interior sound simulator that generates and outputs the generated order component and a random component obtained based on the random component.
A variable delay unit that delays the sound signal by a delay amount corresponding to the control signal, and a signal subtractor that subtracts the sound signal delayed by the variable delay unit from the sound signal before being delayed by the variable delay unit. And the delay amount of the variable delay unit is controlled based on the rotational speed of the rotating body to a delay amount that can obtain a random component from which the order component due to the rotation of the rotating body is removed by the signal subtracting unit. A vehicle interior sound simulating apparatus comprising: an order component removal filter having a delay control unit that generates a control signal and transmits the control signal to the variable delay unit.

In an order component removal filter that inputs a signal including an order component corresponding to the rotational speed of a predetermined rotating body and a random component independent of the rotational speed and extracts the random component by removing the order component. ,
A variable delay unit that delays the signal by a delay amount according to a control signal; a signal subtracting unit that subtracts a signal delayed by the variable delay unit from a signal before being delayed by the variable delay unit; Based on the rotational speed of the rotating body, a control signal for controlling the delay amount of the variable delay unit to a delay amount that obtains a signal from which the order component due to the rotation of the rotating body is removed by the signal subtracting unit. An order component removal filter comprising: a delay control unit that generates and transmits the delay control unit to the variable delay unit.