(57)【要約】デジタル変調信号が双方向チャネル（２１８）を含む信号経路において伝送されるデジタル通信システムにおいて、多段エコーキャンセラ（２０７）が、入力信号（２２０）から出力信号の線形エコー信号および非線形エコー信号を除去するためにシーケンスにおいて接続された２つ以上の段階を有している。結果は、出力信号周波数帯域と入力信号周波数帯域が重なり合うシステムと、帯域が分離されているシステム（ここでは、非線形エコー信号にのみかかわる）の双方にとって重要である。２つ以上の各段階はそれぞれの関数を特徴としており、シーケンスは２つ以上の段階の関数が線形関数と非線形関数の順に交互に並ぶことを特徴とする。
2006505211.tif 000002
【特許請求の範囲】
【請求項１】
デジタル伝送システムにおける、分散性媒体から受け取られたデジタル変調信号から、前記分散性媒体を介して伝送すべきデジタル変調信号のエコー信号を消去するための組み合わせであって、
デジタル変調信号源と、
デジタル変調信号を受け取るように結合された入力側を有するデジタル／アナログ変換器（ＤＡＣ）と、
前記ＤＡＣにより生成されるデジタル変調アナログ信号を受け取り分散性チャネルにおける伝送のために前記デジタル変調アナログ信号を提供するように結合された入力側と、前記分散性チャネルから受け取られたデジタル変調アナログ信号を提供するための出力側とを有するバイブリッドと、
前記分散性チャネルから受け取られた前記デジタル変調信号を受け取るように結合され、出力側とを有するアナログ／デジタル変換器（ＡＤＣ）と、
前記伝送デジタル変調信号源に結合された入力側と、出力側とを有する多段エコーキャンセラと、
ここで、前記多段エコーキャンセラは、前記伝送デジタル変調信号のエコー信号における線形ひずみをモデル化する第１結果を生成する第１関数を特徴とする第１段階と、前記伝送デジタル変調信号の前記エコー信号における非線形ひずみをモデル化する第２結果を前記第１結果から生成する第２関数を特徴とする少なくとも１つの第２段階とを有しており、
前記多段エコーキャンセラの前記出力側に接続された第１入力側と、前記ＡＤＣの前記出力側に接続された第２入力側と、出力側とを備える減算器とを具備する組み合わせ。
【請求項２】
前記多段エコーキャンセラは、
３つ以上の段階を具備しており、
前記各段階は入力側と出力側とを有しており、これら順次の段階が成すシーケンスにおいて別の段階に後続する任意の段階の前記入力側は前記別の段階の前記出力側に接続されており、
前記各段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記段階の前記それぞれの関数が線形関数と非線形関数の順に交互に並ぶものである、請求項１に記載の組み合わせ。
【請求項３】
前記シーケンスにおいて第３段階の後に１つ以上の段階が続いており、前記１つ以上の段階の前記それぞれの関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項２に記載の組み合わせ。
【請求項４】
前記線形関数は以下の関数を含んでおり、
【数１】
2006505211.tif 000004
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項２に記載の組み合わせ。
【請求項５】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項２に記載の組み合わせ。
【請求項６】
デジタル変調信号を受け取るように結合された入力側を有するデジタル／アナログ変換器（ＤＡＣ）と、
前記ＤＡＣにより生成されたデジタル変調アナログ信号をフィルタ処理するために前記ＤＡＣに結合されたフィルタと、
前記ＤＡＣにより生成され前記フィルタによりフィルタ処理された前記デジタル変調アナログ信号を、伝送のために分散性チャネルに提供し、前記分散性チャネルからデジタル変調アナログ信号を受け取るために、前記フィルタに結合されたハイブリッドと、
前記分散性チャネルから受け取られた前記デジタル変調アナログ信号を受け取るように結合され、変換されたデジタル変調信号を提供するための出力側を有するアナログ／デジタル変換器（ＡＤＣ）と、
デジタル変調信号源に結合された入力側と、出力側とを有する多段エコーキャンセラと、
ここで、前記多段エコーキャンセラは、前記デジタル変調信号のエコー信号における第１線形ひずみをモデル化する第１結果を生成する第１関数を特徴とする第１段階と、前記第１段階に直列に結合され、前記デジタル変調信号の前記エコー信号における非線形ひずみをモデル化する第２結果を前記第１結果から生成する第２関数を特徴とする第２段階と、前記デジタル変調信号の前記エコー信号における第２線形ひずみをモデル化する第３結果を前記第２結果から生成する第３関数を特徴とする少なくとも１つの第３段階とを有しており、
線形ひずみと非線形ひずみが補正された、前記伝送デジタル変調信号の前記モデル化されたエコー信号を前記変換された受信デジタル変調信号から減算するために、前記ＡＤＣの前記出力側と前記多段エコーキャンセラの前記出力側とに結合された減算器とを具備する、デジタル変調伝送システム用の送信／受信装置。
【請求項７】
前記多段エコーキャンセラは、
シーケンスを成すように順次に接続された３つ以上の段階を具備しており、
前記各段階は入力側と出力側とを有しており、前記シーケンスにおいて別の段階に後続する任意の段階の前記入力側は前記別の段階の前記出力側に接続されており、
前記各段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記段階のそれぞれの前記関数が線形関数と非線形関数の順に交互に並ぶものである、請求項６に記載の送信／受信装置。
【請求項８】
前記シーケンスにおいて第３段階の後に１つ以上の段階が続いており、前記１つ以上の段階の前記それぞれの関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項7に記載の送信／受信装置。
【請求項９】
前記線形関数が以下の関数を含んでおり、
【数２】
2006505211.tif 000005
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項６に記載の送信／受信装置。
【請求項１０】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項６に記載の送信／受信装置。
【請求項１１】
デジタル伝送システムにおける、分散性媒体から受け取られたデジタル変調信号から、前記分散性媒体を介して伝送すべきデジタル変調信号のエコー信号をモデル化および消去する方法であって、
第１デジタル変調信号を提供するステップと、
前記第１デジタル変調信号を第１デジタル変調アナログ信号に変換するステップと、
前記第１デジタル変調アナログ信号を分散性チャネルを介した伝送のためにハイブリットに提供するステップと、
前記第１デジタル変調アナログ信号のエコー信号を有する前記分散性チャネルから受け取られた第２デジタル変調アナログ信号をハイブリッドから提供するステップと、
前記第２デジタル変調アナログ信号を第２デジタル変調信号に変換するステップと、
前記第１デジタル変調信号のひずんだエコー信号をモデル化するステップと、
第１関数を用いて線形ひずみに関して前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルを補正し、次いで少なくとも、第２関数を用いて非線形ひずみに関して前記第１デジタル変調信号の前記エコー信号の前記モデルを補正することで、前記第１デジタル変調信号の補正されたエコー信号を生成するステップと、
前記補正されたエコー信号を前記第２デジタル変調信号から減算するステップとを具備する方法。
【請求項１２】
補正されたエコー信号を生成する前記ステップが、
シーケンスを成す関数により前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルを補正するステップを含んでおり、
前記各関数は入力信号と出力信号とを有しており、前記シーケンスにおいて別の関数に後続する任意の関数の前記入力信号は前記別の段階の前記出力信号を受け取り、
前記各関数は、出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する、入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．への応答を特徴とし、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記関数が線形関数と非線形関数の順に交互に並ぶものである、請求項１１に記載の方法。
【請求項１３】
前記線形関数が以下の関数を含んでおり、
【数３】
2006505211.tif 000006
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項１２に記載の方法。
【請求項１４】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項１２に記載の方法。
【請求項１５】
第１デジタル変調信号をアナログ形式に変換するステップと、
前記アナログ形式の前記第１デジタル変調信号をフィルタ処理するステップと、
前記フィルタ処理されたアナログ形式の前記第１デジタル変調信号を、伝送のために分散性チャネルの一端に結合するステップと、
前記分散性チャネルの前記一端から第２デジタル変調アナログ信号を受け取るステップと、
前記第１デジタル変調信号のひずんだエコー信号をモデル化するステップと、
第１関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける第１線形ひずみを補正する第１結果を生成することと、第２関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける非線形ひずみを補正する第２結果を前記第１結果から生成することと、少なくとも１つの第３関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける第２線形ひずみを補正する第３結果を前記第２結果から生成することとにより、補正されたエコー信号を生成するステップと、前記第２デジタル変調信号から前記補正されたエコー信号を減算するステップとを具備する、送信機／受信機のエコーキャンセル方法。
【請求項１６】
補正されたエコー信号を生成する前記ステップが、
関数のシーケンスにおいて前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルを補正するステップを含んでおり、
前記各関数は入力信号と出力信号とを有しており、前記シーケンスにおいて別の関数に後続する任意の関数の前記入力信号は前記別の関数の前記出力信号を受け取り、
前記各関数は、出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する、入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．への応答を特徴とし、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記関数が線形関数と非線形関数の順に交互に並ぶものである、請求項１５に記載の方法。
【請求項１７】
前記シーケンスにおいて第３関数の後に１つ以上の関数が続いており、前記１つ以上の関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項16に記載の方法。
【請求項１８】
前記線形関数が以下の関数を含んでおり、
【数４】
2006505211.tif 000007
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項１６に記載の方法。
【請求項１９】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項１６に記載の方法。
【請求項２０】
デジタル変調信号を受け取るように結合された入力側を有するデジタル／アナログ変換器（ＤＡＣ）と、
前記ＤＡＣにより生成されたデジタル変調アナログ信号をフィルタ処理するために前記ＤＡＣに結合されたフィルタと、
前記ＤＡＣにより生成され前記フィルタによりフィルタ処理された前記デジタル変調アナログ信号を、伝送のためにＤＳＬチャネルに提供し、前記ＤＳＬチャネルからデジタル変調アナログ信号を受け取るために、前記フィルタに結合されたハイブリッドと、
前記ＤＳＬチャネルから受け取られた前記デジタル変調アナログ信号を受け取るように結合され、変換されたデジタル変調信号を提供するための出力側を有するアナログ／デジタル変換器（ＡＤＣ）と、
デジタル変調信号源に結合された入力側と、出力側とを有する多段エコーキャンセラと、
ここで、前記多段エコーキャンセラは、前記デジタル変調信号のエコー信号における第１線形ひずみをモデル化する第１結果を生成する第１関数を特徴とする第１段階と、前記第１段階に直列に結合され、前記デジタル変調信号の前記エコー信号における非線形ひずみをモデル化する第２結果を前記第１結果から生成する第２関数を特徴とする第２段階と、前記デジタル変調信号の前記エコー信号における第２線形ひずみをモデル化する第３結果を前記第２結果から生成する第３関数を特徴とする少なくとも１つの第３段階とを有しており、
線形ひずみと非線形ひずみを補正された、前記伝送デジタル変調信号の前記モデル化されたエコー信号を前記変換された受信デジタル変調信号から減算するために、前記ＡＤＣの前記出力側と前記多段エコーキャンセラの前記出力側とに結合された減算器とを具備する、デジタル加入者回線（ＤＳＬ）システム用の送信／受信装置。
【請求項２１】
前記多段エコーキャンセラは、
シーケンスを成すように順次に接続された３つ以上の段階を具備しており、
前記各段階は入力側と出力側とを有しており、前記シーケンスにおいて別の段階に後続する任意の段階の前記入力側は前記別の段階の前記出力側に接続されており、
前記各段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記段階の前記それぞれの関数が線形関数と非線形関数の順に交互に並ぶものである、請求項２０に記載の送信／受信装置。
【請求項２２】
前記シーケンスにおいて第３段階の後に１つ以上の段階が続いており、前記１つ以上の段階の前記それぞれの関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項２１に記載の送信／受信装置。
【請求項２３】
前記線形関数が以下の関数を含んでおり、
【数５】
2006505211.tif 000008
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項２１に記載の送信／受信装置。
【請求項２４】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項２１に記載の送信／受信装置。
【請求項２５】
第１デジタル変調信号をアナログ形式に変換するステップと、
前記アナログ形式の前記第１デジタル変調信号をフィルタ処理するステップと、
前記フィルタ処理されたアナログ形式の前記第１デジタル変調信号を、伝送のためにＤＳＬチャネルの一端に結合するステップと、
前記ＤＳＬチャネルから第２デジタル変調アナログ信号を受け取るステップと、
前記第２デジタル変調アナログ信号を第２デジタル変調信号に変換するステップと、
前記第１デジタル変調信号のひずんだエコー信号をモデル化するステップと、
第１関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける第１線形ひずみを補正する第１結果を生成することと、第２関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける非線形ひずみを補正する第２結果を前記第１結果から生成することと、少なくとも１つの第３関数を用いて、前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルにおける第２線形ひずみを補正する第３結果を前記第２結果から生成することとにより、補正されたエコー信号を生成するステップと、
前記補正されたエコー信号を前記第２デジタル変調信号から減算するステップとを具備する、デジタル加入者回線（ＤＳＬ）システムにおける送信機／受信機のエコーキャンセル方法。
【請求項２６】
補正されたエコー信号を生成する前記ステップが、
シーケンスを成す順次の関数により前記第１デジタル変調信号の前記ひずんだエコー信号の前記モデルを補正するステップを含んでおり、
前記各関数は入力信号と出力信号とを有しており、前記シーケンスにおいて別の関数に後続する任意の関数の前記入力信号は前記別の段階の前記出力信号を受け取り、
前記各関数は、出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する、入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．への応答を特徴とし、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、
前記シーケンスにおいて前記関数が線形関数と非線形関数の順に交互に並ぶものである、請求項２５に記載の方法。
【請求項２７】
前記シーケンスにおいて第３関数の後に１つ以上の関数が続いており、前記１つ以上の関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項２６に記載の方法。
【請求項２８】
前記線形関数が以下の関数を含んでおり、
【数６】
2006505211.tif 000009
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項２６に記載の方法。
【請求項２９】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項２６に記載の方法。
【請求項３０】
デジタル伝送システムにおける、分散性媒体から受け取られたデジタル変調信号から、前記分散性媒体を介して伝送すべきデジタル変調信号のエコー信号を消去するための組み合わせであって、
デジタル変調信号源と、
デジタル変調信号を受け取るように結合された入力側を有するデジタル／アナログ変換器（ＤＡＣ）と、
前記ＤＡＣにより生成されるデジタル変調アナログ信号を受け取り分散性チャネルにおける伝送のために前記デジタル変調アナログ信号を提供するように結合された入力側と、前記分散性チャネルから受け取られたデジタル変調アナログ信号を提供するための出力側とを有するバイブリッドと、
前記分散性チャネルから受け取られた前記デジタル変調信号を受け取るように結合され、入力側を有するアナログ／デジタル変換器（ＡＤＣ）と、
前記伝送デジタル変調信号源に結合された入力側と、出力側とを有する多段エコーキャンセラと、
ここで、前記多段エコーキャンセラは、前記伝送デジタル変調信号のエコー信号におけるひずみをモデル化する第１結果を生成する第１関数を特徴とする第１段階と、前記伝送デジタル変調信号の前記エコー信号におけるひずみをモデル化する第２結果を前記第１結果から生成する第２関数を特徴とする少なくとも１つの第２段階とを含む複数の段階をシーケンスにおいて有しており、
前記シーケンスは少なくとも１つの線形段階と１つの非線形段階とを含んでおり、
前記多段エコーキャンセラの前記出力側に接続された第１入力側と、前記ＡＤＣの前記出力側に接続された第２入力側と、出力側とを備える減算器とを具備する組み合わせ。
【請求項３１】
前記多段エコーキャンセラは、
３つ以上の段階を具備しており、
前記各段階は入力側と出力側とを有しており、シーケンスを成す順次の前記段階において別の段階に後続する任意の段階の前記入力側は前記別の段階の前記出力側に接続されており、
前記各段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものである、請求項３０に記載の組み合わせ。
【請求項３２】
前記シーケンスにおいて第３段階の後に１つ以上の段階が続いており、前記１つ以上の段階の前記それぞれの関数が交互に並ぶ（非線形、線形、非線形、．．．）ものである、請求項３１に記載の組み合わせ。
【請求項３３】
前記線形関数が以下の関数を含んでおり、
【数７】
2006505211.tif 000010
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項３１に記載の組み合わせ。
【請求項３４】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項３１に記載の組み合わせ。
【請求項３５】
１つ以上の前記線形関数が周波数領域に実装される、請求項３１に記載の組み合わせ。
【請求項３６】
１つの前記段階の入力側と同様に接続された入力側と、出力側とを有する線形エコーキャンセラと、
前記線形エコーキャンセラの前記出力側に接続された第１入力側と、前記複数段階の１つの段階の前記出力側に接続された第２入力側とを有する、前記複数段階の別の段階への入力信号または前記多段エコーキャンセラの前記出力信号として組み合わされる出力信号を提供する加算器とをさらに含む、請求項３０に記載の組み合わせ。
【請求項３７】
前記複数段階が最後の段階をさらに含み、
前記線形エコーキャンセラの前記入力側は前記第１段階の前記入力側と同様に接続されており、
前記線形エコーキャンセラの前記出力側は前記加算器の前記第１入力側に接続されており、前記最後の段階の前記出力側は前記加算器の前記第２入力側に接続されている、請求項３６に記載の組み合わせ。
【請求項３８】
前記多段エコーキャンセラが、シーケンスを成す３以上の段階を含んでおり、前記線形エコーキャンセラが単一の段階を成しており、
前記各段階は入力側と出力側とを有しており、
前記各段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものである、請求項３６に記載の組み合わせ。
【請求項３９】
前記線形関数が以下の関数を含んでおり、
【数８】
2006505211.tif 000011
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項３８に記載の組み合わせ。
【請求項４０】
１つ以上の前記非線形関数が、１次元べき級数関数と、既知の非線形関数の逆関数と、既知の区分的線形関数の逆関数と、次数Ｐの１次元差分べき級数関数と、ｋパラメータに依存する一次元関数と、次数Ｐの二次元べき級数と、ｋパラメータに依存するＤ次元関数とを含むグループにおける関数である、請求項３８に記載の組み合わせ。
【請求項４１】
１つ以上の前記線形関数が周波数領域に実装される、請求項３８に記載の組み合わせ。
【請求項４２】
デジタル伝送システムにおける、分散性媒体から受け取られたデジタル変調信号から、前記分散性媒体を介して伝送すべきデジタル変調信号のエコー信号を消去し、前記伝送信号と前記受信信号との非線形結合の効果を除去し、前記受信信号を等化するための装置であって、
デジタル変調信号源と、
デジタル変調信号を受け取るように結合された入力側を有するデジタル／アナログ変換器（ＤＡＣ）と、
前記ＤＡＣにより生成されるデジタル変調アナログ信号を受け取り分散性チャネルにおける伝送のために前記デジタル変調アナログ信号を提供するように結合された入力側と、前記分散性チャネルから受け取られたデジタル変調アナログ信号を提供するための出力側とを有するバイブリッドと、
前記分散性チャネルから受け取られた前記デジタル変調信号を受け取るように結合され、出力側を有するアナログ／デジタル変換器（ＡＤＣ）と、
前記伝送デジタル変調信号源に結合された第１入力側と、前記ＡＤＣの出力側に結合された第２入力側と、出力側とを有する多段組み合わせとを具備しており、
前記多段組み合わせは、前記第１入力側に接続された入力側と出力側とを有するゼロ以上の段階を含む第１シーケンスと、前記第２入力側に接続された入力側と出力側とを有するゼロ以上の段階を含む第２シーケンスと、２つの入力側と１つの出力側とを有し、前記２つの入力側のうちの前記第１入力側は前記第１シーケンスの前記出力側に接続されており、前記第２入力側は前記第２シーケンスの前記出力側に接続されている係合段階と、出力側と入力側とを有するゼロ以上の段階を含んでおり、前記入力側は前記結合段階の前記出力側に接続されている第３シーケンスとを有しており、
前記第１シーケンスは前記伝送信号におけるひずみをモデル化する第１結果を生成するものであり、前記第２シーケンスは前記受信信号におけるひずみを補正する第２結果を生成するものであり、前記結合段階は伝送信号と受信信号との間の結合を補正しエコー信号を消去する第３結果を生成するものであり、前記第３シーケンスは前記受信信号におけるひずみを補正する第４結果を生成するものであり、
前記多段組み合わせは少なくとも１つの非線形段階を含むものである、装置。
【請求項４３】
前記第１シーケンス、前記第２シーケンス、および前記最終シーケンスそれぞれが単一の線形段階を含んでおり、前記各線形段階は入力側と出力側とを有しており、前記各線形段階は入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．に対応して出力信号デジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する各関数を特徴としており、前記出力信号時系列の任意の要素の値は前記入力信号系列の１つ以上の要素の値に依存するものであり、前記結合段階は、入力信号デジタル時系列ｘ1，ｘ2，ｘ3，．．．およびｙ1，ｙ2，ｙ3，．．．に対応して出力信号デジタル時系列ｚ1，ｚ2，ｚ3，Ｋを生成する関数を特徴としており、前記結合段階により生成される前記出力信号時系列の任意の要素の値は前記入力信号系列それぞれの１つ以上の要素の値に依存するものである、請求項４２に記載の装置。
【請求項４４】
前記線形関数が以下の関数を含んでおり、
【数９】
2006505211.tif 000012
ここで、Ａは定数であり、αkは前記関数のパラメータであり、前記パラメータは設定可能な値を有する、請求項４３に記載の装置。
【請求項４５】
前記結合段階の特徴は、各入力側からの１つの要素に依存する２次元べき級数関数と、ｋパラメータと前記２つの入力側からの１つの要素とに依存する２次元関数と、ｋパラメータと前記第１入力側からのＡ要素と前記第２入力側からのＢ要素、ここでＡ＋Ｂ＝Ｄであり、とに依存するＤ次元関数とを含むグループからの非線形関数である、請求項４３に記載の装置。
【請求項４６】
１つ以上の前記線形関数が周波数領域に実装される、請求項４３に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本出願は、２００２年１０月３０日出願の米国仮出願第６０／４２２，６５５号に対して優先権を主張するものである。本出願は、引用することにより本明細書の一部を成すものとする次の２つの係属中の米国特許出願の一部継続出願として優先権を主張するものでもある。
【０００２】
２００１年１０月１日出願第０９／９６８，４６９号の「デジタル変調信号における最大／平均出力信号比の低減（Ｐｅａｋ−ｔｏ−Ａｖｅｒａｇｅ  Ｐｏｗｅｒ  Ｒａｔｉｏ  Ｒｅｄｕｃｔｉｏｎ  ｉｎ  ａ  Ｄｉｇｉｔａｌｌｙ  Ｍｏｄｕｌａｔｅｄ  Ｓｉｇｎａｌ）」
２００１年１０月１日出願第０９／９６８，０６３号の「デジタル変調信号における線形ひずみおよび非線形ひずみを補正する多段等化器（Ａ  Ｍｕｌｔｉｓｔａｇｅ  Ｅｑｕａｌｉｚｅｒ  ｔｈａｔ  Ｃｏｒｒｅｃｔｓ  ｆｏｒ  Ｌｉｎｅａｒ  ａｎｄ  Ｎｏｎｌｉｎｅａｒ  Ｄｉｓｔｏｒｔｉｏｎ  ｉｎ  ａ  Ｄｉｇｉｔａｌｌｙ  Ｍｏｄｕｌａｔｅｄ  Ｓｉｇｎａｌ）」
本発明は、線形ひずみおよび非線形ひずみが存在する場合の双方向チャネルを介したデジタル変調手段による情報伝送に関する。より詳細には、本発明は、特にエコー信号が線形ひずみと非線形ひずみの両方の影響を受ける場合における、伝送信号の反響すなわちエコー効果に対する受信デジタル変調信号の補正に関する。
【背景技術】
【０００３】
デジタル変調とは、デジタル符号を用いて１つ以上の搬送波の１つ以上の特性を変化させて、このようにして形成された変形搬送波に情報を載せることを指す。このようにして、変調された搬送波が情報を「搬送」する。変調されていない搬送波は周波数がゼロである、すなわち、電圧など一定レベルである場合や、正弦波のように時間とともに変化する場合がある。デジタル変調は、搬送波の振幅、位相、および周波数のうちの１つ以上を変調する。デジタル変調の目的は、通信チャネルやデータ記憶チャネルなどで、１つ以上の変調信号を介して情報を伝送することである。
【０００４】
明確にするために、デジタル変調信号の「伝送」とは、チャネルと、信号がチャネルに出入りするために通過しなければならない、いずれかのチャネル端における任意の他の要素とを含む信号経路をデジタル変調信号が通過することを指すものとする。「チャネル」という用語は、信号の伝導または記憶に用いられる物理媒体を意味する。チャネルの例としては、ツイストペア線、同軸ケーブル、光ファイバ、宇宙空間の電磁波などが挙げられる。チャネルに加え、信号経路には、チャネルへデジタル変調信号を提供したりチャネルからデジタル変調信号を受け取るためにいずれかのチャネル端に結合された構成要素または要素が含まれる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
デジタル通信またはデータ記憶のために設計されたシステムにおけるデジタル変調信号の伝送では、信号が線形ひずみおよび非線形ひずみの影響を受けることが多い。このようなひずみにより信号が劣化するため、信号から情報を確実に抽出するためには、信号の受信時に補正措置を講じる必要がある。
【０００６】
線形ひずみは、伝送時の信号の状態を変化させる。これに関連して、信号が伝送されるチャネルは、信号成分の周波数に依存する不均等なレベルまで信号成分の振幅と位相を分散させる。その結果、受信信号が損なわれ、符号間干渉が生じることがある。このようなチャネルは、「分散性チャネル」と呼ばれる。入力信号または一部の入力信号成分の変化に直接比例して出力信号が変化するようなチャネルは、「線形チャネル」とみなされることがある。ただしこのようなチャネルでは、異なる周波数の成分がチャネルを異なる速度で通過し、異なる要因で減衰する。線形ひずみのこれらの効果は、受信信号を等化することで改善することができる。線形等化器は、受信信号の成分を調整して、チャネルを介した伝送が原因で生じた成分変化を補正することで、線形ひずみの効果を除去または低減する。線形効果は、ひずみのない信号の一部ではなかった周波数を新たに生成することはない。
【０００７】
非線形ひずみは、信号がひずむ比例性または線形性がある程度乱れた場合に発生する。このような非線形効果は通常、信号経路全体に広がることはなく、特定の場所に集中する。非線形ひずみの例として、（１）信号振幅または振幅の導関数（スルーレート）に依存するある程度の非線形性を示す、チャネルまたは中間チャネル中継器への入力側におけるドライバ、（２）ある程度の非線形（非オーム）特性を有するチャネルにおける腐食した接点、（３）磁心における磁気ヒステリシスに関連する可能性がある、大きな非線形性を示すチャネルにおけるトランスが挙げられる。一般に、非線形効果は、ひずみのない元の信号の一部ではなかった周波数を新たに生成する。
【０００８】
単一チャネルが、２つの信号経路の逆方向伝送機能を備えることがある。単一の共有チャネルを介した双方向伝送では、チャネルの両端で入力信号から出力信号を分離する手段をチャネル内に設ける必要がある。また、逆方向信号を分離し、チャネルの両端の中間で逆方向信号を再び組み合わせることができる中継器手段をチャネル内に設けなければならない場合もある。
【０００９】
この分離プロセスは完全なものではないため、一方のチャネル端において伝送される信号の一部が反響すなわち「エコー」する場合がある。このエコー信号が入力信号に重なって、あたかも雑音であるかのようにみなされることがある。一部の通信システムでは、入力信号と出力信号が、同じ周波数帯域の全部または一部を共有する。このような場合には、線形エコーキャンセラを用いてエコー信号レベルを低減することが多い。これは、エコー信号の予測値を生成するためにエコー信号に対する様々なシステム構成要素とチャネルの効果をモデル化することで実現される。これを入力信号から減算すると、エコー信号が大幅に低減され、信号対雑音比が向上する。他のシステムでは、２つの方向用に帯域が分離されている。これは、周波数分割二重化すなわちＦＤＤとして知られている方式である。１つの信号のエコー信号の線形部分は適切な帯域に残るが、これは非線形エコー信号には当てはまらない。その理由は、非線形エコー信号には、この帯域の範囲をはるかに超えて広がり反対の帯域に重なり合うことがある新しい周波数成分が含まれることが多いからである。したがって、本明細書の主題である多段非線形エコーキャンセラは、ＦＤＤ方式と非ＦＤＤ方式の両方に利益をもたらすことができる。
【００１０】
場合によっては、伝送すべき信号に低レベルの帯域外成分が含まれることがある。たとえば、限られたサンプリングレートを用いた場合に望ましくない効果が生じたり、隣接するデータシンボル間の推移においてデータにグリッチが生じることがある。これらの成分は、伝送前に信号を慎重にフィルタ処理することで除去することができる。また、後で、その効果を線形エコーキャンセルにより補正することもできる。伝送前に帯域外成分をすべて除去した場合、線形エコーキャンセラはＦＤＤ方式にとって意味をなさない。
【００１１】
エコーキャンセル問題の解決策の一部では、非線形エコーキャンセルと線形エコーキャンセルの両方が必要とされている。たとえば、米国特許第５，６２７，８８５号では、線形エコーキャンセラと非線形エコーキャンセラを個別に設ける方法が検討されている。ただし、これらのエコーキャンセラは直列ではなく並列に接続される。
【課題を解決するための手段】
【００１２】
本発明による多段非線形エコーキャンセラ（ＭＮＥＣ）は、エコー信号の非線形ひずみが線形信号経路における１つ以上の場所で避けられない場合に特に、線形エコーおよび非線形エコー効果を除去するための信号補正問題に対する効果的な解決策を提供するものである。調整可能パラメータを備えるＭＮＥＣは一連の様々な状況に適合できるため、非線形効果（１つまたは複数）の場所および特性を事前に知る必要はない。ＭＮＥＣには複数の段階が含まれる。各段階は、伝送信号のエコー信号におけるひずみをモデル化する関数を特徴とする。各段階は、各段階の関数が線形関数と非線形関数の順に交互に並ぶようなシーケンスにおいて接続される。このシーケンスには複数の段階が含まれる。大きな非線形効果が存在しない信号経路の１つ以上の区間におけるエコー信号に対する線形ひずみを除去するために、１つ以上の線形段階が設けられる。局所的な非線形効果によるひずみを補正するために、１つ以上の非線形段階が設けられる。非線形ひずみが局所的である限り、比較的少数の調整可能パラメータを有する関数を特徴とする非線形段階によるモデル化と補正の精度を高めることができる。
【００１３】
エコー信号経路のいずれかの場所に局所化された非線形ひずみの主要な発生源が１つである場合、多段等化器の好ましい実施形態は、線形、非線形、線形というシーケンスの３つの段階を有する。この場合ＭＮＥＣには、伝送される信号のコピーに対して直接作用して、エコー信号経路の始点と非線形ひずみ発生源の場所との間のエコー信号経路部分で発生する線形ひずみをモデル化する第１結果を生成する関数を特徴とする第１線形段階が含まれる。たとえば、この線形ひずみは、伝送前に帯域外成分を低減するために信号に加えられる一部のフィルタ処理において発生することがある。ＭＮＥＣにはさらに、第１段階に結合される第２非線形段階が含まれる。この第２段階は、第１結果に対して作用して、非線形ひずみ発生源により生成された非線形ひずみをモデル化する第２結果を生成する第２関数を特徴とする。たとえば、このような非線形ひずみは、デジタル／アナログ変換器、ラインドライバ、線路結合トランス、線路における接続の腐食、またはこれらの組み合わせにおいて発生することがある。第２段階に結合される第３線形段階は、第２結果に対して作用して、非線形ひずみ発生源の場所とエコー信号経路の終点との間の信号経路部分における線形ひずみをモデル化する第３結果を生成する第３関数を特徴とするたとえば、この線形ひずみは、ハイブリッド回路と線路との間のインピーダンスの不一致、線路における線番号の変更、線路に取り付けられたブリッジタップ、線路における不良接合、帯域外成分を低減するための受信機入力側のフィルタ、またはこれらの組み合わせにおいて発生することがある。第３結果がＭＮＥＣの出力信号となる。従来のエコーキャンセラでは、エコー信号の全部または一部を消去するために、この出力信号が受信信号から減算される。ＭＮＥＣのこの構成は、非線形ひずみが発生する場所がエコー信号経路に複数存在する場合にも容易に一般化することができる。このような場合には、多段等化器の段階が４つ以上になることがある。これらの場合には、第１段階と最後の段階が両方とも線形となることが多いが、非線形性が正確にエコー信号経路の始点と終点に特定される場合は例外である（この場合は、線形段階を省略してもよい）。
【００１４】
段階を特徴付ける関数は、調整可能パラメータを含んでもよい。特定のエラーレベルを最低限に抑える値に、これらのパラメータを設定するために、エコーキャンセラコントローラが線形段階および非線形段階に結合される。一般に、関数パラメータの調整では、エラーレベルを最低限に抑えるための最適化プロセスが活用される。
【発明を実施するための最良の形態】
【００１５】
本発明を上記の１つ以上の図面において示し、以下の説明において詳細に開示する。ここでは、「接続」される要素を図示し説明しているが、これは、これらの要素を順序立てて、これらの要素が互いに協調して働く仕方を基本的に説明するためのものである。したがって、ここに図示または説明していない他の要素を、図示および説明した要素間の接続に配置することは本発明の範囲内である。本発明は、双方向通信チャネルを介したデジタル通信に適用される。このチャネルは、複数媒体のいずれにおいて具体化してもよい。
【００１６】
従来のエコーキャンセラとは異なり、本発明による多段非線形エコーキャンセラ（ＭＮＥＣ）は、上りデータと下りデータが、分離された周波数帯域を占める場合、すなわち、周波数分割二重化すなわちＦＤＤとして知られるプロセスにおいて使用することができる。図１には、このような分離された上り帯域（１１０）と下り帯域（１２０）が示されている。非線形ひずみのプロセスを通じて、上り信号が、上り帯域よりも広く下り帯域と重なり合うことがある帯域（１３０）にある一連のひずみ係数付き信号要素を生成することがある。同様に、下り信号が、上り帯域（図示していない）と重なり合うことがある一連のひずみ係数付き信号要素が成す帯域を生成することがある。これらの一連のひずみ係数付き信号要素のエコー信号が、その帯域における真の信号と干渉し、最大データレートを下げる。データレートの低下率は、一連のひずみ係数付き信号要素が原因で高くなる実効「ノイズフロアー」が、線路に信号がない場合のバックグラウンドレベルをどの程度上回るかに関連している。ＭＮＥＣの機能は、ひずみとエコー効果の両プロセスをシミュレートし、干渉信号のモデルを作成し、このモデルを入力信号から減算して、これらの望ましくない成分を除去し、最大データレートを高めることである。上り帯域と下り帯域が重なり合う（すなわち、ＦＤＤでない）場合、ＭＮＥＣは標準的な線形エコーキャンセラに置き換わりエコー信号の一連の非線形ひずみ係数付き信号要素の除去を試みるという二重の機能を実行する。
【００１７】
図２はデジタル通信システムの一方の端を示すブロック図であるが、このシステムでは、伝送先に伝送すべき入力データ２０２が、符号化および変調回路２０４に提供される。（デジタル情報は、ハードウェアまたはソフトウェアのいずれにおいても処理できることに注意されたい。これは、信号がアナログ形式である参照番号２１０〜２２４が付された部品を除く図２の全部品に当てはまる。）回路２０４は、入力データ２０２をデジタル符号に対応付ける。この符号化されたデータは、一連のシンボルに分割される。各シンボルは、一定ビット数のデジタルデータを表す。次にこれらのシンボルを用いて、１つ以上の振幅、周波数、および位相で搬送波または搬送波セットが変調される。許容されるあらゆるシンボルについて、これらの搬送波パラメータに固有の設定が存在する。これらのパラメータは、次のシンボルを表すパラメータに切り替わる前の一定期間、固定されたままになる。デジタル変調信号２０６が回路２０４により生成される。これらの信号２０６は、伝送すべき変調された搬送波をデジタル形式で表す。
【００１８】
デジタル変調信号２０６が、デジタル／アナログ変換器（ＤＡＣ）２０８とＭＮＥＣ２０７に提供される。ＤＡＣ２０８は、デジタル変調信号をアナログ形式２１０に変換する。ＤＡＣ２０８から出力された信号２１０がアナログフィルタ２１２に結合される。このアナログフィルタの目的は信号２１０の帯域外成分のレベルを下げることであるが、副作用として、帯域内成分に位相ひずみと振幅ひずみをもたらす場合がある。このフィルタ処理された出力信号は、双方向通信チャネル２１８を介した出力信号と入力信号を分離する働きをするハイブリッド回路２１６に接続される。ハイブリッド回路２１６（明確には図示していない）には、出力側ラインドライバ、入力側前置増幅器、線路結合トランス、および設計によって異なることがあるその他の構成要素が含まれる。
【００１９】
ハイブリッド２１６が線路上で伝送信号と受信信号を完全に分離できれば理想的である。残念ながら、実際には完全に分離されない。通信チャネル２１８は、様々な作用により、出力信号の一部分をエコーすなわち反響し、あたかも出力信号のこの一部分が受信信号２２０の一部分であるかのように出力信号のこの一部分がハイブリッド２１６から出るようにすることがある。ツイストペア電話線の場合、これらの効果としては、ハイブリッド２１６と線路２１８との間のインピーダンスの不整合、線路における線番号の変更、線路に取り付けられたブリッジタップ、および線路における不良接合などが挙げられる。一般にエコー信号は、ＤＡＣ２０８、ハイブリッド２１６の一部、また線路２１８における不良接続や不良接合などで発生する可能性がある一部の非線形ひずみも被る。その結果、エコー信号に一連の非線形ひずみ係数付き信号要素が生じるという問題を解決しようとするのが、本発明によるＭＮＥＣである。受信信号２２０は、アナログフィルタ２２２に接続される。このアナログフィルタの目的は帯域外成分の強度を低減することであるが、副作用として、帯域内成分に位相ひずみと振幅ひずみをもたらす場合がある。アナログフィルタ２２２の出力側は、アナログ／デジタル変換器（ＡＤＣ）２２６に接続される。ＡＤＣを抜けた受信信号は、デジタル形式となる。この時点で、ＭＮＥＣ２０７の出力信号が、受信信号から２３０で減算される。本発明によるＭＮＥＣ２０７の機能は、入力信号に含まれる実際のエコー信号を消去して（あるいは部分的に消去して）信号対雑音比を向上するために入力信号から減算すべき一連のひずみ係数付き信号要素を含む予想エコー信号のモデルを生成することである。２３０で減算した後、信号は、受信信号に対するチャネルの効果を補正するために等化器２３４に接続される。最後に、信号が復調および復号化２３８されて、出力データ２４０が生成される。
【００２０】
本発明による多段非線形エコーキャンセラ（ＭＮＥＣ）２３７の構成、機能、および動作は、図３ａ〜図３ｄ、図４、図５ａ〜図５ｇ、および図６ａ〜図６ｂを参照すれば理解できるであろう。ＭＮＥＣは、少なくとも２つの段階を含むシーケンスから成る。各段階は入力信号としてデジタル時系列ｘ1，ｘ2，ｘ3，．．．を受け取り、出力信号として別のデジタル時系列ｙ1，ｙ2，ｙ3，．．．を生成する。出力信号系列の特定要素が入力信号系列の複数の要素の値に依存することがある（たとえば、ｙnがｘn-2、ｘn-1、ｘn、ｘn+1、およびｘn+2に依存することがある）。段階は、必要に応じて複数の設定可能パラメータに依存することがある各関数を特徴とする。設定可能パラメータの数は、各段階で同じである必要はない。関数は、関数を特徴とする段階の構成と働きを完全に記述するものであるため、以下では、段階に関する説明に代えて関数について説明する。段階を特徴付ける関数の設定可能パラメータにより、ＭＮＥＣは、未知または変わりやすい特性を有するチャネルを含む信号経路に適応することができる。各段階は次のような順序で接続される。すなわち、ＭＮＥＣへの入力側が第１段階への入力側となり、第１段階の出力側が第２段階の入力側に接続され、第２段階の出力側が第３段階の入力側に接続され、最後の段階までこれが繰り返されて、最後の段階の出力側がＭＮＥＣの出力側となるように接続される。段階は、「線形」と「非線形」の二種類に分けられる。いずれの種類であるかは、段階の出力信号時系列の要素が常に（および、すべてのパラメータ設定において）、その段階の入力信号時系列の要素に加えられる変更に対して線形的に（すなわち、直接比例して）変化するかどうかによって決まる。ＭＮＥＣは常に、少なくとも１つの線形段階と１つの非線形段階とを含むシーケンスに従う。線形段階の出力信号値は通常、チャネル分散性の補正に必要な多数の入力信号値に依存する。これに対し、非線形段階の出力信号要素は通常、１つまたは少数の入力信号要素にのみ依存する。一般に、段階のシーケンスにおいては、二種類の段階が交互に並ぶ。ただし、この規則は、特別な場合には順守されないことがある。
【００２１】
特定のＭＮＥＣ構成を選択する際には、以下の指針が適用される。エコー信号が生成されるのは、信号が、線形ひずみおよび非線形ひずみの様々な発生源にさらされやすい「エコー経路」を辿るときである。再び図２を参照すればわかるように、この経路はＤＡＣの入力側２０６からＡＤＣの出力側２２８までである。この経路には、ＤＡＣ２０８、アナログフィルタ２１２、ハイブリッド２１６、チャネル２１８、アナログフィルタ２２２、およびＡＤＣ２２６が含まれる。場合によっては、他の装置がチャネルに接続され、エコー信号と相互作用することもある（たとえば、デジタル通信と同じツイストペア線を共有する電話器）。短いチャネルの場合、エコー信号が、チャネル２１８の反対端（図示していない）に接続されたモデムまたはその他の装置と相互作用することもある。
【００２２】
エコー信号経路には、１つ以上の局所的な非線形ひずみ発生源が含まれると考えられる。各非線形効果は、エコー信号経路の特定箇所において、信号の振幅またはその一次導関数など、１次元（または、少なくとも低次元）の力学変数に従って作用すると考えられる。多段アプローチの基礎となる設計手法は、エコー信号経路をこれらの局所的な非線形効果により分離される一連の線形区間とみなすことができるという考えを前提としている。これらの効果の詳細とその場所を事前に知る必要はない。これらの効果の１つは、エコー信号経路の始点付近で生じる。エコーキャンセルプロセスの最初のステップは、線形段階を用いて、伝送される（未知の）信号から該当する力学変数（１つまたは複数）を再現することである。次に、線形段階に続く非線形段階により、この力学変数（１つまたは複数）から非線形ひずみを除去する。多くの場合、この非線形段階は、わずかな数の項のみ有する単一変数のべき級数展開を特徴としてもよい。この非線形段階の出力信号は、非線形効果を受けた直後のエコー信号を表すのに用いられる。エコー信号経路に顕著な非線形効果が１つのみ存在する場合、この出力信号を別の線形段階に渡してもよい。ＭＮＥＣはこの線形段階を用いて、エコー信号経路の終点までのエコー信号の追加線形ひずみを補正する。以上が、図３ａに示した線形、非線形、線形というシーケンスの３段階構成の説明である。別の方法では、非線形性がすべて補正されるまで、非線形段階と線形段階が順に交互に並ぶようなシーケンスを追加することができる。図３ｂは、２つの異なる場所で発生する非線形効果を含むエコー信号経路をモデル化するのに用いることができる５段階のシーケンスを示している。通常、ＭＮＥＣの最初と最後は線形段階となることに注意されたい。ただし、非線形ひずみ場所の１つが一方のエコー信号経路端にあるかまたはその近くにある場合は、この場所とエコー信号経路端との間に重大な線形ひずみ効果が残らないように、例外を設けてもよい。
【００２３】
システム内の複数の場所でエコー信号が発生することがあるため、その作用を適切にモデル化するために複数の「エコー信号経路」を考慮しなければならない場合がある。これを実現するには、複数のエコーキャンセラを並列に配置すればよい。この場合、すべてを順次ＭＮＥＣタイプにしてもよいし、一部を線形エコーキャンセラにしてもよい。図３ｃに示した典型例では、３段階ＭＮＥＣが１段階線形エコーキャンセラと並列に配置されている。図からわかるように、両方のエコーキャンセラは同じ入力データを受け取るが、両方の出力信号は加算される。追加するエコーキャンセラは、ＭＮＥＣの一部の段階とのみ並行に配置することができることに注意されたい。たとえば、図３ｂに示したＭＮＥＣの最後の３段階と並行に線形段階を配置することができる。さらに複雑な構成も可能である。図３ｄは、その中でも特に有用な構成を示している。
【００２４】
図３ｄに示した構成は、多段等化器の働きと多段非線形エコーキャンセラの働きを組み合わせた多段装置である。この装置では、２つの線形段階が、チャネルに入力される伝送データを受け取る入力側と、チャネルから接続される受信データを受け取る入力側とをそれぞれ有しており、また、２つの入力側と少なくとも１つの出力側とを有する非線形段階に提供される出力側をそれぞれ有している。このような装置は、入力信号（受信データ）と出力信号（伝送データ）のエコー効果の双方が同じ非線形ひずみの発生源と相互作用する場合に有用であることがある。この典型例として、ツイストペア線における不完全な接合の非線形特性が挙げられる。非線形段階の２つの入力側にそれぞれの信号を提供する２つの線形段階の出力信号は必ずしも加算されるとは限らず、むしろ、非線形段階を特徴付ける一部の非線形関数への個別の入力側の役割を果たすことができる。これにより、出力信号と入力信号との結合が可能になる。このような結合が重要とみなされない場合は、必要に応じて、非線形多段等化器とＭＮＥＣの両方を結合しないで使用することができる。
【００２５】
図４は、線形段階の好ましい実施形態を示している。ここでは、この線形段階が有限インパルス応答（ＦＩＲ）フィルタとして具体化されている。当業者であれば、このような要素を、その要素が実装される関数で特徴付けるまたは記述することができることを理解できるであろう。この場合、この実施形態は次の関数４１０により実装される。
【００２６】
【数１０】
2006505211.tif 000003
図４に示した関数４１０を図３ａに示したＭＮＥＣの第１段階として用いる場合、入力信号はデジタル値ｕ0，ｕ1，ｕ2，ｕ3．．．の時間系列となる。この関数では、連続する各デジタル値が、デジタル値と組み合わされて（乗算されて）積を算出する値を有するパラメータ（この場合は係数αk）に関連付けられる。指数ｋの範囲には通常、選択された開始値ｋ0と終了値ｋ1との間のすべての整数値が含まれる。これらの値は、正数、負数、またはゼロのいずれであってもよい。どの値を用いるかは、特定チャネルの分散性およびその他の特性によって決まる。必要であれば、図に示したように、信号レベルの変化を補正するために定数パラメータＡを含めてもよい。すべての積が合計され、Ａに加算されて、デジタル値ｘ0，ｘ1，ｘ2，ｘ3．．．の出力信号時間系列の単一要素が生成される。各パラメータとデジタル値の関連付けを時間系列において１つ前に進めることで、出力信号時間系列の次の要素が生成される。この時間系列が、図３ａに示したＭＮＥＣの第１段階により生成される第１結果を構成する。以下に詳細に説明する方法で、エコーキャンセラコントローラにより係数の値が設定・変更される。
【００２７】
図３ａのＭＮＥＣの第２段階（ここでは、非線形段階と想定する）は入力信号時間系列として第１結果を受け取り、出力信号時間系列ｙ0，ｙ1，ｙ2，ｙ3．．．として第２結果を生成する。図５ａ〜図５ｇを参照すれば、複数の非線形段階の実施形態を理解できるであろう。図５ａでは、基本的なべき級数関数５１４により非線形段階が実装されている。これはほとんどの場合に好ましい実施形態であり、特に振幅ベースの非線形性についてはそうである。これは、１次元（１−Ｄ）段階である。ここで１次元というのは、出力信号ｙnが単一の出力信号ｘnにのみ依存するという意味である。次元が低いため、この段階は、数十または数百という次元または数の「タップ」が実際に存在する場合が多い線形段階（図４）に比べ単純である。図５ｂに示した関数５１６は、チャネルにおいてまたはその近辺において非線形関数の正確な式ｆ（ｘ）が既知であり、この場合に非線形段階をも特徴付けるｆ（ｘ）により具体化される場合に、非線形段階を特徴付けることができる。図５ｃは、非線形段階を特徴付けることができる差分ベースの１−Ｄ非線形関数５１８を示している。.この非線形段階の実施形態は、非線形性にかかわるスルーレートを補正するために関数５１４に代わるものとして用いてもよい。パラメータの数に応じて、この非線形段階を図５ｄに示した一般的な１−Ｄ関数５２０に従って実装してもよい。このことが当てはまるのは、特定の用途にとって好都合である他の種類の級数展開や、非線形関数の形状に直接近似させるスプラインまたはその他の方法である。図５ｅに示した非線形段階の実施形態においては、この段階を特徴付ける関数５２２を、振幅およびスルーレートに同時に依存するような１−Ｄに下げることができない非線形性の２−Ｄべき級数として実装してもよい。これは、関数が２つの連続する入力信号値に依存するような実施形態となることが多い。図５ｆの非線形段階の実施形態は、Ｄ次元の一般的関数５２４を特徴とする。Ｄが増えるほど、この実施形態は望ましくなくなる。図５ｇは、図３ｄに示したような二重入力側非線形段階の一実施形態を示している。この段階を特徴付ける非線形関数を、２つの各入力信号シーケンスからの１つの要素を利用するべき級数として実装してもよい。ｘとｙの両方の様々な累乗を含む級数に交差項が存在する場合に、上り／下りの結合を補正することができる。これを、第１入力信号からのＡ要素と第２入力信号からのＢ要素に依存する任意のＤ次元非線形関数に一般化することができる（ここでＡ＋Ｂ＝Ｄである）ことは明白である。ＡとＢは少なくとも１でなければならないため、Ｄは少なくとも２でなければならない。
【００２８】
図３ｂに示したように、必要に応じて段階を追加してもよい。この場合、一般に線形段階と非線形段階が交互に並ぶ。
【００２９】
ＭＮＥＣを線形段階と非線形段階が交互に並ぶような構成にすることで、エコー信号のモデル化に必要な調整可能パラメータの数を減らすことができる。たとえば、２つの各線形段階が５０個の入力変数（ｋ0からｋ1まで５０個の値）を持つ関数を特徴とし、非線形段階が５乗までの１つの変数のべき級数を用いる関数を特徴とするような、３つの段階から成るＭＮＥＣ構成を検討してみよう。任意の定数項が各線形段階に含まれる場合、この構成では、合計で１０６個の調整可能パラメータが必要となる。これらの段階を組み合わせてＭＮＥＣを形成する場合、各出力信号値は９９個の入力信号値に依存する。このため、単一段階の場合と同じ結果を得るには、非線形エコーキャンセラに９９個の変数を処理する非線形関数が必要となる。適切な非線形関数は、９９個の変数で可能なすべての５乗の形を使って作成することができるが、この場合、９１，９６２，５２０個の調整可能パラメータが必要である。この数は、上記の例で必要とされる調整可能パラメータの個数１０６を大幅に上回るため、望ましくない。
【００３０】
図６ａに示した高速フーリエ変換（ＦＦＴ）６１０などのフーリエ変換を線形段階の前に配置することで、線形段階を周波数領域に実装してもよい。この場合、線形段階は図６ｂに示した形式の関数６２０を特徴とする。この関数は加算がないため、図４の関数４１０よりも単純である。あらゆる線形段階をこの方法で実装することができる。その場合、各線形段階の前にＦＦＴが配置され、各線形段階の後に逆ＦＦＴが配置される。ＭＮＥＣに周波数領域における入力信号および／または出力信号を持たせたい場合は、最初のＦＦＴおよび／または最後の逆ＦＦＴを取り除くことができる。これは、復調プロセスでフーリエ変換が用いられる場合に有用である。このことは、ＡＤＳＬ（非対称型デジタル加入者回線）通信方式でよく利用されているＤＭＴ（離散マルチトーン）変調に当てはまる。周波数領域において必要とされる係数の数（ＤＭＴ信号の各チャネルで２つ）が、最適化するには望ましくないほど多くなることがある。このような場合、これらの一連のパラメータが１つ以上の滑らかな曲線上にあると想定することで、係数の数を減らしてもよい。これらの曲線は、最適化に用いられるまばらなマーカーポイントを特徴とすることができる。完全なパラメータセットは、三次スプラインまたはその他の補間法を用いてこれらのマーカーポイントから取得することができる。
【００３１】
図７ａは、本発明のＭＮＥＣを用いて達成できることを示す例である。図２のシステムは、ハードウェアとソフトウェアにそれぞれ部分的に実装された。既知のデータで変調された標準ＤＭＴ（離散マルチトーン）信号が１２，０００フィートの実際のツイストペア線を通じて伝送され、そのエコー信号が、同じ（伝送）端で処理するために受信された。非線形ひずみの発生源となり得るのは、ラインドライバ（ハイブリットの一部）と線路結合トランスである（ただし、これらに限定されない）。このグラフは、約４ｋＨｚ間隔のＤＭＴチャネルの番号に対応する残留雑音レベルを均等目盛で示したものである。上部の曲線７１０は、線形エコーキャンセルで達成できた最高の結果を示している。中央の曲線７２０は、線形エコーキャンセラ結果７１０と同じくらい複雑な構成を用いた、本発明のＭＮＥＣにより向上した結果を示している。雑音が１／２未満に（ほぼ８ｄＢ）低減されている。一番下の曲線７３０は、信号が伝送されていないときのバックグラウンドレベルを示している。チャネル１４付近で上向き突出部がみられるのは、バックグラウンドの一部である未確認の雑音源が原因である。
【００３２】
図７ｂは、分離された上り帯域と下り帯域を用いた際に、非線形効果のせいで大きな帯域外エコー信号が存在する場合の実験結果を示している。伝送信号（上り帯域）は、チャネル６からチャネル３１まで（図示していない）に限定されている。図に示したチャネルはすべて、下り帯域に含まれる。非線形性がなければ、これらのチャネルに上りエコー信号は存在しないであろう。しかし、線形エコーキャンセルの影響を全く受けない曲線７４０を見れば、大きなエコー信号が存在することがわかる。このエコー信号を補正しなければ、これらの下りチャネルに含まれる入力データの回復において相応のエラーが生じる。本発明によるＭＮＥＣを適用すると、かなりの改善を示す曲線７５０のような結果が得られる。これにより、線路が対応できるデータレートが高くなる。下部の曲線７６０は、伝送信号がないときのバックグラウンドレベルを示している。この特定の例における非線形エコー信号の発生源は、線路における接点不良であることがわかっている。
【００３３】
図８では、パラメータを最適化するために、ＭＮＥＣまたはそのコピーが異なる構成で接続されている。テストデータセットが生成され、線路上で伝送される。この信号のエコー信号が集められて記憶される。エコーキャンセラコントローラ８２４を用いて、最適なパラメータが特定される。エコーキャンセラコントローラが、一連の試験的パラメータをＭＮＥＣ８０６の各段階（８０８、８１０、および８１２）に提供する。次に、テストデータ８０２がＭＮＥＣにより処理され、捕捉されたエコーデータ８１６から８１４で減算される。その結果がコントローラ８２４により処理されて、パラメータ値精緻化プロセスで用いられる。
【００３４】
図９は、本発明によるＭＮＥＣを動作させるための一連の初期構成ステップを示している。最適化技術を利用しているため、図９のプロセスを最適化プロセスと呼ぶこともできる。この図は、誤差測定値に応じてエコーキャンセラコントローラが実行する処理を具体的に示した流れ図である。この流れ図における説明では、具体的には図４および図５ａ〜図５ｆに示した関数を特徴とするＭＮＥＣの各段階の構成が標準的であることを前提としている。このような構成は、汎用または専用のデジタルプロセッサ、プログラマブルロジック、特定用途向け集積回路（ＡＳＩＣ）、デジタル信号プロセッサ（ＤＳＰ）、およびこれらに相当するあらゆる手段により実行されるソフトウェアルーチンを含む（ただし、これに限定されない）、任意の様々な形式で実装してもよい。通常、これらの関数は、順次接続された記憶セル配列に実装される。これらのセルに記憶された値を取得して係数と組み合わせる（一般には乗算）ために、タップからこれらのセルにアクセスすることができる。係数の値は、これらの値の設定および変更を目的としてアクセスできる場所に記憶される。
【００３５】
図９では、本発明によるＭＮＥＣがステップ９４０で初期化される。初期化とは、各関数の全パラメータの値を、予め決められた初期値に設定することを意味する。設定の一選択肢として、入力信号がほとんどまたは全く変わらないまま出力信号に転換できるようにする値に、パラメータを初期設定することができる。たとえば、図４の線形段階を初期化する際に、パラメータαkを１に設定し、定数Ａを含めその他のすべてのパラメータを０に設定することができる。また、図５ａの非線形段階を初期化する際に、すべてのパラメータを０に設定することができる。別の選択肢として、以前に試みた最適化で得られたパラメータの最良値を初期値として用いることができる。次に、ステップ９４２チャネルを介して伝送され受信・記憶されたデータの擬似ランダム系列を通常は含んでいる既知の信号がＭＮＥＣにより処理される。９４４で得られる誤差測定値、通常は減算（図８の８２０）により得られる最小自乗誤差値と特定の最適化ルーチン（たとえば、修正Ｐｏｗｅｌｌルーチン）を用いて、関数パラメータの新しい値セットがステップ９４６で選択される。９４８で、ルーチンがステップ９４２、ステップ９４４、ステップ９４６を繰り返しループする。ループのたびに、記憶されたデータセットが各段階全体で新しいパラメータセットで再処理され、これが誤差測定値が最小値に収束するまで続けられる。次にパラメータが固定され、実際の（未知の）データを得るためにエコー信号をモデル化するのに用いることができるようになる。９４４で誤差値を定義する方法は複数考えられる。たとえば、減算（図８の８２０）後の残留信号をフーリエ変換で処理すれば、対象とする特定の周波数チャネルに残っている誤差の大きさを調べることができる。
【００３６】
図３ｃおよび図３ｄに示したような変更が加えられたＭＮＥＣ構成では、パラメータを最適化するためにステップを追加する必要がある。図３ｃに示した並列線形エコーキャンセラの場合には、少なくとも２つの手法が存在る。
【００３７】
第１の手法では、すべてのパラメータを同等に扱い、上記の方法を用いてすべてのパラメータを同時に最適化する。第２の手法では、当技術分野で周知の標準的方法を用いて線形エコーキャンセラ（図３ｃの段階１Ｂ）のパラメータをまず設定し、残りの段階（図３ｃの段階１Ａ、２Ａ、および３Ａ）をオフにしておく。次に、すべての段階をオンにすることができる。次に、線形エコーキャンセラのパラメータを固定したまま、上記の方法を用いて残りのパラメータを最適化する。図３ｄに示した等化器／エコーキャンセラの多段組み合わせ（ＭＣＯＭ）の場合、最適化プロセスは少し複雑である。この場合、線路を介して同時に伝送・受信される既知の信号を持つ入力データを集める必要がある。次に、この受信データセットと既知の伝送データがＭＣＯＭの適切な入力側に送り込まれ、ＭＣＯＭにより処理される。次に、予測される結果（線路の反対端から伝送された既知の信号）から出力信号を減算することができる。その差を用いて、最適化ルーチンがＭＣＯＭの全パラメータにとって最良の選択値を繰り返し見出すことができるようにする誤差測定値が生成される。
【００３８】
この最適化プロセスを効果的に改良する方法が考えられる。二重ＭＮＥＣ（またはＭＣＯＭ）を配置すれば、最適化プロセス中に実際のデータの流れを中断しないように現在のパラメータセットを使用し続けながら、新しいパラメータセットの最適化を行うことができる。最適化にかかわるパラメータの数が多い場合、最適化の進行が遅くなることがあるため、この方法を使えば、長時間にわたってデータを処理できないような事態を避けることができる。現在、満足できるパラメータセットがない場合、好ましい方法は、最適化のレベルが高いパラメータセットを計算しながら使用できる、最適化のレベルが低いが機能的であるパラメータセットを迅速に計算することである。別の改良方法は、最適化プロセスを２つの（または、それ以上の）ループに分割することである。これは、最終段階が線形段階である場合に特に効果的である。この場合、最終段階のパラメータに対してのみ十分な最適化を実行し、その間それ以外のパラメータを固定したままにしておくことができる。最終段階は線形であるため、速い内部ループにおいてはこれを極めて短時間に実行することができる。その他のパラメータは、異なる最適化ルーチンを用いて、遅い外部ループにより調整される。その他のパラメータがすべて、遅いループを通過する間に、最終段階のパラメータは速いループにより再び十分に最適化される。速いループの代わりに、特異値分解などの手法を用いて最終段階にとって最適なパラメータを直接計算することもできる。
【００３９】
このＭＮＥＣにおいて使用できるような非線形関数では、最小限のエラーをもたらす係数セットが見つかると、真の「大域的最小値」ほど小さくはない「局所的最小値」が生成されることがある。このような曖昧さを解消する周知の手法が存在する。多様な初期条件から最適化プロセスを開始するか、あるいは、大域的最小値を見つけるための幾分遅い最適化ルーチンを使用しなければならないことがある。有用な手法としては、修正パウエル法、共役勾配法、多次元ディザリング、シミュレーティッドアニーリング法などがある（ただし、これらに限定されない）。最適化と最小化のためのこれらおよびその他の手段と方法は、たとえば、『Ｎｕｍｅｒｉｃａｌ  Ｒｅｃｉｐｅｓ：Ｃ言語で記述される科学的計算方法（Ｔｈｅ  Ａｒｔ  ｏｆ  Ｓｃｉｅｎｔｉｆｉｃ  Ｃｏｍｐｕｔｉｎｇ  ｗｉｔｈ  Ｌｉｓｔｉｎｇｓ  ｉｎ  Ｃ）』Ｗｉｌｌｉａｍ  Ｐｒｅｓｓ,  et al.,  Ｃａｍｂｒｉｄｇｅ  Ｕｎｉｖｅｒｓｉｔｙ  Ｐｒｅｓｓ,１９８９に記載されている。
【図面の簡単な説明】
【００４０】
【図１】周波数分割二重化またはＦＤＤを備えるシステムの周波数帯域を示している。上り帯域と下り帯域は分離しているが、上り帯域の非線形ひずみ帯域は下り帯域と大きく重なり合うことがあることに注意されたい。
【図２】本発明による多段非線形エコーキャンセラを組み込んだデジタル通信システムの一方の端の要素を示すブロック図である。
【図３ａ】本発明の多段非線形エコーキャンセラ（ＭＮＥＣ）の２つの実施形態を示している。
【図３ｂ】本発明の多段非線形エコーキャンセラ（ＭＮＥＣ）の２つの実施形態を示している。
【図３ｃ】２つのエコーキャンセラを並列に組み合わせた一実施形態を示している。
【図３ｄ】ＭＮＥＣを多段非線形等化器と組み合わせる方法を示している。
【図４】本発明の多段非線形エコーキャンセラの線形段階の一実施形態を示している。
【図５ａ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｂ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｃ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｄ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｅ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｆ】多段非線形エコーキャンセラの非線形段階の各実施形態を示している。
【図５ｇ】図３ｄに示した多段等化器とエコーキャンセラの組み合わせで用いることができる二重入力側非線形段階の一実施形態を示している。
【図６ａ】線形段階のフーリエ変換と離散的な周波数領域の実施形態をそれぞれ示している。
【図６ｂ】線形段階のフーリエ変換と離散的な周波数領域の実施形態をそれぞれ示している。
【図７ａ】帯域内エコー信号と帯域外エコー信号それぞれについて、本発明による多段非線形エコーキャンセラを用いた場合に、伝送信号の一連の非線形ひずみ係数付き信号要素の残留エコー信号が改善されていることを示している。
【図７ｂ】帯域内エコー信号と帯域外エコー信号それぞれについて、本発明による多段非線形エコーキャンセラを用いた場合に、伝送信号の一連の非線形ひずみ係数付き信号要素の残留エコー信号が改善されていることを示している。
【図８】コントローラを用いて、多段非線形エコーキャンセラの性能を制御するパラメータを最適化する方法を示している。
【図９】多段非線形エコーキャンセラの段階を特徴付ける関数のパラメータ値を初期化および設定するための、本発明による最適化手順を示す一般的な流れ図である。
【符号の説明】
【００４１】
１１０：上り帯域
１２０：下り帯域
１３０：帯域
２０２：入力データ
２０４：符号化および変調回路、回路
２０６：デジタル変調信号、信号、入力信号
２０８：デジタル／アナログ変換器、ＤＡＣ
２０７：ＭＮＥＣ、多段非線形エコーキャンセラ（ＭＮＥＣ）、多段エコーキャンセラ
２１０：アナログ形式、信号
２１２：アナログフィルタ
２１６：ハイブリッド回路、ハイブリッド
２１８：双方向通信チャネル、通信チャネル、線路、チャネル、双方向チャネル
２２０：受信信号、入力信号
２２２：アナログフィルタ
２２６：アナログ／デジタル変換器（ＡＤＣ）、ＡＤＣ
２２８：出力側
２３０：減算
２３４：等化器
２３８：復調および復号化
２４０：出力データ
４１０：関数
５１４：べき級数関数、関数
５１６：関数
５１８：１−Ｄ非線形関数
５２０：１−Ｄ関数
５２２：関数
５２４：関数
６１０：高速フーリエ変換（ＦＦＴ）
６２０：関数
７１０：曲線、線形エコーキャンセラ結果
７２０：曲線
７３０：曲線
７４０：曲線
７５０：曲線
７６０：曲線
８０２：テストデータ
８０６：ＭＮＥＣ
８０８：段階
８１０：段階
８１２：段階
８１４：減算
８１６：エコーデータ
８２０：減算
８２４：エコーキャンセラコントローラ、コントローラ
９４０：ステップ
９４２：ステップ
９４４：エラーレベル、ステップ
９４６：ステップ
９４８：
【図１】
2006505211.tif 000013
【図２】
2006505211.tif 000014
【図３ａ】
2006505211.tif 000015
【図３ｂ】
2006505211.tif 000016
【図３ｃ】
2006505211.tif 000017
【図３ｄ】
2006505211.tif 000018
【図４】
2006505211.tif 000019
【図５ａ】
2006505211.tif 000020
【図５ｂ】
2006505211.tif 000021
【図５ｃ】
2006505211.tif 000022
【図５ｄ】
2006505211.tif 000023
【図５ｅ】
2006505211.tif 000024
【図５ｆ】
2006505211.tif 000025
【図５ｇ】
2006505211.tif 000026
【図６ａ】
2006505211.tif 000027
【図６ｂ】
2006505211.tif 000028
【図７ａ】
2006505211.tif 000029
【図７ｂ】
2006505211.tif 000030
【図８】
2006505211.tif 000031
【図９】
2006505211.tif 000032
【国際調査報告】
2006505211.tif 000033
2006505211.tif 000034