(57)【特許請求の範囲】
【請求項１】
互いに独立な入射波が混信して得られた受信信号から分離信号を出力する独立成分分析処理部と、
前記分離信号の時間を時間シフト量パラメータだけずらすことにより、時間シフトした分離信号を得る時間シフト部と、
各分離信号と時間シフトした分離信号の時間シフト相関係数値を計算する相関係数計算部とを備え、
各分離信号の時間シフト相関係数値を用いて入射波数を推定する入射波数推定装置であって、
前記相関係数計算部は、雑音のみが存在する際に前記時間シフト量パラメータを変化させて雑音の時間シフト相関係数を計算し、計算した前記時間シフト相関係数があらかじめ設定した許容相関係数値よりも小さい値となる時間シフト量パラメータを設定する
ことを特徴とする入射波数推定装置。
【請求項２】
前記相関係数計算部により計算された時間シフト相関係数値が所定の閾値を超えるものの総数を入射波数と推定する閾値判定部をさらに備える
ことを特徴とする請求項１記載の入射波数推定装置。
【請求項３】
前記相関係数計算部は、入射波が存在せず、雑音のみが存在する場合の受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した全ての値より大きい値を前記所定の閾値に設定する
ことを特徴とする請求項２記載の入射波数推定装置。
【請求項４】
前記相関係数計算部は、入射波が存在する状況で、受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した値の最大値より小さい値を前記所定の閾値に設定する
ことを特徴とする請求項２記載の入射波数推定装置。
【請求項５】
前記相関係数計算部は、入射波が存在する状況で、受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した値の平均値を前記所定の閾値に設定する
ことを特徴とする請求項２記載の入射波数推定装置。
【請求項６】
前記相関係数計算部は、フィルタの特性が分からないときには、雑音のみが存在する際に前記時間シフト量パラメータを変化させて雑音の時間シフト相関係数を計算し、計算した前記時間シフト相関係数があらかじめ設定した許容相関係数値よりも小さい値となる時間シフト量パラメータを設定し、フィルタの特性が分かるときには、前記時間シフト量パラメータをフィルタ長に基づき設定する
ことを特徴とする請求項１又は２記載の入射波数推定装置。
【請求項７】
互いに独立な入射波が混信して得られた受信信号から分離信号を出力する独立成分分析処理ステップと、
前記分離信号の時間を時間シフト量パラメータだけずらすことにより、時間シフトした分離信号を得る時間シフトステップと、
各分離信号と時間シフトした分離信号の時間シフト相関係数値を計算する相関係数計算ステップとを含み、
各分離信号の時間シフト相関係数値を用いて入射波数を推定する入射波数推定方法であって、
前記相関係数計算ステップは、雑音のみが存在する際に前記時間シフト量パラメータを変化させて雑音の時間シフト相関係数を計算し、計算した前記時間シフト相関係数があらかじめ設定した許容相関係数値よりも小さい値となる時間シフト量パラメータを設定する
ことを特徴とする入射波数推定方法。
【請求項８】
前記相関係数計算ステップにより計算された時間シフト相関係数値が所定の閾値を超えるものの総数を入射波数と推定する閾値判定ステップをさらに含む
ことを特徴とする請求項７記載の入射波数推定方法。
【請求項９】
前記相関係数計算ステップは、入射波が存在せず、雑音のみが存在する場合の受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した全ての値より大きい値を前記所定の閾値に設定する
ことを特徴とする請求項８記載の入射波数推定方法。
【請求項１０】
前記相関係数計算ステップは、入射波が存在する状況で、受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した値の最大値より小さい値を前記所定の閾値に設定する
ことを特徴とする請求項８記載の入射波数推定方法。
【請求項１１】
前記相関係数計算ステップは、入射波が存在する状況で、受信信号を独立成分分析処理後、得られた各分離信号の時間シフト相関係数値を計算し、この計算した値の平均値を前記所定の閾値に設定する
ことを特徴とする請求項８記載の入射波数推定方法。
【請求項１２】
前記相関係数計算ステップにおいて、フィルタの特性が分からないときには、雑音のみが存在する際に前記時間シフト量パラメータを変化させて雑音の時間シフト相関係数を計算し、計算した前記時間シフト相関係数があらかじめ設定した許容相関係数値よりも小さい値となる時間シフト量パラメータを設定し、フィルタの特性が分かるときには、前記時間シフト量パラメータをフィルタ長に基づき設定する
ことを特徴とする請求項７又は８記載の入射波数推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、独立成分分析ＩＣＡ（Independent Component Analysis）による信号分離処理で、入射波数設定パラメータを事前に設定する必要がある問題に対し、時間シフト相関係数値に基づき入射波数を推定する入射波数推定装置及び入射波数推定方法に関するものである。
【背景技術】
【０００２】
独立成分分析ＩＣＡは、複数の入射波が混信して得られる受信信号から、各入射波が互いに統計的に独立であるという性質のみを用いて、受信信号から各入射波を分離する方式である。この方式は、複数話者の音声の混信分離、通信等での混信電波の分離等、種々の分野で適用されている。信号の統計的独立性のみを用いて分離を行うため、各入射波に関する事前情報が不要であるという利点がある。
【０００３】
独立成分分析ＩＣＡには様々な方式が提案されているが、代表的な方式が、非特許文献１に示されている。方式の詳細は非特許文献１に譲るが、ここで重要なことは、独立成分分析は、設定パラメータとして入射波数をあらかじめ設定しておく必要があることである。入射波数を誤ると、全入射波を正しく分離することができない。
【０００４】
このため、混信した入射波の入射波数を推定する場合、一般には、受信信号の共分散行列の固有値の情報から入射波数を推定する。その中でも、代表的な方式として、ＭＤＬ（Minimum Description Length）と呼ばれる方式（非特許文献２参照）がある。方式の詳細は非特許文献２に譲るが、この方式は、受信信号の共分散行列を固有値分解して得られた固有値を用いて入射波数を判定する。
【０００５】
ここで、後述の実施の形態の説明で必要な要点を中心に、独立成分分析ＩＣＡの概要を説明する。詳細は非特許文献１に譲る。なお、非特許文献１は複素信号を分離する複素ＩＣＡであるが、実数信号を分離するＩＣＡでも同じである。
【０００６】
図５は、独立成分分析処理の概要を示す図である。また、図６は、入射波数パラメータを真値と異なった値に設定したときの分離信号例を示す図である。
【０００７】
独立成分分析ＩＣＡは、上述したように、複数の入射波が混信して得られる受信信号から、入射波が互いに統計的に独立である性質のみを用いて各入射波を分離する方式である。ＩＣＡは、以下の混信モデルにおいて、互いに独立な入射波が混信して得られた受信信号から各入射波を混信分離するアルゴリズムである。
【０００８】
【数１】
0005328210.tif 000002
【０００９】
ここで、ｓ＝［ｓ１，…，ｓＪ］Ｔは各入射波ｓｊ（ｊ＝１，…，Ｊ）をベクトル状に格納した入射波ベクトル、Ａは混信を表す混合行列、ｎ＝［ｎ１，…，ｎＩ］Ｔは各受信アンテナ♯ｉ（ｉ＝１，…，Ｉ）の受信機で受信される受信機雑音をベクトル状に格納した受信機雑音ベクトル、ｘ＝［ｘ１，…，ｘＩ］Ｔは、各受信アンテナ♯ｉで受信した受信信号ｘｉをベクトル状に格納した受信信号ベクトルである。入射波数はＪ、受信信号数はＩとする。混合行列の各要素は、混信を意味する係数であり、信号が入射してくる方位や、信号の周波数、アンテナ１の配置、ゲイン等の諸条件から決まる値である。
【００１０】
例えば、受信アンテナ♯１の受信信号ｘ１は、ｘ１＝Ａ１１ｓ１＋，…，＋Ａ１ＪｓＪ＋ｎ１のような形となり、受信信号ｘ１はｓ１，…，ｓＪが混合行列の係数Ａｉｊが乗算されてから線形的に加算された形で混信しており、受信信号そのままでは、各到来波ｓｊを得られないことが分かる。
【００１１】
なお、これらのベクトルは、Ａ／Ｄ変換器３でサンプルされたタイミングで時系列（ｋ＝１，…，Ｋ）に得られる、すなわち、ｘ＝｛ｘ（１），…，ｘ（Ｋ）｝、ｓ＝｛ｓ（１），…，ｓ（Ｋ）｝、ｎ＝｛ｎ（１），…，ｎ（Ｋ）｝である。なお、図５は入射波数が３波の例（Ｊ＝３）を示している。
【００１２】
独立成分分析ＩＣＡは、あらかじめ入射波の数をパラメータとして設定しておく必要がある。以下、入射波数パラメータＪＰを正しく設定した（ＪＰ＝Ｊ）と仮定し、ＩＣＡの処理の概要を説明する。ＩＣＡでは、前処理としてＰＣＡ（Principal Component Analysis）により受信信号ｘを白色化処理する。白色化とは、信号間で互いに無相関で、同一電力であることを意味する。受信信号ｘの共分散行列ハットＲを次の式（２）で計算する。
【００１３】
【数２】
0005328210.tif 000003
【００１４】
ここで、Ｅ｛・｝は期待値、Ｈは複素共役転置を意味する。上記の行列の固有値分解で得られる固有値及び固有ベクトルを、固有値の大きさが降順になるように並べ、
【００１５】
【数３】
0005328210.tif 000004
【００１６】
とする。上記ＪＰ個の固有値とそれに対応する固有ベクトルから、白色化行列Ｍ∈ＣＪ×Ｉを、
【００１７】
【数４】
0005328210.tif 000005
【００１８】
と定義する。ここで、ｄｉａｇ｛１／√λ１，…，１／√λＪ｝は、１／√λ１，…，１／√λＪを対角要素にもつ正方行列である。Ｍを受信信号に乗算する。
【００１９】
【数５】
0005328210.tif 000006
【００２０】
この場合、
【００２１】
【数６】
0005328210.tif 000007
【００２２】
となる。ＩＪ×ＪはＪ×Ｊの単位行列である。単位行列の非対角項は零である。すなわち、これは、受信信号ベクトルｘは互いに無相関なＪ次のベクトルチルダｘに変換されたことを意味する。無相関とは、独立であるための必要条件であるから、実質、これでＩＣＡ処理の半分が終わったと言える。あとは、互いに無相関な信号チルダｘから、最も独立な信号を見つけ出せばよい。
【００２３】
次に、ＩＣＡでは、チルダｘから最も独立な信号を見つけるため、ネゲントロピーと呼ばれる指標値を最大化する。このネゲントロピーとは、ガウス性の変数に対しては０となるよう定義された尺度である。実際には、ネゲントロピーの計算には各入射波の確率密度関数が必要であるが、それらは一般には未知であるので、何らかの評価関数Ｇを導入し、各分離信号のネゲントロピー近似値の和を最大化する問題に帰着させる。具体的には、次の式の各分離信号のネゲントロピー値Ｅｊの和が最大となるような直交変換行列Ｗを見つける処理である。
【００２４】
【数７】
0005328210.tif 000008
【００２５】
ここで、ｙｇａｕｓｓは、分離信号ハットｓｊ＝ＷｊＨチルダｘと同じ分散をもつガウス性の信号、Ｇは何らかの滑らかな（微分可能な）関数、Ｅ｛・｝は期待値であり、期待値の計算にはサンプル平均値を用いる。関数Ｇ（ｙ）としては、例えば、
【００２６】
【数８】
0005328210.tif 000009
【００２７】
などがあり、ａは設定パラメータである。式（８）の最適化問題の解法は、ここでは重要でないので、省略する。
【００２８】
そして、最後に、式（８）の最適化問題で得られた直交変換行列Ｗの複素共役転置をチルダｘに乗算することにより、分離信号ハットｓを抽出する。
【００２９】
【数９】
0005328210.tif 000010
【００３０】
最終的に得られる混合行列推定値ハットＡの一般化逆行列は、
【００３１】
【数１０】
0005328210.tif 000011
【００３２】
である。このとき、
【００３３】
【数１１】
0005328210.tif 000012
【００３４】
である。仮に、混合行列が誤差なく推定できた、すなわち、ハットＡ＝Ａであるとする。このとき、ハットＡ†Ａは単位行列となるので、分離信号ハットｓは次の式となる。
【００３５】
【数１２】
0005328210.tif 000013
【００３６】
式（１４）では、雑音項Ａ†ｎは残るものの、ｓの項には如何なる混信係数も存在しないため、ｓが正しく得られることが分かる。以上が、入射波数パラメータＪＰを正しく設定（ＪＰ＝Ｊ）したときの処理である。
【００３７】
しかしながら、入射波数パラメータＪＰを真値Ｊより大きく設定（ＪＰ＞Ｊ）した場合、分離信号ハットｓは次の式で得えられる。
【００３８】
【数１３】
0005328210.tif 000014
【００３９】
ここで、Ｂはある行列である。式（１４）と比較すると、入射波成分ｓ＋Ａ†ｎに加え、雑音ｎの線形和信号Ｂｎが出力される。すなわち、分離信号のうちＪ波は入射波成分、残りのＪＰ−Ｊ波は雑音ｎの線形和信号Ｂｎが得られる。ガウス信号の線形和はガウス信号であるので、ｎがガウス雑音の場合、Ｂｎもガウス雑音である。なお、式（１５）は、ベクトル中の上Ｊ波を入射波成分、下ＪＰ−Ｊを雑音成分として表しているが、実際のＩＣＡでは、順番の不確定性の性質があるため、各分離信号の順番は未知である。
【００４０】
独立成分分析ＩＣＡに誤った入射波数パラメータを真値より大きく設定したときの動作例を説明する。図６（ａ）に入射波数を真値より小さく設定した時の分離信号例（正解３波を２波と設定⇒誤った出力）、図６（ｂ）に入射波数を正しく設定したときの分離信号例（全３波を正しく出力）、図６（ｃ）に入射波数を真値より大きく設定した時の分離信号例（正解３波を５波と設定⇒３波以外で雑音を出力）を示す。なお、ここでは、入射波として若干周波数の異なる３波の正弦波とし、観測雑音にはガウス性雑音を付加している。
【００４１】
図６より、（ｂ）では波数を正しく設定しているために全入射波を正しく分離できているのに対し、入射波数パラメータを真値より小さく設定すると、正しく信号を分離抽出できないことが分かる。一方、入射波数を真値より大きく設定すると、全信号を正しく分離抽出し、残りの信号成分として、雑音信号を出力することが分かる。
【００４２】
次に、入射波数パラメータが真値より少なく設定（ＪＰ＜Ｊ）した場合、式（４）で信号空間より狭い空間への射影となるため、全入射波を正しく分離できない。
【００４３】
【非特許文献１】Bingham E. and Hyvarinen A.,“A fast fixed-point algorithm for independent component analysis of complex valued signals,” International Journal of Neural Systems, Vol. 10, No.1, pp.1-8, Feb. 2000.
【非特許文献２】Mati Wax and Thomas Kailath,"Detection of Signals by Information Theoretic Criteria“，IEEE Trans. on Acoustics, Speech and Signal Processing, Vol. ASSP-33, No. 2, Apr. 1985.
【発明の開示】
【発明が解決しようとする課題】
【００４４】
独立成分分析ＩＣＡでは、事前に入射波数パラメータを設定する必要があるが、実際には未知の複数の混信波を分離する場合、入射波数を事前に知ることは困難である。更に、入射波数パラメータが真の値と異なると、正しく全入射波のみを分離抽出できないという問題点があった。また、ＭＤＬなどの受信信号の共分散行列の固有値に基づく方式は、実環境では、雑音の有色性などが原因で入射波数を誤るという問題点があった。
【００４５】
この発明は、上述のような課題を解決するためになされたもので、その目的は、固有値に基づかない新たな方法、すなわち各分離信号の時間シフト相関係数値に基づいて入射波数を推定することができる入射波数推定装置及び入射波数推定方法を得るものである。
【課題を解決するための手段】
【００４６】
この発明に係る入射波数推定装置は、互いに独立な入射波が混信して得られた受信信号から分離信号を出力する独立成分分析処理部と、前記分離信号の時間を時間シフト量パラメータだけずらすことにより、時間シフトした分離信号を得る時間シフト部と、各分離信号と時間シフトした分離信号の時間シフト相関係数値を計算する相関係数計算部とを設け、各分離信号の時間シフト相関係数値を用いて入射波数を推定する入射波数推定装置であって、前記相関係数計算部は、雑音のみが存在する際に前記時間シフト量パラメータを変化させて雑音の時間シフト相関係数を計算し、計算した前記時間シフト相関係数があらかじめ設定した許容相関係数値よりも小さい値となる時間シフト量パラメータを設定するものである。
【発明の効果】
【００４７】
この発明に係る入射波数推定装置は、固有値に基づかない新たな方法、すなわち各分離信号の時間シフト相関係数値に基づいて入射波数を推定することができるという効果を奏する。
【発明を実施するための最良の形態】
【００４８】
実施の形態１．
この発明の実施の形態１に係る入射波数推定装置について図１から図３までを参照しながら説明する。図１は、この発明の実施の形態１に係る入射波数推定装置を含む受信機の構成を示す図である。なお、各図中、同一符号は同一又は相当部分を示す。
【００４９】
図１において、受信機は、複数（Ｉ本）の受信アンテナ１と、加算器２と、Ａ／Ｄ変換器３と、入射波数推定装置４とが設けられている。
【００５０】
また、入射波数推定装置４は、独立成分分析（ＩＣＡ）処理部４１と、時間シフト部４２と、相関係数計算部４３と、閾値判定部４４とが設けられている。
【００５１】
つぎに、この実施の形態１に係る入射波数推定装置の動作について図面を参照しながら説明する。図２は、この発明の実施の形態１に係る入射波数推定装置の動作を示すフローチャートである。また、図３は、この発明の実施の形態１に係る入射波数推定装置の時間シフト相関係数値例を示す図である。
【００５２】
ここでは、各分離信号の時間シフト相関係数値に基づく入射波数推定装置を説明する。
【００５３】
本入射波数推定装置は、以下の考え方に基づいている。
１）一般に、雑音はランダム性が強いため、時間をずらして相関係数値を計算すると、理想的にはその値は零となる。
２）一般に、入射波は、あるキャリア周波数及び何らかの変調方式が加わった信号であるため、雑音よりランダム性が弱い。このため、時間をずらしてもその相関係数値は零とならない。
３）入射波数を真値に大きく設定すると、入射波の他に、雑音が出力される。
【００５４】
これらの性質を利用し、入射波数を大きく設定してＩＣＡ処理後、得られた各分離信号とそれらを時間シフトしたものの相関係数値を計算し、その値の大小から、入射波か雑音かを判定すれば、入射波と判定された総数により入射波数の推定を行うことができる。
【００５５】
まず、未知の入射波数の入射波ｓ＝［ｓ１ｓ２ｓ３］Ｔ（図１では一例とし３波）が、混信してＩ本の受信アンテナ１で受信される。
【００５６】
次に、ＩＣＡ処理部４１（ステップ１０１）では、大きめの入射波数パラメータ値ＪＰを設定し、混信受信信号から信号分離を行う。このとき、ハットｓ１，…，ハットｓＪＰの計ＪＰ波の分離信号が得られる。
【００５７】
次に、時間シフト部４２（ステップ１０２）では、分離信号の時間をｓｈだけずらすことにより、ハットｓ１（ｋ−ｓｈ），…，ハットｓＪＰ（ｋ−ｓｈ）を得る。
【００５８】
次に、相関係数計算部４３（ステップ１０３）では、各分離信号と時間シフトした信号の相関係数を計算する。時間シフト相関係数値Ｃｊは、次の式で計算する。
【００５９】
【数１４】
0005328210.tif 000015
【００６０】
ここで、ｓｈは時間シフト量パラメータである。なお、ＩＣＡでは分離信号ハットｓｊの分散値が１に規格化されているので、式（１６）の分母は１であり，実際の計算では分子のみを計算すれば良い。
【００６１】
最後に、閾値判定部４４（ステップ１０４）により、式（１６）で計算した値に閾値判定を行うことにより、入射波数を推定する。閾値を超えるものの総数が入射波数推定結果となる。
【００６２】
【数１５】
0005328210.tif 000016
【００６３】
ここで、Ｃｔｈは閾値であり、０≦Ｃｔｈ≦１の値をとる。
【００６４】
図３に時間シフト相関係数値例を示す。この図３では，正解数が３波の場合、設定数を６波としてＩＣＡを実行後、得られた各分離信号の時間シフト相関係数値を示している。分かりやすいように、大きい順にならべている。モンテカルロ１０試行分を示している。このように、３つの時間シフト相関係数値が閾値を超えており、閾値判定により３波と判定され、正解である。
【００６５】
実施の形態２．
この発明の実施の形態２に係る入射波数推定装置について図４を参照しながら説明する。図４は、この発明の実施の形態２に係る入射波数推定装置の閾値の設定方法を示す図である。なお、この発明の実施の形態２に係る入射波数推定装置の構成は、上記の実施の形態１と同様である。
【００６６】
閾値判定部４４の閾値の設定には、オペレータが、時間シフト相関係数値を見ながら手動で決定することも可能であるが、自動化するには、例えば以下の方法が考えられる。
１）入射波が存在せず、雑音のみが存在する場合の受信信号をＩＣＡ処理後、得られた各分離信号の時間シフト相関係数値を計算し、その値よりある程度大きい値に設定しておく。そして、実際に入射波が進入してきた場合にその閾値を用いる。こうすれば、閾値を低く設定しすぎて雑音を信号と誤判定することが無い。
２）入射波が存在する状況で、受信信号をＩＣＡ処理後、得られた各分離信号の時間シフト相関係数値を計算し、その最大値からある程度小さい値に設定する。こうすれば、閾値を高く設定しすぎて信号を雑音と誤判定することが無い。
３）入射波が存在する状況で、受信信号をＩＣＡ処理後、得られた各分離信号の時間シフト相関係数値を計算し、それらの平均値を閾値とする。こうすれば、１）及び２）のトレードオフ的な値を設定することができる。
【００６７】
図４（ａ）、（ｂ）、（ｃ）に、上記の１）〜３）の各方法の概念を示す。この図４では、入射波３波で、ＪＰを６波と設定した場合の例を示している。
【００６８】
実施の形態３．
この発明の実施の形態３に係る入射波数推定装置について説明する。なお、この発明の実施の形態３に係る入射波数推定装置の構成は、上記の実施の形態１と同様である。
【００６９】
すなわち、この発明の実施の形態３に係る入射波数推定装置４は、独立成分分析（ＩＣＡ）処理部４１と、時間シフト部４２と、相関係数計算部４３と、閾値判定部４４とが設けられている。相関係数計算部４３は、時間シフト量パラメータｓｈをフィルタ長Ｎに基づき設定する。
【００７０】
時間シフト相関係数値は、時間シフト量パラメータｓｈに依存する。雑音が完全にランダムな場合、ｓｈが１以上であれば、理論的には零となる。しかしながら、例えば、受信機雑音が、受信機のＡ／Ｄ変換器３とＩＣＡ処理部４１の間に挿入された帯域制限フィルタなどにより有色性を帯び、時間的相関性が強くなってしまっている場合には、ｓｈが小さいと、時間シフト相関係数値が大きくなってしまう。この場合、入射波の時間シフト相関係数値との差が小さくなるため、場合によっては入射波と雑音の判別が難しくなる。このような場合、ｓｈを大きめに設定する方法が考えられる。例えば帯域制限フィルタがＦＩＲである場合、フィルタ長をＮとし、フィルタ係数をｈ（０），…，ｈ（Ｎ−１）とすると、フィルタ後の雑音ｎｆは、
【００７１】
【数１６】
0005328210.tif 000017
【００７２】
時間シフト相関は、
【００７３】
【数１７】
0005328210.tif 000018
【００７４】
であるので、ｓｈをフィルタ長Ｎより大きい値に設定しておけば、式（１９）は零、すなわち、雑音の時間シフト相関係数は零となる。但し、無駄に大きくしてしまうと、相関係数を計算する際に必要なハットｓｊ（ｋ）とハットｓｊ（ｋ−ｓｈ）が重複する時間幅を十分にとるのに必要な信号サンプル数が増えるので、この点を考慮して設定すればよい。
【００７５】
また、フィルタの特性が分からなくても、雑音のみが存在する際にｓｈを変化させて雑音の時間シフト相関係数を計算し、値が十分低くなるようなｓｈを設定する。
【図面の簡単な説明】
【００７６】
【図１】この発明の実施の形態１に係る入射波数推定装置を含む受信機の構成を示す図である。
【図２】この発明の実施の形態１に係る入射波数推定装置の動作を示すフローチャートである。
【図３】この発明の実施の形態１に係る入射波数推定装置の時間シフト相関係数値例を示す図である。
【図４】この発明の実施の形態２に係る入射波数推定装置の閾値の設定方法を示す図である。
【図５】独立成分分析処理の概要を示す図である。
【図６】入射波数パラメータを真値と異なった値に設定したときの分離信号例を示す図である。
【符号の説明】
【００７７】
１  受信アンテナ、２  加算器、３  Ａ／Ｄ変換器、４  入射波数推定装置、４１  独立成分分析処理部、４２  時間シフト部、４３  相関係数計算部、４４  閾値判定部。
【図１】
0005328210.tif 000019
【図２】
0005328210.tif 000020
【図３】
0005328210.tif 000021
【図４】
0005328210.tif 000022
【図５】
0005328210.tif 000023
【図６】
0005328210.tif 000024