(57)【要約】【課題】所望の音源から発せられる音のみを入力する。【解決手段】死点形成手段２から出力される音圧Ｏ(f,k)には死点に存在する音源から発せられる音以外の音、すなわち、雑音の音圧のみが含まれ、一方、集音センサ手段１で検出される音圧Ｍ(f,k)には死点に存在する音源から発せられる音の音圧と雑音の音圧の双方が含まれている。そして、目的話者音声抽出手段３によって死点に存在する音源から発せられる音の音圧Ｓ(f,k)のみを抽出するので、雑音が死点の方向から到来する場合においても死点に存在する音源から発せられる音のみを入力することができる。【選択図】図１
2008304555.tif 000002
【特許請求の範囲】
【請求項１】
音圧、当該音圧の時間微分値、当該音圧を二次元直交座標系の各軸方向に微分した空間微分値をそれぞれ検出する集音センサ手段と、集音センサ手段で検出される音圧、時間微分値、空間微分値に対して所定の係数ベクトルとの荷重和及び低域通過フィルタ処理を行うことにより集音感度が最小となる死点を予め設定した目的話者の位置に形成する死点形成手段と、集音センサ手段で検出される音圧と死点形成手段から出力される音圧を用いて前記目的話者から発せられる音声の音圧のみを抽出する目的話者音声抽出手段とを備えたことを特徴とする音響入力装置。
【請求項２】
目的話者音声抽出手段は、スペクトル・サブトラクション法によって前記音源から発せられる音の音圧を抽出することを特徴とする請求項１記載の音響入力装置。
【請求項３】
目的話者音声抽出手段は、独立成分分析によって前記音源から発せられる音の音圧を抽出することを特徴とする請求項１記載の音響入力装置。
【請求項４】
目的話者音声抽出手段は、独立成分分析を行う前に主成分分析を行うことを特徴とする請求項３記載の音響入力装置。
【請求項５】
集音センサ手段は、前記二次元直交座標系の各軸と直交する向きに複数が配置されたマイクロホンを具備することを特徴とする請求項１〜４の何れか１項に記載の音響入力装置。
【請求項６】
集音センサ手段は、振動板が中央の１点で支持された２軸直交型のジンバル構造を有するマイクロホンを具備することを特徴とする請求項１〜４の何れか１項に記載の音響入力装置。
【請求項７】
前記音源位置を変更する音源位置変更手段を備え、音源位置変更手段は、集音センサ手段で検出される音圧の瞬時パワーを死点形成手段から出力される音圧の瞬時パワーで除した値が所定のしきい値以上となるときに集音センサ手段で検出される音圧、時間微分値、空間微分値に基づいて推定される位置に前記音源位置を変更することを特徴とする請求項１〜６の何れか１項に記載の音響入力装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音源から発せられる音を入力する音響入力装置に関するものである。
【背景技術】
【０００２】
周囲騒音や残響の存在する環境下で特定の音源から発せられる音、例えば、人の発する音声（話者音声）のみを入力する音響入力装置が従来より種々提供されている。例えば、特許文献１に記載されている従来例は、２つのマイクロホンと、各マイクロホンで集音した音響信号毎に周波数分析して２つのチャンネル別の周波数成分を得るとともに、各チャネルの周波数成分について適応ビームフォーマ処理を行うことにより、目的音（話者音声）の方向以外の感度を低くして周囲雑音等が抑圧された音響信号を取得し、同様に周囲雑音等の方向以外の感度を低くして目的音が抑圧された音響信号を取得し、適応ビームフォーマ処理で用いられるフィルタ係数から目的音方向と周囲雑音方向を推定して逐次修正し、さらに、スペクトル・サブトラクション処理によって前者の音響信号と後者の音響信号から周囲雑音成分を除去することによって目的音（話者音声）のみを入力するものである。
【特許文献１】特開２０００−４７６９９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、上記従来例においては、例えばマイクロホンからみて話者と雑音源が同じ方向に存在する場合には目的音と目的音以外の雑音とを分離することができないという問題がある。
【０００４】
本発明は上記事情に鑑みて為されたものであり、その目的は、所望の音源から発せられる音のみを入力することができる音響入力装置を提供することにある。
【課題を解決するための手段】
【０００５】
請求項１の発明は、上記目的を達成するために、音圧、当該音圧の時間微分値、当該音圧を二次元直交座標系の各軸方向に微分した空間微分値をそれぞれ検出する集音センサ手段と、集音センサ手段で検出される音圧、時間微分値、空間微分値に対して所定の係数ベクトルとの荷重和及び低域通過フィルタ処理を行うことにより集音感度が最小となる死点を予め設定した目的話者の位置に形成する死点形成手段と、集音センサ手段で検出される音圧と死点形成手段から出力される音圧を用いて前記目的話者から発せられる音声の音圧のみを抽出する目的話者音声抽出手段とを備えたことを特徴とする。
【０００６】
請求項１の発明によれば、死点形成手段から出力される音圧には死点に存在する音源から発せられる音以外の音、すなわち、雑音の音圧のみが含まれ、一方、集音センサ手段で検出される音圧には死点に存在する音源から発せられる音の音圧と雑音の音圧の双方が含まれており、目的話者音声抽出手段によって死点に存在する音源から発せられる音の音圧のみを抽出するので、雑音が死点の方向から到来する場合においても死点に存在する音源から発せられる音のみを入力することができる。
【０００７】
請求項２の発明は、請求項１の発明において、目的話者音声抽出手段は、スペクトル・サブトラクション法によって前記音源から発せられる音の音圧を抽出することを特徴とする。
【０００８】
請求項３の発明は、請求項１の発明において、目的話者音声抽出手段は、独立成分分析によって前記音源から発せられる音の音圧を抽出することを特徴とする。
【０００９】
請求項４の発明は、請求項３の発明において、目的話者音声抽出手段は、独立成分分析を行う前に主成分分析を行うことを特徴とする。
【００１０】
請求項５の発明は、請求項１〜４の何れか１項の発明において、集音センサ手段は、前記二次元直交座標系の各軸と直交する向きに複数が配置されたマイクロホンを具備することを特徴とする。
【００１１】
請求項６の発明は、請求項１〜４の何れか１項の発明において、集音センサ手段は、振動板が中央の１点で支持された２軸直交型のジンバル構造を有するマイクロホンを具備することを特徴とする。
【００１２】
請求項７の発明は、請求項１〜６の何れか１項の発明において、前記音源位置を変更する音源位置変更手段を備え、音源位置変更手段は、集音センサ手段で検出される音圧の瞬時パワーを死点形成手段から出力される音圧の瞬時パワーで除した値が所定のしきい値以上となるときに集音センサ手段で検出される音圧、時間微分値、空間微分値に基づいて推定される位置に前記音源位置を変更することを特徴とする。
【００１３】
請求項７の発明によれば、音源の位置が変動した場合でも音源位置変更手段によって音源位置が変更されるために音源位置と死点がずれることがなく、その結果、音源が移動する場合においても当該音源から発せられる音のみを入力することができる。
【発明の効果】
【００１４】
本発明によれば、雑音が死点の方向から到来する場合においても死点に存在する音源から発せられる音のみを入力することができる。
【発明を実施するための最良の形態】
【００１５】
以下、図面を参照して本発明の実施形態を詳細に説明する。尚、実施形態の説明において参照する文献については、参考文献１，参考文献２，…のように表記し、それぞれの参考文献１，２，…の一覧を末尾に掲載する。
【００１６】
（実施形態１）
本実施形態の音響入力装置は、図１に示すように音圧、当該音圧の時間微分値、当該音圧を二次元直交座標系の各軸方向に微分した空間微分値をそれぞれ検出する集音センサ手段１と、集音センサ手段１で検出される音圧、時間微分値、空間微分値に対して所定の係数ベクトルとの荷重和及び低域通過フィルタ処理を行うことにより集音感度が最小となる死点を予め設定した目的話者の位置に形成する死点形成手段２と、集音センサ手段１で検出される音圧と死点形成手段２から出力される音圧を用いて前記死点に存在する音源（話者）から発せられる音の音圧のみを抽出する目的話者音声抽出手段３とを備えている。
【００１７】
集音センサ手段１は、図１に示すように三次元直交座標系のｘ軸並びにｙ軸と直交する向き（ｚ軸の正の向き）に複数（図示例では４本）が配置された無指向性のマイクロホン１０Ａ，１０Ｂ，１０Ｃ，１０Ｄと、マイクロホン１０Ａ〜１０Ｄの出力信号ｆA(t)，ｆB(t)，ｆC(t)，ｆD(t)に対して時空間勾配測定処理を行う時空間勾配測定処理部１１とを具備する。時空間勾配測定処理部１１では、各マイクロホン１０Ａ〜１０Ｄの出力信号ｆA(t)，ｆB(t)，ｆC(t)，ｆD(t)から音圧の同相成分Ｍ(t)、時間微分値（時間勾配成分）Ｍt(t)、ｘ軸方向空間微分値（ｘ軸方向空間勾配成分）Ｍx(t)、ｙ軸方向空間微分値（ｙ軸方向空間勾配成分）Ｍy(t)をそれぞれ下式より求める。
【００１８】
Ｍ(t)=ｆA(t)+ｆB(t)+ｆC(t)+ｆD(t)
Ｍt(t)=dｆA(t)/dt+dｆB(t)/dt+dｆC(t)/dt+dｆD(t)/dt
Ｍx(t)=ｆA(t)+ｆB(t)−ｆC(t)−ｆD(t)
Ｍy(t)=ｆA(t)−ｆB(t)+ｆC(t)−ｆD(t)
死点形成手段２は、集音センサ手段１から出力される同相成分Ｍ(t)、時間微分値Ｍt(t)、ｘ軸方向空間微分値Ｍx(t)、ｙ軸方向空間微分値Ｍy(t)を用いて、時空間勾配法を応用することで死点を形成するものである。ここで、死点形成手段２による死点形成処理を説明するに当たって、初めに時空間勾配法について詳しく説明する。
【００１９】
時空間勾配法とは、そもそも動画像中の見かけの速度場であるオプティカルフローを決定する手法の一つとして提案されたものである（参考文献１参照）。動画像中の濃淡パターンの特徴を表す画像関数ｆ(x,y,t)が、運動に際し不変に保たれるとの仮定（ｆ(x,y,t)=ｆ(x+δx,y+δy,t+δt)）より、ある点(x,y)におけるオプティカルフローの速度と、動画像の濃淡分布の空間勾配および時間勾配を関係付ける式をもとにした解析手法である。以下、この手法について詳しく解説する。
【００２０】
時刻t+δtにおいて、座標(x+δx,y+δy)での濃淡パターンｆ(x+δx,y+δy,t+δt)を(x,y,t)のまわりでテーラー展開すると、
ｆ(x+δx,y+δy,t+δt)=ｆ(x,y,t)+ｆxδx+ｆyδy+ｆtδt+Ｏ(δx+δy+δt）    …（１）
となる。ここで、Ｏ(δx+δy+δt）はδx,δy,δtの２次以上の項であるが、微小量であるために以降では無視する。この時、時刻ｔにおいて座標(x,y)にある濃淡パターンが、δt時刻経過した後に座標(x+δx,y+δy)にその濃度値分布を一定に保ったまま移動した時、その対応付けから次式が成り立つ。
【００２１】
ｆ(x,y,t)=ｆ(x+δx,y+δy,t+δt)
=ｆ(x,y,t)+ｆxδx+ｆyδy+ｆtδt    …（２）
ｆxδx+ｆyδy+ｆtδt=０                 …（３）
式（３）の両辺をδtで割ると、
ｆxδx/δt+ｆyδy/δt+ｆt=０              …（４）
を得る。ここで、δtが無限小であると仮定して、δt→０とすると次式を得る。
【００２２】
ｆxdx/dt+ｆydy/dt+ｆt=０                  …（５）
オプティカルフロー速度v=(u，v)=(dx/dt，dy/dt)を用いると、式（５）は、
uｆx+vｆy+ｆt=０                        …（６）
となり、式（６）は動画像の濃淡値の時間、空間に関する勾配とオプティカルフロー速度vとを関係付ける式である。
【００２３】
次に、「ある着目点の近傍領域Γにおいて速度場はほぼ一定であると近似できる」という仮定を行う。この時、領域Γ内のいたるところで式（６）が成立しなければならない。そこで、式（６）の左辺の２乗積分（下記の式（７））を用いて評価し、最小自乗法によって速度場を求める。
【００２４】
【数１】
2008304555.tif 000003
【００２５】
式（７）をu，vに関して微分し、０とおくと、
uSxx+vSxy+Sxt=0，uSxy+vSyy+Syt=0          …（８）
【００２６】
【数２】
2008304555.tif 000004
【００２７】
が得られる。式（８）を解くと速度ベクトル(u，v)は
u=(SytSxy-SxtSyy)/(SxxSyy-S2xy)，v=(SxtSxy-SytSxx)/(SxxSyy-S2xy)  …（１０）
のように求められる。
【００２８】
次に、上述の動画像中のオプティカルフロー速度を求める時空間勾配法を応用して、音源が空間中に作る音場のある１点における音圧とその時空間勾配の間に成り立つ線形関係に基づいて、音源位置を定位する手法について説明する（参考文献２参照）。
【００２９】
図２に示すように観測点を原点とする三次元直交座標系を取り、その前方（z>0）に互いに無相関な点音源が複数個あるとする。音速をc、i番目の音源の座標を(xi,yi,zi)、音源と観測点との距離をRi=(x2i+y2i+z2i)1/2、音源音をgi(t)、各音源が観測点に形成する音場をfi(t)とすると、観測点に形成される合成音場ｆはこれらからの球面波の和として、
【００３０】
【数３】
2008304555.tif 000005
【００３１】
と表される。これを偏微分することにより、観測点での音場のx,y微分、時間微分は下記の式(12),(13),(14)で表される。
【００３２】
【数４】
2008304555.tif 000006
【００３３】
ここで、
ξix=xi/R2i,ξiy=yi/R2i    …(15)
は強度勾配と呼ばれ、
τix=xi/cRi,τiy=yi/cRi    …(16)
はx,y方向時間勾配と呼ばれる。
【００３４】
次に簡単のため、１音源の場合の音源定位手法について述べる。１音源の場合、式(１２)，（１３）は
fx=-ξxf-τxft，fy=-ξyf-τyft        …(17)
となり、式（１）と同様に最小自乗法を適用してτx,τy,ξx,ξyを求める。短時間の時間窓Γにおいて評価関数を
J=∫Γ｛(fx+ξxf+τxft)2+(fy+ξyf+τyft)2｝dt        …(18)
とする。式(18)をτx,τy,ξx,ξyに関して偏微分し、０とおくと下式が得られる。
【００３５】
∂J/∂τx=∫Γ2(fx+ξxf+τxft)・ftdt=0，∂J/∂τy=∫Γ2(fy+ξyf+τyft)・ftdt=0        …(19)
∂J/∂ξx=∫Γ2(fx+ξxf+τxft)・ftdt=0，∂J/∂ξy=∫Γ2(fy+ξyf+τyft)・ftdt=0        …(20)
ここで、観測窓Γから推定される共分散行列を
【００３６】
【数５】
2008304555.tif 000007
【００３７】
とおくと、式(19)，(20)は
Sxt+ξxSt+τxStt=0,Syt+ξySt+τyStt=0        …(22)
Sx+ξxS+τxSt=0,Sy+ξyS+τySt=0          …(23)
と書き直される。式(22)，(23)を解くことにより、τx,τy,ξx,ξyが次式のように求められる。
【００３８】
τx=(SxSt-SSxt)/(SStt-S2t)，τy=(SySt-SSyt)/(SStt-S2t)    …（24）
ξx=(SxtSt-SxStt)/(SStt-S2t)，ξy=(SytSt-SyStt)/(SStt-S2t)  …（25）
音源の方位角(x/R,y/R)=(cτx,cτy)については式(21),(24)から求められる。音源までの距離Ｒについては、式(15)，(16)から最小自乗法を適用することにより求められる。評価関数を
【００３９】
【数６】
2008304555.tif 000008
【００４０】
とし、これを1/Rで偏微分して0とおくと
【００４１】
【数７】
2008304555.tif 000009
【００４２】
となる。これを解くと
R=c(τ2x+τ2y)/(τxξx+τyξy)          …(28)
のように音源までの距離が求められる。
【００４３】
次に、音場の時空間勾配を利用して、指向性制御を行う手法について解説する(参考文献３〜５参照)。今、１音源の場合を仮定すると、観測点における音圧信号ｆ(t)のx,y方向の空間勾配は式(12)，(13)より
【００４４】
【数８】
2008304555.tif 000010
【００４５】
となる。この式を音源から観測点に向かうベクトルｒ=(x,y,z)を用いて書き直すと
【００４６】
【数９】
2008304555.tif 000011
【００４７】
となる。次にｆ(t),ｆt(t),∇ｆ(t)が観測される時、これらの荷重和は
【００４８】
【数１０】
2008304555.tif 000012
【００４９】
と表される。ただし、u,utは実数定数、w=(wx,wy,0)は観測点を原点とし、任意の方向を向いている単位ベクトルである。式(30)を式(31)に代入すると、
【００５０】
【数１１】
2008304555.tif 000013
【００５１】
となる。よって時空間勾配の荷重和は、ｆ(t),ｆt(t)に対してそれぞれ異なる指向特性Ｈ(ｒ)，Ｈt(ｒ)をもつフィルタの和として表される。Ｈ(ｒ)=αのとき、式(33)は
【００５２】
【数１２】
2008304555.tif 000014
【００５３】
と変形できる。ここで、２つのベクトルａ，ｂの成す角をθとすると以下の公式が成り立つ。
【００５４】
【数１３】
2008304555.tif 000015
【００５５】
式(38)の公式を用いると式(36)は次式のように書き換えられる。
【００５６】
【数１４】
2008304555.tif 000016
【００５７】
ここで、｜w｜=１より、
【００５８】
【数１５】
2008304555.tif 000017
【００５９】
という球の方程式で表される。u+α=0の場合には、式(35)は
ｒ・w=０      …(42)
となる。また、Ｈt(r)=αの時には式(34)は
【００６０】
【数１６】
2008304555.tif 000018
【００６１】
となるので、ベクトルｒとｗの成す角をθ(r)とすると｜w｜=１より
【００６２】
【数１７】
2008304555.tif 000019
【００６３】
となる。よって、式(43)は
【００６４】
【数１８】
2008304555.tif 000020
【００６５】
となる。
【００６６】
式(41)、(42)、(45)より、Ｈ(r)，Ｈt(r)について次のような性質を得る。
１）２つの指向特性Ｈ(r)，Ｈt(r)はｗを軸とする回転対称体をもつ
２）Ｈ(r)=０の時、ｒの分布は直径1/u(u≠0)の球面または平面(u=0)を成す
３）Ｈt(r)=０の時、ｒの分布は頂角2cut(ut≠0）の円錐面または平面(ut=0)を成す
４）Ｈ(r)=０とＨt(r)=０の時のｒの分布の交わりは円または平面を成す
式(32)を周波数領域に変換すると、
【００６７】
【数１９】
2008304555.tif 000021
【００６８】
を得る。よって音源ｒからs(t)への周波数応答Ｔ(r,w)は、
Ｔ(r,w)=H(r)+jωHt(r)    …(47)
となり、Ｈ(r)，Ｈt(r)が実数であればＴ(r,w)=0となる場合には
Ｈ(r)=0，Ｈt(r)=0    …(48)
となる。故に、式(47)からＳ(ω)=0となる零点分布は、周波数ωに依存せず、音源位置ｒのみに依存することが分かる。したがって、観測点における音圧の時間勾配とx,y方向の空間勾配が得られる時に、零感度領域（死点）を形成するには、ある瞬間においてｆ，ｆt，ｆx，ｆyの荷重和を取り、補償フィルタ処理（低域通過フィルタ処理）を施すだけでよい。
【００６９】
而して、本実施形態における死点形成手段２においては、集音センサ手段１から出力される同相成分Ｍ(t)、時間微分値Ｍt(t)、ｘ軸方向空間微分値Ｍx(t)、ｙ軸方向空間微分値Ｍy(t)をそれぞれ上述のｆ，ｆt，ｆx，ｆyに置き換え、これらの値を要素とするベクトルＭ＝(Ｍ(t)  Ｍt(t)  Ｍx(t)  Ｍy(t))Tを定義し、これらに対する荷重を要素とする係数ベクトルＷ＝(Ｗ  Ｗt  Ｗx  Ｗy)Tとの荷重和を演算した後、低域通過フィルタ処理を施すことによって、予め決められた任意位置に死点を形成する。具体的には、上述の指向特性Ｈ(r)，Ｈt(r)をそれぞれＨ1(ri)，Ｈ2(ri)と置き換えて下記のように定義し（但し、ｒiは音源ｉの位置ベクトル、ｎix，ｎiyはそれぞれｒi/|ｒi|のｘ成分とy成分である。）、
【００７０】
【数２０】
2008304555.tif 000022
【００７１】
さらに、図１において(p+jωq)-1で表した低域通過フィルタ２０のフィルタ係数p=WHH1(ri)，q=WHH2(ri)が何れもゼロとなるような係数ベクトルＷを選択することで音源ｉの位置ｒiに、周波数に依存しない死点が形成できる。そして、死点形成手段２の出力Ｏ(t)には死点に存在する音源から発せられる音の音圧が含まれない、言い換えると、死点に居る話者の音声を除く周囲騒音や残響音（以下、雑音という。）のみが含まれていることになる。
【００７２】
目的話者音声抽出手段３は、集音センサ手段１から出力される音圧（目的音＜死点に存在する話者の音声＞と雑音を含む音圧）Ｍ(t)と、死点形成手段２から出力される雑音成分Ｏ(t)とから、従来周知のスペクトル・サブトラクション法によって目的音Ｓ(t)を抽出する処理を行う（参考文献６参照）。まず、目的話者音声抽出手段３では同相成分Ｍ(t)及び雑音成分Ｏ(t)をフレーム分割部３０にて単位時間（フレーム時間）毎に分割し、分割された音圧Ｍ(t,k)及び雑音成分Ｏ(t,k)を高速フーリエ変換部３１で時間領域から周波数領域に変換する（但し、ｋはフレーム番号を示す）。そして、雑音成分Ｏ(f,k)の平均振幅μ(=E{|O(f,k|})を雑音平均振幅算出部３２で算出し、振幅算出部３３で算出した音圧Ｍ(f,k)の振幅値｜Ｍ(f,k)｜から雑音成分Ｏ(f,k)の平均振幅μを減算するとともに、減算した値（｜Ｍ(f,k)｜−μ）に、位相算出部３４で算出した音圧Ｍ(f,k)の位相（=exp{j∠M(f,k)}）を乗算することで雑音が含まれていない出力Ｓ(f,k)=(|M(f,k)|−μ)・exp{j∠M(f,k)}を取り出し、この出力Ｓ(f,k)を高速フーリエ逆変換して周波数領域から時間領域に戻すことで雑音が含まれない目的音Ｓ(t)のみを得ることができる。
【００７３】
このように本実施形態の音響入力装置によれば、死点形成手段２から出力される音圧（同相成分）Ｏ(f,k)には死点に存在する音源から発せられる音以外の音、すなわち、雑音の音圧のみが含まれ、一方、集音センサ手段１で検出される音圧（同相成分）Ｍ(f,k)には死点に存在する音源から発せられる音の音圧と雑音の音圧の双方が含まれており、目的話者音声抽出手段３によって死点に存在する音源から発せられる音の音圧Ｓ(f,k)のみを抽出するので、雑音が死点の方向から到来する場合においても死点に存在する音源から発せられる音のみを入力することができる。図３は目的音と雑音を含む音圧Ｍ(f)、目的音を含まない雑音のみの音圧Ｏ(f)、目的話者音声抽出手段３で抽出される目的音の音圧Ｓ(f)の周波数特性の一例を示しており、音圧Ｓ(f)に含まれる雑音成分が十分に抑圧されていることが判る。
【００７４】
ここで、本実施形態の音響入力装置をインターホン装置（ドアホン子器）に搭載すれば、周囲騒音の大きい環境下においても話者（来訪者）の音声のみを抽出して通話することができる。
【００７５】
ところで、本実施形態における集音センサ手段１では４本のマイクロホン１０Ａ〜１０Ｄをｘｙ平面上に配置しているが、図４に示すように平面視円形の振動板１３が中央の１点（中心）で支持された２軸直交型のジンバル構造を有するマイクロホン１２をマイクロホン１０Ａ〜１０Ｄの代わりに用いても構わない。振動板１３は全体が薄い円盤状であって、その中央部分における同心円上に二重の溝１４，１５が形成され、内側の溝１４を仕切る一対のビーム１４ａ，１４ａと、外側の溝１５を仕切る一対のビーム１５ａ，１５ａとを有し、支持棒１６に支持された点（中心）を支点として各ビーム１４ａ，１５ａのねじれによってｘ軸及びｙ軸の回りに回転可能となっている（参考文献７，８参照）。したがって、振動板１３上に４つの観測点Ａ，Ｂ，Ｃ，Ｄを設定し、各観測点における変位量をマイクロホン１０Ａ〜１０Ｄの出力に置き換えれば、同相成分Ｍ(t)、時間微分値Ｍt(t)、ｘ軸方向空間微分値Ｍx(t)、ｙ軸方向空間微分値Ｍy(t)を検出することが可能である。
【００７６】
（実施形態２）
本実施形態は、目的話者音声抽出手段３における抽出処理としてスペクトル・サブトラクション法の代わりに独立成分分析の手法を利用する点に特徴があり、その他の構成並びに動作は実施形態１と共通であるから、共通の構成要素には同一の符号を付して図示並びに説明を省略する。
【００７７】
独立成分分析(ICA:Independent Component Analysis)の目的は、複数の観測される変数を統計的に独立な変数の線形結合として表現することであり、観測変数から計算で求められる独立な変数が独立成分である。例えば、観測変数ベクトルをＸとし、この観測変数ベクトルＸが未知の独立変数ベクトルＳの線形結合で与えられると仮定すると、未知の混合行列をＡとして、Ｘ＝ＡＳの関係が成立する。そして、独立成分分析とは、独立成分及び混合行列に関する知識を一切利用せずに観測データのみから分離行列Ｗを用いてＹ＝ＷＸで求められる復元データベクトルＹの各成分が独立となるような分離行列Ｗを求める手法であり、理想的には分離行列Ｗが混合行列Ａの逆行列（Ｗ＝Ａ-1）となればよい。ここで、観測データが二次元の場合の独立成分分析のモデルを図５に示す。
【００７８】
ここで、目的音と雑音を含む音圧Ｍ(f,k)と、死点形成手段２より出力される雑音の音圧Ｏ(f,k)を成分に持つ観測行列をＸ(f,k)=[Ｍ(f,k) Ｏ(f,k)]Tとし、目的音の音圧Ｓ(f,k)と騒音源の音圧Ｎ(f,k)を成分に持つ行列をＳ(f,k)=[Ｓ(f,k) Ｎ(f,k)]Tとし、空間伝達行列をＡとすると、
【００７９】
【数２１】
2008304555.tif 000023
【００８０】
と表される。そして、目的話者音声抽出手段３においては、Ｓ(f,k)＝Ｗi+1(f,k)Ｘ(f,k)を満足する分離行列Ｗを適応的に同定することで目的音と雑音を含む音圧Ｏ(f,k)から目的音Ｓ(f,k)のみを分離して抽出する。この分離行列Ｗi+1(f)は、
【００８１】
【数２２】
2008304555.tif 000024
【００８２】
と表される（参考文献７参照）。ここで、diagは対角行列を示し、Φ(Y)はφ(y)=(1+exp(-y))-1で表されるシグモイド関数、若しくはφ(y)=3/4y11+25/4y9+14/4y7+47/4y5+29/4y3などの多項式で近似される非線形ベクトル関数である。例えば、上記式をシグモイド関数によって近似した場合、
【００８３】
【数２３】
2008304555.tif 000025
【００８４】
となる。
【００８５】
ところで、観測変数の数が独立変数の数よりも多いならば、観測変数は線形従属であり、分離行列は低次元化を行う行列となる。また、変数が互いに独立であれば、それらは無相関であるため、分離行列は変数を無相関化する行列でもある。無相関化とそれに伴う低次元化を同時に行う統計的手法に主成分分析(PCA:Principal Component Analysis)があり、独立成分分析の前処理として主成分分析が利用されることがある。
【００８６】
そこで本実施形態においても、雑音として反射や残響がある場合に、観測行列Ｘ(f,k)を低次元化するために独立成分分析の前処理として主成分分析を行うようにしてもよい。なお、変数が互いに独立であれば無相関でもあるため、主成分分析によって無相関化と低次元化を同時に行うことができる。
【００８７】
例えば、音源の数がｒである場合に、ｍ次元の観測行列Ｘ(f,k)の特異値分解は下記のようになる。
【００８８】
【数２４】
2008304555.tif 000026
【００８９】
また、各音源スペクトルS1,S2,…,Srとしたとき、
【００９０】
【数２５】
2008304555.tif 000027
【００９１】
ここで、
【００９２】
【数２６】
2008304555.tif 000028
【００９３】
という変換を考えると、その分散は
【００９４】
【数２７】
2008304555.tif 000029
【００９５】
となる。分散行列が対角行列となることから、変換後の変数は互いに無相関となる。また、ｍ次元の観測行列Ｘをｒ次元に圧縮することができ、後段の独立成分分析における処理量を低減することが可能となる。
【００９６】
（実施形態３）
実施形態１，２では目的音の音源（例えば、話者）の位置が既知であることを前提として、当該位置に死点を形成することで音源から発せられる目的音のみを集音している。しかしながら、ドアホン子器のように目的音の音源（来訪者）の位置が一意に定まらない場合も多い。一方、雑音（周囲騒音並びに残響音）が非常に少ない環境下においては、既に説明した時空間勾配法による音源定位の技術を用いて音源の位置を推定することができ、音源の位置が一意に定まらない場合においても、音源の位置を推定して当該位置に死点を形成することで音源から発せられる目的音のみを入力することが可能である。
【００９７】
そのために本実施形態では、集音センサ手段１で検出される同相成分Ｍ(t)の瞬時パワーＰM(t)を死点形成手段２から出力される同相成分Ｏ(t)の瞬時パワーＰO(t)で除した値（＝ＰM(t)／ＰO(t)）が所定のしきい値δ以上となるときに集音センサ手段１で検出される音圧（同相成分Ｍ）、時間微分値Ｍt、空間微分値Ｍx,Ｍyに基づいて推定される位置に前記音源位置を変更する音源位置変更手段を備えている。この音源位置変更手段は死点形成手段２とほぼ共通の処理を行うものであるから死点形成手段２で兼用することも可能である。そして、上述の条件が満たされるときに死点形成手段２で音源位置を推定し、推定された音源位置がそのときの死点の位置と異なっていれば、推定された音源位置に死点を変更するのである。
【００９８】
このように本実施形態によれば、音源の位置が変動した場合でも音源位置変更手段によって音源位置が変更されるために音源位置と死点がずれることがなく、その結果、音源が移動する場合においても当該音源から発せられる音のみを入力することができる。
＜参考文献一覧＞
参考文献１：安藤  繁  「画像の時空間微分算法を用いた速度ベクトル分布計測システム」  計測自動制御学会論文集  22-12，1330/1336(1986)
参考文献２：安藤  繁・篠田  裕之・小川  勝也・光山  訓  「時空間勾配法に基づく３次元音源定位センサシステム」  計測自動制御学会論文集  第２９巻第５号，p520~528，1993
参考文献３：N. Ono, T. Arita, Y. Senjo, and S. Ando, “Directivity steering principle for biomimicry silicon microphone”, Proc. Int. Conf. Solid State Sensors, Actuators, and Microsystems (Transducers '05), pp. 792-795, 2005.
参考文献４：小野, 安藤, “音場の計測と指向性制御, 第22回センシングフォーラム資料, pp. 305-310, 2005. 
参考文献５：小野, 有田, 千條, 安藤, “時空間勾配計測に基づく指向性制御と音源分離の理論, 日本音響学会2005年春季研究発表会講演論文集, 2-6-13, pp. 607-608, 2005.
参考文献６：S.F.Boll "Suppression of Acoustic Noise in Speech. using Spectral Subtraction" IEEE Trans.on.Acoustics,Speech and Signal Processing Vol.ASSP-27,No.2,pp.113-1,1979
参考文献７：小野 順貴,斎藤 章人,安藤 繁「ヤドリバエを模倣した超小型音源定位センサの理論と実験（第２報）」,第１９回センシングフォーラム,pp.379-382,2002
参考文献８：小野 順貴,斎藤 章人,安藤 繁「ヤドリバエを模倣した微分検出型音源定位センサの理論と実験」,聴覚研究会資料,pp.187-192,2002
【図面の簡単な説明】
【００９９】
【図１】本発明の実施形態１を示すブロック図である。
【図２】同上における時空間勾配法を説明するための説明図である。
【図３】同上の説明図である。
【図４】同上におけるジンバル構造型のマイクロホンを示し、（ａ）は振動板の平面図、（ｂ）は断面図である。
【図５】本発明の実施形態２における目的話者音声抽出手段の説明図である。
【符号の説明】
【０１００】
１  集音センサ手段
２  死点形成手段
３  目的話者音声抽出手段
【図１】
2008304555.tif 000030
【図２】
2008304555.tif 000031
【図３】
2008304555.tif 000032
【図４】
2008304555.tif 000033
【図５】
2008304555.tif 000034