(57)【要約】【課題】少ない演算量で、かつ所望の周波数特性を確保しながら、信号を高域と低域とに分割することのできる帯域分割フィルターを提供する。【解決手段】入力信号を１サンプルずつ遅延させる遅延器１−１〜１−ｍと、遅延器１−１〜１−ｍのそれぞれの入出力に得られる遅延量の異なる信号にそれぞれ係数を乗算する乗算器２−０〜２−ｍと、乗算器２−０〜２−ｍの出力を加算して第１のフィルター出力とする加算器３、４と、遅延器１−１〜１−ｍのそれぞれの入出力のうち遅延量が中央値となる信号に、対応する乗算器２−（ｍ／２）の係数ａｍ／２を１から減算した値を係数として乗算した値と、中央値となる信号以外の遅延器１−１〜１−ｍのそれぞれの入出力に、対応する乗算器の係数の逆符号の値を係数として乗算した値との和を生成して第２のフィルター出力とする乗算器５、６および加算器７とを有する。【選択図】図１
2010252010.tif 000002
【特許請求の範囲】
【請求項１】
入力信号を１サンプルずつ遅延させる複数の遅延器と、
この複数の遅延器のそれぞれの入出力に得られる遅延量の異なる信号にそれぞれ係数を乗算する複数の乗算器と、
この複数の乗算器の出力を加算して第１のフィルター出力とする加算手段と、
上記複数の遅延器のそれぞれの入出力のうち遅延量が中央値となる信号に、上記複数の乗算器のうち対応する乗算器の係数を１から減算した値を係数として乗算した値と、上記中央値となる信号以外の上記複数の遅延器のそれぞれの入出力に上記複数の乗算器のうちそれぞれ対応する乗算器の係数の逆符号の値を係数として乗算した値との和を生成して第２のフィルター出力とする回路手段と
を有することを特徴とする帯域分割フィルター。
【請求項２】
請求項１記載の帯域分割フィルターにおいて、前記複数の乗算器がそれぞれ乗算する係数は、前記第１のフィルター出力が高域通過特性または低域通過特性となる係数に設定され、前記第２フィルター出力は低域通過特性または高域通過特性となることを特徴とする帯域分割フィルター。
【請求項３】
請求項１または２記載の帯域分割フィルターにおいて、
前記加算手段は、前記中央値となる信号以外の信号に関して前記複数の乗算器の対応する出力を加算する第１の加算器と、この第１の加算器の出力に前記中央値となる信号に関して前記複数の乗算器の対応する出力を加算する第２の加算器とを有し、
前記回手段は、前記中央値となる信号に、前記複数の乗算器のうち対応する乗算器の係数を１から減算した値を係数として乗算する減算値乗算器と、上記第１の加算器の出力を逆符号にする逆符号乗算器と、上記減算値乗算器の出力と上記逆符号乗算器の出力とを加算する第３の加算器とを有する
ことを特徴とする帯域分割フィルター。
【請求項４】
コンピュータに、
入力信号を１サンプルずつ遅延させる複数の遅延器、
この複数の遅延器のそれぞれの入出力に得られる遅延量の異なる信号にそれぞれ係数を乗算する複数の乗算器、
この複数の乗算器の出力を加算して第１のフィルター出力とする第１の加算手段、
この複数の乗算器の出力を加算して第１のフィルター出力とする加算手段、および
上記複数の遅延器のそれぞれの入出力のうち遅延量が中央値となる信号に、上記複数の乗算器のうち対応する乗算器の係数を１から減算した値を係数として乗算した値と、上記中央値となる信号以外の上記複数の遅延器のそれぞれの入出力に上記複数の乗算器のうちそれぞれ対応する乗算器の係数の逆符号の値を係数として乗算した値との和を生成して第２のフィルター出力とする回路手段
の各機能を実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、帯域分割フィルターおよびプログラムに関する。
【背景技術】
【０００２】
信号帯域を高域と低域に分割する帯域分割フィルターとして、ハイパスフィルター（ＨＰＦ）およびローパスフィルター（ＬＰＦ）が組み合わせて用いられる。以下では、可聴信号、特に音声信号を、高域専用再生スピーカーと低域専用再生スピーカーとに振り分ける場合を例に説明する。
【０００３】
図２は、帯域分割フィルターを有限インパルス応答（ＦＩＲ）フィルターで実現する構成例を示す図である。この帯域分割フィルターは、１チャンネルあたり、ＨＰＦとＬＰＦのふたつのフィルターで構成される。
【０００４】
すなわち、ＨＰＦは、遅延器１１−１〜１１−ｍ、乗算器２１−０〜２１−ｍおよび加算器３１からなるＦＩＲフィルターにより構成される。ＬＰＦは、遅延器１２−１〜１２−ｍ、乗算器２２−０〜２２−ｍおよび加算器３２からなるＦＩＲフィルターにより構成される。
【０００５】
遅延器１１−１〜１１−ｍはそれぞれ、入力信号を１サンプルずつ遅延させる。乗算器２１−０〜２１−ｍは、遅延器１１−１〜１１−ｍの入出力（タップ）に得られる遅延量の異なる信号に、それぞれ係数ａ０〜ａｍを乗算する。加算器３１は、乗算器２１−０〜２１−ｍの出力を加算し、加算器３２は、乗算器２２−０〜２２−ｍの出力を加算する。係数ａ０〜ａｍは、加算器３１の出力に低域通過特性が得られるように設定される。
【０００６】
遅延器１２−１〜１２−ｍはそれぞれ、入力信号を１サンプルずつ遅延させる。乗算器２２−０〜２２−ｍは、遅延器１２−１〜１２−ｍのそれぞれの入出力（タップ）に得られる遅延量の異なる信号に、それぞれ係数ｂ０〜ｂｍを乗算する。加算器３２は、乗算器２２−０〜２２−ｍの出力を加算する。係数ｂ０〜ｂｍは、加算器３２の出力に高域通過特性が得られるように設定される。
【０００７】
図２に示す個々のＦＩＲをｍ＝２０４８構成（タップ数２０４９）で実現すると、１チャンネルあたり、（２０４８＋１）×２＝４０９８回の乗算が必要となる。また、ステレオ音声ではその２倍、さらに６チャンネル出力だと、４０９８×６＝２４５８８回の乗算が必要となる。
【０００８】
図２に示す帯域分割フィルターは、実際には、信号処理プロセッサー（ＤＳＰ）により実現される。しかし、２４５８８回の乗算を行う能力のあるＤＳＰは高価であり、他のイコライザーなどの処理と同時に使用することが、処理能力的は難しい。ＦＩＲのタップ数を減らすことで演算回数は減るが、意図する周波数特性を確保することができなくなる。また、遅延器１１−１〜１１−ｍ、１２−１〜１２−ｍとしてのバッファが大量に必要であり、ＤＳＰに外付けＲＡＭ（ランダムアクセスメモリ）を付加したり、内蔵ＲＡＭの多いＤＳＰを使用する必要がある。
【０００９】
演算量の少ない帯域分割フィルターとして、特許文献１には、元データとＬＰＦの出力との差分をＬＰＦの出力とすることが開示されている。また、特許文献２には、標本化周波数ｆＳの半分で演算する構成が開示されている。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２００４−２０１１０９号公報
【特許文献２】特開平０５−０６３５１０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１１】
特許文献１に開示の技術は、平坦な周波数特性からＬＰＦの周波数を差し引くことで、ＨＰＦの周波数特性を得るものである。ＬＰＦにより処理された周波数帯は、演算精度によって特性が劣化し、ＨＰＦの通過周波数帯ではその特性がさらに劣化してしまう。特許文献２に記載の技術は、標本化周波数ｆＳを半分にしてしまうので、高域の表現力が無くなってしまう。
【００１２】
本発明は、このような課題を解決し、少ない演算量で、かつ所望の周波数特性を確保しながら、信号を高域と低域とに分割することのできる帯域分割フィルターを提供することを目的とする。
【課題を解決するための手段】
【００１３】
本発明の帯域分割フィルターは、入力信号を１サンプルずつ遅延させる複数の遅延器と、この複数の遅延器のそれぞれの入出力に得られる遅延量の異なる信号にそれぞれ係数を乗算する複数の乗算器と、この複数の乗算器の出力を加算して第１のフィルター出力とする加算手段と、複数の遅延器のそれぞれの入出力のうち遅延量が中央値となる信号に、複数の乗算器のうち対応する乗算器の係数を１から減算した値を係数として乗算した値と、中央値となる信号以外の複数の遅延器のそれぞれの入出力に複数の乗算器のうちそれぞれ対応する乗算器の係数の逆符号の値を係数として乗算した値との和を生成して第２のフィルター出力とする回路手段と有することを特徴とする。
【発明の効果】
【００１４】
本発明によれば、少ない演算量で、かつ所望の周波数特性を確保しながら、信号を高域と低域とに分割することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態に係る帯域分割フィルターの構成を示す図である。
【図２】公知の帯域分割フィルターの構成例を示す図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施の形態について、図面を参照して説明する。
【００１７】
[構成]
図１は本発明の実施の形態に係る帯域分割フィルターの構成を示す図である。この帯域分割フィルターは、遅延器１−１〜１−ｍ（ｍは偶数）、乗算器２−０〜２−ｍ、５、６および加算器３、４、７を有する。遅延器１−１〜１−ｍは、入力信号を１サンプルずつ遅延させる。乗算器２−０〜２−ｍは、遅延器１−１〜１−ｍのそれぞれの入出力に得られる遅延量の異なる信号に、それぞれタップ係数ａ０〜ａｍを乗算する。加算器３、４は、加算手段を構成し、乗算器２−０〜２−ｍの出力を加算して第１のフィルター出力とする。
【００１８】
乗算器５、６および加算器７は、複数の遅延器１−１〜１−ｍのそれぞれの入出力のうち遅延量が中央値となる信号、すなわち遅延器１−（ｍ／２）と遅延器１−（ｍ／２＋１）との間の信号に、対応する乗算器２−（ｍ／２）のタップ係数ａｍ／２を１から減算した値１−ａｍ／２を係数として乗算した値と、遅延量が中央値となる信号以外の複数の遅延器１−１〜１−ｍのそれぞれの入出力にそれぞれ対応する乗算器２−０〜２−（ｍ／２−１）、２−（ｍ／２＋１）〜２−ｍのタップ係数ａ０〜ａｍ／２−１、ａｍ／２−１ａｍの逆符号の値を係数として乗算した値との和を生成して第２のフィルター出力とする回路手段として動作する。
【００１９】
上述したように、乗算器２−０〜２−ｍの出力を加算して第１のフィルター出力とする加算手段として、この実施の形態では、加算器３、４を有する。加算器３は、遅延値が中央値となる信号以外の信号に関して、対応する乗算器２−０〜２−（ｍ／２−１)、２−（ｍ／２＋１)〜２−ｍの出力を加算する。加算器４は、加算器３の出力に、遅延値が中央値となる信号に関して、対応する乗算器２−（ｍ／２）の出力を加算する。一方、乗算器５は、減算値乗算器として動作し、遅延器１−（ｍ／２）と遅延器１−（ｍ／２＋１）との間の信号に、乗算器２−（ｍ／２）のタップ係数ａｍ／２を１から減算した値１−ａｍ／２を、タップ係数として乗算する。乗算器６は、逆符号乗算器として動作し、加算器４の逆符号にする。加算器７は、乗算器５、６の出力を加算する。
【００２０】
ここで、乗算器２−０〜２−ｍがそれぞれ乗算する係数ａ０〜ａｍは、加算器４の出力が低域通過特性となる係数に設定される。このとき、加算器７の出力は、高域通過特性となる。加算器４の出力が高域通過特性となるようにタップ係数ａ０〜ａｍを設定することもでき、その場合には、加算器７の出力は、低域通過特性となる。
【００２１】
［ＦＩＲフィルターの設計］
ここで、本発明に係るＦＩＲフィルター設計の原理を説明する。ここでは、窓関数法を使用し、ＬＰＦとＨＰＦとでカットオフ周波数ＦＣを同じ周波数とする。
【００２２】
フィルターの周波数特性をＧ（ｅｊωｔ）とし、標本化角周波数をωＳとし、遮断周波数をωＣとして、フーリエ級数展開を行うと、
【数１】
2010252010.tif 000003
となる。これをｇｎで解くと、
【数２】
2010252010.tif 000004
となる。
【００２３】
有限サンプルで打ち切るために、カイザー窓などを使用し、ＦＩＲのタップ数に収める。窓の幅をＯｒｄｅｒとすると、有限サンプルのｇｎは、
【数３】
2010252010.tif 000005
と表される。窓関数をｗｉｎｄｏｗとし、数２の式に掛け合わせると、ｇｎは、
【数４】
2010252010.tif 000006
となる。
【００２４】
ところで、以下のようにカットオフ周波数（ｆc）とｇｎを加工することで、ＬＰＦとＨＰＦとの変換を行うことができる。
【数５】
2010252010.tif 000007
ここで、
【数６】
2010252010.tif 000008
と置き換え、ｇｎとの値の関連を見つけると、ｇ０以外では逆符号となる。また、
【数７】
2010252010.tif 000009
と置き換えることで、ＬＰＦまたはＨＰＦにおけるｇ０が、ＨＰＦまたはＬＰＦにおける１−ｇ０となることがわかる。
【００２５】
以上の窓関数法でＬＰＦまたはＨＰＦに関して求めたｇｎに対して、図１におけるタップ係数ａ０〜ａｍを設定する。これにより、加算器４の出力がＨＰＦ出力またはＬＰＦ出力となる。一方、ＬＰＦまたはＨＰＦのタップ係数は、ｎ＝０に相当するタップ係数が１−ａｍ／２、ｎ≠０に相当するタップ係数が−ａ０〜−ａｍ／２−１、−ａｍ／２＋１〜−ａｍとなる。したがって、加算器７の出力が、ＬＰＦ出力またはＨＰＦ出力となる。
【００２６】
［実施の形態の効果］
ＦＩＲフィルターを用いた帯域分割処理は、演算負荷が大きいので、これまでは、高価なＤＳＰでしか実現できなかった。これに対して以上説明した実施の形態によれば、１個のフィルターで、ふたつのフィルターの演算を行うことができる。これにより、安価なＤＳＰでも、タップ数を減らすことなく、所望の周波数特性を確保しながら、ＦＩＲフィルターを用いた帯域分割が可能となる。
【００２７】
これを具体的な数値により説明する。図１に示す波長分割フィルターをｍ＝２０４８で実現する場合、１チャンネルあたりの乗算処理は、２０４８+３＝２０５１回で実現できる。６チャンネル出力だと、２０５１×６＝１２３０６回の乗算で済む。これは、図２を参照して説明した構成の場合のほぼ半分の演算量である。また、必要な遅延器およびフィルター係数も１組で済み、メモリの消費量も、図２を参照して説明した構成の半分となる。
【００２８】
［他の実施の形態］
乗算器５、６および加算器７による演算は、図１に示した実施の形態とは異なる順序で行うこともできる。たとえば、乗算器５ではタップ係数ａｍ／２−１を乗算することとし、乗算器６を加算器７の出力側に配置することもできる。すなわち、加算器３の出力と乗算器５の出力とを加算器７で加算した後に、加算器７の出力を逆符号とする。
【００２９】
また、ふたつのフィルターで加算器３を共有するのではなく、タップ係数ａ０〜ａｍ／２−１、ａｍ／２−１ａｍが乗算された信号の加算とは別に、それぞれの信号を逆符号にしてから加算することもできる。この場合、演算量は増えるが、逆符号にする処理は通常の乗算よりは処理負荷が小さいので、図２に示す波長分割フィルターよりは演算量は少なくて済む。
【００３０】
以上の説明では、可聴信号を高域専用再生スピーカーと低域専用再生スピーカーとに振り分ける場合を例に説明したが、本発明は、可聴信号以外の信号の帯域分割にも利用することができる。また、図１に示した波長分割フィルターを多段に接続し、帯域をさらに分割することもできる。
【符号の説明】
【００３１】
１−１〜１−ｍ、１１−１〜１１−ｍ、１２−１〜１２−ｍ  遅延器
２−０〜２−ｍ、２１−０〜２１−ｍ、２２−０〜２２−ｍ  乗算器
３、４  加算器（加算手段）
５  乗算器（回路手段の一部、減算値乗算器）
６  乗算器（回路手段の一部、逆符号乗算器）
７  加算器（回路手段の一部、第３の加算器）
３１、３２  加算器
【図１】
2010252010.tif 000010
【図２】
2010252010.tif 000011