"スピーカの位置とその総数Nと、制御点の位置とその総数Mと、上記制御点における複素振幅の値を指定する制御量Cmと、を入力とし、フィルタ係数Hnωを要素とする上記N行×1列の複素行列H、スピーカnから上記制御点mまでの伝達関数Gmnωを要素とする上記M行×上記N列の複素行列G、その複素行列Gの複素共役転置行列をG†、上記制御量Cmを要素とするM行×1列の複素行列をCとしたときに、次の線形方程式を上記複素行列Hについて解くことでフィルタ係数を算出するフィルタ係数生成方法であって、上記スピーカの総数Nと上記制御点の総数Mの大小関係によって、上記複素行列Hを次式、で求めることを特徴とするフィルタ係数生成方法。請求項1に記載したフィルタ係数生成方法において、さらに上記制御点における複素振幅の値を拘束する拘束条件を入力とし、上記制御点内の当該拘束条件を設けた制御点における制御量をC1=1,0,0,...,0Tと表し、当該制御点において成り立つ線形方程式を次式、上記拘束条件を設けた制御点以外の上記制御点において成り立つ線形方程式を次式、とした時に、上記拘束制御点における上記拘束条件が成り立つ条件の下で、複素行列Hの複素共役転置行列をH†、空間相関行列KをK=G2G2Tとして、次式に示す最小二乗解を、フィルタ係数Hnωに関して解いたフィルタ係数Hnωを、求めることを特徴とするフィルタ係数生成方法。入力されたディジタル音響信号に、任意の位置に配置される2個以上のスピーカa1,...,aNただし、Nは2以上の正数。各スピーカをanとする。にそれぞれ対応するディジタルフィルタ係数を時間領域で畳み込み、上記各スピーカに供給する音声信号を生成する音声信号生成過程を含む局所再生方法であって、上記スピーカa1,...,aNと、任意の位置における複素振幅の値を指定する制御点である上記複素振幅の値を0とする上記スピーカanからの距離がそれぞれ異なる位置に配置される3個以上の抑圧制御点と、上記スピーカと上記抑圧制御点との間の位置に上記複素振幅の値を1とする制御点である再生制御点とが設けられ、上記ディジタルフィルタ係数は請求項1又は2に記載したフィルタ係数生成方法で生成された値であることを特徴とする局所再生方法。請求項3に記載した局所再生方法において、上記スピーカa1,...,aNは、各スピーカanを直線で結ぶと平面が形成されるように配置された3個以上のスピーカであって、任意の位置における複素振幅の値を指定する制御点である上記複素振幅の値を0とする上記スピーカanからの距離がそれぞれ異なる位置に配置される4個以上の抑圧制御点と、上記スピーカanと上記抑圧制御点との間の位置に上記複素振幅の値を1とする制御点である再生制御点とが設けられ、上記ディジタルフィルタ係数は請求項1又は2に記載したフィルタ係数生成方法で生成された値であることを特徴とする局所再生方法。請求項4に記載した局所再生方法において、上記3個以上のスピーカa1,...,aNは、2次元平面上に配置され、上記制御点も当該2次元平面上に設けられることを特徴とする局所再生方法。スピーカの位置とその総数Nと、制御点の位置とその総数Mと、上記制御点における複素振幅の値を指定する制御量Cmと、を入力とし、フィルタ係数Hnωを要素とする上記N行×1列の複素行列H、スピーカnから上記制御点mまでの伝達関数Gmnωを要素とする上記M行×上記N列の複素行列G、その複素行列Gの複素共役転置行列をG†、上記制御量Cmを要素とするM行×1列の複素行列をCとしたときに、次の線形方程式を上記複素行列Hについて解くことでフィルタ係数を算出するフィルタ係数生成装置であって、上記スピーカの総数Nと上記制御点の総数Mの大小関係によって、上記複素行列Hを次式、で求めることを特徴とするフィルタ係数生成装置。任意の位置に配置される2個以上のスピーカと、入力されたディジタル音響信号に上記各スピーカに対応するディジタルフィルタ係数を時間領域で畳み込み上記スピーカに供給する供給信号を生成する信号供給部と、を備えた局所再生装置であって、上記スピーカと、任意の位置における複素振幅の値を指定する制御点であり上記複素振幅の値を0とする上記スピーカからの距離がそれぞれ異なる位置に設けられる3個以上の抑圧制御点と、上記スピーカと上記抑圧制御点との間の位置に上記複素振幅の値を1とする制御点である再生制御点とが設けられ、上記ディジタルフィルタ係数は請求項1又は2に記載したフィルタ係数生成装置で生成された値であることを特徴とする局所再生装置。請求項1又は2に記載したフィルタ係数生成方法、又は請求項3乃至5の何れかに記載した局所再生方法をコンピュータに実行させるためのプログラム。しかし、従来の方法では、複数のスピーカを直線状に配置して構成したスピーカアレーを対向させて、その中間点中心を再生領域と指定しているので各スピーカからの音の指向特性はその中間点を中心とした同心円状に制御される。そのため、再生領域を任意の場所に設定することが出来ない課題があった。この発明は、この課題に鑑みてなされたものであり、再生領域を任意の場所に設定することが可能な局所再生装置とフィルタ係数生成装置と、それらの方法とプログラムを提供することを目的とする。この発明のフィルタ係数生成装置は、スピーカの位置とその総数Nと、制御点の位置とその総数Mと、上記制御点における複素振幅の値を指定する制御量Cmと、を入力とし、フィルタ係数Hnωを要素とするN行×1列の複素行列H、スピーカnから制御点mまでの伝達関数Gmnωを要素とするM行×N列の複素行列G、その複素行列Gの複素共役転置行列をG†、制御量Cmを要素とするM行×1列の複素行列をCとしたときに、次の線形方程式を複素行列Hについて解くことでフィルタ係数を算出するフィルタ係数生成装置であって、スピーカの総数Nと制御点の総数Mの大小関係によって、複素行列Hを次式、で求める。また、この発明の局所再生装置は、2個以上のスピーカと、信号供給部とを備える。2個以上のスピーカは任意の位置に配置される。そして、任意の位置における複素振幅の値を0とする各スピーカからの距離がそれぞれ異なる位置に設けられる3個以上の抑圧制御点と、各スピーカと抑圧制御点との間の位置に複素振幅の値を1とする再生制御点とが設けられ、信号供給部内のディジタルフィルタのフィルタ係数は、上記したフィルタ係数生成装置で生成された値である。この発明のフィルタ係数生成装置によれば、スピーカの総数Nと制御点の総数Mの大小関係に基づいて、任意に配置されたスピーカにより任意の場所に再生領域を設定することができるフィルタ係数を生成することが出来る。また、この発明の局所再生装置は、この発明のフィルタ係数生成装置で求めたフィルタ係数を用いることで、2個以上のスピーカと3個以上の抑圧制御点との間の任意の位置に再生制御点を設けることが出来る。以下、この発明の実施の形態を図面を参照して説明する。複数の図面中同一のものには同じ参照符号を付し、説明は繰り返さない。フィルタ係数生成装置図1に、この発明のフィルタ係数生成装置100の機能構成例を示す。フィルタ係数生成装置100は、後述する局所再生装置が具備するディジタルフィルタのフィルタ係数を生成するものであり、スピーカの位置とその総数Nと、制御点の位置とその総数Mと、各制御点における複素振幅の値を指定する制御量Cmと、制御点における複素振幅の値を拘束する拘束条件の有無と、を入力として、スピーカの総数Nと制御点の総数Mに対応したフィルタ係数を自動的に生成する点で新しい。フィルタ係数生成装置100は、伝達関数算出部10と、フィルタ係数算出部20と、で構成される。伝達関数算出部10は、スピーカの位置Xn=xn,yn,znと、制御点の位置Xm=xm,ym,zmを入力として、スピーカから各制御点までの伝達関数Gmnωを次式によって求める。ここで、∥Xm-Xn∥はスピーカnから制御点mまでのユークリッド距離を表す。cは音速である。ωは角周波数であり周波数fとの関係はω=2πfとして表される。フィルタ係数算出部20は、その伝達関数Gmnωと、局所再生装置のスピーカの総数Nと、制御点の総数Mと、制御点における複素振幅の値を指定する制御量Cmと、拘束条件の有無を入力とする。拘束条件とは、制御点における観測音圧複素振幅が必ず制御量となる条件のことであり、再生音を必ず再生したい点と、音圧を抑圧したい点がある場合は、この拘束条件で設定する。フィルタ係数算出部20の動作を、拘束条件が無い場合で説明する。フィルタ係数算出部20では、スピーカの位置とその総数Nと、制御点の位置とその総数Mと、上記制御点における複素振幅の値を指定する制御量Cmと、制御点における複素振幅の値を拘束する拘束条件の有無と、を入力とし、ディジタルフィルタに与えるフィルタ係数Hnωを算出して出力する。制御量Cmは以下のように決定する。拘束条件が無い場合ステップS21のNO、フィルタ係数Hnωを要素とするN行×1列の複素行列H、スピーカnから制御点mまでの伝達関数Gmnωを要素とするM行×N列の複素行列G、制御量Cmを要素とするM行×1列の複素行列をCとしたときに、式5の線形方程式を複素行列Hについて解くことでフィルタ係数Hnωを算出する。そこで更に、スピーカの総数Nと制御点の総数Mの大小関係によって、複素行列Hを次式で求める。図2に示すフィルタ係数決定装置100の動作フローを参照してその動作を説明する。伝達関数算出部10は、スピーカの位置Xn=xn,yn,znと、制御点の位置Xm=xm,ym,zmを入力として、各スピーカnと制御点mの位置情報から、スピーカと制御点間の距離∥Xm-Xn∥を算出するステップS11。そして、スピーカから各制御点までの伝達関数Gmnωを式3によって求めるステップS12。そして拘束条件の無い場合ステップS21のNO、フィルタ係数算出部20は、スピーカの総数Nと減衰制御点の総数Mとの数が等しい時ステップS22のYES、フィルタ係数Hnωを、Hnω=G-1C式61として算出するステップS23。減衰制御点の総数Mが、スピーカの総数Nよりも多い場合ステップS24のYES、フィルタ係数Hnωは、最小二乗法による近似解として求められる。上記した式5の右辺と左辺の誤差の2乗をε式7とし、この誤差の2乗εが最小となるようなHを求める。式8は∂ε∂H=0とすることで解くことができる式62ステップS25。制御点の総数Mが、スピーカの総数Nよりも少ない場合ステップS24のNO、フィルタ係数Hnωは、自身のノムル∥H∥が最小となるような解として求められる。上記した式4の解がH=G-Cと表されるとき、G-はGの一般逆行列として定義される。ノムル∥H∥が最小となるとき、Gは以下の式を満たすステップS26。以上より、したがって、以上の処理は、全てのフィルタ係数Hnωの算出が終了するまでステップS27のYES繰り返される。拘束条件がある場合ステップS21のYES、制御点内の当該拘束条件を設けた制御点における制御量をC1=1,0,0,...,0Tと表し、当該制御点において成り立つ線形方程式を次式、この拘束条件を設けた制御点以外の制御点において成り立つ線形方程式を次式、とした時に、拘束制御点における上記拘束条件が成り立つ条件の下で、複素行列Hの複素共役転置行列をH†、空間相関行列KをK=G2G2Tとして、次式に示す最小二乗解を、フィルタ係数Hnωに関して解いたフィルタ係数Hnωを求める。このような条件付き最小化問題は、ラグランジェの未定定数法を用いることで解くことができる。そこで、この方法に基づいて最小化すべき評価関数QHを次のように定義する。ここでΛはM個の未定係数からなる列ベクトルである。評価関数をHで微分した結果が零に等しいと置くことで、フィルタ係数を求めることができ。Λを決定するために、G1H=C1に代入すると、式16を式15に代入することで、フィルタ係数が以下のように与えられるステップS28。なお、拘束条件がない場合に、スピーカの総数Nと制御点の数Mの大小関係によって、フィルタ係数Hnωを要素とする複素行列Hを求める計算式を変更する例を説明したが、拘束条件がない場合、フィルタ係数を生成しないようにしても良いステップS21のNO~終了の破線。また、拘束条件がある場合にフィルタ係数を生成しないようにしたフィルタ係数生成装置としても良いステップS21のYES~終了の破線。つまり、拘束条件の有無にそれぞれ対応する専用の装置として、この発明のフィルタ係数生成装置を構成しても良い。局所再生装置図3に、局所再生装置200の機能構成例を示す。局所再生装置200は、任意の位置に配置される放音方向も任意である2個以上のスピーカa1,...,aNただし、Nは2以上の正数。各スピーカをanとする。と、入力されたディジタル音響信号に各スピーカanに対応するディジタルフィルタ係数を時間領域に畳み込み各スピーカanに供給する供給信号を生成する信号供給部40とを備える。そして、任意の位置における複素振幅の値を指定する制御点であるその複素振幅の値を0とする各スピーカanからの距離がそれぞれ異なる位置に配置される3個以上の抑圧制御点×図中の抑圧制御点は○付きと、各スピーカanと抑圧制御点との間の任意の位置に複素振幅の値を1とする制御点である再生制御点◎とが設けられ、上記ディジタルフィルタ係数は、上記したフィルタ係数生成装置100で生成したフィルタ係数である。抑圧制御点×は、スピーカanと再生制御点◎を取り囲むように設けても良い。スピーカanは、2次元平面上の任意の位置に配置しても良いし、3次元空間上の任意の位置に配置しても良い。なお、スピーカanの数が2個の場合は、抑圧制御点の数は3個以上、スピーカanの数がN個の場合でも抑圧制御点は3個以上数を設ければ良い。信号供給部40は、各スピーカanに、それぞれ接続されるAD変換器41nとディジタルフィルタ42nとDA変換器43nとAMP44nの直列構成が、スピーカanの個数分並列に用意されて構成される。そして、複数のAD変換器41nに一つの音響信号が入力される。ディジタルフィルタ42nのそれぞれには、上記したフィルタ係数生成装置100で算出されたディジタルフィルタ係数Hnωが設定される。各ディジタルフィルタ42nは、入力されるディジタル音響信号にフィルタ係数Hnωを畳み込むことで、各スピーカanを駆動するためのディジタル駆動信号Dnωを生成し、この駆動信号に対応するDA変換器43nに出力する。なお、駆動信号Dnωは次式で表せる。ディジタルフィルタ係数Hnωは、Q個の角周波数領域ωq毎に生成される。周波数の帯域分割数をQとし、q=0,1,...,Q-1とした場合、角周波数帯域は次式で表せる。ここでfsは、AD変換器41nのサンプリング周波数である。フィルタ係数生成装置100は、角周波数帯域ωqにおけるフィルタ係数Hnωqを算出し、それを逆フーリエ変換したものを各ディジタルフィルタ42nにおけるフィルタ係数として設定する。なお、周波数の帯域分割数Qの値を大きくとればフィルタとしての精度は高くなるが計算コストが大きくなる。分割数Qは任意に設定することができ、例えば数百程度の値とする。また、2のべき乗数に設定すると演算が高速化できる。なお、AD変換器41nは、入力される音響信号がディジタル信号であれば無くても良い。また、周波数帯域は角周波数でなくても良い。図4に、局所再生装置300の機能構成例を示す。3個以上のスピーカa1,...,aNただし、Nは3以上の正数。各スピーカをanとする。は、直線で結ぶと平面が形成されるように配置される。3個以上のスピーカanが配置される座標系の原点αから所定の半径Rの円周上に4個以上の抑圧制御点×が設けられる。そして、3個以上のスピーカanと抑圧制御点との間に再生制御点◎が設けられる。4個以上の抑圧制御点×は、3次元空間上に配置される。抑圧制御点×は、スピーカanと再生制御点◎を取り囲むように設けられる。3個のスピーカanは、例えば図5に示すように3次元座標系のxy平面の角に配置するようにしても良い。その場合の抑圧制御点×は、同一xy平面上の原点αから所定半径r上の任意の位置に設けても良い。また、再生制御点◎は同一xy平面上に設けられる。この例では、同一xy平面上の再生制御点の音圧の制御精度を高めることが可能である。シミュレーション結果この発明の局所再生装置の性能を確認する目的で、コンピュータシュミレーションを行った。実験条件を表1と図6に示す。1組が12個のスピーカで構成されるスピーカユニットの中心は、xy座標系の0.40、-0.40、00.4、0-0.4の位置に放音方向を原点に向けた向きに配置される。そして、再生制御点◎は、0.10.25の位置に設けた。図7に、y=025上の音圧を測定した結果を示す。横軸はx軸上の位置、縦軸は音圧dBである。再生制御点◎のx座標:01上の音圧が0dBと最も高く制御されていることが分かる。このように、この発明の局所再生装置によれば音響信号の再生領域を特定エリア内に限定することが可能である。なお、フィルタ係数生成装置100及び特定エリア音響再生装置200,300,400は、例えばROM、RAM、CPU等で構成されるコンピュータに所定のプログラムが読み込まれて、CPUがそのプログラムを実行することで実現されるように構成してもよい。その場合、その処理内容を記述したプログラムは、コンピュータで読み取り可能な任意の記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリがある。より具体的には、例えば、磁気記録装置として、ハードディスク装置、フレキシブルディスク、磁気テープ等を、光ディスクとして、DVDDigitalVersatileDisc、DVDRAMRandomAccessMemory、CDROMCompactDiscReadOnlyMemory、CDRRecordableRWReWritable等を、光磁気記録媒体として、MOMagnetoOpticaldisc等を、半導体メモリとしてEEPROMElectronicallyErasableandProgrammableReadOnlyMemory等を用いることができる。また、このプログラムの流通は、例えば、そのプログラムを記録したDVD、CDROM等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記録装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。また、各手段は、コンピュータ上で所定のプログラムを実行させることにより構成することにしてもよいし、これらの処理内容の少なくとも一部をハードウェアとして実現することとしてもよい。"