(57)【要約】【課題】  入力信号の信号レベルに依存することなく演算処理のみでハウリングを検出し、ハウリングを聴覚的に確認できるようになる前に抑制する。【解決手段】  ハウリングを検出する対象となる入力信号が入力される入力端６と、入力端６から入力されたアナログの入力信号をデジタルへ変換する入力部２と、入力部２から出力された信号のうち有限長の区間をサンプリングする信号区分化部３と、サンプリングした信号の複素周波数を算出する複素周波数算出部４と、複素周波数の虚部の符号に基づいてハウリングの有無を検出するハウリング判定部５と、ハウリングの有無の結果を出力する出力端７と、でもってハウリング検出装置１’’を構成する。虚部の符号が正ならばハウリングとして検出する。検出したハウリングはハウリング抑制部９にて抑制される。【選択図】      図５
2006311356.tif 000002
【特許請求の範囲】
【請求項１】
アナログの入力信号をデジタルへ変換する入力部と、
デジタルに変換された入力信号の時系列から有限長の区間をサンプリングする信号区分化部と、
信号区分化部によってサンプリングされた時系列の区間ごとに、この区間の時系列信号を特徴づける１つ以上の複素周波数Ｆ＝ｆ＋ｉｂを算出する複素周波数算出部と、
算出した複素周波数Ｆ＝ｆ＋ｉｂの虚部ｂの符号に基づいてハウリングの発生の有無を検出するハウリング検出部と、
を備えることを特徴とするハウリング検出装置。
【請求項２】
前記信号区分化部において、
出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数１】
2006311356.tif 000018
がランダムな時系列になる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行することを特徴とする請求項１に記載のハウリング検出装置。
【請求項３】
前記信号区分化部において、
出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数２】
2006311356.tif 000019
の自己相関が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行することを特徴とする請求項１に記載のハウリング検出装置。
【請求項４】
前記信号区分化部において、
出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数３】
2006311356.tif 000020
が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行することを特徴とする請求項１記載のハウリング検出装置。
【請求項５】
前記信号区分化部が出力する時系列信号ｓ（ｔ）から線形予測法を用いてＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkを抽出することを特徴とする請求項１記載のハウリング検出装置。
【請求項６】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数が１つ以上存在する場合をハウリングの発生として検出することを特徴とする請求項１記載のハウリング検出装置。
【請求項７】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数においてハウリングが発生していると判定することを特徴とする請求項１記載のハウリング検出装置。
【請求項８】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数に対して抑圧特性を持つノッチフィルタを用いてハウリングを抑制するハウリング抑制部を備えることを特徴とする請求項１記載のハウリング検出装置。
【請求項９】
入力信号ｓ（ｔ）を区分的にＫ個（Ｋは１以上とする）の複素係数ａkと複素周波数Ｆk＝ｆk＋ｉｂkとを用いて
【数４】
2006311356.tif 000021
なる波形で近似し、各区分ごとに虚部ｂkの符号に基づいてハウリングの発生を検出することを特徴とするハウリング検出方法。
【請求項１０】
前記入力信号ｓ（ｔ）を最小自乗法を用いて
【数５】
2006311356.tif 000022
なる波形で近似することを特徴とする請求項９に記載のハウリング検出方法。
【請求項１１】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正であるものが１つ以上存在する場合にハウリングが発生していると判定することを特徴とする請求項９に記載のハウリング検出方法。
【請求項１２】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数においてハウリングが発生していると判定することを特徴とする請求項９に記載のハウリング検出方法。
【請求項１３】
前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数に対して抑圧特性を持つノッチフィルタを用いてハウリングを抑制するハウリング抑制部を有することを特徴とする請求項９に記載のハウリング検出方法。
【請求項１４】
前記複素周波数Ｆk＝ｆk＋ｉｂkを算出するのに用いた時系列の区間に対するハウリング抑制に対し、当該区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkに加え、これに先行する０個以上の区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkをも前記ノッチフィルタのパラメータとして用いてハウリングを抑制することを特徴とする請求項１３に記載のハウリング検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、アナログ信号を扱う電気回路におけるハウリングを検出する技術に関する。
【背景技術】
【０００２】
従来より音声や音楽などの物理的な振幅を電気的な信号に変換して採集し、採集した電気信号を増幅して再び物理的な振幅として出力する構成の電気回路が用いられている。こうした電気回路のうちアナログ信号を扱う電気回路において不要なハウリングを検出する技術が提案されている。
【０００３】
そのハウリング検出技術として、たとえば入力信号におけるハウリングの発生に伴い増大する信号レベルをハウリングによる音量の増大として聴覚的に検出するか、あるいは増大する信号レベルをハウリング成分として電気的に検出している。
【０００４】
また他のハウリング検出技術として、ハウリング成分を含まない正常な状態の周波数特性とハウリング状態にあってハウリング成分を含んだ周波数特性とを弁別するために、パワー周波数特性を周波数軸上において微分して求めた値をもってハウリング検出に用いる技術が開示されている（特許文献１参照）。
【特許文献１】特開平７−１４７７００号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の技術によれば、入力信号中のハウリングの信号レベルの増大に依存してハウリングを聴覚的に確認しハウリング検出するか、あるいはハウリングが発生した信号レベルでの電気信号のパワースペクトルやパワーエンベロープなどを解析してハウリング検出をする必要があり、このためハウリングが聴覚的に不快感を生じる強度にまでハウリングの信号レベルが増大してからでなければハウリング検出をすることができなかった。
【０００６】
本発明は、上記の課題に鑑みてなされたもので、その目的は、入力信号中のハウリングの信号レベルに依存することなくハウリングが聴覚的に確認できる信号レベルになる前にハウリングを検出することにある。
【課題を解決するための手段】
【０００７】
課題を解決するために、請求項１に記載の本発明は、アナログの入力信号をデジタルへ変換する入力部と、デジタルに変換された入力信号の時系列から有限長の区間をサンプリングする信号区分化部と、信号区分化部によってサンプリングされた時系列の区間ごとに、この区間の時系列信号を特徴づける１つ以上の複素周波数Ｆ＝ｆ＋ｉｂを算出する複素周波数算出部と、算出した複素周波数Ｆ＝ｆ＋ｉｂの虚部ｂの符号に基づいてハウリングの発生の有無を検出するハウリング検出部と、を備える。
【０００８】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、時系列からサンプリングした時系列信号について複素周波数算出部により複素周波数Ｆ＝ｆ＋iｂの算出処理のみを行い、算出結果の虚部ｂの符号に基づいてハウリングが聴覚的に確認できる信号レベルになる前にハウリングの発生を検出している。
【０００９】
請求項２に記載の本発明は、請求項１において、前記信号区分化部において、出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数１】
2006311356.tif 000003
【００１０】
がランダムな時系列になる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行する。
【００１１】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、複素周波数Ｆ＝ｆ＋ｉｂの算出において複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする請求項２に記載の演算子を演算して得られる時系列がランダムな時系列になる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を行い、虚部ｂの符号に基づいてハウリングが聴覚的に確認できる信号レベルになる前にハウリングの発生を検出する。
【００１２】
請求項３に記載の本発明は、請求項１において、前記信号区分化部において、出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数２】
2006311356.tif 000004
【００１３】
の自己相関が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行する。
【００１４】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、複素周波数Ｆ＝ｆ＋ｉｂの算出において複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする請求項３に記載の演算子を演算して得られる時系列の自己相関が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を行い、虚部ｂの符号に基づいてハウリングが聴覚的に確認できる信号レベルになる前にハウリングの発生を検出する。
【００１５】
請求項４に記載の本発明は、請求項１において、前記信号区分化部において、出力した時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数３】
2006311356.tif 000005
【００１６】
が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行する。
【００１７】
本発明においては、入力信号の信号レベルに依存することなく、複素周波数Ｆ＝ｆ＋ｉｂの算出において複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする請求項４に記載の演算子を演算して得られる時系列の二次ノルムが最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を行い、虚部ｂの符号に基づいてハウリングが聴覚的に確認できる信号レベルになる前にハウリングの発生を検出する。
【００１８】
請求項５に記載の本発明は、請求項１において、前記信号区分化部が出力する時系列信号ｓ（ｔ）から線形予測法を用いてＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkを抽出する。
【００１９】
本発明においては、Ｋ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkを抽出するのに線形予測法を用いている。
【００２０】
請求項６に記載の本発明は、請求項１において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数が１つ以上存在する場合をハウリングの発生として検出する。
【００２１】
本発明においては、ハウリングの発生の検出基準として、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数が１つ以上存在する場合で定めている。
【００２２】
請求項７に記載の本発明は、請求項１において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数においてハウリングが発生していると判定する。
【００２３】
本発明においては、ハウリングの発生の検出を前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数において行っている。
【００２４】
請求項８に記載の本発明は、請求項１において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数に対して抑圧特性を持つノッチフィルタを用いてハウリングを抑制するハウリング抑制部を備える。
【００２５】
本発明においては、ハウリングが聴覚的に確認できる信号レベルになる前にノッチフィルタを用いて発生したハウリングを抑制している。
【００２６】
請求項９に記載の本発明は、入力信号ｓ（ｔ）を区分的にＫ個（Ｋは１以上とする）の複素係数ａkと複素周波数Ｆk＝ｆk＋ｉｂkとを用いて
【数４】
2006311356.tif 000006
【００２７】
なる波形で近似し、各区分ごとに虚部ｂkの符号に基づいてハウリングの発生を検出する。
【００２８】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、入力信号の時系列からサンプリングした時系列信号について複素周波数算出部により複素周波数Ｆ＝ｆ＋iｂを算出し、虚部ｂの符号に基づいてハウリングが聴覚的に確認できる信号レベルになる前にハウリングの発生を検出する。
【００２９】
請求項１０に記載の本発明は、請求項９において、前記入力信号ｓ（ｔ）を最小自乗法を用いて
【数５】
2006311356.tif 000007
【００３０】
なる波形で近似する。
【００３１】
本発明においては、入力信号ｓ（ｔ）の近似に最小自乗法を用いている。
【００３２】
請求項１１に記載の本発明は、請求項９において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正であるものが１つ以上存在する場合にハウリングが発生していると判定する。
【００３３】
本発明においては、ハウリングが聴覚的に確認できる信号レベルになる前に検出するために、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正であるものが１つ以上存在する場合をもってハウリングの発生の検出を行っている。
【００３４】
請求項１２に記載の本発明は、請求項９において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数においてハウリングが発生していると判定する。
【００３５】
本発明においては、ハウリングが聴覚的に確認できる信号レベルになる前に検出するために、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数においてハウリングの発生の検出を行っている。
【００３６】
請求項１３に記載の本発明は、請求項９において、前記複素周波数Ｆk＝ｆk＋ｉｂkの虚部ｂkの符号が正である複素周波数の実部ｆkに対応する周波数に対して抑圧特性を持つノッチフィルタを用いてハウリングを抑制するハウリング抑制部を有する。
【００３７】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、ノッチフィルタを用いて発生したハウリングを聴覚的に確認できる信号レベルになる前に抑制している。
【００３８】
請求項１４に記載の本発明は、請求項１３において、前記複素周波数Ｆk＝ｆk＋ｉｂkを算出するのに用いた時系列の区間に対するハウリング抑制に対し、当該区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkに加え、これに先行する０個以上の区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkをも前記ノッチフィルタのパラメータとして用いてハウリングを抑制する。
【００３９】
本発明においては、入力信号中のハウリングの信号レベルに依存することなく、前記複素周波数Ｆk＝ｆk＋ｉｂkを算出するのに用いた時系列の区間に対するハウリング抑制に対し、当該区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkに加え、これに先行する０個以上の区間で算出された虚部ｂｋの符号が正な複素周波数の実部ｆkをノッチフィルタのパラメータとして用いて、発生したハウリングを聴覚的に確認できる信号レベルになる前に抑制している。
【発明の効果】
【００４０】
本発明によれば、入力信号中のハウリングの信号レベルに依存することなく演算処理のみでハウリングを検出することができる。
【００４１】
また本発明の別の効果は、ハウリングを聴覚的に確認できるようになる前に抑制することができる。
【発明を実施するための最良の形態】
【００４２】
＜第１の実施の形態＞
図１は第１の実施の形態のハウリング検出装置の構成図を示している。この図１にはハウリング検出装置１と、このハウリング検出装置１が備える、ハウリングを検出する対象となる入力信号が入力される入力端６と、入力端６から入力されたアナログの入力信号をデジタルへ変換する入力部２と、入力部２から出力された信号のうち有限長の区間をサンプリングする信号区分化部３と、サンプリングした信号の複素周波数を算出して虚部を求める複素周波数算出部４と、虚部の符号に基づいてハウリングの発生を検出するハウリング判定部５と、ハウリングの有無の結果を出力する出力端７と、が示されている。
【００４３】
このような構成において、入力部２はＡ／Ｄコンバータなどの信号変換回路を用いることができる。そして信号区分化部３と、複素周波数算出部４と、ハウリング判定部５の構成は、たとえばハードウェアによる論理回路の構成により実現することができる。
【００４４】
次に、こうした構成のハウリング検出装置１によるハウリングの検出のための信号処理の流れを図２のフロー図に示して説明する。
【００４５】
まず入力端６へ入力信号が入力される。この入力信号はハウリングの発生の発生を検出する対象となる信号であり、アナログ形式の信号である。入力された入力信号は入力部２においてＡ／Ｄ変換の信号処理がされてデジタル形式の信号となり時系列ｓ（ｔ）を得る（ステップ１：Ｓ１）。そして入力部２で実行されるＡ／Ｄ変換のサンプリング周波数は、例えば２０ｋＨｚまでのハウリング検出に対しては４０ｋＨｚ以上が選択され、５００ＭＨｚまでのハウリング検出に対しては１ＧＨｚ以上が選択される。
【００４６】
つぎに信号区分化部３で時系列ｓ（ｔ）のｎ番目の区分として、時刻ｔnから時刻ｔn+1の間の時系列ｓ（ｔ）を切り出す（ステップ２：Ｓ２）。信号区分化部３で実行する区分化には、例えばｔn+1−ｔnが一定になるようにとる等の方法を用いる。
【００４７】
つぎに複素周波数算出部４において信号区分化部３で切り出された時系列ｓ（ｔ）に対して演算子
【数６】
2006311356.tif 000008
【００４８】
がランダムな時系列になるような複素周波数Ｆk＝ｆk＋ｉｂkを算出する（ステップ３：Ｓ３）。
【００４９】
なお、複素周波数Ｆk＝ｆk＋ｉｂkを算出する方法としては、時系列
【数７】
2006311356.tif 000009
【００５０】
が最も小さくなるような複素周波数Ｆk＝ｆk＋ｉｂkを算出する方法や、時系列
【数８】
2006311356.tif 000010
【００５１】
の自己相関が最も小さくなるような複素周波数Ｆk＝ｆk＋ｉｂkを算出する等の方法を用いてもよい。
【００５２】
また、入力信号ｓ（ｔ）を区分的にＫ個（Ｋは１以上とする）の複素係数ａkと複素周波数Ｆk＝ｆk＋ｉｂkとを用いて
【数４】
2006311356.tif 000011
【００５３】
なる波形で近似する近似方法を用いて複素周波数Ｆk＝ｆk＋ｉｂkを算出してもよい。
【００５４】
上記いずれの算出方法を用いても、入力された入力信号に発生したハウリング信号レベルの増大に依存することなく、演算による信号処理のみで同様に抽出できる。
【００５５】
つぎにハウリング判定部５において複素周波数算出部４で算出された複素周波数Ｆk＝ｆk＋ｉｂkのうち、一つでも虚部ｂkの符号が正となるものが存在した場合は、これをもってハウリングを検出し、出力端７へ検出結果を出力する（ステップ４：Ｓ４）。
【００５６】
以上説明した第１の実施の形態によれば、入力信号中のハウリングの信号レベルに依存することなく演算処理のみでハウリングを検出することができる。
【００５７】
＜第２の実施の形態＞
図３に第２の実施の形態のハウリング検出装置の構成図を示す。第２の実施の形態は既に説明した第１の実施の形態に対して、ハウリング検出装置１’に線形予測係数算出部８とハウリング抑制部９を備える点で異なっている。第１の実施の形態では入力信号にハウリングが含まれているか否かを検出し、その検出結果を出力するのみであったが、第２の実施の形態においては発生したハウリングを抑制し、ハウリングを除去した信号を出力している。
【００５８】
図３に示した構成において実行されるハウリング検出とその抑制の信号処理を説明するフロー図を図４に示す。
【００５９】
まず入力端６に入力された入力信号は入力部２でＡ／Ｄ変換されてデジタル形式の信号となり時系列ｓ（ｔ）を得る（ステップ１０：Ｓ１０）。入力部２におけるＡ／Ｄ変換のサンプリング周波数は第１の実施の形態と同じく、たとえば２０ｋＨｚまでのハウリング検出に対しては４０ｋＨｚ以上が選択され、５０ＭＨｚまでのハウリング検出に対しては１ＧＨｚ以上が選択される。
【００６０】
つぎに信号区分化部２で時系列ｓ（ｔ）のｎ番目の区分として、時刻ｔnから時刻ｔn+1の間の時系列ｓ（ｔ）および同内容の時系列ｓ（ｔ）’を切り出す（ステップ１１：Ｓ１１）。
【００６１】
つぎに線形予測係数算出部８において信号区分化部３で切り出された時系列ｓ（ｔ）に対して線形予測法を用いてＫ個（Ｋは１以上とする）の予測係数ａkを算出する（ステップ１２：Ｓ１２）。
【００６２】
つぎにハウリング判定部５において線形予測係数算出部８で算出された予測係数ａkを係数とする方程式
【数９】
2006311356.tif 000012
【００６３】
をｘについて解いて得られるＫ個の解ｃkのうち、一つでも｜ｃk｜＞１なる解ｃkが存在するときにハウリングとして検出する（ステップ１３：Ｓ１３）。
【００６４】
つぎに複素周波数算出部４において｜ｃk｜＞１なる解ｃkの対数を取る（Ｆｋ＝−ｉｌｎｃｋ）ことによって、ハウリングに対応する複素周波数Ｆk＝ｆk＋ｉｂkを算出する（ステップ１４：Ｓ１４）。
【００６５】
つぎにハウリング抑制部９において、ハウリングに対応する複素周波数Ｆk＝ｆk＋ｉｂkの実部ｆkに対応する周波数に減衰を与える図示しないノッチフィルタを構成し、信号区分化部３で切り出された時系列ｓ（ｔ）'について算出してハウリングを抑制する（ステップ１５：Ｓ１５）。なお、ノッチフィルタはハードウェアにより構成された濾波回路により実現することができる。
【００６６】
つぎに出力部１０においてハウリング抑制部９でハウリング抑制された入力信号をＤ／Ａ変換して再び元のアナログ形式の信号へ戻した後に出力端７より出力する。
【００６７】
以上説明した第２の実施の形態によれば、入力信号中のハウリングの信号レベルに依存することなく演算処理のみでハウリングを検出することができる。
【００６８】
また別の効果として、ハウリングを聴覚的に確認できるようになる前に抑制することができる。
【００６９】
＜第３の実施の形態＞
図５に第３の実施の形態のハウリング検出装置の構成図を示す。第３の実施の形態は既に説明した第１の実施の形態に対して、ハウリング検出装置１’’にハウリング抑制部９を備える点で異なっている。第１の実施の形態では入力信号にハウリングが含まれているか否かを検出し、その検出結果を出力するのみであったが、第３の実施の形態においては第２の実施の形態と同様に発生したハウリングを抑制し、ハウリングを除去した信号を出力している。
【００７０】
図５に示した構成において実行されるハウリング検出とその抑制の信号処理を説明するフロー図を図６に示す。
【００７１】
まず入力端６に入力された入力信号は入力部２でＡ／Ｄ変換されてデジタル形式の信号となり時系列ｓ（ｔ）を得る（ステップ２０：Ｓ２０）。入力部２におけるＡ／Ｄ変換のサンプリング周波数は第１および第２の実施の形態と同じく、たとえば２０ｋＨｚまでのハウリング検出に対しては４０ｋＨｚ以上が選択され、５０ＭＨｚまでのハウリング検出に対しては１ＧＨｚ以上が選択される。
【００７２】
つぎに信号区分化部２で時系列ｓ（ｔ）のｎ番目の区分として、時刻ｔnから時刻ｔn+1の間の時系列ｓ（ｔ）および同内容の時系列ｓ（ｔ）’を切り出す（ステップ２１：Ｓ２１）。
【００７３】
つぎに複素周波数算出部４で、時系列ｓ（ｔ）の最小自乗法を用いて
【数１０】
2006311356.tif 000013
【００７４】
なる波形で近似することによって複素周波数Ｆk＝ｆk＋ｉｂkを算出する（ステップ２２：Ｓ２２）。
【００７５】
なお複素周波数Ｆk＝ｆk＋ｉｂkの算出は、例えば、近似の残差Ｅを与える式
【数１１】
2006311356.tif 000014
【００７６】
となるようにする等の方法をとればよい。
【００７７】
なお、複素周波数Ｆk＝ｆk＋ｉｂkの算出には時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数１】
2006311356.tif 000015
【００７８】
がランダムな時系列になる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行してもよい。
【００７９】
または時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数２】
2006311356.tif 000016
【００８０】
の自己相関が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkの算出を実行してもよい。
【００８１】
あるいは時系列信号ｓ（ｔ）にＫ個（Ｋは１以上とする）の複素周波数Ｆk＝ｆk＋ｉｂkをパラメータとする演算子
【数３】
2006311356.tif 000017
【００８２】
が最も小さくなる前記複素周波数Ｆk＝ｆk＋ｉｂkを求める算出を実行してもよい。これらのいずれの算出方法を用いても、入力された入力信号に発生したハウリング信号レベルの増大に依存することなく、演算による信号処理のみで同様に抽出できる。
【００８３】
つぎにハウリング判定部５において、複素周波数算出部４で算出されたＦk＝ｆk＋ｉｂkのうち一つでも虚部ｂkの符号が正となるものが存在するときにハウリングが存在しているものと検出する（ステップ２３：Ｓ２３）。
【００８４】
つぎにハウリング抑制部９において、信号区分化部３で切り出された本区間とこれに先行するＬ個（Ｌは０以上とする）の区間におけるハウリングに対応する複素周波数Ｆk＝ｆk＋ｉｂkの実部ｆkに対応する周波数に減衰を与える図示しないノッチフィルタを構成し、信号区分化部３で切り出された時系列ｓ（ｔ）’について算出してハウリングを抑制する（ステップ２４：Ｓ２４）。
【００８５】
つぎに出力部１０においてハウリング抑制部９でハウリング抑制された信号をＤ／Ａ変換して再び元のアナログ形式の入力信号へ戻し、出力端７から出力する。
【００８６】
以上説明した第３の実施の形態によれば、入力信号中のハウリングの信号レベルに依存することなく演算処理のみでハウリングを検出することができる。
【００８７】
また別の効果として、ハウリングを聴覚的に確認できるようになる前に抑制することができる。
【図面の簡単な説明】
【００８８】
【図１】第１の実施の形態のハウリング検出装置の構成図を示す。
【図２】第１の実施の形態におけるハウリング検出処理のフロー図を示す。
【図３】第２の実施の形態のハウリング検出装置の構成図を示す。
【図４】第２の実施の形態におけるハウリング抑制処理のフロー図を示す。
【図５】第３の実施の形態のハウリング検出装置の構成図を示す。
【図６】第３の実施の形態におけるハウリング抑制処理のフロー図を示す。
【符号の説明】
【００８９】
１、１’、１’’  ハウリング抑制装置
２  入力部
３  信号区分化部7
４  複素周波数算出部
５  ハウリング判定部
６  入力端
７  出力端
８  線形予測係数算出部
９  ハウリング抑制部
１０  出力部
【図１】
2006311356.tif 000023
【図２】
2006311356.tif 000024
【図３】
2006311356.tif 000025
【図４】
2006311356.tif 000026
【図５】
2006311356.tif 000027
【図６】
2006311356.tif 000028