JPS61212996A - 多チャンネル制御装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 「産業上の利用分野」 こ□の発明はスピーカ素子より放射された音波に対し、
これと同一の任意の３次元空間に設定さ杆た１つ、また
は２つ以」二の点の音圧波形を上記スピーカ素子へ供給
する入力信号を制御することにより、所望の音圧分布を
つくりだす多点放射形音圧制御装置に関するものである
。

［従来の技術−１ 近年「音像定位」を対象とする研究から、複数のスピー
カ素子を用いて複数個の点にお＆）る音圧をに二乗誤差
最小」なる評価尺度に基づいて制御する装置を構成する
手法の概念〔柳田他：Ｊ。

八ｃｏｕｓｔ、ｓｏｃ、Ｊｐｎ（Ｅ）、４，２，１０７
−１０９（１９８３）　）力く報告されている。ここで
は、この方法の概念を１時間軸−Ｌ」で展開して説明を
行う。以下第１Ｉ図に示す簡単なモデルに基づいて説明
する。第１１図は説明を簡単にするため、１点の音圧を
制御する従来法を示しているが、２点以上に対する手法
も全く同じ原理に基づいている。

第１１図において、室内反射音の存在する音場内に置か
れているマイクロホン素子１１．で受音される仮想スピ
ーカ素子Ｓ°の音圧を、仮想スピーカ素子Ｓ”を用いず
に、スピーカ素子１．を使って再現することができれば
、スピーカ素子１１が音響信号を放射しているにもかか
わらず、あたかも仮想スピーカ素子Ｓ゛から音響信号が
放射されているかのような状況をつくり出すことができ
る・。この明細書では、このような状況をつくり出す為
のスピーカ制御を「音圧制御」と呼ふことにする。上記
状況をつくり出す為には、仮想スピーカ素子Ｓ゛とマイ
クロホン素子１１１　　との間の信号伝達特性、即ちイ
ンパルス応答と、スピーカ素子１．及びマイクロホン素
子１１３間のインパルス応答とが等しくなるように、Ｆ
ＩＲフィルタ（トランスバーザルフィルタ）２１□の係
数を適切に定め、そのＦＩＲフィルタ２１１　を通じて
信号をスピーカ相素子１．へ供給すればよい。

さて、スピーカ素子Ｓ′とマイクロホン素子１１゜との
間のインパルス応答をｒ（ｋ）　（ｋ　＝　１．２．−
＝−、ｎ）とし、スピーカ素子１１とマイクロホン素子
１１．との間のインパルス応答をｇ　（ｋ）（ｋ　＝　
１　、２　、−、　ｍ）とすれば、次式を満たすフィル
タ係数ｈ（ｋ）（ｋ　−１、２、−−−−−、Ｔ、、’
）をＦＩＲフィルタ２１１　に与えることにより所望の
音場制御が実現される。

ｒ（ｋ）　＝　ｇ　（ｋ）■ｈ　（ｋ　）　　−−−−
−−−−−−−−−−−−−−−−−−−−−（１，１
但し、■は離散的畳み込み演算を表す。

（１）式において、ｒ（ｋ）を（ｒ、ｒｚ　−”−ｒｎ
　）”（Ｔは転置を表す）、ｇ（ｋ）をＣｇ＋ｇｚ−ｇ
ｍ　）　”　、ｈ（ｋ）　を（ｈ＋　　ｈｚ−ｈｔ　）
　Ｔとそれぞれ離散値で表現されるベクトルで表せば、
（１）式％式％ と書き換えられる。ところが（Ｉ）゛　式からフィルタ
係数ｈＩ　　＋　　２　　＋−’−＋　ｈｌ　を求める
場合、未知数り の数りに比べて方程式の数ｎの方が多い為、（１）式は
通常不能方程式となり、この式を満足する解り、、ｈ２
．−、ｈ、、を求めることは不可能である。

そこで従来においては「二乗誤差最小」を評価尺度とし
て、 Σ１ｅ（ｋ）１２−Σ巨（ｋ）　−ｇ　（ｋ）■ｈ（ｋ
）ビ　−＋２１で表される二乗誤差を最小とするフィル
タ係数ｈ（ｋ）を求め、これを用いて所望のインパルス
応答ｒ（ｋ）を近似的に実現するような音圧制御が行わ
れていた。このｈ（ｋ）は次式によって求められる。

しかしこの手法では、 （ａｌ　　誤差が最小になるようにするものであり、（
１）式を満足させることができないことによる誤差を原
理的に含んでいる。

（ｂｌ　　ｔａ＋の誤差は、各スピーカ素子、各被制御
点及び周囲環境の位置関係に大きく影響される。

という欠点があった。

この発明の目的は、任意の３次元空間内に設定された１
つあるいはそれ以上の数の被制御点において、各被制御
点ごとに定められる任意の信号伝達特性、即ち、インパ
ルス応答を、上記従来法のような近似誤差を伴わずに実
現することにより、所望の音圧分布を正確につくり出す
ことができる音圧制御装置を提供することにある。

「問題点を解決するための手段」 この発明によれば空間の複数の点にスピーカ素子が配置
され、これらスピーカ素子にそれぞれ各別のフィルタを
通じて同一人力信号が供給され、これらスピーカ素子と
同一空間に、そのスピーカ素子数より少なくとも１つは
少ない数の点（以後、被制御点と呼ぶ）についてそれぞ
れ定められた所望の信号伝達特性と、上記各スピーカ素
子及び上記各被制御点間の信号伝達特性とを用いて上記
各フィルタの係数が算出され、その係数を対応フィルタ
に設定することにより、上記人力信号が、上記各被制御
点において、その被制御点に与えられた上記所望の信号
伝達特性を満足するように制御される。従来においては
スピーカ素子数と被制御点とが同数又は後者が大とされ
たが、この発明では被制御点の数がスピーカ素子数より
必ず小とされる。

発明の原理 まず第１図を用いてこの発明の詳細な説明する。

第１図において、８個（Ｓ≧２）のスピーカ素子１１　
＋　１２　＋−−−−−＋　Ｉ　Ｓが空間に配される。

またその空間のＮ個の被制御点１１．、−、１１Ｎが設
けられる。

第１図では被制御点１１゜、−、１１Ｎにそれぞれマイ
クロホン素子が配された場合でこれら被制御点ＬＬ　、
−−−−−１１ｕをマイクロホン素子と呼ぶこともある
。ＮはＮ≦Ｓ−１とされる。入力端子１２から同一の入
力信号がそれぞれＦＩＲフィルタ２１１　、２１２゜−
−−、２１３を通じてスピーカ素子１１．１□、−１ｓ
に印加される。これらＦＩＲフィルタ２１＋　、　２１
□、−２１゜の各フィルタ係数を制御することによりそ
れぞれスピーカ素子１１，１□、　−’−１ｓに印加さ
れる入力信号を各別に制御することができる。

はじめに次の条件を仮定する。これは、通常の音響的環
境ではおおむね満たされるものである。

（条件−１） ８個のスピーカ素子１１＋　１ｚ＋−３１ｓ　とＮ個の
被制御点１１．、−、１１Ｎ　との間の各インパルス応
答はすべて、同じ有限長さｍを持つ離散的時系列（以下
、このような離散的時系列をヘクトルと呼ぶ）で表すこ
とができる。

（条件−２） （条件−１）で述べた各インパルス応答をＺ変換して得
られる伝達関数は互いに同じ零点を持たない。

さて、各被制御点り、−、ＩＮについてそれぞれ設定さ
れるＮ個の所望のインパルス応答をｒ、（ｋ）　　Ｕ＝
１．２．−　Ｎ）とし、８個のスピーカ素子１．、−、
　　Ｉ　Ｓ及びＮ個のマイクロホン素子１１＋　、　’
−−−’ＩＩＮ間のＳ−Ｎ個のインパルス応答をｇｉｊ
（ｋ）　（ｉ＝　１．２．−３　；　ｊ　＝　１　、２
．−、　Ｎ）とし、８個のスピーカ素子のそれぞれの前
段に設置するフィルタ２１．、−、２１．の各フィルタ
係数をｈ＋　　（ｋ）（ｉ　＝　１．２．−−−、Ｓ）
としたとき、ｒＪ（ｋ）　　−ΣｇｔＪ（ｋ）■ｈ　ｉ
　（ｋ）　−−−−−ｍ＝−−−１４１を満足するフィ
ルタ係数ｈｉ（ｋ）を、正確に求めることができれば、
従来技術の欠点を解消することができる。ここで、ｒＪ
（ｋ）　１ｇ　＋’；（ｋ）　　、　　ｈ　ｒ（ｋ）を
それぞれ（ｒｉ＋　ＴＪ□’−’ｒｊｎ）　Ｔ＋　　’
（ｇ　ｉｊ＋ｇ　ｉｊ□−ｇ　ｔａｌｌ）　”　、（ｈ
ｉｔ　ｈｔｚ　−−−−−−−−ｈｔｔ）　”（ｉ　＝
１．２．−−−、　Ｓ　ｉ　ｊ＝１．２．−、　Ｎ）な
るベクトルで表すことにすれば（４）式は、ｎ＝ｍ＋１
−−１ となる。以後の説明を容易にする為（４Ｆ式の左辺を１
（Ｎ−ｎｘｌのベクトル）、右辺の行列をＧ（Ｎ−ｎＸ
ｓ−Ｌの行列）、右辺のベクトルを■（Ｓ−ＬＸＩのベ
クトル）とおけば、（４）式は（４ビ式となる。

１１！＝Ｇ−１１１−−−−−−−−＋４１”（４ビ式
右辺の行列のに着目して、列の数が行の数より多くなる
条件、つまり未知数が方程式の数と同−又はこれより多
くなり、解が得られる条件、Ｓ−Ｌ≧Ｎ−ｎ＝Ｎ　・（
ｍ＋！−−１）４°、Ｌ≧Ｎ　・（ｍ−１）　／　（Ｓ
　−Ｎ　）　−−−−−−−［５１を満たすようなフィ
ルタ長りを定めれば、（条件−２）によって、行列Ｇの
階数がＮ・　（ｍ→−Ｌ−１）となる。従って（４１式
は ｆａｌ　　Ｌ　＝　Ｎ　・（ｍ−１）　／　（Ｓ　−Ｎ
　）のとき正規方程式 ％式％ 方程式 とすることができ、次のように解くことができる。

（正規方程式の場合） 行列Ｇの逆行列が存在するから（６）式を演算すればよ
い。

ＩＴ　＝　Ｇｌ；刊・Ｉ’Ｒ−−−−−−−−−−−−
−−−−〜−−−−−−−−−（６）（不定方程式の場
合） 行列Ｇの逆行列をそのまま求めることはできない。また
（４１　”式を満足するベクトルＢＩは唯一ではないが
、例えば１ノルム最小型一般逆行列」　（柳井他：“射
影行列、一般逆行列、特異値分解”。

東大出版会、　６５−６８（１９８３）　）を用いて、
１１（−［１；Ｔ　　・　（（Ｇ　　−（１；丁　）−
１・　Ｉｌ！−−−−−−−−−（６）’なる解を求め
ることができる。

即ちこの発明の最も主要な特徴は、音圧被制御点の数Ｎ
より多い数の制御用スピーカ素子及びフィルタを用いる
ことにより、音圧制御の為の基本式（４）式または（４
）式を満足するフィルタ係数が求められ、その結果各被
制御点それぞれについて定められるインパルス応答を、
８個のスピーカによって近似ではなく正確に再現できる
ことである。

従ってこの発明は、原理的には、「各被制御点において
近似的なインパルス応答しか再現できない」従来法の欠
点を解消するものである。

「実施例−１」 第２図はこの発明の実施例を示し、全系をディジタル信
号処理系で構成した例である。第２図において、多点放
射形音圧制御装置（この発明の装置）への入力信号Ｕ（
ｔ）は、入力端子１２よりＡ／Ｄ変換器３０に供給され
、Ａ／Ｄ変換器３０は入力信号Ｕ（ｔ）を離散化した信
号Ｕ　（ｋ）を出力し、これらは８個の同じ信号Ｕ＋　
（ｋ）（＋　＝　１　、　２、−、　Ｓ）に分岐して８
個のＦｌｌ？フィルタ２１．　、２１□、−、２１゜に
それぞれ供給され、これらフィルタ２１．　、２１□。

−−−−、２１３の出力はそれぞれＤ／Ａ変換器５０Ｉ
。

５０ｇ　、　−−−−、５０ｓを通じて８個のスピーカ
素子１、。

１□、−・−１，（Ｓ≧２）に供給される。

１つ又はそれ以上の数Ｎ個（’Ｎ≦５−１）の被制御点
１１．、−、　ＩＩＮ　として、この例ではＮ個のマイ
クロホン素子を用いており、これらマイクロホン素子１
１．、−１１．の出力はそれぞれ開閉器２００を経て、
８個のスピーカ素子１８．１□＋　’−”＋　　Ｉ　Ｓ
　とＮ個の被制御点１１１、−、１１．との間のＳ−Ｎ
個のインパルス応答ベクトルｇ、□ｊ　（ｉ　＝Ｌ　２
、−。

Ｓ；ｊ＝１．２．・−、Ｎ）を記憶する波形記録器６０
に供給される。

波形記録器６０の出力と、各被制御点１１．、−。

１１Ｎごとに定められる所望のＮ個のインパルス応答ベ
クトル（以後、目標インパルス応答ヘクトルと呼ぶ）を
記憶する波形記録器６１の出力とは、それぞれ波形操作
器７０に供給される。波形操作器７０の出力は、演算設
定器８０に供給される。演算設定器８０の出力はフィル
タ係数決定器１００に供給され、この出力は８個のＦＩ
Ｒフィルタ２１＋　、　２１ｚ　、−’。

２１、の対応する設定入力側に供給される。

次にこの実施例の動作を説明する。この実施例において
、８個のＦＩＲフィルタ２１．１　、２１ｚ　、−。

２１、は、最初入力信号Ｕ、（ｋ）（＋　＝　１　、　
２　、’　　−・。

Ｓ）を未処理のまま直ちに出力するような状態に設定さ
れているものとする。波形記録器６０には、スピーカ素
子ＩＩ、■□、−，１ｓ　と被制御点１１□−’　＋１
１Ｎ　との間のインパルス応答ベクトル（以後、室内イ
ンパルス応答ベクトルと呼ぶ）が、必要に応じて開閉器
２００を閉じることにより記憶され、又波形記録器６１
には各被制御点ＩＬ　、−、ｌｌｏごとに定められるＮ
個の目標インパルス応答ベクトルが、あらかじめ記憶さ
れているものとする。

さて波形操作器７０は、波形器記録６０及び６１にそれ
ぞれ記憶されている室内インパルス応答ペクトル及び目
標インパルス応答ベクトルを用いて、以下にｊホベる［
１］〜［４］の処理を行う。

［１］　室内インパルス応答ベクトル（ｇｉｊ（を−Ｌ
２．−、Ｓ　；　ｊ　＝　１．２．−　、Ｎ）。

ｇ　ｒＪ−［ｇ　ｒＪＩｇ　ｉ、＋□−ｇ　＝；、］　
”　　（ｍ”２ｍ）とする）のそれぞれについて、あら
かじめ決められた基準値■（■≧０）と室内インパルス
応答ベクトルｇｉｊの構成要素とを先頭から比較して、
ｇｒｊｐ　＞　Ｖ　（Ｐ　＝　１　、　２、−）　　　
　−−−−−１７１を満たす値ｇ、４２が最初にでるま
でを雑音と見做し、先頭からＰ−１までの要素を０と置
き換えるとともに、この間を群遅延Ｄ　ｉ　ｊとしてＤ
　、　、　＝　Ｐ　−１−−−−−−−−−−−−−−
・−−−−−−−−−−−−−−−−−−＋８１にセソ
１へする。

［２］　目標インパルス応答ベクトル（ｆＪ　（ｊ−１
，２，−、Ｎ）　、ｔｉ　＝　［ｒｓ＋＋　　ｒｊ２＋
　　１ｒｊｎ’　　］”　　（ｎ’　　≧ｎ）とする）
のそれぞれについて、［１］と同じ処理を行い、群遅延
Ｄ′４を求める。目標インパルス応答ベクトルはあらか
じめ計算で与えられる場合のみならず、例えばある音源
からの音響信号を打消したい場合にその音源からの各被
制御点との間のインパルス応答が目標インパルス応答ベ
クトルとなり、これを求める。

［３］　　　［１］で求めた群遅延り、□Ｊについて（
９）式の処理を行う。

Ｄ　ｓ＝　ＭＩＮ（Ｄ　Ｉｊ＋　　Ｄ　ｚｊ、−・−、
Ｄ　ｓＪ）　−−−−一−−−−−（９１ＭＩＮ（ａｔ
　ｌ　　ｚ　、−１ａｑ　）　　ｉ　ａＱ　（ｑ−１，
２＋一−−−−，＞のうち最小値となるａ９を選び出す
関数［４］　第ｊ（ｊ＝　１．２．−、Ｎ）番目の被制
御点に定められている目標インパルス応答ベクトルｆ、
と、室内インパルス応答ベクトルｇ内について、それぞ
れのインパルス応答ベクトルの先頭からり１個の要素を
削除する。つまり各インパルス応答ベクトルについて共
通となる同一の遅延量を除去する。但し、Ｄ゛、≧Ｄ、
を満だしていなければ、後述する制御において正確な目
標インパルス応答の再現はできない。それは［時間的に
遅れた信号でこれより進んだ信号を制御する」ことはで
きないという物理的な理由による。

以上に述べた［１］〜［４］の処理は、前述した（４）
°″式の階数が行列Ｇの行数で決まるようにする為の工
夫である。

次に波形操作器７０は、演算設定器８０に上記［１］〜
［４］の処理を受けた目標インパルス応答ベクトル及び
室内インパルス応答ベクトルを出力する。

ここで室内インパルス応答は、前に述べた２つの条件 （条件−１） 室内インパルス応答はすべて、長さｍのベクトルで表す
ことができる。

（条件−２） 室内インパルス応答をＺ変換して得られる伝達関数は、
互いに同し零点を持たない。

を満たしているものとする。

演算設定器８０は、波形操作器７０からの入力を用いて
、次の［６１〜［９］の処理を行う。

［６］　あらかじめ定められている室内インパルス応答
ベクトルの長さｍ、目標インパルス応答へクトルの長さ
ｎ、スピーカ素子の数Ｓ、被制御点の数Ｎを用い、以下
の処理に従って、Ｆｉｌ？フィルタ２１．　、２１□、
−、２１，の必要タップ長りを決める。

Ｌ≧Ｎ　・（ｍ　−１）　／　（Ｓ　−Ｎ）　　−−一
−−−−−−−−−−４０）ｎ　＝　ｍ　＋　Ｌ　−１
−−−−−−−−−−αυを同時に満足するタップ長り
を求める。但し、００１式において、 ｎ＜　ｍ　十Ｌ　−１−−−−−−’−’−−−−−−
’−’−−−−−−−αυ”となるししか求められない
場合には、目標インパルス応答ベクトルに適当な零づめ
を行い００式を満足するような長さのベクトルにつくり
変え、改めてこれを目標インパルス応答ヘクトルとする
。

［７１Ｎ個の目標インパルス応答ベクトルを用いて、（
４）”式左辺に示されているようなＮ・　（ｍ＋Ｌ−１
）の長さを持つ目標インパルス応答ベクトルＲをつくる
。

［８］　　Ｓ−Ｎ個の室内インパルス応答ベクトルを用
いて、（４）”式右辺にしめされているような（Ｎ・　
（ｍ＋Ｌ−１）、５−Ｌ）の行列のをっくる。

［９］　　　［６］で定めたＬについて、（Ｌ＝Ｎ−（
ｍ−１）／　（Ｓ−Ｎ）のとき）（Ｎ　・（ｍ十Ｌ　−
１）　、　Ｎ　・（ｍ＋Ｌ　−１））の単位行列■をつ
くる。

（Ｌ　＞Ｎ　・（ｍ　−１）　／　（Ｓ−Ｎ）のとき）
Ｇの転置行列ｃｒ、Ｔをつくる。

以上［６１〜［９］の処理を行った後、演算設定器８０
は、Ｉ’Ｒ，ＣＧ、　ＧＩＦＴ　　（又は■）をフィル
タ係数決定器１００に出力する。

フィルタ係数決定器１ｏｏは第３図に示すように、３つ
の行列乗算器１０１　、１０２　、１０３と、１つの逆
行列演算器１１０　と、係数分配器１２０　とで構成さ
れている。以下、第３図を参照してフィルタ係数決定器
１００の動作を説明する。

［１０］　　行列乗算器１０１　（７）入力端子１０１
−１１ｍ［１；、１０２−２にＣＴ　（又は■）が入力
され、これらの積が計算されてｃ−、ｃＴ　（又はＧ）
が出力される。

［１１］　　逆行列演算器１１０により、行列乗算器１
０１の出力から（（［；−（［；Ｔ）−１（又はＧ　−
１）が得られ、行列乗算器１０２へ出力される。

［１２コ　行列乗算器１０２の入力端子１０２−１にＧ
Ｔ（又は■）が入力され、出力としてｃＴ　　・　（Ｇ
・（ＣＴ）−１（又は（Ｇ　−１）かえられ、これは行
列乗算器１０３へ出力される。

［１３］　　行列乗算器１０３の入力端子１０３−１に
Ｒが入力され、その結果衣のようなフィルタ係数■が求
められる。

（＋　＝　１　、　２．−−−−−、　　Ｓ）［１４］
　　係数分配器１２０が、８個のフィルタ２１．。

２１□、−、２１３に対してＱ２）式で求めたフィルタ
係数Ｉｈ、　、　Ｉｈｚ　、−、Ｉｈｓを分配する。

以上述べた［１］〜［１４］の処理を行うことにより、
ＦＩＲフィルタ２１．　、２１□、−、２１Ｓの設定終
了後、この装置の入力端子１２に入力信号Ｕ（ｔ）を加
えれば、スピーカ素子１．．１□、−、１ｓから放射さ
れる音響信号は、Ｎ個の被制御点１１１゜−−−−、１
１，において、これらの点に対して定められた目標イン
パルス応答の持つ特性を完全に再現するものとなる。

ところでフィルタ係数の演算には、前述した［１］〜［
１４コの処理による以外に、種々の逐次近似法を用いる
こともできる。逐次近似法を用いる場合には、アルゴリ
ズムの収束性などに注意する必要もあるが、上述したよ
うな逆行列あるいは、ノルム最小型一般逆行列［０４式
で得られるｚ？　　・（（１；−（ＣＴ）−’のこと］
を直接求める方法に比べて、計算量やメモリ量などの面
で有利となる。

逐次近似法による場合は、第２図におけるフィルタ係数
決定器１００を、第４図に示すような逐次型フィルタ係
数決定器１３０に取り換えればよい。

この逐次型フィルタ係数決定器１３０は、逐次計算を受
は持つ逐次演算回路１３１　と係数分配器１２０　とよ
り構成されている。逐次演算回路１３１では、例えば、
次のような最急降下アルゴリズムで〔ファジェーエフ他
：線形代数の計算法、産業図書。

（１９７６）　）フィルタ係数を求めることができる。

（正規方程式又は、不定方程式の場合）ξ；アルゴリズ
ムの繰り返し回数を表すパラメータ β（ξ）；収束係数 （正規方程式となる場合） α（ξ）５収束係数 更にこの逐次型フィルタ係数決定器としては、第５図に
示すような、８個の逐次型フィルタ係数決定回路１４１
□１４１２．−、１４１ｓ　と誤差演算回路Ｌ０　とで
構成することも可能である。誤差演算回路Ｌ０は次式に
従って誤差ヘクトルを算出する。

ξ；アルゴリズムの繰り返し回数を表パラメータ［ｉｉ
１番目のスピーカとＮ個の制御点との間のインパルス応
答ベクトル ｇｉｊ（ｊ＝１．２．−、Ｎ）を用いて作る次のような
行列。

１ｈ、（ξ）；第１番目のスピーカ素子に接続されてい
るＦｌｌ？　フィルタの係数ベクトルそして、８個の逐
次型フィルタ係数決定回路１４１．。

−−−、１４１，は、次のようなアルゴリズムでフィル
タ係数を求める。

（正規方程式、不定方程式の場合） １ｈ、（ξ→−１）＝ｌｈ、（ξ）＋β（ξ）・ｇ＋Ｔ
−Ｅ（ξ）　　　　　　　　〜−−−−−−−−−−０
５１（正規方程式の場合に用いる。） Ｉｈ１（ξ＋１）−１ｈ、（ξ） ＋α（ξ）・ｓｔ　　（ξ＞　　　　−−−−−−−−
−−−−−−−−−−０５）゛（１４１式で求められる
Ｅ（ξ）を［ｅｌ　（ξ）　　ｌｌ！２（ξ）　−’−
ｅ　Ｎ、　Ｉｍ＋Ｌ−１１（ξ）］７と表せば、Ｎ−（
ｍ十Ｌ−１）　−３−Ｌであるから、Ｅ　’（ξ）を先
頭からＬ個づつ区切ることができる。

これを［ｆｆｉ、Ｔ（ξ）＠、”（ξ＞　−リＴ　（ξ
）］丁とすれば、ｅｉ　　（ξ）はこのベクトルの第１
番目の要素ベクトルである。

この方法は、８個の演算回路１４１．、−、１４１３が
並列に計算を行う為、フィルタ係数の決定に要する計算
量が、前記Ｑ３１．Ｑ３）”式を計算する場合に比べて
、見かけ上１／Ｓとなる。従って、フィルタ係数の算出
に要する時間を大幅に短縮するこができる。

〔実施例−２〕 既述した実施例−１は、８個（Ｓ≧２）のスピーカ素子
から放射される音響信号を、Ｎ個（Ｎ≦５−１）の被制
御点において、それぞれの被制御点ごとに定められた任
意の信号伝達特性（目標イン六ルス応答）を満足するよ
うに制御することにより、任意の音圧分布を実現するも
のであった。

しかし「Ｎ個の被制御点すべての音圧を零にする」とい
う制御を行う場合には、α２）式のＲが零ベクトルとな
ることからも分かるように、８個のＦＴ１１フィルタ２
１＋　、−、２１，のタップ係数がすべて零となる、実
用上無意味な制御をすることになる。

ところが、Ｎ個の被制御点の音圧はすべて零に保ち、か
つ被制御点以外の場所では音響信号が存在するような制
御は、特に、拡声会議システム等において、スピーカ素
子とマイクロホン素子との音響結合により引き起こされ
る「ハウリング」を防止する為の手段として、極めて要
求の高いものである。

そこで、実施例−１に示した装置に若干の改良を加える
ことによって上記の制御を可能にするようなこの発明の
別の構成例を第６図に示す。第６図において第２図と対
応する箇所には同一の番号を付けて重複説明を省略する
。第６図において、入力端子１２からの入力信号Ｕ（ｔ
）は、遅延回路１６０にも供給され、その出力はスピー
カ素子９に供給されている。

遅延回路１６０ば、スピーカ素子９からＮ個の被制御点
１１．、−、１１Ｎにそれぞれ放射される音響信号が、
８個のスピーカ素子１１．１□、−、Ｉｓから各被制御
点１１．、−、　ＩＩＮにそれぞれ放射される音響信号
のいずれよりも、遅れて到達するように適切な遅延量を
与える為のものである。これは言い換えれば、スピーカ
素子９とＮ個の被制御点１１１　、−−−−−、１１．
との距離が、スピーカ素子１．．１□。

−、ｉ　、と各被制御点１１゜、−、１１Ｎ　との各距
離より遠くなるように設置してあれば遅延回路１６０は
不要であることを意味する。遅延回路１６０を使用する
のは、「時間的に遅れた信号で、これより進んだ信号を
制御する」というような矛盾した状況を避ける為である
。

又波形記録器６２が追加されており、波形記録器６１と
ともに減算器１７０に接続され、減算器１７０の出力は
波形操作器７０に供給される。

次にこの実施例の動作説明をする。この実施例において
、波形記録器６２には、既にスピーカ素子９と各被制御
点１．１．、−、’　ＩＩＮ　との間のＮ個のインパル
ス応答ベクトルｒｓ　　（ｊ＝１．２．−・−９Ｎ）が
記憶されているものとする。又他の部分は実施例−１と
同様に初期設定されているものとする。

まず、波形記録器６１．６２に記憶されているインパル
ス応答ベクトルｒｊ　、ｒｊ　　（ｊ＝１．２．−。

Ｎ）は減算器１７０により、 シ・、＝（’ｊ−ｒｊ−−−−−−−−・−−−一〜−
一−−−−−０６１の操作をうけ、目標インパルス応答
ベクトルとして、Ｃ４が波形操作器７０に人力される。

０６）式において、ｒｊがすべて零ベクトルであるとす
れば、目標インパルス応答ベクトルは、 となり、これはちょうどスピーカ素子９から各被制御点
］、１１．−．１１ＮまでのＮ個のインパルス応答ベク
トルｒｊ　の位相を逆転したものとなっている。つまり
この処理によりインパルス応答へクトル７４の位相を逆
転したものを得る。

これ以後は実施例−１に述べた［１］〜［１４］とまっ
たく同じ処理を行って、８個のＦＩＲフィルタ２１．、
−、２１．の係数を決定すればよい。

ＦＩＲフィルタ２１．、−、２１．に対するフィルタ係
数の設定が完了した後、スピーカ素子１１，１□。

−，１，に入力端子１２より入力信号を与えれば、Ｎ個
の被制御点１１．、−、１１．における音響信号は面式
で与えられる目標インパルス応答を満足し、ちょうど、
スピーカ素子り、−，１，ｓがスピーカ素子９に対し逆
相で鳴っているかのような状況がつくりだされる。従っ
てスピーカ素子９とスピーカ素子１．、■□、−，１ｓ
に同し入力信号を与えれば、Ｎ個の被制御点ＩＬ　、−
、ＩＩＮにおける音圧は、ことごとく零となり、かつ、
被制御点１１．。

−−−、ＩＩＮ以外の場所には音響信号が存在すること
になる。

即ち実施例−１で述べた装置に常時有限振幅の音響信号
を放射できるように構成されたスピーカ素子を１つ又は
それ以上つけ加えることにより、音場内のＮ個の被制御
点ＩＬ　、−、１１Ｎの音圧だけを零にし、他の点では
十分な音量の音響信号を放射するような制御が可能とな
る。このことは、この発明装置がハウリングを防止する
手段の１つとしても利用できることを示している。

〔実施例−３〕 実施例−２では常時、有限振幅の音響信号を放射するよ
うに構成されたスピーカ素子９を１つ以上加えることに
より、この発明装置がノ＼ウリングを防止する１手法と
して有効であることを説明した。

実施例−３では、Ｎ個の被制御点に設置する各マイクロ
ホン素子の出力信号を零にすることによ　　　−リ、ハ
ウリングを抑圧するこの発明の他の実施例を説明する。

第７図にこの実施例を示す。第７図においで、第２図と
対応する箇所には同一の番号を付けである。

第７図では、新たに目標インパルス応答畳み込み回路１
８０がつけ加えられている。この回路１８０は、Ｎ個の
被制御点、即ちマイクロホン素子ごとに与えられる目標
インパルス応答ベクトルをＩ’ＩＲフィルタ等を用いて
表したものであり、実時間で時系列との畳み込み演算を
実行するものである。

又目標インパルス応答畳み込み回路１．８０の入力側に
は入力信号１１　（ｋ）が供給されるようになっており
、その出力は各マイクロホン１１１、−、１１Ｎの出力
とともに減算器１９０に供給されている。

次にこの実施例の動作説明を行う。まず、実施例−１に
示した［１］〜［］４Ｊの処理を行い８個のＦｉｌｌ　
フィルタ２１＋　、−、２１ｓに対するフィルタ係数の
設定がなされているものとする。このとき目標インパル
ス応答畳み込み回路１８０と、スピーカ素子１１，１□
、　−”、　　１　ｓ　とに同じ入力信号を加えれば、
各マイクロホン素子１１１　、−、　ｌＩＮの出力信号
は、目標インパルス応答を再現するものとなり、減算器
１９０には全く同じ信号が入力されることになる。従っ
て、減算器１９０の出力信号は零となり、実施例−２で
述べた発明装置と同様にハウリングを抑圧する効果が得
られることになる。

加えてこの実施例は、Ｎ個の被制御点１１１、−。

１１、における音圧を零にする必要は全くないから、空
間内（被制御点１１１、−、　ＩＩＮ　）にいる人に対
してはスピーカ素子り、−、１、からの音響信号を十分
な明瞭性で与え、即ち人の聞きよい環境となるように各
被制御点１１．、−、１１．を制御しながらハウリング
をも抑圧できるという極めて優れた特徴を有するもので
ある。このことは、この発明を用いた装置でなくてはで
きない大きな長所である。

なお上述の各側において、被制御点１１．、−。

１１Ｎ　とスピーカ素子１．、−、　１　、との間の各
インパルス応答は１度測定しておけば、これらの配置関
係１周囲の音響環境に変化がなければ、必ずしもマイク
ロホン素子を設置しておくことなく、被制御点に所望の
インパルス応答特性を持つ音響信号を得ることができる
。

夫１湛１ 次にこの発明装置の有効性を確認する為、第８図に示す
ようなスピーカ素子とマイクロホン素子の配置の下に行
った実験の結果を述べる。実験は、（実験諸元） スピーカ素子；直径１０　ｃｍのフルレンジスピーカ３
個、スピーカ間隔；１１ｃｍ、マイクロホン素子；１個
、スピーカ素子〜マイクロホン素子間最短距離；９０ｃ
ｍ、室容積；１．２．Ｍ、室内平均吸音率；０．８＋残
響時間；７５ｍ５ｅｃ　、目標インパルス応答ヘクトル
長；４００点、室内インパルス応答ヘクトル長；４００
点、インパルス応答測定方法、１０’ＯＯ回の同期加算
１周波数帯域、　３１．５〜３Ｌ０Ｈｚ　、サンプリン
グ周期；８ＫＨｚ、使用アルゴリズム；０５）式２００
回繰り返しＦＩＲフィルタ；３９９　タップ。

このような状況下でスピーカ素子９から放射された音を
スピーカ素子１１．１□を用いて消去するというもので
ある。従来のこの種の装置では音を消す為に用いるスピ
ーカ素子、即ち、音圧制御用のスピーカ素子は１つ用い
るのが一般的であったが、この実験では１点の音圧制御
の為に２つの制御用スピーカ素子１１．１□を用いてい
る。これは（５）式を満たすような装置を構成する為で
あり、この発明の実施例として必要な事項である。又第
９図に、スピーカ素子１１及びマイクロホン素子１１間
、スピーカ素子１□及びマイクロホン素子１１間、スピ
ーカ素子９及びマイクロホン素子１１間の、それぞれの
インパルス応答の周波数特性をそれぞれ曲線Ａ、Ｂ。

Ｃで示す。ここでは、スピーカ素子９とマイクロホン素
子ＩＩとの間のインパルス応答（曲線Ｃ）を目標インパ
ルス応答とした。

第１０図の周波数特性は、スピーカ素子ＩＩとスピーカ
素子１２とを用いて合成したインパルス応答と目標イン
パルス応答との違いを表す誤差波形（曲線Ｅ１□）を、
目標インパルス応答の周波数特性（曲線Ｃ）を基準に描
いたものである。又この図には、（２）式に示した従来
の手法をスピーカ素子１．。

スピーカ素子１２につきそれぞれ用いた実験で得た誤差
波形の周波数特性（それぞれ１０００タツプのＦＩＲフ
ィルタを用いた。）も曲線Ｅｌｌ　Ｂｇとして同時に描
いである。それぞれの誤差波形のパワーは目標インパル
ス応答のパワーに対し、従来法を用いた場合には、−１
１ｄＢ（Ｅｌ）、−１９ｄＢ（Ｅ２）、この発明の装置
を用いた場合は一２５６Ｂ（Ｂｒｚ）であった。

以上の結果よりこの発明装置によれば、マイクロホン素
子位置に対する音圧制御が従来法を上回る十分なパフォ
ーマンス（性能）を持って実現されることが実験によっ
ても確認された。特に第１０図はそのままハウリングマ
ージンをも示しており、この装置をハウリング防止に適
用することができ、かつ極めて大きなハウリングマージ
ンを得ることができることも明らかになった。又０５）
式のアルゴリズムが確認されたことによりフィルタ係数
算出にかかる計算量が、フィルタ係数決定器を並列に働
かせる方法により見かけ上の計算量を節約できることも
確認された。

「発明の効果」 以上説明したようにこの発明装置は、複数のスピーカ素
子それぞれの入力信号に対して、適切な制御を行うフィ
ルタの係数を厳密に求め、上記スピーカ素子数未満の被
制御点に与える音圧を任意の目標値に設定することを可
能とするものであり、この装置が実用上十分な性能を発
揮することも実験によって確認された。又この発明装置
は、目標インパルス応答にあたる信号を自由に設定でき
る為、 （１）　　この装置の設置されている三次元空間内にあ
る他の音響装置による音圧分布以外に、その空間とは異
なった空間（例えば他の部屋）における音圧分布も再現
できる、 （２）　　全く新しい音圧分布を創造できる、（３）　
　この装置と同一空間にあり、目標インパルス応答を与
えるスピーカ素子に供給する信号をこの装置の入力信号
とは逆相の信号とすることにより、任意の３次元空間に
設定する１つ以上のマイクロホン素子の音圧をほぼ零に
保つことができる。従って、有効にハウリングを防止す
る一手法を提供することができる、 という特徴を有するものである。

尚、この発明の装置を構成するにあたり、スピーカ素子
として複数のスピーカユニット（例えば、ウーハ−とツ
イータ）を用いること、被制御点として用いるマイクロ
ホン素子として、複数のマイクロホンユニットの組み合
わせ（例えば、指向性マイクロホン）を用いることは、
この発明の装置の性能向上に有用であると考えられる。

又、この発明の信号処理系は、以上に述べた構成に基づ
く演算、遅延手順などを行うＤＳＰあるいは、マイクロ
プロセッサ等で構成することも可能である。加えてこの
明細書では、フィルタとしてＦＩＲフィルタを用いるこ
とを前提として説明したが、ＦＩＲフィルタに替えてＩ
ＩＲフィルタ等を用いることも十分に可能である。

【図面の簡単な説明】

第１図はこの発明の装置の原理を説明するための図、第
２図はこの発明装置の一実施例を示すブロック図、第３
図はフィルタ係数決定器１００の一実施例を示すブロッ
ク図、第４図は逐次近似法でフィルタ係数を計算する系
の一実施例を示すブロック図、第５図は複数のフィルタ
係数決定回路によって構成することにより、第４図のフ
ィルタ係数決定器より見かけ上の計算量を減らした、逐
次近似法でフィルタ係数を算出する系の一実施例を示す
ブロック図、第６図は被制御点の音圧をすべて零とし、
かつ被制御点以外の場所には十分な音量の音響信号を放
射するようにしたこの発明装置の一実施例を示すブロッ
ク図、第７図は被制御点におかれた受音素子の出力信号
はすべて零にし、被制御点の位置に対しては任意の音圧
分布を与えることを可能としたこの発明の一実施例を示
すブロック図、第８図は実験におけるスピーカ素子とマ
イクロホン素子との位置関係を示す図、第９図は実験に
用いた各スピーカ素子とマイクロホン素子との間のイン
パルス応答の周波数特性を示す図、第１０図はこの発明
に属する装置及び従来技術を用いた場合の音圧制御性能
の周波数特性を表す実験結果の図、第１１図は従来の音
圧制御装置の原理を説明するための図である。 １、　、１２．−、　　Ｉ　Ｓ　　：　２つ以上のスピ
ーカ素子、９：有限振幅信号を放射するように構成され
たスピーカ素子、１１．、−、１１８　　：　１つ又は
複数個の被制御点、２１．　、２１□、−、２１，：　
ＦＩＲフィルタ、３（１：Ａ／Ｄ変換器、５０：Ｄ／Ａ
変換器、６０、６１．６２．　　：波形記録器、７０：
波形操作器、８０：演算設定器、１００：フィルタ係数
決定器、１０１　、１０２　、１０３　　：行列乗算器
、１１０：逆行列演算器、１２０：係数分配器、１３０
：逐次型演算回路、１３１：逐次演算回路、１４１．、
１４１ｇ、−。 １４１、　　：逐次型フィルタ係数決定回路、Ｌ０：誤
差演算回路、１６０：遅延回路、１７０　、１９０　　
：減算器、１８０：目標インパルス応答前に込み回路、
２００：開閉器。 特許出願人　　日本電信電話公社 代　　理　　人　　　　草　　野　　　　　　　卓手続
補正書（自発） 昭和６０年９月２６日 １事件の表示　　特願昭６０−５３８８６２発明の名称
　多点放射形音圧制御装置３補正をする者 事件との関係　　特許出願人 日本電信電話株式会社 ４代　理　人　　東京都新宿区新宿４−２−２１　　相
撲ビル５補正の対象　　明細暑中発明の名称、特許請求
の範囲の欄およ６補正の内容 （１）　　発明の名称「ン点放射形音圧制御装置」を「
多チヤンネル制御装置」と訂正する。 （２、特許請求の範囲を別紙のとおり訂正する。 （３）　　明細書３頁Ｌ〜２０行「この発明はスピーカ
素子より・・・・・・・関するものである。」を下記の
とおり訂正する。 「二の発明は有限ど見なすことの出来るインパルス応答
を持ち、原理的に多入力多出力である線形系に対し、複
数の出力素子に加えられる入力信号を制御することによ
り、七の線形系が有する１つ、または２つ以」二の被制
御点における出力信号の１云達特性を所望の特性とする
多チヤンネル制御装置に関するものである。 上記線形系の例としては、複数のスピーカ素子と１つ以
上の受音点から成る音響系や、複数の送信用アンテナと
１つ以上の受信用アンテナとで構成される電波１云搬系
などがある。 上記音響系にこの発明を適用する場合には、上記複数の
出力素子が複数のスピーカ素子に対応し、上記被制御点
が受音点に対応する多チャンイ・ル制御装置となり、そ
の受音点（以下、被制御点と呼ぶ）において、任意の音
圧分布を実現する「音圧制御」を行うものとなる。 以下この明細書では」１記音圧制御を例としてこの発明
の説明を行うが、曲の系についても」１記複数の入力端
及び被制御点との対応関係をこの音響系ど同嘩に定める
ことにより、全く同じ説明が出来る。」 （４）　　同書５頁９行「スピーカ（・目素子１１へ」
を「スピーカ素子１□へ」と訂正する。 （５）　　同書５頁１７行「音場制御が実現される。」
を「音圧制御が実現される。」と訂正する。 （６）　　同書７頁１９行「二の発明の目的は、任意の
」を１この発明の目的は、音響系を対象とした場合で述
べると、任意の」と訂正する。 （７）　　同書８頁５行「音圧制御装置を提供すること
にある。」を「多チヤンネル制御装置を提供することに
ある。」と訂正する。 （８）同壱８頁７〜８行「この発明によれば・・・・そ
れぞれ各」を「二の発明は、音響系を例として述べると
、空間の複数の点にスピーカ素子が配置され、これらス
ピーカ素子にそれぞれ個」と訂正する。 （９）同書８頁Ｌ行「ヒの係数を対応フィ」を１その係
数を対応するフィ」と訂正する。 （１０）　　同書９頁５行「のスピーカ素子」を「の出
力素子であるスピーカ素子」と訂正する。 （１１）　　同書１３頁６〜１１行「音圧被制御点の数
・・・・・・８個のスピーカによ」を「被制御点の数Ｎ
より多い数の出力素子（この場合はスピーカ素子）及び
フィルタを用いることにより、出力制御の為の基本式（
４）式または（４）′式を満足す”るフィルタ係数が求
められ、その結果各被制御点それぞれについて定められ
るインパルス応答を、８個の出力素子からの信号によ」
と訂正する。 （１２）　　同書１３頁１９行「多点放射形音圧制御装
置」を「多チヤンネル制御装置」と訂正する。 （１３）　　同書１４頁７行「８個のスピーカ素子１１
．」を１Ｓ個の出力素子どしてのスピーカ素子１１．」
と訂正する。 （Ｌ）　　同書１７頁３〜４行「ベクトルとなり、これ
を求める。」を「くクトルとなる。」と訂正する。 （１６）　　同書１８頁１８行、１９行、２０頁２行、
８行［’　６　’ｌ］Ｊをｒ［５〕Ｊとそれぞれ訂正す
る。 （１７）　　同１１：１８頁１８行、２０頁２行、８行
「

〔９〕」をｒ［８］１とそれぞれ訂正する。 （１８）同書１９頁１４行［７１Ｊを［６’：］Ｊと訂
正する。 （１９）同書１９頁１８行「〔８〕」を「〔７〕」と訂
正する。 （２０）同書２０頁１６行［１０１Ｊを「

〔９〕」と訂
正する。 （２１）同書２０頁１９行［〔１１〕」を［１０］Ｊと
訂正する。 （２２）同書２１頁６行「〔１２〕」をｒｃｌｌｌＪと
訂正する。 （２３）同書２１頁６行［１：Ｄ、ｊをｒ、［］ｌ）Ｊ
と訂正する。 （２４）同書２１頁１３行、１ｇ行、２２頁５行、２８
頁１７行、３１頁４行［１，４］Ｊを［〔１３，１）Ｊ
とそれぞれ訂正する。 （２５）　　同書３１頁１９行「明瞭性で与え、」を１
明瞭性を確保して与え、」と訂正する。 （２６）　　同書３５頁９〜１７行「嵐」−説明したよ
うに・・・・・・・・・自由に設定できる為、」を下記
のとおり訂正する。 １以上説明したようにこの発明装置は、多入力多出力線
形系の複数の出力素子（即ち実施例におけるスピーカ素
子）のそれぞれの入力信号に対して、適切な制御を行う
フィルタの係数を厳密に求め、上記出力素子数未満の被
制御点からの出力信号の１云達特性（実施例における音
圧）を任意の目標値に設定することを可能とするもので
あり、この装置が実用上十分な性能を発揮することも音
圧制御を行う実験によって確認された。又この発明装置
は、目標インパルス応答にあたる信号を自由に設定でき
る為、特に音圧制御系に用いる場合には、」 （２７）　　同書３６頁１１行［尚、この発明の装置を
構成するにあたり、」を「尚、この発明の装置を音圧制
御に用いる場合に、」と訂正する。 （２８）　　同書３６頁１７行［有用であると考えられ
る。又、この」を「有用であると考えられる。 父、この」と訂正する。 （２９）　　同書３７頁４行「分に可能である。」の次
に下記を加入する。 「：の明細Ｈ＃では、音圧制御を例にこの発明の説明を
行った。この発明は有限と見なすことの出来るインパル
ス応答を有し、原理的に多入力＆出力である線形系を制
御する為のものであり、既に説明した音圧制御系におい
ては複数の出力素子が複数のスピーカ素子で、さらに１
つ以」二の被制御点が受音素子として具現化されていた
。 従ってこれと同様の（也の多入力多出力線形系について
も、出力素子及び被制御点を適当に具体化することによ
り、音圧制御系について述べた事柄と全く同様の説明が
出来る。即ちこの発明によれば上記線形系に対して、出
力素子数未満の被制御点における信号の特性を、任意に
かつ厳密に制御することが可能である。」 以上 特許請求の範囲 （１）　　多入力多出力である線形系に接続する複数の
出力素子と、これら出力素子のそれぞれに出力を印加し
、同一の入力信号が供給される複数のフィルタと、 」１記線形系に接続する上記複数の出力素子数未満の数
の素子（以後、被制御点と呼ぶ）についてそれぞれ定め
られる所望の信号伝達特性と、上記各出力素子と上記各
被制御点との間の信号伝達特性とを用いた演算により、
上記各フィルタの係数を算出して設定する手段とを備え
、上記同一の入力信号が、上記各被制御点において、各
被制御点に与えられた上記所望の信号伝達特性を満足す
るように制御する多チヤンネル制御装置。 （２）　上記複数の出力素子及び上記複数のフィルタを
それぞれＳ個とし、上記被制御点の数をＮ個とし、かつ
上記出力素子と被制御点との間の各信号伝達特性のそれ
ぞれをｍ個の離散信号として取り扱うように構成され、
上記Ｓ個のフィルタがそれぞれ、 Ｌ≧Ｎ・（ｍ−１）／（Ｓ−Ｎ） により与えられるＬ個のタップを持つ離散形ＦＩＲフィ
ルタにより構成されていることを特徴とする特許請求の
範囲第１項記載の多テヤンイ・ル制御装置。 （３）上記複数のフィルタの出力素子どしてＳ個のスピ
ーカ素子を用い、Ｓ個のスピーカ素子とＮ個の被制御点
から成る音響系を制御することを特徴とする特許請求の
範囲第１項記載の多チヤンネル制御装置。

Claims (4)

【特許請求の範囲】
1. （１）空間の複数の点に配置されたスピーカ素子と、こ
   れらスピーカ素子のそれぞれに出力を印加し、同一の入
   力信号が供給される複数のフィルタと、 上記複数のスピーカ素子と同一空間に設定され、上記ス
   ピーカ素子数未満の数の点（以後、被制御点と呼ぶ）に
   ついてそれぞれ定められる所望の信号伝達特性と、上記
   各スピーカ素子及び上記各被制御点間の信号伝達特性と
   を用いた演算により、上記各フィルタの係数を算出して
   設定する手段とを備え、 上記同一の入力信号が、上記各被制御点において、各被
   制御点に与えられた上記所望の信号伝達特性を満足する
   ように制御する多点放射形音圧制御装置。
2. （２）上記複数のスピーカ素子及び上記フィルタをそれ
   ぞれＳ個とし、上記被制御点の数をＮ個とし、かつ上記
   スピーカ素子及び上記被制御点間の各信号伝達特性のそ
   れぞれをｍ個の離散信号として取り扱うように構成され
   、上記Ｓ個のフィルタはそれぞれ、 Ｌ≧Ｎ・（ｍ−１）／（Ｓ−Ｎ） により与えられるＬ個のタップを持つ、離散形ＦＩＲフ
   ィルタにより構成されていることを特徴とする特許請求
   の範囲第１項記載の多点放射形音圧制御装置。
3. （３）Ｌ＝Ｎ・（ｍ−１）／（Ｓ−Ｎ） なるＬ個のタップをもつ上記離散形ＦＩＲフィルタに対
   し、「正方行列」の逆行列演算を用い、Ｌ＝〔Ｎ・（ｍ
   −１）／（Ｓ−Ｎ）〕＋ｉ （ｉ；１、２、・・・・・） なるＬ個のタップをもつ上記離散形ＦＩＲフィルタに対
   し、「ノルム最小型一般逆行列」等を用いて、上記各フ
   ィルタ係数の算出として上記各タップ係数を決定するこ
   とを特徴とする特許請求の範囲第２項記載の多点放射形
   音圧制御装置。
4. （４）上記各被制御点にマイクロホン素子が用いられ、 上記同一の入力信号が供給され、有限振幅信号を放射す
   るように構成された少なくとも１つのスピーカ素子が更
   に設けられ、 上記複数のスピーカ素子及び上記有限振幅信号を放射す
   るように構成されたスピーカ素子から放射され、上記マ
   イクロホン素子にそれぞれ受音される音響信号を、各マ
   イクロホン素子ごとに定められる任意の信号伝達特性を
   満足するように制御することを特徴とする特許請求の範
   囲第１項記載の多点放射形音圧制御装置。