(57)【特許請求の範囲】
【請求項１】
マルチパスや自局の送信部の送信信号の回り込み波を受信する受信部、送信部、及び、受信部からの信号が＋端子に入力され、適応フィルタの信号が−端子に入力される減算器と、前記適応フィルタとを備え、マルチパスや回り込み波をキャンセルする回り込みキャンセラ、を有する中継装置おいて、
前記回り込みキャンセラは、
減算器出力信号の周波数特性を算出し、算出した周波数特性と周波数特性の逆数を逆フーリエ変換して時間信号に変換する手段と、
時間信号から時間的に先行した成分を主波成分として抽出する手段と、
時間信号及び主波成分に基づいて、マルチパスや自局の送信信号の回り込み波をキャンセルするように前記適応フィルタのフィルタ係数を算出する手段を備えたことを特徴とする中継装置。
【請求項２】
請求項１において、算出した周波数特性に含まれるＦＦＴ窓位置の影響をサンプル時間未満の精度で除去する手段を備えたことを特徴とする中継装置。
【請求項３】
請求項１又は請求項２において、周波数特性と周波数特性の逆数をＮ倍（Ｎ≧２）の動作周波数に変換すると共に時分割し、
上記時分割信号を時分割演算処理にて逆フーリエ変換して時間信号に変換する手段を備えたことを特徴とする中継装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、直交周波数分割多重変調方式の中継装置に関わり、特にＳＦＮ（ Single Frequency Network：単一周波数ネットワーク）における中継装置の送受信アンテナ間での信号の回り込みをキャンセルする回り込みキャンセラに関する。
【背景技術】
【０００２】
地上デジタル放送等の直交周波数分割多重変調（ＯＦＤＭ：Orthogonal Frequency Division Multiplexing ）方式による中継方式には、オンエアの放送波を受信し、受信周波数と同一の周波数で再送伝するＳＦＮ（ Single Frequency Network ）中継方式がある。ＳＦＮ中継方式では、送受信の周波数が同一であるため、送信信号が受信アンテナに回り込み、信号の劣化や発振を引き起こす可能性がある。そのため、このような信号の回り込み（以下、単に「回り込み」と言う）をキャンセルする回り込みキャンセラを備えた中継方式が採用されており、例えば、特許文献１にその原理が記載されている。
また、親局から中継局までの伝搬路で生じたマルチパスのキャンセル方式に関しても、例えば、特許文献２に記載されている。
【０００３】
特許文献１若しくは特許文献２等に記載のキャンセル方式では、マルチパスや回り込みを受け、キャンセルした後の周波数特性Ｘ（ω）（即ち、親局から到来する希望波の周波数特性）を算出し、周波数特性Ｘ（ω）から主波Ｄを抽出し、式（１）に示す回り込み波のキャンセル残差成分ＥＣ（ω）を算出する。
【０００４】
【数１】
0005016444.tif 000002
また、マルチパスの等化に適したキャンセル残差成分ＥＭ（ω）を式（２）によって算出する。
【０００５】
【数２】
0005016444.tif 000003
ここで、主波Ｄは、暫定的に、周波数特性Ｘ（ω）を周波数領域（すべてのサブキャリア領域）で平均化した値としている。仮の算出方式を式（３）に示す。
【０００６】
【数３】
0005016444.tif 000004
キャンセル残差ＥＣ（ω）及びＥＭ（ω）は、式（４）または式（５）で示すＩＦＦＴ（逆フーリエ変換：Inverse Fast Fourier Transform 、又はＩＤＦＦＴ：Inverse Discrete Fast Fourier Transform ）処理により、それぞれ、キャンセル残差成分の時間信号ｅＭ（ｔ）及びｅＣ（ｔ）に変換される。そして、式（６）に示すフィルタ係数更新式に適用して、キャンセル動作を実現する。
【０００７】
【数４】
0005016444.tif 000005
【０００８】
【数５】
0005016444.tif 000006
【０００９】
【数６】
0005016444.tif 000007
式（６）において、ｅ（ｔ）はｅＣ（ｔ）、ｅＭ（ｔ）に基づき算出している。その詳細は、特許文献２に記載されている。また、μは０＜μ＜１の範囲で設定する。
【００１０】
また、受信したＯＦＤＭ信号から有効シンボル期間に相当する期間の信号をＦＦＴ（フーリエ変換：Fast Fourier Transform 、又はＤＦＦＴ：Discrete Fast Fourier Transform ）窓を抽出し、抽出されたＯＦＤＭ信号の有効シンボル期間（ＦＦＴ窓期間）をＦＦＴ処理することによって伝送回路特性を推定する回り込みキャンセラが、特許文献３に記載されている。
【００１１】
【特許文献１】特開平１１−３５５１６０号公報
【特許文献２】特開２００２−１５２０６５号公報
【特許文献３】特開２００４−３２０６７７号公報
【発明の開示】
【発明が解決しようとする課題】
【００１２】
従来技術において、周波数特性Ｘ（ω）の算出にはＯＦＤＭ信号に既知信号として挿入されているパイロットキャリアを用いることが多い。算出した周波数特性Ｘ（ω）は、ＯＦＤＭ時間信号に設けるＦＦＴ窓位置に起因する回転（位相回転）が生じる。
ＦＦＴ窓位置が、シンボル（シンボル：ＯＦＤＭ信号の１回の変調で送信する１ビット若しくは複数ビットのデータ）の有効シンボルタイミングに一致する時、周波数特性Ｘ（ω）に回転は発生しない。しかし、ＦＦＴ窓位置が有効シンボルタイミングに一致せず、時間ずれ（ＦＦＴ窓位置ずれ）が発生すると、有効シンボルからの時間ずれのずれ量τに応じて式（７）に示す回転が発生する（回転した周波数特性をＸ’（ω）とする）。
【００１３】
【数７】
0005016444.tif 000008
【００１４】
前述のように、主波Ｄは、周波数特性Ｘ（ω）に含まれる直流成分としている。しかし、式（７）に示すような回転が生じている場合には、周波数特性Ｘ’（ω）を用いて式（３）に示す主波抽出処理を行っても、主波Ｄを正しく算出することができない。
そのため、上述の特許文献３には、ＦＦＴ窓の位置ずれの影響を除去する方式が記載されている。しかし、低ＤＵ（ Desired to Undesired ）比のマルチパスが存在する場合には、マルチパスの影響によりＦＦＴ窓位置ずれを完全に除去することができない。また、主波よりもマルチパスのレベルの方が大きい場合には、マルチパスを主波と誤検出してしまう。
このように、回転が生じている周波数特性Ｘ’（ω）から算出したキャンセル残差Ｅ（ω）を時間信号に変換して式（６）によりキャンセル係数を算出した場合には、正しいタップ位置に係数を適用することかできず、信号品質劣化や発振などを引き起こしてしまうことがある。
【００１５】
本発明の目的は、上記のような問題を解決し、低ＤＵ比のマルチパスが存在する場合でも高精度な回り込みキャンセル動作を行うことができる回り込みキャンセラを備えた中継装置を提供することにある。
本発明の目的は、また、上記の中継装置において、算出した周波数特性に含まれるＦＦＴ窓位置の影響をサンプル時間未満の精度で除去することができる回り込みキャンセラを備えた中継装置を提供する。
【課題を解決するための手段】
【００１６】
上記の目的を達成するために、本発明の中継装置は、マルチパスや自局の送信部の送信信号の回り込み波を受信する受信部、送信部、及び、受信部からの信号が＋端子に入力され、適応フィルタの信号が−端子に入力される減算器と、前記適応フィルタとを備え、マルチパスや回り込み波をキャンセルする回り込みキャンセラ、を有する中継装置おいて、前記回り込みキャンセラは、減算器出力信号の周波数特性を算出し、算出した周波数特性と周波数特性の逆数を逆フーリエ変換して時間信号に変換する手段と、時間信号から時間的に先行した成分を主波成分として抽出する手段と、時間信号及び主波成分に基づいて、マルチパスや自局の送信信号の回り込み波をキャンセルするように前記適応フィルタのフィルタ係数を算出する手段を備えたものである。
また好ましくは、上記中継装置は、算出した周波数特性に含まれるＦＦＴ窓位置の影響をサンプル時間未満の精度で除去する手段を備えたものである。
また好ましくは、上記中継装置は、算出した周波数特性に含まれるＦＦＴ窓位置の影響をサンプル時間未満の精度で除去する手段を備え、周波数特性と周波数特性の逆数をＮ倍（Ｎ≧２）の動作周波数に変換すると共に時分割し、時分割信号を時分割演算処理にて逆フーリエ変換して時間信号に変換する手段を備えたものである。
【００１７】
また、本発明の中継装置は、受信装置からの信号が減算器の＋端子に供給され、減算器の−端子には減算器出力が適応フィルタを経由した信号が供給され、適応フィルタはマルチパスや自局の送信信号の回り込み波をキャンセルするように機能する中継装置おいて、減算器出力信号の周波数特性を算出し、算出した周波数特性と周波数特性の逆数を逆フーリエ変換して時間信号に変換し、逆フーリエ変換により算出された時間信号から時間的に先行した成分を主波成分として抽出し、時間信号及び主波成分に基づいて、マルチパスや自局の送信信号の回り込み波をキャンセルするように適応フィルタのフィルタ係数を算出する手段を備えたものである。
【発明の効果】
【００１８】
本発明によれば、低ＤＵ比のマルチパスが存在する場合でも高精度な回り込みキャンセル動作を行うことができる回り込みキャンセラを備えた中継装置が実現できる。
【発明を実施するための最良の形態】
【００１９】
本発明の第一の実施例を、図１を用いて説明する。図１は、本発明の回り込みキャンセラを備えた中継装置の構成を示すブロック図である。１０１は回り込みキャンセラ、１０２は受信アンテナ、１０３は受信変換部、１０４は送信変換部、１０５は送信アンテナ、１１は減算器、１２は回り込み信号の複製を作るトランスバーサルフィルタ、１３は周波数特性算出部、１４は回転補正部、１５は除算部、１６ａと１６ｂはＩＦＦＴ処理部、１７ａと１７ｂは主波抽出部、１８はキャンセル残差算出部、１９はＳｉｎｃ信号発生部、１ａはイメージ除去部、１ｂは係数更新部、１ｃは主波位置制御部である。また、トランスバーサルフィルタ１２は、例えば、ＦＩＲ（ Finite Impulse Response ）フィルタである。
なお、中継装置の受信アンテナ１０２と受信変換部１０３とを、受信部と称し、送信アンテナ１０５と送信変換部１０４とを、送信部と称する。
【００２０】
図１において、受信部側では、受信アンテナ１０２は、親局信号Ｘ（ω）及びそのマルチパスによる信号、並びに自局回り込み波（即ち、自局の送信部側の送信アンテナ１０５から出力される放送波の回り込み波）の合成された信号（受信信号）Ｘ’（ω）を受信アンテナで受信し、受信変換部１０３に出力する。受信変換部１０３は、受信アンテナ１０２から入力された信号Ｘ’（ω）に対して、フィルタ処理、周波数変換処理、等を行いベースバンド帯域に変換した後、その信号を回り込みキャンセラ１０１の減算器１１の被減算端子（＋端子）に出力する。
【００２１】
減算器１１には、受信変換部１０３から＋端子に入力された信号（Ｓａ）の他に、その減算端子（−端子）にトランスバーサルフィルタ１２からの出力信号（Ｓｂ）が入力され、出力端子からは、＋端子に入力された信号が−端子から入力された信号で減算された信号（Ｓｃ＝Ｓａ−Ｓｂ）が出力される。
この出力信号は、回り込みキャンセラの出力信号として送信部側の送信変換部１０４に出力されると共に、トランスバーサルフィルタ１２、及び周波数特性算出部１３に出力される。
トランスバーサルフィルタ１２は、入力された信号を、係数更新部１ｂから入力されるフィルタ係数を適用してフィルタリングした信号を減算器１１の−端子に出力する。
また、送信部側の送信変換部１０４は、入力された出力信号について、周波数変換、フィルタ処理、等を行って送信部側の送信アンテナ１０５から送信する。
【００２２】
周波数特性算出部１３は、入力された信号に挿入されているＯＦＤＭ信号特有のパイロットキャリアなどを用いて、周波数特性Ｘ’（ω）を推定する。
ここで、周波数特性算出部１３で算出（推定）した周波数特性Ｘ’（ω）は、式（７）に示したように、ＦＦＴ窓位置ずれにより理想的な周波数特性ＸＤ（ω）に対して位相回転が生じている。周波数特性算出部１３の出力信号Ｘ’（ω）は、回転補正部１４に出力される。
回転補正部１４は、主波位置制御部１ｃから入力される制御量τＤに基づいて、式（８）に示すように上記のＦＦＴ窓位置ずれによる位相回転成分を除去する。
【００２３】
【数８】
0005016444.tif 000009
【００２４】
式（７）に示したＦＦＴ窓位置ずれが生じている場合の周波数特性と、式（８）により、回転補正部１４への制御量τＤがＦＦＴ窓位置ずれのずれ量τと一致した時、周波数特性Ｘ（ω）は理想周波数特性ＸＤ（ω）に一致させることができる。
回転補正部１４は、入力された信号の周波数特性を理想周波数特性ＸＤ（ω）に一致させた信号を、除算部１５及びＩＦＦＴ処理部１６ｂに出力する。
この制御量τＤとＦＦＴ窓位置ずれのずれ量τを一致させる制御は、主波位置制御部１ｃにおいて行うが、主波位置制御部ｌｃの動作については後述する。
【００２５】
除算部１５からＩＦＦＴ処理部１６ａを通って主波抽出部１７ａに至る経路は、式（１）と式（４）で示した回り込みキャンセル残差に相当する信号を算出する。  また、ＩＦＦＴ処理部１６ｂから主波抽出部１７ｂまでの経路は、式（２）と式（５）に相当する信号を算出する。
従来技術では、低ＤＵのマルチパス環境下においては正しい主波を抽出することできず、誤ったキャンセル残差信号を生成してしまうという問題が挙げられた。これは、周波数領域で主波を抽出しているために発生する問題であるが、本発明では時間領域で主波を高精度に抽出し、時間領域信号からキャンセル残差信号を生成する。
【００２６】
以下に、主波抽出処理とキャンセル残差算出処理について図４に示すマルチパスと回り込みが混在するモデルを用いて説明する。図４は、親局信号のマルチパスと、自局（中継装置）の送信部側の送信アンテナからの送信波の回り込みが混在する伝播路モデルを模式的に示した図である。４１は伝播路特性算出部、４２は加算器、４３は回り込み伝播路特性算出部である。
図４において、主波Ｄが伝播路特性算出部４１に入力し、伝播路特性算出部４１が伝播路特性を算出して加算器４２の一方の端子に入力する。加算器４２は、他方の入力端子に入力される回り込み伝播路特性算出部４３からの出力信号と、一方の端子から入力される伝播路特性算出部４１の出力信号とを加算して、出力端子からキャンセル動作を行わない場合の周波数特性Ｘ（ω）の信号を出力すると共に、回り込み伝播路特性算出部４３に出力する。回り込み伝播路特性算出部４３は、入力された信号について回り込み伝播路特性を算出して、加算器４２の他方の入力端子に入力する。
図４のような伝播路モデルにおいて、親局から中継局（本発明の中継装置）までの伝搬路特性を｛１十Ｍ（ω）｝とし、回り込み伝搬路特性Ｃ（ω）とすると、キャンセル動作を行わない場合の周波数特性Ｘ（ω）は式（９）で表される。
【００２７】
【数９】
0005016444.tif 000010
ここで、Ｄは親局信号である主波とし、Ｘ（ω）にはＦＦＴ窓位置ずれ（ずれ量τ）が生じているものとする。除算部１５は、検出した周波数特性Ｘ（ω）に対して、式（１０）に示す複素除算演算を行う。
【００２８】
【数１０】
0005016444.tif 000011
式（１０）に式（９）に示すモデルを代入すると、式（１１）となる。
【００２９】
【数１１】
0005016444.tif 000012
なお、式（１０）においてデジタル演算精度を向上させるため、所定係数Ｋ（Ｋ＞１）を乗じても良い。
【００３０】
次に図１に戻って、除算部１５からの出力信号Ｙ（ω）は、ＩＦＦＴ処理部１６ａに出力される。ＩＦＦＴ処理部１６ａとＩＦＦＴ処理部１６ｂは同一の構成であり、出力信号Ｙ（ω）とＸ（ω）の信号を時間軸領域の時間信号ｙ（ｔ）、ｘ（ｔ）に変換する。
式（９）の時間領域変換を考えると、マルチパス特性Ｍ（ω）と回り込み特性Ｃ（ω）は、それぞれ遅延要素を含んでいる。このため、時間的に最も先行した成分は、主波ＤをＩＦＦＴ処理した信号となっている。
また同様に、式（１１）の時間領域変換を考えると、先行波成分は、１／ＤのＩＦＦＴ処理結果となっている。
【００３１】
この時間信号ｙ（ｔ）とｘ（ｔ）の一例を図５に示す。図５は、時間を横軸にし、Ｑ成分を縦軸に、Ｉ成分を斜軸にした時間信号の一例を示す図である。図５(a) は時間信号ｙ（ｔ）を示す図、図５(b) は時間信号ｘ（ｔ）を示す図である。
図５に示すように、先行波成分以外の成分について、時間信号ｙ（ｔ）に関しては回り込み成分、マルチパス成分の他にマルチパスのイメージ成分が生じている。また、時間信号ｘ（ｔ）にも、同様に回り込み波のイメージ成分が生じている。これらのイメージ成分は、後述するイメージ除去部１ａにて除去する。
【００３２】
ＩＦＦＴ処理部１６ａは、出力信号Ｙ（ω）の信号を時間信号ｙ（ｔ）に変換し、主波抽出部１７ａとキャンセル残差算出部１８に出力する。同様に、ＩＦＦＴ処理部１６ｂは、出力信号Ｘ（ω）の信号を時間信号ｘ（ｔ）に変換し、主波抽出部１７ｂとキャンセル残差算出部１８に出力する。
主波抽出部１７ａは、時間信号ｙ（ｔ）から主波成分ｙＤ（ｔ）を抽出し、キャンセル残差算出部１８に出力する。同様に、主波抽出部１７ｂは、時間信号ｘ（ｔ）から主波成分ｘＤ（ｔ）を抽出し、キャンセル残差算出部１８に出力する。
【００３３】
主波抽出部１７ａ、１７ｂの構成を図３に示し、主波抽出部の動作を図６を用いて説明する。図３は、本発明の主波抽出部の一実施例の構成を示すブロック図である。１７１は絶対値器、１７２は比較器、１７３は先行波位置検出部、１７４は保持器である。また、図６は、本発明の主波抽出部の動作の一例を説明するための図で、時間を横軸にし、信号成分の絶対値を縦軸にした時間信号の一例を示す図である。図６(a) は時間信号ｙ（ｔ）の絶対値｜ｙ（ｔ）｜を示す図、図６(b) は時間信号ｘ（ｔ）の絶対値｜ｘ（ｔ）｜を示す図である。
【００３４】
図３において、主波抽出部１７ａ若しくは１７ｂに入力された時間信号ｙ（ｔ）若しくはｘ（ｔ）は、絶対値器１７１と保持部１７４に入力される。絶対値器１７１は、入力された時間信号ｙ（ｔ）若しくはｘ（ｔ）の絶対値｜ｙ（ｔ）｜若しくは｜ｘ（ｔ）｜を算出し、比較器１７２に出力する。比較器１７２は、入力された時間信号の絶対値について、各時間毎に、雑音成分を抽出し、抽出された雑音成分の信号の値を０に置き換え処理を行う。比較器１７２は、処理後の信号を先行波位置検出部１７３に出力する。雑音成分の抽出は、例えば図６に示すように、信号成分と雑音成分を区別するための閾値を予め設定しておき、閾値以下の信号は雑音と判断することで実行する。
【００３５】
先行波位置検出部１７３は、入力された雑音成分を除去された時間信号から、時間的に最も先行しているピーク（先行ピーク）を検出し、検出した先行ピーク波を主波成分によるものであると判断し、先行ピーク波の先行ピーク位置信号ＰＯＳを保持器１７４に出力する。
先行ピーク位置信号ＰＯＳの一例を、図６の矢印の位置で示す。
先行ピークの具体的な検出方法としては、例えば、絶対値が０以外の値となる最も時間的に早いタイミングを検出し、検出したタイミングから数サンプル程度の窓を設け、窓内の最大値を先行ピークとする方法がある。
【００３６】
保持器１７４は、主波成分を表わす先行ピーク位置信号ＰＯＳと、時間信号ｙ（ｔ）若しくはｘ（ｔ）が入力される。保持器１７４は、先行ピーク位置信号ＰＯＳのタイミングの時間信号ｙ（ＰＯＳｙ）、ｘ（ＰＯＳｘ）の値を保持し、より真に近い主波Ｄの成分として、それぞれ、キャンセル残差算出部１８に出力する。
主波抽出部１７ａからの出力をｙＤ、主波抽出部１７ｂからの出力をｘＤとし、それぞれの出力を算出する式（１２）（１３）を以下に示す。
【００３７】
【数１２】
0005016444.tif 000013
【００３８】
【数１３】
0005016444.tif 000014
【００３９】
これら主波成分ｙＤ及びｘＤが入力されるキャンセル残差算出部１８には、更に、ＩＦＦＴ処理部１６ａ及び１６ｂの出力である時間信号ｙ（ｔ）及びｘ（ｔ）、並びに、Ｓｉｎｃ信号発生部１９の出力信号ｓｉｎｃ（ｔ）が入力される。
キャンセル残差算出部１８は、式（１）（２）で説明したキャンセル残差信号と等化な演算を行う。即ち、式（１）（２）を、式（４）（５）に代入し、変形を行い、キャンセル残差成分の時間信号ｅＭ（ｔ）及びｅＣ（ｔ）を算出し、イメージ除去部１ａに出力する。
キャンセル残差成分の時間信号ｅＭ（ｔ）及びｅＣ（ｔ）の算出式を式（１４）（１５）に示す。また、係数αについては後述する。
【００４０】
【数１４】
0005016444.tif 000015
【００４１】
【数１５】
0005016444.tif 000016
【００４２】
ところで、ＩＦＦＴ処理部１６ａ、１６ｂにてＩＦＦＴ処理する際、ＩＦＦＴの処理ポイント数より出力信号Ｙ（ω）及びＸ（ω）のポイント数が少ない場合が多い。このため、出力信号Ｙ（ω）及びＸ（ω）に、０を付加してポイント数を一致させた後にＩＦＦＴ処理する。このような処理方法は、例えば、特開２００４−１５３７９９号公報に記載されている。
【００４３】
このような処理では、ＩＦＦＴ結果に、時間的な拡がりが生じることとなる。式（１４）、（１５）におけるＩＦＦＴ［１］も出力信号Ｙ（ω）及びＸ（ω）のポイント数分のみが１であり、それ以外は出力信号Ｙ（ω）及びＸ（ω）のＩＦＦＴ処理と同様に０とする。従って、式（１４）、（１５）におけるＩＦＦＴ［１］はｓｉｎｃ関数で表されることになる。
このｓｉｎｃ信号はｓｉｎｃ信号発生器１９にて生成する。ｓｉｎｃ信号の生成方式としては、例えば、予め計算させておいたｓｉｎｃ信号のデータをＲＯＭ（ Read Only Memory ）等のメモリに記憶させることもできる。これにより、直接演算する構成よりハードウェアの規模を削城することができる。
Ｓｉｎｃ信号発生器１９の出力信号ｓｉｎｃ（ｔ）を用いて式（１２）（１３）を表すと式（１６）（１７）となる。
【００４４】
【数１６】
0005016444.tif 000017
【００４５】
【数１７】
0005016444.tif 000018
【００４６】
係数αは、ＩＦＦＴのポイント数と出力信号Ｙ（ω）及びＸ（ω）のポイント数の差異を補償する所定係数を示している。キャンセル残差算出部１８では、式（１６）及び（１７）の演算を行う。以上の処理により、低ＤＵマルチパス環境下であってもキャンセル残差信号ｅＭ（ｔ）及びｅＣ（ｔ）を高精度に算出することができる。
【００４７】
キャンセル残差信号ｅＭ（ｔ）及びｅＣ（ｔ）が入力されたイメージ除去部１ａは、イメージ成分を除去した信号ｅ（ｔ）を係数更新部１ｂに出力する。イメージ成分の除去に関しては、特許文献２の記載の如く行う。従って、ここでは説明を省略する。
キャンセル残差信号ｅ（ｔ）を入力された係数更新部１ｂは、フィルタ係数を算出し、算出したフィルタ係数をトランスバーサルフィルタ１２に出力する。これによって、トランスバーサルフィルタ１２は、最適な係数に更新されて、適用される。
【００４８】
上記の説明では、主波Ｄの検出精度は１サンプル単位である。しかし、高精度なキャンセル動作を行うためにはサンプル単位未満の検出精度が必要となる。これは、前述した回転補正部１４において、サンプル単位未満の回転補正を行うことで実現できる。
ＦＦＴ窓位置ずれのずれ量τと回転補正部１４の制御量τＤが一致した時、周波数特性Ｘ（ω）は、理想周波数特性ＸＤ（ω）に一致させることができる。そして、主波抽出部１７ａ及び１７ｂにおいて算出する主波位置ＰＯＳは、０となる。
【００４９】
主波位置制御部１ｃは、τ＝τＤとなるように制御量τＤを算出して位相回転量を制御する。即ち、主波位置制御部１ｃには、主波抽出部１７ａ及び１７ｂからの信号が入力され、この信号に基づいて制御を行う。
具体的には、主波抽出部１７ａ及び１７ｂは、主波位置となる先行ピーク位置信号ＰＯＳを算出すると共に、先行ピーク位置信号ＰＯＳの前後１サンプルの時間信号ｙ（ｔ）及びｘ（ｔ）を主波位置制御部１ｃに出力する。主演位置制御部１ｃは、制御量τＤの算出方式として、理想的な制御量τとの誤差ｅｒｒをサンプル未満の精度で算出し、誤差ｅｒｒが０若しくは所定値になるように制御すれば良い。誤差ｅｒｒは、式（１８）及び（１９）で定義する。
【００５０】
【数１８】
0005016444.tif 000019
【００５１】
【数１９】
0005016444.tif 000020
【００５２】
式（１８）及び（１９）は、誤差ｅｒｒが０に収束することを前提とした式である。式（１８）及び（１９）は、時間信号ｙ（ｔ）、ｘ（ｔ）において、サンプル単位のＦＦＴ窓位置誤差として式（１８）及び（１９）の第一項を主波位置である先行ピーク位置信号ＰＯＳとしている。また、サンプル単位未満の誤差にとしては、第二項で示した先行ピーク位置信号ＰＯＳの前後サンプルの絶対値の差分を主波レベルで正規化することにより算出している。
また、式（１９）の第一項において、ｙ（ＰＯＳｙ）にマイナスの符号が付加されている。これは、出力信号Ｙ（ω）が出力信号Ｘ（ω）の逆数であるため、ＦＦＴ窓位置の影響による時間シフト量が、時間信号ｘ（ｔ）と比較して符号反転しているためである。
第二項に関しては、ＩＦＦＴのポイント数と出力信号Ｙ（ω）及びＸ（ω）の周波数特性のポイント数が異なること起因する波形の時間拡がりを利用している。  τ＝τＤである場合には、主波成分により生成された先行ピークの前後１サンプルは同値になるが、τ≠τＤの場合には前後サンプルの絶対値に差が発生するため、この特性を利用してサンプル未満の位置誤差ｅｒｒを算出している。
【００５３】
図７は、ＦＦＴ窓位置誤差τに対する検出誤差ｅｒｒを示している。図７は、ＦＦＴ窓位置誤差と検出誤差ｅｒｒの関係の一例を示す図である。横軸がＦＦＴ窓位置誤差τ、縦軸が検出誤差ｅｒｒを示す。
図７から分かる様に、ほぼ線形的な特性を有していることか分かる。
この誤差ｅｒｒを用いて、誤差ｅｒｒが０若しくは所定値になるように制御を行う。制御方式としては、一般的なＰＩＤ（ Proportional Integral Derivative ）制御やＰＩ（ Proportional Integral ）制御等のフィードバック制御方式を用いれば良い。
【００５４】
また、ＦＦＴ窓位置ずれの補償方式は、上記で説明した方式であっても良いが、Ａ／Ｄコンバータをサンプリングするサンプリングクロックの位相を制御する構成であっても良い。
【００５５】
以上説明した本発明の第一の実施例により低ＤＵのマルチパス環境下であっても、高精度なキャンセル動作を実現できる。
【００５６】
次に本発明の第二の実施例について図２を用いて詳細に説明する。図２は、本発明の中継装置における回り込みキャンセラの一実施例の構成を示すブロック図である。第二の実施例は、第一の実施例における図１の中継装置の回り込みキャンセラー１０１の別の実施例である。図２の回り込みキャンセラー２０１は、図２の回り込みキャンセラー１０１において、ＩＦＦＴ処理部１６ａと１６ｂを一系統のＩＦＦＴ処理部１６とし、主波抽出部１７ａと１７ｂを一系統の主波抽出部１７とした他、除算部１５とＩＦＦＴ処理部１６の間に、時分割ＭＵＸ（多重：Multiplex ）処理部１ｄを付加し、主波抽出部１７とキャンセル残差算出部１８の間に時分割ＤＥＭＵＸ（多重分離：De Multiplex ）処理部１ｅを付加した構成である。それ以外の同一の符号を付した構成部（処理部）は、第一の実施例と同等の機能を有する。
【００５７】
第一の実施例では、２系統のＩＦＦＴ処理部１６ａ及び１６ｂ、並びに、主波抽出部１７ａ及び１７ｂを用いて構成していた。しかし、ＩＦＦＴ処理では多くの演算が必要であり、ＩＦＦＴ処理を２系統有することはハードウェア規模か大きくなってしまう。
従って、本発明の第二の実施例では、高速の演算を行うことにより、１系統のハードウェアを用いて第一の実施例と等価な処理を実現する。
【００５８】
第二の実施例の動作について、図２と図８を参照して説明する。図８は、本発明の中継装置の、回り込みキャンセラーの一実施例のタイミングチャートである。横軸は時間で、左から右に時間が推移している。図８（ａ）は、時分割ＭＵＸ処理部１ｄに入力される出力信号Ｙ（ω）とＸ（ω）のタイミングを示す図、図８（ｂ）は、ＩＦＦＴ処理部１６に入力される信号のタイミングを示す図で、時分割ＭＵＸ処理部１ｄで時分割多重された信号である。図８（ｃ）は、ＩＦＦＴ処理部１６で処理された信号が主波抽出部１７に入力する信号のタイミングを示す図で、図８（ｄ）は、時分割ＤＥＭＵＸ処理部１ｅに入力される主波成分の信号のタイミングを示す図で、主波抽出部１７が、時間信号ｙ（ｔ）及びｘ（ｔ）からそれぞれ抽出した主波成分ｙＤ（ｔ）及びｘＤ（ｔ）である。また、図８（ｅ）は、キャンセル残差算出部１８に入力される信号のタイミングを示す図である。
【００５９】
図２において、時分割ＭＵＸ処理部１ｄには、出力信号Ｙ（ω）とＸ（ω）が入力される（図８（ａ））。時分割ＭＵＸ処理部ｌｄは、後段のＩＦＦＴ処理部１６で時分割処理させるため、出力信号Ｙ（ω）とＸ（ω）のクロック周波数をＮ倍（Ｎ≧２）すると共に、時分割多重してＩＦＦＴ処理部１６に出力する（図８（ｂ））。
ＩＦＦＴ処理部１６は、Ｎ倍のクロック周波数で演算を行い、第一の実施例と比較して１／Ｎの時間で演算を終了する（図８（ｃ））。
主波抽出部１７でも同様に、高速演算を行い、主波成分ｙＤの算出と主波成分ｘＤの算出を行うための演算回路を共用しながら、時分割処理を行う（図８（ｄ））。
このようにして得られた時間信号ｙ（ｔ）及びｘ（ｔ）、並びに、主波成分ｙＤ及びｘＤは、時分割ＤＥＭＵＸ処理部１ｅに出力される。
時分割ＤＥＭＵＸ処理部１ｅは、時分割ＭＵＸ処理部１ｄと逆に、動作周波数を元に戻す処理を行い、２系統の信号を算出する（図８（ｅ））。
キャンセル残差算出部１８及びそれ以降の処理は、第一の実施例と同様であり、説明は省略する。以上の処理により、第一の実施例に対してハードウェア規模を大幅に削減させることができる。
【００６０】
上記実施例では、干渉波をキャンセルするための減算器出力信号の周波数特性を算出し、算出した周波数特性と、その逆数をそれぞれ逆フーリエ変換して時間信号を算出する。更に、算出した時間信号から時間的に最も先行した成分を主波成分と判断して抽出する。そして、時間信号から主波成分を算出することにより低ＤＵのマルチパス環境下でも高精度に主波成分を抽出することかできる。
この時間信号と主波成分に基づいてキャンセル残差成分を算出し、キャンセル残差成分は適応フィルタの係数として使用し、適応フィルタでは干渉波のレプリカを生成する。生成したレプリカ信号を、干渉波が混入している入力信号から減算することにより、低ＤＵのマルチパス環境下でも干渉波を除去することができる。
【図面の簡単な説明】
【００６１】
【図１】本発明の中継装置の一実施例の構成を示すブロック図。
【図２】本発明の中継装置の回り込みキャンセラの一実施例の構成を示すブロック図。
【図３】本発明の主波抽出部の一実施例の構成を示すブロック図。
【図４】マルチパスと回り込みが混在する伝播路モデルを模式的に示した図。
【図５】時間を横軸にした時間信号の一例を示す図。
【図６】本発明の主波抽出部の動作の一例を説明するための図。
【図７】ＦＦＴ窓位置誤差と検出誤差ｅｒｒの関係の一例を示す図。
【図８】本発明の一実施例のタイミングチャート。
【符号の説明】
【００６２】
１ａ：イメージ除去部、  １ｂ：係数更新部、  １ｃ：主波位置制御部、  １１：減算器、  １ｄ：時分割ＭＵＸ処理部、  １ｅ：時分割ＤＥＭＵＸ処理部、  １２：トランスバーサルフィルタ、  １３：周波数特性算出部、  １４：回転補正部、  １５：除算部、  １６ａ，１６ｂ：ＩＦＦＴ処理部、  １７ａ，１７ｂ：主波抽出部、  １８：キャンセル残差算出部、  １９：Ｓｉｎｃ信号発生部、  １０１：回り込みキャンセラ、  １０２：受信アンテナ、  １０３：受信変換部、  １０４：送信変換部、  １０５：送信アンテナ、  １７１：絶対値器、  １７２：比較器、  １７３：先行波位置検出部、  １７４：保持器。
【図１】
0005016444.tif 000021
【図２】
0005016444.tif 000022
【図３】
0005016444.tif 000023
【図４】
0005016444.tif 000024
【図５】
0005016444.tif 000025
【図６】
0005016444.tif 000026
【図７】
0005016444.tif 000027
【図８】
0005016444.tif 000028