JPH0654883B2 - 多入力形制御装置 - Google Patents
多入力形制御装置Info
- Publication number
- JPH0654883B2 JPH0654883B2 JP3265986A JP3265986A JPH0654883B2 JP H0654883 B2 JPH0654883 B2 JP H0654883B2 JP 3265986 A JP3265986 A JP 3265986A JP 3265986 A JP3265986 A JP 3265986A JP H0654883 B2 JPH0654883 B2 JP H0654883B2
- Authority
- JP
- Japan
- Prior art keywords
- output
- input
- filter
- fir
- impulse response
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Description
【発明の詳細な説明】 「産業上の利用分野」 この発明はｍ入力点とｎ出力点（ｍ＜ｎ）との間の各イ
ンパルス応答が有限とみなせる線形系（ｍ入力ｎ出力線
形ＦＩＲ系と記す）を対象として、 (a)その線形ＦＩＲ系の入力信号を推定する (b)その線形ＦＩＲ系のその入力信号に所望の特性を付
与して又は付与した後取り出す などの動作で代表される逆フィルタ処理を実現する多入
力形制御装置に関するものである。
ンパルス応答が有限とみなせる線形系（ｍ入力ｎ出力線
形ＦＩＲ系と記す）を対象として、 (a)その線形ＦＩＲ系の入力信号を推定する (b)その線形ＦＩＲ系のその入力信号に所望の特性を付
与して又は付与した後取り出す などの動作で代表される逆フィルタ処理を実現する多入
力形制御装置に関するものである。
この発明の多入力形制御装置は、上記の逆フィルタ処理
の特徴を生かすことに依り、次に述べる分野への応用を
可能とする。
の特徴を生かすことに依り、次に述べる分野への応用を
可能とする。
○音響系への応用 この発明装置における受信素子としてマイクロホン素子
を用い、 ｍ個の音源（入力点）から任意の１つの信号を選択する
ことに依る雑音除去（選択受聴）を可能とし、また 音源から放射される直接音のみを受音することに依る残
響除去を可能とする。
を用い、 ｍ個の音源（入力点）から任意の１つの信号を選択する
ことに依る雑音除去（選択受聴）を可能とし、また 音源から放射される直接音のみを受音することに依る残
響除去を可能とする。
○電波系への応用 この発明装置における受信素子として受信用アンテナ素
子を用い、 希望局から放射される信号のみを選択することに依る妨
害波除去を可能とし、 希望局からの直接波のみを受信することに依るゴースト
受信防止、又はマルチパス受信防止を可能とする。
子を用い、 希望局から放射される信号のみを選択することに依る妨
害波除去を可能とし、 希望局からの直接波のみを受信することに依るゴースト
受信防止、又はマルチパス受信防止を可能とする。
「従来の技術」 １入力１出力線形ＦＩＲ系の逆フィルタ処理を例とし
て、第１２図を参照して従来技術の説明を行う。対象と
する１入力１出力線形ＦＩＲ系１０１の出力点に受信素
子１１が設けられ、受信素子１１の出力は逆フィルタ処
理を行う為のフィルタ11０へ供給される。以下の説明で
は離散化した信号を用いる。尚ｍ入力ｎ出力線形ＦＩＲ
系に対しても以下と全く同様の説明が行える。
て、第１２図を参照して従来技術の説明を行う。対象と
する１入力１出力線形ＦＩＲ系１０１の出力点に受信素
子１１が設けられ、受信素子１１の出力は逆フィルタ処
理を行う為のフィルタ11０へ供給される。以下の説明で
は離散化した信号を用いる。尚ｍ入力ｎ出力線形ＦＩＲ
系に対しても以下と全く同様の説明が行える。
線形ＦＩＲ系１０１への入力信号x(k)(k＝1,2…）とフ
ィルタ11０の出力y(k)を等しくするような逆フィルタ処
理（入力信号推定）を考えることにする。説明を簡単に
する為、線形ＦＩＲ系１０１には遅延（応答の遅れ）は
無いものとする。
ィルタ11０の出力y(k)を等しくするような逆フィルタ処
理（入力信号推定）を考えることにする。説明を簡単に
する為、線形ＦＩＲ系１０１には遅延（応答の遅れ）は
無いものとする。
さて上記逆フィルタ処理を実現するフィルタ11０の係数
は、 g1(k)：線形ＦＩＲ系１０１のインパルス応答 h1(k)：フィルタ11０の係数 を満足するh1(k)で与えられる。例えば線形ＦＩＲ系１
０１が一般室内音場である場合に、室内の１点に置かれ
た信号源（スピーカ、人の口等、線形ＦＩＲ系１０１の
入力点に相当する）から音響信号（入力信号）ｘ（ｋ）
が放射されれば、該室内の１点（出力点）に配されたマ
イクロホン素子（受信素子）１１は上記信号源から該マ
イクロホン素子への直接音の他に、反射による残響音も
受音してしまう。この時室内インパルス応答（線形ＦＩ
Ｒ系１０１のインパルス応答）g1(k)の影響を打ち消す
ように、即ち がｋ＝１で１、ｋ＝２，３，…で０となるようにフィル
タ係数h1(x)が決められると、フィルタ11０の出力信号y
(k)は入力信号x(k)と等しくなり、残響信号が抑圧され
た信号が得られる。ところが、線形ＦＩＲ系１０１が非
最小位相推移系、つまりこの線形系の伝達関数G1(z)の
零点がｚ平面上の単位円外にも存在する場合、（反射の
存在する一般室内音場等は典型的な例である）には(2)
式をｚ変換して得られる次の関係式 H(z)＝１／G(z) ……(3) を満たすH(z)が不安定となり、上記逆フィルタ処理は実
現不可能であることが知られている(S.T.Neely and J.
B.Allen,“Invertibility of a Room Impulse Respons
e”J.Acoust.Soc.Am,66(1),pp.163-169,July,1979)。
は、 g1(k)：線形ＦＩＲ系１０１のインパルス応答 h1(k)：フィルタ11０の係数 を満足するh1(k)で与えられる。例えば線形ＦＩＲ系１
０１が一般室内音場である場合に、室内の１点に置かれ
た信号源（スピーカ、人の口等、線形ＦＩＲ系１０１の
入力点に相当する）から音響信号（入力信号）ｘ（ｋ）
が放射されれば、該室内の１点（出力点）に配されたマ
イクロホン素子（受信素子）１１は上記信号源から該マ
イクロホン素子への直接音の他に、反射による残響音も
受音してしまう。この時室内インパルス応答（線形ＦＩ
Ｒ系１０１のインパルス応答）g1(k)の影響を打ち消す
ように、即ち がｋ＝１で１、ｋ＝２，３，…で０となるようにフィル
タ係数h1(x)が決められると、フィルタ11０の出力信号y
(k)は入力信号x(k)と等しくなり、残響信号が抑圧され
た信号が得られる。ところが、線形ＦＩＲ系１０１が非
最小位相推移系、つまりこの線形系の伝達関数G1(z)の
零点がｚ平面上の単位円外にも存在する場合、（反射の
存在する一般室内音場等は典型的な例である）には(2)
式をｚ変換して得られる次の関係式 H(z)＝１／G(z) ……(3) を満たすH(z)が不安定となり、上記逆フィルタ処理は実
現不可能であることが知られている(S.T.Neely and J.
B.Allen,“Invertibility of a Room Impulse Respons
e”J.Acoust.Soc.Am,66(1),pp.163-169,July,1979)。
そこで従来の技術では次の評価量 を最小にする係数 を持つフィルタを安定かつ構成の簡単なＦＩＲ(Finite
Impulse Response)型で実現する方法がとられていた。
しかしこの手法は以下に示す欠点を原理的に有してい
る。
Impulse Response)型で実現する方法がとられていた。
しかしこの手法は以下に示す欠点を原理的に有してい
る。
(a)先に述べたように(2)式を満たす安定なフィルタは実
現出来ない。従って前記従来の手法で求まるＦＩＲフィ
ルタは を零にすることはできず、正確な逆フィルタ処理はでき
ない。
現出来ない。従って前記従来の手法で求まるＦＩＲフィ
ルタは を零にすることはできず、正確な逆フィルタ処理はでき
ない。
(b) の大きさやe(k)の特性は、線形ＦＩＲ系のインパルス応
答g1(k)に依存する為、上記従来技術で達成される性能
は対象とする線形ＦＩＲ系毎に大きく変動する。
答g1(k)に依存する為、上記従来技術で達成される性能
は対象とする線形ＦＩＲ系毎に大きく変動する。
この発明の目的は、ｍ入力ｎ出力線形ＦＩＲ系に加えら
れる任意の入力信号を正確に推定すること、 そのｍ入力ｎ出力線形ＦＩＲ系に加えられる入力信号に
所望の特性（インパルス応答）を与えた後に取り出すこ
と、 の逆フィルタ処理を実現する多入力形制御装置を提供す
ることにある。
れる任意の入力信号を正確に推定すること、 そのｍ入力ｎ出力線形ＦＩＲ系に加えられる入力信号に
所望の特性（インパルス応答）を与えた後に取り出すこ
と、 の逆フィルタ処理を実現する多入力形制御装置を提供す
ることにある。
この発明によれば入力点より出力点が多いｍ入力ｎ出力
線形ＦＩＲ系（ｍ，ｎは１以上の整数ｎ＞ｍ）を対象と
し、そのｎ個の出力点に受信素子がそれぞれ配され、こ
れらｎ個の受信素子にＦＩＲフィルタが接続され、その
ｎ個のＦＩＲフィルタの出力は加算手段で加算される。
ｍ個の入力点の第ｉ番目（ｉ＝１，２…，ｍ）とｎ個の
出力点の第ｊ番目（ｊ＝１，２…，ｎ）との間のインパ
ルス応答gij(k)をWij個の離散信号で表わし、上記第ｉ
番目の入力点と上記加算手段の出力との間の上記所望の
伝達特性ri(k)をPi個の離散信号で表すと、上記ｊ番目
の出力点に接続された第ｊ番目のＦＩＲフィルタは 但し、ｉ＝１，２，…，ｍ ｊ＝１，２，…，ｎ を満足するタップ数ｌｊを有し、かつその第ｊ番目のＦ
ＩＲフィルタの係数hj(k)が 但し、ｉ＝１，２，…，ｍ を満足しているものである。
線形ＦＩＲ系（ｍ，ｎは１以上の整数ｎ＞ｍ）を対象と
し、そのｎ個の出力点に受信素子がそれぞれ配され、こ
れらｎ個の受信素子にＦＩＲフィルタが接続され、その
ｎ個のＦＩＲフィルタの出力は加算手段で加算される。
ｍ個の入力点の第ｉ番目（ｉ＝１，２…，ｍ）とｎ個の
出力点の第ｊ番目（ｊ＝１，２…，ｎ）との間のインパ
ルス応答gij(k)をWij個の離散信号で表わし、上記第ｉ
番目の入力点と上記加算手段の出力との間の上記所望の
伝達特性ri(k)をPi個の離散信号で表すと、上記ｊ番目
の出力点に接続された第ｊ番目のＦＩＲフィルタは 但し、ｉ＝１，２，…，ｍ ｊ＝１，２，…，ｎ を満足するタップ数ｌｊを有し、かつその第ｊ番目のＦ
ＩＲフィルタの係数hj(k)が 但し、ｉ＝１，２，…，ｍ を満足しているものである。
発明の基本原理 第１図を用いてこの発明の基本原理を説明する。第１図
は第１２図に示した１入力１出力線形ＦＩＲ系１０１に
この発明を適用した場合でこの線形ＦＩＲ系１０１に信
号伝達経路１１１を並列接続して１入力２出力線形ＦＩ
Ｒ系１２１とし、この線形ＦＩＲ系１２１の２つの出力
点に受信素子１１，１２がそれぞれ配され、これら受信
素子１１，１２はＦＩＲフィルタ11１，11２にそれぞれ
接続され、これらＦＩＲフィルタ11１，11２の各出力は
加算器２２で加算される。この第１図はこの発明装置の
最小構成を表わしている。
は第１２図に示した１入力１出力線形ＦＩＲ系１０１に
この発明を適用した場合でこの線形ＦＩＲ系１０１に信
号伝達経路１１１を並列接続して１入力２出力線形ＦＩ
Ｒ系１２１とし、この線形ＦＩＲ系１２１の２つの出力
点に受信素子１１，１２がそれぞれ配され、これら受信
素子１１，１２はＦＩＲフィルタ11１，11２にそれぞれ
接続され、これらＦＩＲフィルタ11１，11２の各出力は
加算器２２で加算される。この第１図はこの発明装置の
最小構成を表わしている。
説明を簡単にする為、線形系１０１及び信号伝達経路１
１１には遅延がないものとする。この第１図において線
形ＦＩＲ系１２１への入力信号x(k)と加算器２２の出力
y(k)を等しくするような逆フィルタ処理を実現させる為
には、 g1j(k)(j＝１，２）：入力点と第ｊ番出力点 との間のインパルス応答 ｈｊ（ｋ）（ｊ＝１，２）：ＦＩＲフィルタ11ｊ （ｊ＝１，２）の係数 を満足する係数h1(k),h2(k)を求めれば良い。(5)式を満
足する係数h1(k),h2(k)を有するＦＩＲフィルタ11１，1
1２は、 （条件１）g11(k),g12(k)の伝達関数G11(z),G12(z)は共
通の零点を持たない。
１１には遅延がないものとする。この第１図において線
形ＦＩＲ系１２１への入力信号x(k)と加算器２２の出力
y(k)を等しくするような逆フィルタ処理を実現させる為
には、 g1j(k)(j＝１，２）：入力点と第ｊ番出力点 との間のインパルス応答 ｈｊ（ｋ）（ｊ＝１，２）：ＦＩＲフィルタ11ｊ （ｊ＝１，２）の係数 を満足する係数h1(k),h2(k)を求めれば良い。(5)式を満
足する係数h1(k),h2(k)を有するＦＩＲフィルタ11１，1
1２は、 （条件１）g11(k),g12(k)の伝達関数G11(z),G12(z)は共
通の零点を持たない。
なる条件下で (a)フィルタ係数h1(k),h2(k)をそれぞれｌ１，ｌ２個の
離散信号で表わすことにすれば、 ｌ１w12−１ ｌ２w11−１ …(6a) 但し、インパルス応答g11(k),g12(k)がそれぞれw11,w12
個の離散信号で表わされているものとする なる関係を満足する(5)式の解h1(k),h2(k)は唯一組存在
する。
離散信号で表わすことにすれば、 ｌ１w12−１ ｌ２w11−１ …(6a) 但し、インパルス応答g11(k),g12(k)がそれぞれw11,w12
個の離散信号で表わされているものとする なる関係を満足する(5)式の解h1(k),h2(k)は唯一組存在
する。
(b)ｌ１＞w12−１ ｌ２＞w11−１ …(6b) なる関係を満足する(5)式の解h1(k),h2(k)は無数に存在
する。
する。
のように実現され得ることを示すことが出来る。即ち
(5)式をｚ変換すれば １＝G11(z)H1(z)+G12(z)H2(z) …(7) G11(z),G12(z),H1(z),H2(z)：ｚ多項式 となる。この式においてG11(z),G12(z)が互いに素、つ
まりG11(z),G12(z)の零点が一致しないならば、ユーク
リッドの互除法を用いて、 である。
(5)式をｚ変換すれば １＝G11(z)H1(z)+G12(z)H2(z) …(7) G11(z),G12(z),H1(z),H2(z)：ｚ多項式 となる。この式においてG11(z),G12(z)が互いに素、つ
まりG11(z),G12(z)の零点が一致しないならば、ユーク
リッドの互除法を用いて、 である。
K(z)：任意のｚ多項式 と表わされる(7)式の一般解H1(z),H2(z)が存在する
（“理工学のための数学ハンドブック”丸善，PP,10-1
1,1983）。
（“理工学のための数学ハンドブック”丸善，PP,10-1
1,1983）。
ここで上記 以外にも deg H1′(z)＜deg G12(z)， deg H2′(z)＜deg G11(z) …(10) を満足する(7)式の解H1′(z),H2′(z)が存在するものと
仮定する。(8)式より、H1′(z),H2′(z)は次のように表
わされる。
仮定する。(8)式より、H1′(z),H2′(z)は次のように表
わされる。
ここでH1′(z),H2′(z)の次数は となり、これは上記の仮定と矛盾する。つまり他に解が
存在するとしたこの仮定は誤ったものであって(9)式を
満足する(7)式の解は唯一組である。
存在するとしたこの仮定は誤ったものであって(9)式を
満足する(7)式の解は唯一組である。
即ち、(a),(b)で述べた様な係数h1(z),h2(z)を持つＦＩ
Ｒフィルタ11１，11２の存在性は保証された。
Ｒフィルタ11１，11２の存在性は保証された。
次にＦＩＲフィルタ11１，11２の実現方法、即ちフィル
タ係数h1(k)h2(k)の求め方について説明する。インパル
ス応答g11(k),g12(k)をベクトル フィルタ係数h1(k),h2(k)をベクトル とすれば(5)式は、 但し、 w11＋ｌ１−１＝w12＋ｌ２−１ と書き換えることが出来る。前記(a),(b)の論議より、
(7)式を満足するフィルタ係数 の存在性は（条件１）の下に保証されているから、(7)
式は以下のように解くことが出来る。
タ係数h1(k)h2(k)の求め方について説明する。インパル
ス応答g11(k),g12(k)をベクトル フィルタ係数h1(k),h2(k)をベクトル とすれば(5)式は、 但し、 w11＋ｌ１−１＝w12＋ｌ２−１ と書き換えることが出来る。前記(a),(b)の論議より、
(7)式を満足するフィルタ係数 の存在性は（条件１）の下に保証されているから、(7)
式は以下のように解くことが出来る。
(a)′フィルタ係数h1(k),h2(k)の個数ｌ１，ｌ２が ｌ１w12−１ ｌ２w11−１ を満足する解を求める場合には、 を定める事で、(13)式右辺の畳み込み行列 が正列行列となるから、 とフィルタ係数を求める事が出来る。
(b)′ｌ１＞w12−１ ｌ２＞w11−１ を満足する解を求める場合には、 ｌ１＋ｌ２＞w1j＋ｌｊ−１＝(j＝1,2) …(15) となる為(7)式右辺の畳み込み行列 は横長行列となる。即ち、(13)式は不定方程式となる
が、例えば の様に解く事が出来る。
が、例えば の様に解く事が出来る。
以上の説明では、１入力２出力線形ＦＩＲ系１２１の入
力信号x(k)を正確に推定する逆フィルタ処理を例とした
が、(13)式左辺の目標ベクトル をＰ１個の離散信号で与えられる任意の伝達特性（イン
パルス応答）r1(k)とする事も可能である。
力信号x(k)を正確に推定する逆フィルタ処理を例とした
が、(13)式左辺の目標ベクトル をＰ１個の離散信号で与えられる任意の伝達特性（イン
パルス応答）r1(k)とする事も可能である。
この場合には、 w1j＋ｌ１−１Ｐ１（ｊ−１，２） …(17) となる様にフィルタ係数h1(k),h2(k)の個数、即ちＦＩ
Ｒフィルタ１１１，１１２のタップ数ｌ１，ｌ２を定め
る必要がある。そして、 (a)″ｌ１＝w12−１ ｌ２＝w11−１ となる場合には、 但し、 Ｑ１＝w1j＋ｌｊ−１−Ｐ１０ （ｊ＝１，２） (b)″ｌ１＞w12−１ ｌ２＞w11−１ となる場合には、 とフィルタ係数を定める事が出来る。
Ｒフィルタ１１１，１１２のタップ数ｌ１，ｌ２を定め
る必要がある。そして、 (a)″ｌ１＝w12−１ ｌ２＝w11−１ となる場合には、 但し、 Ｑ１＝w1j＋ｌｊ−１−Ｐ１０ （ｊ＝１，２） (b)″ｌ１＞w12−１ ｌ２＞w11−１ となる場合には、 とフィルタ係数を定める事が出来る。
尚インパルス応答g11(k),g12(k)に次のような遅延D11,D
12（D11≧D12とする） がある場合にはg11(k),g12(k)に共通する遅延D12を次の
ように削除することに依り、上記(a)′，(b)′の手続き
を用いることが出来る。
12（D11≧D12とする） がある場合にはg11(k),g12(k)に共通する遅延D12を次の
ように削除することに依り、上記(a)′，(b)′の手続き
を用いることが出来る。
以上のようにこの発明に依れば、フィルタ係数h1(k),h2
(k)を正しく求めることができ従来技術では実現不可能
だった逆フィルタ処理を正確に実現することが出来る。
(k)を正しく求めることができ従来技術では実現不可能
だった逆フィルタ処理を正確に実現することが出来る。
さらに上記基本原理は、多入力多出力線形ＦＩＲ系に１
つ或いはそれ以上の信号伝達経路を並列接続することに
依って得られる第２図に示すｍ入力ｎ出力（ｍ＜ｎ）線
形ＦＩＲ系１２２を対象とする場合へと容易に拡張出来
る。
つ或いはそれ以上の信号伝達経路を並列接続することに
依って得られる第２図に示すｍ入力ｎ出力（ｍ＜ｎ）線
形ＦＩＲ系１２２を対象とする場合へと容易に拡張出来
る。
ｍ入力ｎ出力線形ＦＩＲ系１２２のインパルス応答をベ
クトル インパルス応答gij(k)はwij個の離散信号で表わされて
いるものとする）、Pi個の離散信号で与えられる所望の
伝達特性（インパルス応答）ri(k)をベクトル Ｑ１＝wij＋ｌｊ−１−P1０ （但し、ｉ＝1,2,…，ｍ；ｊ＝1,2,…，ｎ）） とすれば、 （ｉ） 但しｊ＝１，２，…，ｎ を満足するフィルタタップ数ｌｊを用いる場合には 但し、 (ii) を満足するフィルタタップ数ｌｊを用いる場合には例え
ば、 Ｔは換置を表わす によりフィルタ係数を求めることにより逆フィルタ処理
を正確に実現するＦＩＲフィルタ111,112,…,11nが求め
られる。
クトル インパルス応答gij(k)はwij個の離散信号で表わされて
いるものとする）、Pi個の離散信号で与えられる所望の
伝達特性（インパルス応答）ri(k)をベクトル Ｑ１＝wij＋ｌｊ−１−P1０ （但し、ｉ＝1,2,…，ｍ；ｊ＝1,2,…，ｎ）） とすれば、 （ｉ） 但しｊ＝１，２，…，ｎ を満足するフィルタタップ数ｌｊを用いる場合には 但し、 (ii) を満足するフィルタタップ数ｌｊを用いる場合には例え
ば、 Ｔは換置を表わす によりフィルタ係数を求めることにより逆フィルタ処理
を正確に実現するＦＩＲフィルタ111,112,…,11nが求め
られる。
この発明は、以上述べたように従来技術では実現不可能
であった、多入力多出力線形ＦＩＲ系を対象とする逆フ
ィルタ処理を正確に実現するものであり、以下の特徴を
有する。
であった、多入力多出力線形ＦＩＲ系を対象とする逆フ
ィルタ処理を正確に実現するものであり、以下の特徴を
有する。
(a)多入力多出力線形ＦＩＲ系に１つ或いは複数の信号
伝達経路（有限のインパルス応答を有する適切な線形Ｆ
ＩＲ系）を付け加えることに依り、全体の系をｍ入力ｎ
出力（ｍ＜ｎ）線形ＦＩＲ系に作り換える。
伝達経路（有限のインパルス応答を有する適切な線形Ｆ
ＩＲ系）を付け加えることに依り、全体の系をｍ入力ｎ
出力（ｍ＜ｎ）線形ＦＩＲ系に作り換える。
(b)(20)式で与えられるタップ数ｌｊ（ｊ＝１，２，…
ｎ）を持つｎ個のＦＩＲフィルタの係数を で与えられる連立方程式の解として求める。
ｎ）を持つｎ個のＦＩＲフィルタの係数を で与えられる連立方程式の解として求める。
さらにこの発明は、 (c)音響系へ適用することに依り、雑音抑圧、選択受
聴、残響除去を実現出来る。
聴、残響除去を実現出来る。
(d)電波系へ適用することに依り、妨害波抑圧、マルチ
パス或いはゴーストの防止を実現出来る、ものである。
従来においては音響系や電波系においては、受信素子
（マイクロホン素子、アンテナ素子）数に従って形成さ
れる指向特性を制御することに依って、(c),(d)に挙げ
た問題に対処していたが、 (a)′反射の存在する一般室内音場や電場内で上記問題
を解消する為には、膨大な受信素子が必要とされる、 (b)′指向性制御の結果、受音或いは受信を希望する信
号まで歪んでしまう、 と云う欠点があった。しかしこの発明では、その原理か
ら明らかなように、前記(a)′，(b)′の欠点を伴わずに
目的信号を受音又は受信する装置として、音響系或いは
電波系の分野へ適用出来る。
パス或いはゴーストの防止を実現出来る、ものである。
従来においては音響系や電波系においては、受信素子
（マイクロホン素子、アンテナ素子）数に従って形成さ
れる指向特性を制御することに依って、(c),(d)に挙げ
た問題に対処していたが、 (a)′反射の存在する一般室内音場や電場内で上記問題
を解消する為には、膨大な受信素子が必要とされる、 (b)′指向性制御の結果、受音或いは受信を希望する信
号まで歪んでしまう、 と云う欠点があった。しかしこの発明では、その原理か
ら明らかなように、前記(a)′，(b)′の欠点を伴わずに
目的信号を受音又は受信する装置として、音響系或いは
電波系の分野へ適用出来る。
「実施例」 第３図はこの発明の一実施例を示し、全系をディジタル
信号処理系で構成したものである。第３図においてｍ入
力ｎ出力線形ＦＩＲ系１２２のｍ個の入力点に配された
ｍ個の信号源S1,S2,…,Smから放射される信号はｍ入力
ｎ出力線形ＦＩＲ系１２２を経てそのｎ個の出力点に配
された受信素子11,12,…,1nでそれぞれ受信され、受信
素子11,12,…,1nの出力u1(t),u2(t)…,un(t)はｎチャン
ネルＡ／Ｄ変換器２１に供給される。Ａ／Ｄ変換器２１
はその入力信号u1(t),u2(t)…,un(t)をそれぞれ離散化
した信号u1(k),u2(k)…,un(k)（ｋ＝１，２，…）とし
て出力し、これらは信号処理部１１を構成するｎ個のＦ
ＩＲフィルタ111,112,…,11nに入力され、これらＦＩＲ
フィルタの出力は加算器２２に依り足し合わされ、この
発明装置の出力y(k)となる。
信号処理系で構成したものである。第３図においてｍ入
力ｎ出力線形ＦＩＲ系１２２のｍ個の入力点に配された
ｍ個の信号源S1,S2,…,Smから放射される信号はｍ入力
ｎ出力線形ＦＩＲ系１２２を経てそのｎ個の出力点に配
された受信素子11,12,…,1nでそれぞれ受信され、受信
素子11,12,…,1nの出力u1(t),u2(t)…,un(t)はｎチャン
ネルＡ／Ｄ変換器２１に供給される。Ａ／Ｄ変換器２１
はその入力信号u1(t),u2(t)…,un(t)をそれぞれ離散化
した信号u1(k),u2(k)…,un(k)（ｋ＝１，２，…）とし
て出力し、これらは信号処理部１１を構成するｎ個のＦ
ＩＲフィルタ111,112,…,11nに入力され、これらＦＩＲ
フィルタの出力は加算器２２に依り足し合わされ、この
発明装置の出力y(k)となる。
波形記録器２３及び目標波形記録器５０の出力は演算設
定器２４に供給され、演算測定器２４の出力はフィルタ
係数決定器２５へ入力され、フィルタ係数決定器２５の
出力は信号処理部１１に出力され、これを構成するｎ個
のＦＩＲフィルタ111,112,…,11nの設定入力側に供給さ
れ、その各フィルタ係数を設定する。
定器２４に供給され、演算測定器２４の出力はフィルタ
係数決定器２５へ入力され、フィルタ係数決定器２５の
出力は信号処理部１１に出力され、これを構成するｎ個
のＦＩＲフィルタ111,112,…,11nの設定入力側に供給さ
れ、その各フィルタ係数を設定する。
次にこの実施例の動作を説明する。この実施例において
ｎ個のＦＩＲフィルタ111,112,…,11nは、入力信号u
1(k),u2(k)…,un(k)を未処理のまま直ちに出力するよう
な状態に初期設定されているものとする。また波形記録
器２３及び目標波形記録器５０にはｍｎ個のインパルス
応答ベクトル インパルス応答長）とｍ個の目標インパルス応答（目標
伝達特性）ベクトル Pi：目標インパルス応答長）があらかじめ記憶されてい
るものとする。
ｎ個のＦＩＲフィルタ111,112,…,11nは、入力信号u
1(k),u2(k)…,un(k)を未処理のまま直ちに出力するよう
な状態に初期設定されているものとする。また波形記録
器２３及び目標波形記録器５０にはｍｎ個のインパルス
応答ベクトル インパルス応答長）とｍ個の目標インパルス応答（目標
伝達特性）ベクトル Pi：目標インパルス応答長）があらかじめ記憶されてい
るものとする。
演算設定器２４は波形記録器２３から入力されるインパ
ルス応答ベクトル 及び目標波形記録器５０から入力される目標インパルス
応答ベクトル を用いて、(20b)式に示した様な畳み込み行列 と目標ベクトル とを作り、それらをフィルタ係数決定器２５へ出力す
る。該演算設定器２４では第４図に示すように以下の処
理が行われる。
ルス応答ベクトル 及び目標波形記録器５０から入力される目標インパルス
応答ベクトル を用いて、(20b)式に示した様な畳み込み行列 と目標ベクトル とを作り、それらをフィルタ係数決定器２５へ出力す
る。該演算設定器２４では第４図に示すように以下の処
理が行われる。
(1)波形記録器２３から入力されるインパルス応答ベク
トル の長さwijを用いた次の連立方程式 但し、ｊ＝１，２，…，ｎ を満足するフィルタタップ数ｌｊを求める。
トル の長さwijを用いた次の連立方程式 但し、ｊ＝１，２，…，ｎ を満足するフィルタタップ数ｌｊを求める。
(2)目標波形記録器５０から入力される目標インパルス
応答ベクトル の長さＰｉとwij＋ｌｊ−１（但しｉ＝１，２，…，
ｍ；ｊは１，２，…，ｎの中の任意の数、第５図ではｊ
＝１としている）の大きさを比較し、 (2-1)Ｐｉ＝wij＋ｌｊ−１ …(24a) ならば処理を行わない。
応答ベクトル の長さＰｉとwij＋ｌｊ−１（但しｉ＝１，２，…，
ｍ；ｊは１，２，…，ｎの中の任意の数、第５図ではｊ
＝１としている）の大きさを比較し、 (2-1)Ｐｉ＝wij＋ｌｊ−１ …(24a) ならば処理を行わない。
(2-2)Ｐｉ＞wij＋ｌｊ−１ …(24b) であれば Δｌ＝Ｐｉ−（wij＋ｌｊ−１） …(24c) で与えられるΔｌを(1)で求めたフィルタタップ数ｌｊ
に加える。即ち ｌｊ′＝Δｌ＋ｌｊ(j=1,2,…,n) …(24d) で与えられるタップ数ｌｊ′を改めてフィルタタップ数
ｌｊとする。
に加える。即ち ｌｊ′＝Δｌ＋ｌｊ(j=1,2,…,n) …(24d) で与えられるタップ数ｌｊ′を改めてフィルタタップ数
ｌｊとする。
(2-3)Ｐｉ＜wij＋ｌｊ−１ …(24e) であれば、各目標インパルス応答ベクトル を次の様に作り変える。
Ｑ１＝wij＋ｌｊ−１−Ｐｉ なる３つの処理をｉ＝１からｉ＝ｍまで繰り返す (3)各インパルス応答ベクトル を用いて、 の畳み込み行列 （(20b)式参照） を作る。
(4)各目標インパルス応答ベクトル を用いて （ｊは１，２，…，ｎのうちの任意の数、（図４ではｊ
＝１とした）の目標ベクトル （(20b)式参照） を作る。
＝１とした）の目標ベクトル （(20b)式参照） を作る。
フィルタ係数決定器２５は演算器２６と係数分配器２７
で構成されており、演算設定器２４から入力される を用いた連立方程式 を解くことに依って得られるフィルタ係数を、ＦＩＲフ
ィルタ111,112,…,11nに分配するものである。演算器２
６では第５図に示すような処理が行われる。
で構成されており、演算設定器２４から入力される を用いた連立方程式 を解くことに依って得られるフィルタ係数を、ＦＩＲフ
ィルタ111,112,…,11nに分配するものである。演算器２
６では第５図に示すような処理が行われる。
(5)畳み込み行列 の行数 （但しｉは１，２，…，ｎの中の任意の数）と列数 を比較し、 (5-1) ならば 但し、 （ｉ＝１，２，…，ｎ） (5-2) ならば、 を計算にフィルタ係数ベクトル を求める。
又(5)と同等の処理は (6) α(q)：収束係数 ｑ：アルゴリズム(32)式の繰り返し回数 のような逐次型アルゴリズム（(32)式は最急降下法に基
づくもの）を行う逐次型演算器２８（第６図）を用いて
も可能であり、この逐次型演算は第７図に示すように を初期設定し、誤差 を求め、この誤差の二乗が設定値E0になるまでαを求
め、更に次の を求めることを繰返す。このような逐次型演算を行え
ば、(28),(30)式のような条件判断は不要となる。
づくもの）を行う逐次型演算器２８（第６図）を用いて
も可能であり、この逐次型演算は第７図に示すように を初期設定し、誤差 を求め、この誤差の二乗が設定値E0になるまでαを求
め、更に次の を求めることを繰返す。このような逐次型演算を行え
ば、(28),(30)式のような条件判断は不要となる。
以上に述べた(1)〜(5)又は(1)〜(4),(6)の手続きで得ら
れるフィルタ係数 に依って設定されたＦＩＲフィルタ111,112,…,11nは、
線形ＦＩＲ系１２２に対する逆フィルタ処理を実現する
ものとなっている。
れるフィルタ係数 に依って設定されたＦＩＲフィルタ111,112,…,11nは、
線形ＦＩＲ系１２２に対する逆フィルタ処理を実現する
ものとなっている。
第８図に示すように、波形記録器２３は磁気ディスクや
半導体メモリ等で構成される記憶装置２９とマイクロプ
ロセッサ等で構成出来る処理装置３０の組み合わせに依
って実現することが可能であり、目標波形記録器も同様
である。また、演算設定器２４やフィルタ係数決定器２
５はマイクロプロセッサ等に依って、ＦＩＲフィルタ11
1,112,…,11mは第９図に示すように各フィルタ係数h
1(1)h1(2)…h1(l1)が与えられる乗算器311,312,…,3
1L1、u1(k)を順次その信号間隔だけ遅延し、出力を乗算
器312…31L1へ供給する遅延素子321,322,…,32L1-1、乗
算器311…31L1の出力を加算する和算器３２を組み合わ
せたり、ＤＳＰを用いることに依って実現することが可
能である。受信素子11,12,…,1nとしては、対象となる
線形ＦＩＲ系の出力点の信号を取出す線形な素子であれ
ば任意のものが使用可能である。例としては、電気的な
接続を行うコネクタの他、マイクロホン素子や受信用ア
ンテナ素子等を用いることも出来る。
半導体メモリ等で構成される記憶装置２９とマイクロプ
ロセッサ等で構成出来る処理装置３０の組み合わせに依
って実現することが可能であり、目標波形記録器も同様
である。また、演算設定器２４やフィルタ係数決定器２
５はマイクロプロセッサ等に依って、ＦＩＲフィルタ11
1,112,…,11mは第９図に示すように各フィルタ係数h
1(1)h1(2)…h1(l1)が与えられる乗算器311,312,…,3
1L1、u1(k)を順次その信号間隔だけ遅延し、出力を乗算
器312…31L1へ供給する遅延素子321,322,…,32L1-1、乗
算器311…31L1の出力を加算する和算器３２を組み合わ
せたり、ＤＳＰを用いることに依って実現することが可
能である。受信素子11,12,…,1nとしては、対象となる
線形ＦＩＲ系の出力点の信号を取出す線形な素子であれ
ば任意のものが使用可能である。例としては、電気的な
接続を行うコネクタの他、マイクロホン素子や受信用ア
ンテナ素子等を用いることも出来る。
なお１度フィルタ係数を設定すればその変更する必要が
ない場合は、波形記録器２３、目標波形記録器５０、演
算設定器２４、フィルタ係数設定器２５を付属させて設
けておく必要はない。
ない場合は、波形記録器２３、目標波形記録器５０、演
算設定器２４、フィルタ係数設定器２５を付属させて設
けておく必要はない。
音響系への応用 第１０図に示すように反射のある一般室内音場に目標信
号源Ｓと雑音源N1,N2,…,Nm-1が存在するものとする。
受信素子11,12,…,1nとしてマイクロホン素子を用い
る。他の部分は第３図において同一番号が付してあるも
のに対応している。
号源Ｓと雑音源N1,N2,…,Nm-1が存在するものとする。
受信素子11,12,…,1nとしてマイクロホン素子を用い
る。他の部分は第３図において同一番号が付してあるも
のに対応している。
この場合には、使用するマイクロホン素子数ｎをｍ＋１
とすることに依り、必要なタップ数ｌｊを 但し、ｊは１，２，…，ｎの中の任意の数 なる関係を満足する数とし、目標信号源Ｓから各マイク
ロホンまでのインパルス応答を とすれば、実施例の(2)(4)に示した手続きで得られる目
標ベクトル を 但し、ｊは１，２，…，ｎの中の任意の数 と定めることに依って、 (a)雑音源N1,N2,…,Nm-1からの雑音は受音しない（雑音
抑圧） (b)目標信号源Ｓからの音は無歪で受音する（残響除
去） ことが、ｍ＋１個のマイクロホン素子で構成される小規
模な装置で実現出来る（一般室内音場では極めて多くの
反射が生じる為、指向性を制御する従来法では膨大な数
のマイクロホン素子が必要でかつ目的信号を無歪で受音
することは不可能だった。
とすることに依り、必要なタップ数ｌｊを 但し、ｊは１，２，…，ｎの中の任意の数 なる関係を満足する数とし、目標信号源Ｓから各マイク
ロホンまでのインパルス応答を とすれば、実施例の(2)(4)に示した手続きで得られる目
標ベクトル を 但し、ｊは１，２，…，ｎの中の任意の数 と定めることに依って、 (a)雑音源N1,N2,…,Nm-1からの雑音は受音しない（雑音
抑圧） (b)目標信号源Ｓからの音は無歪で受音する（残響除
去） ことが、ｍ＋１個のマイクロホン素子で構成される小規
模な装置で実現出来る（一般室内音場では極めて多くの
反射が生じる為、指向性を制御する従来法では膨大な数
のマイクロホン素子が必要でかつ目的信号を無歪で受音
することは不可能だった。
電波系への応用 第１１図に示すように例えばビルディング４１などの反
射体が存在する一般的な電場内にS1,S2,S3の信号源（送
信用アンテナ）とN1,N2,…,Nm-3の妨害波源があるもの
とする。この場合、この発明装置の受信素子11,12,…,1
nとして受信用アンテナ素子を用いる。他の部分は第３
図と対応する部分に同一番号を付けてある。
射体が存在する一般的な電場内にS1,S2,S3の信号源（送
信用アンテナ）とN1,N2,…,Nm-3の妨害波源があるもの
とする。この場合、この発明装置の受信素子11,12,…,1
nとして受信用アンテナ素子を用いる。他の部分は第３
図と対応する部分に同一番号を付けてある。
使用するアンテナ素子数ｎ及び目標ベクトル を(33),(34)式に従って定めれば、 (a)妨害波源N1,N2,…,Nm-1からの信号は受信しない（妨
害波抑圧） (b)信号源S1,S2,S3からS1を選択する（選択受信） (c)S1からの信号を無歪で受信する（マルチパス或いは
ゴーストの除去） を、ｍ＋１個の受信用アンテナ素子で構成される小規模
な装置で実現出来る。この受信装置は、指向性を制御す
る従来法の持つ欠点（音響系における従来技術の欠点と
同じ）を解消するものであり、特にビルの谷間や家の立
て込んでいる場所で有効となる。
害波抑圧） (b)信号源S1,S2,S3からS1を選択する（選択受信） (c)S1からの信号を無歪で受信する（マルチパス或いは
ゴーストの除去） を、ｍ＋１個の受信用アンテナ素子で構成される小規模
な装置で実現出来る。この受信装置は、指向性を制御す
る従来法の持つ欠点（音響系における従来技術の欠点と
同じ）を解消するものであり、特にビルの谷間や家の立
て込んでいる場所で有効となる。
「発明の効果」 以上説明したようにこの発明の多入力形制御装置は、ｍ
入力ｎ出力（ｍ＜ｎ）と見做すことの出来る線形ＦＩＲ
系を対象として (a)この線形ＦＩＲ系の入力信号を推定する (b)この線形ＦＩＲ系の入力信号に所望の特性を付与し
て取り出す と云う逆フィルタ処理を正確に実現するものであり、特
に後述する応用において、 (c)目的信号を無歪で受信出来る（上記(a),(b)の逆フィ
ルタ処理に相当する） (d)受信素子数が信号源数＋１でも良い為、少素子数構
成の受信装置が実現できる という効果をもっている。
入力ｎ出力（ｍ＜ｎ）と見做すことの出来る線形ＦＩＲ
系を対象として (a)この線形ＦＩＲ系の入力信号を推定する (b)この線形ＦＩＲ系の入力信号に所望の特性を付与し
て取り出す と云う逆フィルタ処理を正確に実現するものであり、特
に後述する応用において、 (c)目的信号を無歪で受信出来る（上記(a),(b)の逆フィ
ルタ処理に相当する） (d)受信素子数が信号源数＋１でも良い為、少素子数構
成の受信装置が実現できる という効果をもっている。
第１図はこの発明装置の原理を説明する為の図、第２図
はこの発明装置を多入力多出力線形ＦＩＲ系へ適応させ
る為第１図を拡張したこの発明装置の原理図、第３図は
この発明装置の一実施例を示すブロック図、第４図は演
算設定器２４の処理内容の例を示す流れ図、第５図はフ
ィルタ係数決定器２５の処理内容の例を示す流れ図、第
６図は逐次型演算器２８を用いた場合のフィルタ係数決
定器２５の一実施例を示すブロック図、第７図は第６図
に示したフィルタ係数決定器２５の処理内容の例を表わ
す流れ図、第８図は波形記録器２３の一実施例を示すブ
ロック図、第９図はＦＩＲフィルタ11の一実施例を示す
ブロック図、第１０図はこの発明を音響系へ応用した例
を示すブロック図、第１１図はこの発明を電波系に応用
した例を示すブロック図、第１２図は従来の逆フィルタ
処理装置の原理を説明するための図である。 11,12,…,1n：受信素子（第１０図ではマイクロホン素
子、第１１図では受信用アンテナ素子を表わしてい
る）、110：従来技術で用いられた逆フィルタ処理用フ
ィルタ、111,112,…,11n：ＦＩＲフィルタ、１１：111
〜11nで構成される信号処理部、２１：ｎチャンネルＡ
／Ｄ変換器、２２：加算器、２３：波形記録器、２４：
演算設定器、２５：フィルタ係数決定器、２６：演算
器、２７：係数分配器、２８：逐次型演算器、２９：記
憶装置、３０：処理装置、311,312,…,31Lj：乗算器、3
21,322,…,32Lj-1：遅延素子、５０：目標波形記録器、
１０１：１入力１出力線形ＦＩＲ系、１１１：１０１に
並列接続する信号伝達経路（１入力１出力の適当な線形
ＦＩＲ系、１２１：１０１と１１１で構成される１入力
２出力線形ＦＩＲ系システム、１２２：ｍ入力ｎ出力
（ｍ＜ｎ）線形ＦＩＲ系、x(k)：１０１，１２１への離
散化された入力信号、y(k)：１１０の離散化された出力
信号、ｙ′(k)：11１と11２の離散化された出力信号の
和、S1,S2,…,Sm：１２２へのｍ個の入力信号源、Ｓ：
音源、N1,N2,…,Nm-1：雑音源、S1,S2,S3：信号源、N1,
N2,…,Nm-3：妨害波源。
はこの発明装置を多入力多出力線形ＦＩＲ系へ適応させ
る為第１図を拡張したこの発明装置の原理図、第３図は
この発明装置の一実施例を示すブロック図、第４図は演
算設定器２４の処理内容の例を示す流れ図、第５図はフ
ィルタ係数決定器２５の処理内容の例を示す流れ図、第
６図は逐次型演算器２８を用いた場合のフィルタ係数決
定器２５の一実施例を示すブロック図、第７図は第６図
に示したフィルタ係数決定器２５の処理内容の例を表わ
す流れ図、第８図は波形記録器２３の一実施例を示すブ
ロック図、第９図はＦＩＲフィルタ11の一実施例を示す
ブロック図、第１０図はこの発明を音響系へ応用した例
を示すブロック図、第１１図はこの発明を電波系に応用
した例を示すブロック図、第１２図は従来の逆フィルタ
処理装置の原理を説明するための図である。 11,12,…,1n：受信素子（第１０図ではマイクロホン素
子、第１１図では受信用アンテナ素子を表わしてい
る）、110：従来技術で用いられた逆フィルタ処理用フ
ィルタ、111,112,…,11n：ＦＩＲフィルタ、１１：111
〜11nで構成される信号処理部、２１：ｎチャンネルＡ
／Ｄ変換器、２２：加算器、２３：波形記録器、２４：
演算設定器、２５：フィルタ係数決定器、２６：演算
器、２７：係数分配器、２８：逐次型演算器、２９：記
憶装置、３０：処理装置、311,312,…,31Lj：乗算器、3
21,322,…,32Lj-1：遅延素子、５０：目標波形記録器、
１０１：１入力１出力線形ＦＩＲ系、１１１：１０１に
並列接続する信号伝達経路（１入力１出力の適当な線形
ＦＩＲ系、１２１：１０１と１１１で構成される１入力
２出力線形ＦＩＲ系システム、１２２：ｍ入力ｎ出力
（ｍ＜ｎ）線形ＦＩＲ系、x(k)：１０１，１２１への離
散化された入力信号、y(k)：１１０の離散化された出力
信号、ｙ′(k)：11１と11２の離散化された出力信号の
和、S1,S2,…,Sm：１２２へのｍ個の入力信号源、Ｓ：
音源、N1,N2,…,Nm-1：雑音源、S1,S2,S3：信号源、N1,
N2,…,Nm-3：妨害波源。
Claims (1)
- 【請求項１】ｍ個の入力点とｎ個の出力点（ｍ、ｎは１
以上の整数、かつｎ＞ｍ）との間に存在する複数の信号
伝達経路それぞれが有するインパルス応答が有限である
と見做せる線形系を対象とする装置であって、 その線形系から出力される信号を受信し、上記ｎ個の出
力点に配されたｎ個の受信素子と、 そのｎ個の受信素子にそれぞれ接続されたｎ個のＦＩＲ
フィルタと、 これらｎ個のＦＩＲフィルタの出力側と接続され、その
ｎ個の出力を加算する加算手段とを備え、 上記ｍ個の入力点の第ｉ番目（ｉ＝１，２，…，ｍ）と
上記ｎ個の出力点の第ｊ番目（ｊ＝１，２，…，ｎ）と
の間の上記線形系のインパルス応答gij（ｋ）（ｋ＝
１，２，…）をWij個の離散信号で表わし、上記第ｉ番
目の入力点と上記加算手段の出力との間の所望の伝達特
性ri（ｋ）をＰｉ個の離散信号で表わすと、上記第ｊ番
目の出力点に接続される第ｊ番目のＦＩＲフィルタは、 但し、ｉ＝１，２，…，ｍ ｊ＝１，２，…，ｎ を満足するタップ数１ｊを有し、かつその第ｊ番目のＦ
ＩＲフィルタの係数ｈｊ（ｋ）が 但し、ｉ＝１，２，…，ｍ は離散畳み込み演算 を満足している多入力形制御装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3265986A JPH0654883B2 (ja) | 1986-02-17 | 1986-02-17 | 多入力形制御装置 |
| US06/839,677 US4683590A (en) | 1985-03-18 | 1986-03-14 | Inverse control system |
| CA000504255A CA1242003A (en) | 1985-03-18 | 1986-03-17 | Inverse control system |
| DE8686103668T DE3686497T2 (de) | 1985-03-18 | 1986-03-18 | Inverse kontrolleinrichtung. |
| EP86103668A EP0195416B1 (en) | 1985-03-18 | 1986-03-18 | Inverse control system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3265986A JPH0654883B2 (ja) | 1986-02-17 | 1986-02-17 | 多入力形制御装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS62190935A JPS62190935A (ja) | 1987-08-21 |
| JPH0654883B2 true JPH0654883B2 (ja) | 1994-07-20 |
Family
ID=12364991
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3265986A Expired - Fee Related JPH0654883B2 (ja) | 1985-03-18 | 1986-02-17 | 多入力形制御装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0654883B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8271277B2 (en) | 2006-03-03 | 2012-09-18 | Nippon Telegraph And Telephone Corporation | Dereverberation apparatus, dereverberation method, dereverberation program, and recording medium |
| JP4787727B2 (ja) * | 2006-12-04 | 2011-10-05 | 日本電信電話株式会社 | 音声収音装置、その方法、そのプログラム、およびその記録媒体 |
-
1986
- 1986-02-17 JP JP3265986A patent/JPH0654883B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPS62190935A (ja) | 1987-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0195416A2 (en) | Inverse control system | |
| JP3216704B2 (ja) | 適応アレイ装置 | |
| US5383164A (en) | Adaptive system for broadband multisignal discrimination in a channel with reverberation | |
| CN108141691B (zh) | 自适应混响消除系统 | |
| US6895093B1 (en) | Acoustic echo-cancellation system | |
| JP3701940B2 (ja) | 目的信号源から雑音環境に放射される信号を処理するシステム及び方法 | |
| Kirkeby et al. | Local sound field reproduction using digital signal processing | |
| JP3397269B2 (ja) | 多チャネル反響消去方法 | |
| KR100723283B1 (ko) | 음향 에코 및 잡음 제거 적응성 필터 | |
| EP2754307B1 (en) | Apparatus and method for listening room equalization using a scalable filtering structure in the wave domain | |
| US6381272B1 (en) | Multi-channel adaptive filtering | |
| JPH04150400A (ja) | 音像定位装置 | |
| JP2008141593A (ja) | 音声収音装置、その方法、そのプログラム、およびその記録媒体 | |
| JPH0654883B2 (ja) | 多入力形制御装置 | |
| JP2558445B2 (ja) | 多チャンネル制御装置 | |
| JP4041770B2 (ja) | 音響エコー消去方法、その装置、プログラム及びその記録媒体 | |
| JP3579508B2 (ja) | オーディオ装置 | |
| WO1994024662A1 (en) | Method of calculating filter weights for compression wave cancellation systems | |
| JP2002261659A (ja) | 多チャネルエコーキャンセル方法、その装置、そのプログラム及び記録媒体 | |
| JP6892111B2 (ja) | 多入力多出力系の直接逆同定方法及び装置及びプログラム及び記憶媒体、多入力多出力逆フィルタ装置及び方法 | |
| JP2548103B2 (ja) | 音響再生装置 | |
| JP3001092B2 (ja) | 適応アレイの制御方法および適応アレイ装置 | |
| Muron et al. | Modelling of reverberations and audioconference rooms | |
| Uncini et al. | Introduction to Adaptive Signal and Array Processing | |
| JPS61195099A (ja) | 音響再生装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |