(57)【要約】【課題】遅延プロファイル測定において、時間分解能を向上させる。【解決手段】フーリエ変換された受信信号Ｓ（ｋ）について、Ｎ→Ｅ・Ｎポイント数増加部１４で離散ポイント数を増加させＳ０（ｋ）、Ｅ・ＮポイントＩＤＦＴ部１６で逆離散フーリエ変換して、複素数列ｓ０（ｎ）を取得する。最大振幅値・遅延時間検出部２６で、各補正ステップｉ（０≦ｉ≦Ｊ−１）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント，０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、およびそのときの離散ポイントｎｉを検出する。減算部２０では、各補正ステップｉについて、その検出結果と補正関数ｂ（ｎ）とを含む所定の補正式ｓｉ＋１（ｎ）＝ｓｉ（ｎ）−Ｐｉ・ｂ（ｎ−ｎｉ）を用いて複素数列ｓｉ（ｎ）を取得する。このｓｉ＋１（ｎ）が次の補正ステップｉ＋１の複素数列ｓｉ＋１（ｎ）となる。上記補正ステップｉを所定回数（１回以上）反復して得られた複素数列ｓｉ（ｎ）を、遅延プロファイルとして取得する。【選択図】　　　　図１
2004023405.tif 000002
【特許請求の範囲】
【請求項１】
ＯＦＤＭ信号をＮポイント（Ｎ：整数）で逆離散フーリエ変換して取得された離散伝達関数（複素インパルス応答）Ｓ（ｋ）（ｋ：離散ポイント（整数）、０≦ｋ≦Ｎ−１）から、次の式（１）
【数１】
2004023405.tif 000020
（ここに、Ｅ：１より大きい有理数）により離散ポイント数をＥ・Ｎに増加した離散伝達関数Ｓ０（ｋ）を取得する離散ポイント数増加部と、
前記離散伝達関数Ｓ０（ｋ）にＥ・Ｎポイントで逆離散フーリエ変換を施し、Ｅ・Ｎポイントの複素数列ｓ０（ｎ）を取得する逆離散フーリエ変換部と、
各補正ステップｉ（０≦ｉ≦Ｊ−１，ｉ，Ｊともに整数）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント（整数），０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、および該最大値｜Ｐｉ｜をとる離散ポイントｎとしての最大値ポイントｎｉを検出する検出部と（ただし、ｎ≠ｎｉｐｒｅｖ，ｎｉｐｒｅｖ：現補正ステップｉより前の補正ステップｉｐｒｅｖ（０≦ｉｐｒｅｖ＜ｉ）での最大値ポイント）、
各補正ステップｉについて、前記最大値｜Ｐｉ｜および最大値ポイントｎｉに基づいて次の式（２）
【数２】
2004023405.tif 000021
（ここに、ｂ（ｎ）：補正関数）により複素数列ｓｉ＋１（ｎ）を取得する補正処理部と、を備え、補正ステップｉ（整数、０≦ｉ≦Ｊ−１）をＪ回デクリメントして取得された複素数列ｓＪ（ｎ）を、ＯＦＤＭ信号の遅延プロファイル（複素インパルス応答）として取得する遅延プロファイル測定装置。
【請求項２】
前記最大値｜Ｐｉ｜に基づいてデクリメント数Ｊを決定することを特徴とする請求項１に記載の遅延プロファイル測定装置。
【請求項３】
前記補正関数ｂ（ｎ）の振幅が最大となる離散ポイントより時間的に前および後の離散ポイント領域を含む離散ポイントについてのみ補正関数ｂ（ｎ）を保持する補正関数記憶部を備えることを特徴とする請求項１または２に記載の遅延プロファイル測定装置。
【請求項４】
請求項１〜３のうちいずれか一つに記載の遅延プロファイル測定装置によって取得された遅延プロファイルを用いて受信信号（再送信信号）に含まれる回り込み波および遅延波成分を除去する回り込みキャンセラ。
【請求項５】
請求項１〜３のうちいずれか一つに記載の遅延プロファイル測定装置によって取得された遅延プロファイルを用いて受信信号に含まれる遅延波成分を除去する伝送路等化器。
【請求項６】
ＯＦＤＭ信号をＮポイント（Ｎ：整数）で逆離散フーリエ変換して取得された離散伝達関数（複素インパルス応答）Ｓ（ｋ）（ｋ：離散ポイント（整数）、０≦ｋ≦Ｎ−１）から、次の式（１）
【数３】
2004023405.tif 000022
（ここに、Ｅ：１より大きい有理数）により離散ポイント数をＥ・Ｎに増加した離散伝達関数Ｓ０（ｋ）を取得するステップと、
前記離散伝達関数Ｓ０（ｋ）にＥ・Ｎポイントで逆離散フーリエ変換を施し、Ｅ・Ｎポイントの複素数列ｓ０（ｎ）を取得するステップと、
各補正ステップｉ（０≦ｉ≦Ｊ−１，ｉ，Ｊともに整数）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント（整数），０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、および該最大値｜Ｐｉ｜をとる離散ポイントｎとしての最大値ポイントｎｉを検出するステップと（ただし、ｎ≠ｎｉｐｒｅｖ，ｎｉｐｒｅｖ：現補正ステップｉより前の補正ステップｉｐｒｅｖ（０≦ｉｐｒｅｖ＜ｉ）での最大値ポイント）、
各補正ステップｉについて、前記最大値｜Ｐｉ｜および最大値ポイントｎｉに基づいて次の式（２）
【数４】
2004023405.tif 000023
（ここに、ｂ（ｎ）：補正関数）により複素数列ｓｉ＋１（ｎ）を取得するステップと、
を含み、補正ステップｉ（整数、０≦ｉ≦Ｊ−１）をＪ回デクリメントして取得された複素数列ｓＪ（ｎ）を、ＯＦＤＭ信号の遅延プロファイル（複素インパルス応答）として取得する遅延プロファイル測定方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号の受信、並びに中継において、遅延波や単一周波数ネットワーク（Ｓｉｎｇｌｅ　Ｆｒｅｑｕｅｎｃｙ　Ｎｅｔｗｏｒｋ：ＳＦＮ）中継局において生じる回り込み波が存在する伝送路の遅延プロファイル（複素インパルス応答）の時間分解能向上方法、および測定装置、並びにその遅延プロファイル推定方法を利用した回りこみキャンセラ、および伝送路等化器に関し、特に、得られた遅延プロファイル（複素インパルス応答）の時間分解能をさらに細かい時間分解能とし、なおかつ従来法に比べて正確な振幅値を得ることで、遅延プロファイル（複素インパルス応答）を推定し、並びにその遅延プロファイル（複素インパルス応答）の時間分解能向上方法を利用した、より精度の高い回り込みの除去を行う回り込みキャンセラ、および伝送路等化を行う伝送路等化器に関する。
【０００２】
【従来の技術および発明が解決しようとする課題】
従来の遅延プロファイル（複素インパルス応答）推定方法、および測定装置では、得られる遅延プロファイル（複素インパルス応答）の時間分解能は制限されてしまう。従来の遅延プロファイル（複素インパルス応答）推定方法としては、相関法と、何等かの方法により伝達関数を得て（例えば、パイロット信号を用いて得る）からこの伝達関数を逆フーリエ変換して遅延プロファイル（複素インパルス応答）を求める方法の二つがあるが、このうち前者（相関法）では、時間信号をサンプリングした際のサンプリング周波数によって、遅延プロファイル（複素インパルス応答）の時間分解能が制限される。また、後者（何等かの方法により伝達関数を得てからこの伝達関数を逆フーリエ変換して遅延プロファイル求める方法）では、伝達関数の帯域幅によって、遅延プロファイル（複素インパルス応答）の時間分解能が制限される。
【０００３】
遅延プロファイル（複素インパルス）の時間分解能が粗いと、時間分解能の逓倍（整数倍）ではない遅延時間の遅延波、もしくは回り込み波が存在する場合、インパルス応答値が低下し、なおかつ広い時間範囲に渡りインパルス応答が広がる。また、時間分解能の１サンプル間隔に数波の遅延波、および回り込み波が存在する場合に、遅延波、および回り込み波が分離、識別できない。
【０００４】
しかしながら、伝送路状況をより詳細に調査したい場合、サンプリング周波数を高くするか、あるいは伝達関数の帯域幅を広げて測定するしか方法がなかった。
【０００５】
ところで、得られた遅延プロファイル（複素インパルス応答）の時間分解能を任意に向上させる方法としては、自己回帰モデル（Ａｕｔｏｒｅｇｒｅｓｓｉｎｖｅ　Ｍｏｄｅｌ：ＡＲモデル）、最大エントロピー法（Ｍａｘｉｍｕｍ　Ｅｎｔｒｏｐｙ　Ｍｅｔｈｏｄ：ＭＥＭ）等を利用する方法がある。しかしながら、これらの方法においては次の問題点を抱えている。
［１］モデルの次数推定決定が困難である。
［２］モデルの次数を的確に求めたとしても振幅値が正確に求まらない。
モデルの次数推定を誤ると得られる遅延プロファイル（複素インパルス応答）結果は間違ったものとなる。そればかりか、モデルの次数推定が正確に行われたとしても得られる遅延プロファイル（複素インパルス応答）結果における振幅値が正確でないため、厳密な意味で時間分解能を向上する方法とは言えない。このようなことから、遅延プロファイル（複素インパルス応答）の時間分解能を向上させ、なおかつ振幅値もより高精度に求めることが可能な方法が希求されている。
【０００６】
一方、遅延プロファイル（複素インパルス応答）推定方法により得られた遅延プロファイル（複素インパルス応答）結果を利用する回り込みキャンセラ、および伝送路等化器においては、遅延プロファイル（複素インパルス応答）の遅延時間の時間分解能により等化精度を含めた等化能力（もしくはキャンセル能力）が不十分な場合があった。このため、特に、遅延波、および回り込み波等の削除しようとする非希望波の遅延時間が時間分解能の逓倍（整数倍）とならないとき、削除すべき遅延波および回り込み波が１波であっても、多数の遅延波による近似により等化を行うこととなるため希望波電力対非希望波電力比レベルが悪くなる（小さくなる）と等化精度が低下し、等化可能な希望波電力対非希望波電力比レベルが劣化し（大きくなり）、等化能力が低下するという課題があった。また、上述の自己回帰モデル、最大エントロピー法を利用して得られた時間分解能を向上させた遅延プロファイル（複素インパルス応答）結果はモデルの次数推定、および振幅値の不確定さから、回り込みキャンセラおよび伝送路等化器において用いるのは適切とは言えなかった。
【０００７】
そこで、本発明の目的は、ＯＦＤＭ方式によるデジタル伝送において、何等かの理由により時間分解能がＴまでの遅延プロファイルしか得られていない場合に、その遅延プロファイル（複素インパルス応答）の時間分解能をさらに向上させ、細かい時間分解能にて遅延プロファイル（複素インパルス応答）を推定する方法、および測定装置を提供せんとするものである。また、本発明を用いて得られるより細かい時間分解能の遅延プロファイルを利用し、従来の回り込みキャンセラ、および伝送路等化器に比べ、伝送路等化精度（非希望波キャンセル精度）が高く、等化可能な希望波電力対非希望波電力比レベルを改善（小さく）した安定な単一周波数ネットワーク中継装置向け回り込みキャンセラ、および伝送路等化器を提供せんとするものである。
【０００８】
【課題を解決するための手段】
本発明にかかるＯＦＤＭ信号をＮポイント（Ｎ：整数）離散フーリエ変換して取得された離散伝達関数Ｓ（ｋ）（ｋ：離散ポイント（整数）、０≦ｋ≦Ｎ−１）から、次の式（１）
【数５】
2004023405.tif 000003
（ここに、Ｅ：１より大きい有理数）により離散ポイント数をＥ・Ｎに増加した離散伝達関数Ｓ０（ｋ）を取得する離散ポイント数増加部と、上記離散伝達関数Ｓ０（ｋ）にＥ・Ｎポイントで逆離散フーリエ変換を施し、Ｅ・Ｎポイントの複素数列ｓ０（ｎ）を取得する逆離散フーリエ変換部と、各補正ステップｉ（０≦ｉ≦Ｊ−１，ｉ，Ｊともに整数）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント（整数），０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、および該最大値｜Ｐｉ｜をとる離散ポイントｎとしての最大値ポイントｎｉを検出する検出部と（ただし、ｎ≠ｎｉｐｒｅｖ，ｎｉｐｒｅｖ：現補正ステップｉより前の補正ステップｉｐｒｅｖ（０≦ｉｐｒｅｖ＜ｉ）での最大値ポイント）、各補正ステップｉについて、上記最大値｜Ｐｉ｜および最大値ポイントｎｉに基づいて次の式（２）
【数６】
2004023405.tif 000004
（ここに、ｂ（ｎ）：補正関数）により複素数列ｓｉ＋１（ｎ）を取得する補正処理部と、を備え、補正ステップｉ（整数、０≦ｉ≦Ｊ−１）をＪ回デクリメントして取得された複素数列ｓＪ（ｎ）を、ＯＦＤＭ信号の遅延プロファイル（複素インパルス応答）として取得する。
【０００９】
従来の遅延プロファイル（複素インパルス応答）推定方法、および測定装置では、前述のようにサンプリング周波数、または得られる伝達関数の帯域幅により得られる遅延プロファイルの時間分解能に制約があった。また、仮に自己回帰モデルや最大エントロピー法を用いて遅延プロファイル（複素インパルス応答）の時間分解能を向上させたとしても、モデルの次数推定、および振幅値の不確定さから、厳密な意味での遅延プロファイル（複素インパルス応答）の時間分解能向上方法とは言えなかった。遅延プロファイル（複素インパルス応答）の時間分解能向上方法は、時間分解能Ｔの遅延プロファイル（複素インパルス応答）を扱う離散フーリエ変換のポイント数Ｎ（Ｎ：２以上の整数、ただし２のべき乗とは限らない）で離散フーリエ変換し、得られたポイント数Ｎの伝達関数をポイント数Ｅ・Ｎ（Ｅ：有理数）の離散フーリエ変換の周波数領域の低域部に配置し、それ以外の周波数成分は零とし、再び逆離散フーリエ変換し、最大値検出、および補正関数を用いた補正処理を繰り返し行うことで時間分解能Ｔ／Ｅの遅延プロファイルを生成することで実現される。また得られた時間分解能Ｔ／Ｅの遅延プロファイル（複素インパルス応答）を用いることで、従来の回り込みキャンセラ、および伝送路等化器では削除不能、もしくは削除可能であっても削除しきれなかった遅延波、および回り込み波を削除することが可能となる。
【００１０】
また、本発明にかかる遅延プロファイル測定装置では、上記最大値｜Ｐｉ｜に基づいてデクリメント数Ｊを決定するのが好適である。
【００１１】
また、本発明にかかる遅延プロファイル測定装置では、上記補正関数ｂ（ｎ）の振幅が最大となる離散ポイントより時間的に前および後のいくつかの離散ポイント領域を含む離散ポイントについてのみ補正関数ｂ（ｎ）を保持する補正関数記憶部を備えるのが好適である。
【００１２】
【発明の実施の形態】
図１は、本発明第一の実施形態にかかる遅延プロファイル（複素インパルス応答）測定装置１０のブロック図である。
【００１３】
まず、Ｎポイント離散フーリエ変換部（ＮポイントＤＦＴ部）１２において、公知の一般的な遅延プロファイル（複素インパルス応答）測定装置（推定方法）によって得られたＮポイントの遅延プロファイル（複素インパルス応答）結果ｓ（ｎ）（ただし、０≦ｎ≦Ｎ−１）が、次の式（３）によりＮポイントの離散伝達関数Ｓ（ｋ）に変換される。
【数７】
2004023405.tif 000005
次いで、Ｎ→Ｅ・Ｎポイント数増加部（離散ポイント数増加部に相当する）１４において、Ｓ（ｋ）（ただし、０≦ｋ≦Ｎ−１）のＥ倍の値をＥ・Ｎポイントの低域部のＮポイントに配置し、Ｅ・Ｎポイントの残りの（Ｅ−１）・Ｎポイントには全て零を設定して、次の式（１）により離散伝達関数Ｓ０（ｋ）（ただし、０≦ｋ≦Ｅ・Ｎ−１）が得られる。
【数８】
2004023405.tif 000006
次いで、Ｅ・Ｎポイント逆離散フーリエ変換部（Ｅ・ＮポイントＩＤＦＴ部）１６において、Ｓ０（ｋ）は次の式（４）により複素数列ｓ０（ｎ）（ただし、０≦ｎ≦Ｅ・Ｎ−１）に変換される。
【数９】
2004023405.tif 000007
次いで、ｓｉ（ｎ）保持部１８において、ｉ＝０ではｓ０（ｎ）を、ｉ＞０では減算部２０で得られたｓｉ＋１（ｎ）を、ｓｉ（ｎ）として保持する。そして減算部（補正処理部に相当する）２０は、次の式（２）により、ｓｉ（ｎ）から複素乗算・可変遅延部２２において得たＰｉ・ｂ（ｎ−ｎｉ）を減算し、
【数１０】
2004023405.tif 000008
を得る。ただし、ｎ−ｎｉ＜０のとき、ｂ（ｎ−ｎｉ）＝ｂ（Ｅ・Ｎ＋ｎ−ｎｉ）である。ｓｉ＋１（ｎ）はｓｉ（ｎ）保持部１８に出力される。このとき、ｉ＋１をｉとする。なお、補正関数ｂ（ｎ）は、補正関数記憶部２４に記憶されており、これが複素乗算・可変遅延部２２に出力される。この補正関数ｂ（ｎ）は、Ｅ・Ｎポイント離散フーリエ変換の周波数領域において、低域部Ｎポイントの値をＥとし、残りの（Ｅ−１）・Ｎポイントの値を零としたスペクトルをＥ・Ｎポイント逆離散フーリエ変換することで得られた結果である。すなわち、補正関数ｂ（ｎ）は、下記の式（５）で表せる。
【数１１】
2004023405.tif 000009
最大振幅値・遅延時間検出部２６は、集合Ａ（初期状態時Ａ＝φ（空集合））に属さない遅延時間（遅延ポイント数）の中からｓｉ（ｎ）の振幅値｜ｓｉ（ｎ）｜の最大値｜Ｐｉ｜、およびその｜Ｐｉ｜を持つＰｉに対応する遅延時間（遅延ポイント数）ｎｉを検出し、複素乗算・可変遅延部２２にそれらＰｉおよびｎｉを出力するとともに、集合Ａ保持部２８には遅延時間（遅延ポイント数）ｎｉを出力する。集合Ａ保持部２８では最大振幅値・遅延時間検出部２６において一度検出された最大振幅値に対応する遅延時間（遅延ポイント数）ｎｉを保持し、常に最大振幅値・遅延時間検出部２６に集合Ａを出力する。ｓｌ（ｎ）出力部３０においては、ｓｉ（ｎ）が入力され、ｉ＝Ｅ・Ｎとなった時点で、ｓｉ（ｎ）すなわちｓＥ・Ｎ（ｎ）を出力する。こうして得られたｓＥ・Ｎ（ｎ）が、時間分解能を向上した遅延プロファイル（複素インパルス応答）ｓｌ（ｎ）となる。
【００１４】
ここで、本実施形態にかかる遅延プロファイル測定装置１０による遅延プロファイルの時間分解能向上の一例について図５〜図１０を参照して説明する。なお、各図において、横軸は遅延時間ｎ（正規化値：時間軸における離散ポイント）であり、縦軸は振幅レベル（デシベル［ｄＢ］）である。
【００１５】
前準備：前準備として、Ｅ・Ｎポイントの離散フーリエ変換の周波数領域において、低域部Ｎポイントの実数部をＥ、虚数部を零とし、残りの（Ｅ−１）・Ｎポイントの実数部、虚数部ともに零とした伝達関数を作成し、Ｅ・Ｎポイントの逆離散フーリエ変換を施し、得られた時間領域のＥ・Ｎポイントの時間領域の複素数列をｂ（ｎ）として予め補正関数記憶部２４に保持しておく。図５にＥ＝４、Ｎ＝６４の場合の｜ｂ（ｎ）｜を示す。
【００１６】
ステップ１：Ｅ・Ｎポイントの離散フーリエ変換の周波数領域において、低域部Ｎポイントに、元の時間分解能Ｔの遅延プロファイル（複素インパルス応答）ｓ（ｎ）の逆離散フーリエ変換である（離散）伝達関数Ｓ（ｋ）のＥ倍の複素数値を与え、残りの（Ｅ−１）・Ｎポイントの実数部、虚数部ともに零とし、（Ｅ・Ｎ）ポイントの逆離散フーリエ変換を施し、得られた（Ｅ・Ｎ）ポイントの時間領域の複素数列をｓ０（ｎ）とし、ｉ＝０とする。図６にＮ＝６４、遅延波が１波存在（振幅レベル＝（遅延波レベル／希望波レベル）＝−１０［ｄＢ］、遅延時間＝２０．２５Ｔ、位相＝０度）する場合の｜ｓ（ｎ）｜、図６に対応するＥ＝４、Ｎ＝６４の場合の｜ｓ０（ｎ）｜を図７に示す。
【００１７】
ステップ２：集合Ａに属する時間ポイントの値を除いた
【数１２】
2004023405.tif 000010
の振幅値
【数１３】
2004023405.tif 000011
の最大値｜Ｐｉ｜およびそのときの離散ポイント（最大値ポイント）ｎｉを検出し、｜Ｐｉ｜における複素数値
【数１４】
2004023405.tif 000012
およびＰｉに対応するポイント数ｎｉを保持し、ｎｉを集合Ａの要素に加える。
【００１８】
ステップ３：ｎ＝ｎｉ以外で、ｓｉ（ｎ）からＰｉ・ｂ（ｎ−ｎｉ）を減じ、得られた結果をｓｉ＋１（ｎ）とし、ｉ＝ｉ＋１とする。ここで、ｉ＝Ｅ・Ｎならばステップ４へ、そうでなければステップ２へ進む。
【数１５】
2004023405.tif 000013
図８にＥ＝４、Ｎ＝６４の場合の図６に対応する｜ｓ１（ｎ）｜を、図９にＥ＝４、Ｎ＝６４の場合の図６に対応する｜ｓ２（ｎ）｜を示す。
【００１９】
ステップ４：ｓｉ（ｎ）を時間分解能Ｔ／Ｅの遅延プロファイル（複素インパルス応答）ｓｌ（ｎ）とする。図１０にＥ＝４、Ｎ＝６４の場合の図６に対応する｜ｓｌ（ｎ）｜を示す。
【００２０】
図２は、本発明の第二の実施形態にかかる遅延プロファイル（複素インパルス応答）測定装置のブロック線図を示している。本実施形態の最大振幅値・遅延時間検出部４２は、最大振幅値｜Ｐｉ｜を出力する点で、第一の実施形態の最大振幅値・遅延時間検出部２６と異なる。この第二の実施形態では、検出された最大振幅値Ｐｉに基づいてデクリメント数ｊを決定している。具体的には、図２の例では、第一の実施形態のステップ２に、新たな終了条件：「最大振幅値｜Ｐｉ｜＜Ｌｔｈ（Ｌｔｈ：閾値，例えばノイズレベル）ならばステップ４へ。」を追加している。これにより、上記第一の実施形態において行われていた無駄な演算（例えばノイズレベルＬｔｈ以下のインパルスに対する演算）を回避でき、計算量が大幅に削減されるとともに、遅延プロファイル（複素インパルス応答）算出結果の精度が高くなる。
【００２１】
＜補正関数ｂ（ｎ）の簡略化＞　さらに、第一の実施形態、第二の実施形態において、補正関数ｂ（ｎ）の保持する値の数をＥ・Ｎポイントとせず、ｂ（ｎ）の最大値の前後の数ポイントの値（例えば遅延プロファイル上のノイズレベルＬｔｈより上回る値）のみを保持し、差し引く（図１２はノイズレベルＬｔｈ＝−３０［ｄＢ］とした場合のｂ（ｎ）の一例である。）ことで、ｂ（ｎ）を保持しておく記憶容量が削減され、ｓｉ（ｎ）からＰｉ・ｂ（ｎ−ｎｉ）を減算する計算量が削減されることに加え、ノイズレベルにおける無駄な演算を回避できる。なお、同じく図１２に示すように、保持する補正関数ｂ（ｎ）を、時間軸（の離散ポイント）により、例えばｂ（ｎ）の最大値（０［ｄＢ］）の前にＩｆポイント分、最大値の後にＩｂポイント分として決定しても、同様の効果が得られる（Ｉｆ、Ｉｂともに正の整数）。
【００２２】
さらに、上記第一の実施形態において、最大値検出、および補正処理演算の繰り返し回数を伝搬路状況に応じてポイント数Ｅ・Ｎよりも少ない一定回数とする（図１においてＥ・ＮをＪ＜Ｅ・Ｎを満たすＪとする（ステップ３においてｉ＝Ｊならばステップ４へ））ことで処理の簡略化を図ることも可能である。ただし、伝搬路状況によらず演算の繰り返し回数が一定となるため、厳密な意味での正確な遅延プロファイル（複素インパルス応答）ではなくなるが、繰り返し回数を伝搬路状況に合わせて誤差が小さくなるように決定すれば遅延プロファイル（複素インパルス応答）結果の精度を高くすることができる。
【００２３】
なお、元の遅延プロファイル（複素インパルス応答）ｓ（ｎ）に対応する伝達関数Ｓ（ｋ）が既知であれば、図１、２のＳ（ｋ）を求める手順（ＮポイントＤＦＴ部１２の処理）を省略し、図３、４に示すように、Ｓ（ｋ）から時間分解能が向上した遅延プロファイル（複素インパルス応答）ｓｌ（ｎ）の離散フーリエ変換ｓｌ（ｋ）を得ることも可能である。
【００２４】
また、本実施形態にかかる遅延プロファイル測定装置による処理の別の例として、遅延波の遅延時間（時間方向の離散ポイント数）を逓倍（整数倍）ではなく２０．２５Ｔとした場合に得られる遅延プロファイル（複素インパルス応答）の結果を図６（時間分解能Ｔの遅延プロファイル結果）に、実際の遅延プロファイル（複素インパルス応答）を図１１（時間分解能０．２５Ｔの遅延プロファイル結果）に、本発明により時間分解能Ｔの遅延プロファイル（複素インパルス応答）結果から遅延プロファイル（複素インパルス応答）の時間分解能向上法によって時間分解能が０．２５Ｔに向上した遅延プロファイル（複素インパルス応答）結果を図１０に示す。結果より、図６ではインパルス応答が広い遅延時間範囲に広がり、実際の遅延プロファイルである図１１と異なった形状となっているのに対し、図１０では図１１と同じ結果が得られている。このように、本発明により、時間分解能Ｔの遅延プロファイル（複素インパルス応答）から時間分解能Ｔ／４＝０．２５Ｔの遅延プロファイル（複素インパルス応答）が求まることがわかる。なお、図６および図１１では、遅延波が１波のみ存在する場合を示したが、遅延波および／または回り込み波が複数波存在したとしても同様の結果が得られる。
【００２５】
上述の方法は既に得られた遅延プロファイルの時間分解能を向上させ、かつ振幅値を正確に求める方法であるが、得られている伝達関数のポイント数が離散フーリエ変換のポイント数に達していないことから遅延プロファイル（複素インパルス応答）が正確に求まらない場合にも適用が可能である。例えば、処理時間の問題から高速離散フーリエ変換を用いなければならないにもかかわらず、得られる（離散）伝達関数のポイント数Ｍが高速離散フーリエ変換を適用する一般的なポイント数である２のべき乗でないとき、伝達関数のポイント数Ｍよりも大きい最小の２のべき乗ポイント数Ｎの高速フーリエ変換を使用し、伝達関数を周波数領域の低域部に配置し、残り（Ｎ−Ｍ）ポイントの値は零としてＮポイント逆高速離散フーリエ変換することで遅延プロファイル（複素インパルス応答）を得るが、このままだと遅延プロファイル（複素インパルス応答）として、精度の悪いものとなってしまう。従って、このような状況下においても、本発明による遅延プロファイル（複素インパルス応答）の時間分解能向上方法を用いると精度の高い遅延プロファイル（複素インパルス応答）が得られる。ただし、前準備において得る補正関数ｂ（ｎ）は、Ｎポイントの離散フーリエ変換の周波数領域において、低域部Ｍポイントの実数部をＮ／Ｍ、虚数部を零とし、残りの（Ｎ−Ｍ）ポイントの実数部、虚数部ともに零とした伝達関数を作成し、Ｎポイントの逆離散フーリエ変換を施し、得られた時間領域のＮポイントの時間領域の複素数列をｂ（ｎ）として補正関数記憶部２４に保持しておく。
【００２６】
次に、上記実施形態にかかる遅延プロファイル（複素インパルス応答）測定装置１１または４１（図３，図４）を、図１３に示す中継装置５０のフィードバック型回り込みキャンセラ６０ａ、および図１４に示す受信装置７０のフィードバック型伝送路適応等化器６０ｂのフィルタタップ係数値算出部５２（図１５）に導入した実施形態について説明する。
【００２７】
図１３に示す中継装置５０では、フィードバック型回り込みキャンセラ６０ａにより、受信アンテナ５１で受信された受信信号から回り込み波および遅延波成分が除去される。その信号が増幅部（送信装置を含む）５５で増幅され送信アンテナ５３から送出される。また、図１４に示す受信装置７０では、フィードバック型伝送路適応等化器６０ｂにより、受信アンテナ５１で受信された受信信号から遅延波成分が除去され、再送信信号が取得される。その再送信信号がＯＦＤＭ復調部７２に入力される。
【００２８】
フィードバック型回り込みキャンセラ６０ａ（図１３）およびフィードバック型伝送路適応等化器６０ｂ（図１４）は、いずれも、Ａ／Ｄ変換部６４、直交復調部６６、複素トランスバーサルフィルタ６２、直交変調部６８、Ｄ／Ａ変換部６９、加算器５７、およびフィルタタップ係数値算出部５２を含む。本実施形態では、複素トランスバーサルフィルタ６２によってキャンセル信号を生成し、加算器５７において受信信号からキャンセル信号を差し引くことで、受信信号に含まれる回り込み波および遅延波成分を除去する。キャンセル信号の特性は、複素トランスバーサルフィルタ６２のフィルタタップ係数値に応じたものとなるが、本実施形態ではこのフィルタタップ係数値の決定に際して上述した時間分解能向上手法を用いている点で従来とは異なる。
【００２９】
図１５に示すように、逆数化された伝達関数推定部５４において、ＯＦＤＭ時間信号（再送信信号）から逆数化された伝送路伝達関数Ｔ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部５６で、時間分解能がＴ／Ｅに向上したｔｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＴ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｔ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。ここで、逆数化された伝送路伝達関数Ｔ（ｋ）は、元の伝送路伝達関数Ｓ（ｋ）と
【数１６】
2004023405.tif 000014
なる関係にあり、Ａは補正係数である（この補正係数Ａは、例えばＳ（ｋ）の平均値としてもよいし、予め定めた所定値としてもよい）。以上の処理により得られたｔｌ（ｎ）を利用して、更新タップ係数値算出部５８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値（複素タップ係数値）ｈ（ｎ）（０≦ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数１７】
2004023405.tif 000015
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図１５に示すフィルタタップ係数値算出部５２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）６２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）６２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。
【００３０】
次に、本発明にかかる上記遅延プロファイル（複素インパルス応答）測定装置１１または４１（図３，図４）を、図１６に示す中継装置８０のフィードバック型回り込みキャンセラ９０ａ、および図１７に示す受信装置１００のフィードバック型伝送路適応等化器９０ｂのフィルタタップ係数値算出部８２（図１８）に導入した実施形態について説明する。
【００３１】
図１６に示す中継装置８０では、フィードバック型回り込みキャンセラ９０ａにより、受信アンテナ５１で受信された受信信号から回り込み波および遅延波成分が除去される。その信号が増幅器（送信装置を含む）５５で増幅され送信アンテナ５３から送出される。また、図１７に示す受信装置１００では、フィードバック型伝送路適応等化器９０ｂにより、受信アンテナ５１で受信された受信信号から遅延波成分が除去され、再送信信号が取得される。その再送信信号がＯＦＤＭ復調部７２に入力される。
【００３２】
フィードバック型回り込みキャンセラ９０ａ（図１６）およびフィードバック型伝送路適応等化器９０ｂ（図１７）は、いずれも、Ａ／Ｄ変換部８１、加算器８３、Ｄ／Ａ変換部８５、直交復調部９４、複素トランスバーサルフィルタ９２、直交変調部９６、およびフィルタタップ係数値算出部８２を含む。本実施形態では、複素トランスバーサルフィルタ９２によってキャンセル信号を生成し、加算器８３においてＡ／Ｄ変換された受信信号からキャンセル信号を差し引くことで、受信信号に含まれる回り込み波および遅延波成分を除去する。キャンセル信号の特性は、複素トランスバーサルフィルタ９２のフィルタタップ係数値に応じたものとなるが、本実施形態ではこのフィルタタップ係数値の決定に際して上述した時間分解能向上手法を用いている点で従来とは異なる。
【００３３】
図１８に示すように、逆数化された伝達関数推定部８４において、Ａ／Ｄ変換後のＯＦＤＭ時間信号（再送信信号または等化処理後の受信信号）から逆数化された伝送路伝達関数Ｔ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部８６で、時間分解能がＴ／Ｅに向上したｔｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＴ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｔ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。ここで、逆数化された伝送路伝達関数Ｔ（ｋ）は、元の伝送路伝達関数Ｓ（ｋ）と
【数１８】
2004023405.tif 000016
なる関係にあり、Ａは補正係数である（この補正係数Ａは、例えばＳ（ｋ）の平均値としてもよいし、予め定めた所定値としてもよい）。以上の処理により得られたｔｌ（ｎ）を利用して、更新タップ係数値算出部８８において更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０≦ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数１９】
2004023405.tif 000017
によって算出する。以上のような方法により、本発明にかかる上記遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図１８に示すフィルタタップ係数値算出部８２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）９２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）９２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。
【００３４】
次に、本発明にかかる上記遅延プロファイル（複素インパルス応答）測定装置１１または４１（図３，図４）を、図１９に示す中継装置１１０のフィードフォワード型回り込みキャンセラ１２０ａ、および図２０に示す受信装置１３０のフィードフォワード型伝送路適応等化器１２０ｂのフィルタタップ係数値算出部１１２（図２１）に導入した実施形態について説明する。
【００３５】
図１９に示す中継装置１１０では、フィードフォワード型回り込みキャンセラ１２０ａにより、受信アンテナ５１で受信された受信信号から回り込み波および遅延波成分が除去される。その信号が増幅器（送信装置を含む）５５で増幅され送信アンテナ５３から送出される。また、図２０に示す受信装置１３０では、フィードフォワード型伝送路適応等化器１２０ｂにより、受信アンテナ５１で受信された受信信号から遅延波成分が除去され、再送信信号が取得される。その再送信信号がＯＦＤＭ復調部７２に入力される。
【００３６】
フィードフォワード型回り込みキャンセラ１２０ａ（図１９）およびフィードフォワード型伝送路適応等化器１２０ｂ（図２０）は、いずれも、Ａ／Ｄ変換部１２９、直交復調部１２８、直交変調部１２６、Ｄ／Ａ変換部１２４、加算器１２１、複素トランスバーサルフィルタ１２２、およびフィルタタップ係数値算出部１１２を含む。本実施形態では、複素トランスバーサルフィルタ１２２によってキャンセル信号を生成し、加算器１２１において受信信号からキャンセル信号を差し引くことで、受信信号に含まれる回り込み波および遅延波成分を除去する。キャンセル信号の特性は、複素トランスバーサルフィルタ１２２のフィルタタップ係数値に応じたものとなるが、本実施形態ではこのフィルタタップ係数値の決定に際して上述した時間分解能向上手法を用いている点で従来と異なる。
【００３７】
図２１に示すように、伝達関数推定部１１４において、ＯＦＤＭ時間信号（再送信信号または等化処理後の受信信号）から逆数化された伝送路伝達関数Ｓ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部１１６で、時間分解能がＴ／Ｅに向上したｓｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＳ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｓ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。このｓｌ（ｎ）を利用して、更新タップ係数値算出部１１８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０≦ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数２０】
2004023405.tif 000018
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図２１に示すフィルタタップ係数値算出部１１２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１２２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１２２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。
【００３８】
次に、本発明にかかる上記遅延プロファイル（複素インパルス応答）測定装置１１または４１（図３，図４）を、図２２に示す中継装置１４０のフィードフォワード型回り込みキャンセラ１５０ａ、および図２３に示す受信装置１６０のフィードフォワード型伝送路適応等化器１５０ｂのフィルタタップ係数値算出部１４２（図２４）に導入した実施形態について説明する。
【００３９】
図２２に示す中継装置１４０では、フィードフォワード型回り込みキャンセラ１５０ａにより、受信アンテナ５１で受信された受信信号から回り込み波および遅延波成分が除去される。その信号が増幅器（送信装置を含む）５５で増幅され送信アンテナ５３から送出される。また、図２３に示す受信装置１６０では、フィードフォワード型伝送路適応等化器１５０ｂにより、受信アンテナ５１で受信された受信信号から遅延波成分が除去され、再送信信号が取得される。その再送信信号がＯＦＤＭ復調部７２に入力される。
【００４０】
フィードフォワード型回り込みキャンセラ１５０ａ（図２２）およびフィードフォワード型伝送路適応等化器１５０ｂ（図２３）は、いずれも、Ａ／Ｄ変換部８１、Ｄ／Ａ変換部８５、直交復調部１５６、直交変調部１５４、加算器１５１、複素トランスバーサルフィルタ１５２、およびフィルタタップ係数値算出部１４２を含む。本実施形態では、複素トランスバーサルフィルタ１５２によってキャンセル信号を生成し、加算器１５１において受信信号からキャンセル信号を差し引くことで、受信信号に含まれる回り込み波および遅延波成分を除去する。キャンセル信号の特性は、複素トランスバーサルフィルタ１５２のフィルタタップ係数値に応じたものとなるが、本実施形態ではこのフィルタタップ係数値の決定に際して上述した時間分解能向上手法を用いている点で従来とは異なる。
【００４１】
図２４に示すように、伝達関数推定部１４４において、Ａ／Ｄ変換された後のＯＦＤＭ時間信号（再送信信号、等化処理後の受信信号）から伝送路伝達関数Ｓ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、次いで時間分解能向上複素インパルス応答（遅延プロファイル）推定部１４６で、時間分解能がＴ／Ｅに向上したｓｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＳ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｓ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。このｓｌ（ｎ）を利用して、更新タップ係数値算出部１４８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０≦ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、次の式
【数２１】
2004023405.tif 000019
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）時間分解能向上方法を図２４に示すフィルタタップ係数値算出部１４２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１５２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善されて小さくなる。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１５２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。
【００４２】
図１〜図４、および図１３〜図２４までを対象とする信号をＯＦＤＭ信号として示したが、遅延プロファイル（複素インパルス応答）または対応する伝達関数を利用した遅延プロファイル（複素インパルス応答）測定装置、回り込みキャンセラ、および伝送路適応等化器であれば全て適用可能な構成であり、信号の形式はＯＦＤＭ信号には限定されない。
【００４３】
図１３から図２４まで、回り込みキャンセラ、伝送路適応等化器について対象とする信号をＯＦＤＭ信号として示したが、遅延プロファイル（複素インパルス応答）または対応する伝達関数を利用した回り込みキャンセラ、および伝送路適応等化器であれば全て適用可能な構成であり、信号の形式はＯＦＤＭ信号に限らない。
【００４４】
【発明の効果】
本発明によれば、時間分解能Ｔの遅延プロファイル（複素インパルス応答）から時間分解能Ｔ／Ｅに向上させた遅延プロファイル（複素インパルス応答）が得られる。また、もともとの遅延プロファイル（複素インパルス応答）上におけるノイズレベル以下の振幅値を持つインパルスに対して補正処理を行わないことにより、処理量、および計算量を削減しつつ、同等の効果を得られる。また、補正処理に使用する補正関数値のポイント数を少なくすることで、記憶容量、および計算量を削減しつつ、同等の効果を得られる。
【００４５】
また、ノイズレベル以下の振幅値を持つインパルスに対して補正処理を行わず、かつ補正処理に使用する補正関数のポイント数を少なくすることを併用することで、記憶容量、および計算量を削減しつつ、同等の効果が得られる。さらに、本発明による時間分解能を向上させた遅延プロファイル（複素インパルス応答）を伝送路等化器に用いると、伝送路等化精度（不要波キャンセル精度）が高く、等化可能な希望波電力対非希望波電力比レベルの改善された（小さくなった）回り込みキャンセラ、および伝送路等化器を実現できる。
【図面の簡単な説明】
【図１】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置のブロック線図である。
【図２】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置のブロック線図である。
【図３】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置のブロック線図である。
【図４】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置のブロック線図である。
【図５】Ｅ＝４、Ｎ＝６４の場合の補正関数ｂ（ｎ）の振幅値｜ｂ（ｎ）｜を示す図である。
【図６】Ｎ＝６４、遅延波が１波存在（（希望波電力／非希望波電力）＝１０［ｄＢ］、遅延時間＝２０．２５Ｔ、位相＝０［ｄｅｇ］）する場合に時間分解能Ｔで得られるｓ（ｎ）の振幅値｜ｓ（ｎ）｜を示す図である。
【図７】図６のｓ（ｎ）に対応するＥ＝４、Ｎ＝６４の場合のｓ０（ｎ）の振幅値｜ｓ０（ｎ）｜を示す図である。
【図８】図６のｓ（ｎ）に対応するＥ＝４、Ｎ＝６４の場合のｓ１（ｎ）の振幅値｜ｓ１（ｎ）｜を示す図である。
【図９】図６のｓ（ｎ）に対応するＥ＝４、Ｎ＝６４の場合のｓ２（ｎ）の振幅値｜ｓ２（ｎ）｜を示す図である。
【図１０】図６のｓ（ｎ）に対応するＥ＝４、Ｎ＝６４の場合のｓｌ（ｎ）の振幅値｜ｓｌ（ｎ）｜を示す図である。
【図１１】Ｎ＝６４、遅延波が１波存在（（希望波電力／非希望波電力）＝１０［ｄＢ］、遅延時間＝２０．２５Ｔ、位相＝０［ｄｅｇ］）する場合に時間分解能Ｔ／４＝０．２５で得られる遅延プロファイルの振幅値を示す図である。
【図１２】Ｌｔｈ＝−３０［ｄＢ］としたときの補正関数ｂ（ｎ）を示す図である。
【図１３】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードバック型回り込みキャンセラを用いた中継装置を示すブロック線図である。
【図１４】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードバック型伝送路適応等化器を用いた受信装置を示すブロック線図である。
【図１５】図１３のフィードバック型回り込みキャンセラ、および図１４のフィードバック型伝送路適応等化器のフィルタタップ係数値算出部を示すブロック線図である。
【図１６】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードバック型回り込みキャンセラを用いた中継装置を示すブロック線図である。
【図１７】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードバック型伝送路適応等化器を用いた受信装置を示すブロック線図である。
【図１８】図１６のフィードバック型回り込みキャンセラ、および図１７のフィードバック型伝送路適応等化器のフィルタタップ係数値算出部を示すブロック線図である。
【図１９】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードフォワード型回り込みキャンセラを用いた中継装置を示すブロック線図である。
【図２０】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードフォワード型伝送路適応等化器を用いた受信装置を示すブロック線図である。
【図２１】図１９のフィードフォワード型回り込みキャンセラ、および図２０のフィードフォワード型伝送路適応等化器のフィルタタップ係数値算出部を用いた受信装置を示すブロック線図である。
【図２２】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードフォワード型回り込みキャンセラを用いた中継装置を示すブロック線図である。
【図２３】本発明の実施形態にかかる時間分解能を向上した遅延プロファイル（複素インパルス応答）測定装置を用いたフィードフォワード型伝送路適応等化器を用いた受信装置を示すブロック線図である。
【図２４】図２２のフィードフォワード型回り込みキャンセラ、および図２３のフィードフォワード型伝送路適応等化器のフィルタタップ係数値算出部を示すブロック線図である。
【符号の説明】
１０，１１，４０，４１　遅延プロファイル測定装置
１２　ＮポイントＤＦＴ部
１４　Ｎ→Ｅ・Ｎポイント数増加部（離散ポイント数増加部）
１６　Ｅ・ＮポイントＩＤＦＴ部（逆離散フーリエ変換部）
１８　ｓｉ（ｎ）保持部
２０　減算部（補正処理部）
２２　複素乗算・可変遅延部
２４　補正関数記憶部
２６，４２　最大振幅値・遅延時間検出部
２８　集合Ａ保持部
３０，４４　ｓｌ（ｎ）出力部
５０，８０，１１０，１４０　中継装置
５２，８２，１１２，１４２　フィルタタップ係数値算出部
５４，８４　逆数化された伝達関数推定部
５６，８６，１１６，１４６　時間分解能向上複素インパルス応答（遅延プロファイル）推定部
５８，８８，１１８，１４８　更新タップ係数値算出部
６０ａ，９０ａ　フィードバック型回り込みキャンセラ
６０ｂ，９０ｂ　フィードバック型伝送路適応等化器
６２，９２，１２２，１５２　複素トランスバーサルフィルタ
７０，１００，１３０，１６０　受信装置
１１４，１４４　伝達関数推定部
１２０ａ，１５０ａ　フィードフォワード型回り込みキャンセラ
１２０ｂ，１５０ｂ　フィードフォワード型伝送路適応等化器。
【図１】
2004023405.tif 000024
【図２】
2004023405.tif 000025
【図３】
2004023405.tif 000026
【図４】
2004023405.tif 000027
【図５】
2004023405.tif 000028
【図６】
2004023405.tif 000029
【図７】
2004023405.tif 000030
【図８】
2004023405.tif 000031
【図９】
2004023405.tif 000032
【図１０】
2004023405.tif 000033
【図１１】
2004023405.tif 000034
【図１２】
2004023405.tif 000035
【図１３】
2004023405.tif 000036
【図１４】
2004023405.tif 000037
【図１５】
2004023405.tif 000038
【図１６】
2004023405.tif 000039
【図１７】
2004023405.tif 000040
【図１８】
2004023405.tif 000041
【図１９】
2004023405.tif 000042
【図２０】
2004023405.tif 000043
【図２１】
2004023405.tif 000044
【図２２】
2004023405.tif 000045
【図２３】
2004023405.tif 000046
【図２４】
2004023405.tif 000047
【手続補正書】
【提出日】平成15年4月8日(2003.4.8)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】請求項１
【補正方法】変更
【補正の内容】
【請求項１】
ＯＦＤＭ信号をＮポイント（Ｎ：整数）で離散フーリエ変換して取得された離散伝達関数（複素インパルス応答）Ｓ（ｋ）（ｋ：離散ポイント（整数）、０≦ｋ≦Ｎ−１）から、次の式（１）
【数１】
2004023405.tif 000048
（ここに、Ｅ：１より大きい有理数）により離散ポイント数をＥ・Ｎに増加した離散伝達関数Ｓ０（ｋ）を取得する離散ポイント数増加部と、
前記離散伝達関数Ｓ０（ｋ）にＥ・Ｎポイントで逆離散フーリエ変換を施し、Ｅ・Ｎポイントの複素数列ｓ０（ｎ）を取得する逆離散フーリエ変換部と、
各補正ステップｉ（０≦ｉ≦Ｊ−１，ｉ，Ｊともに整数）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント（整数），０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、および該最大値｜Ｐｉ｜をとる離散ポイントｎとしての最大値ポイントｎｉを検出する検出部と（ただし、ｎ≠ｎｉｐｒｅｖ，ｎｉｐｒｅｖ：現補正ステップｉより前の補正ステップｉｐｒｅｖ（０≦ｉｐｒｅｖ＜ｉ）での最大値ポイント）、
各補正ステップｉについて、前記最大値｜Ｐｉ｜および最大値ポイントｎｉに基づいて次の式（２）
【数２】
2004023405.tif 000049
（ここに、ｂ（ｎ）：補正関数）により複素数列ｓｉ＋１（ｎ）を取得する補正処理部と、を備え、補正ステップｉ（整数、０≦ｉ≦Ｊ−１）をＪ回デクリメントして取得された複素数列ｓＪ（ｎ）を、ＯＦＤＭ信号の遅延プロファイル（複素インパルス応答）として取得する遅延プロファイル測定装置。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】請求項６
【補正方法】変更
【補正の内容】
【請求項６】
ＯＦＤＭ信号をＮポイント（Ｎ：整数）で離散フーリエ変換して取得された離散伝達関数（複素インパルス応答）Ｓ（ｋ）（ｋ：離散ポイント（整数）、０≦ｋ≦Ｎ−１）から、次の式（１）
【数３】
2004023405.tif 000050
（ここに、Ｅ：１より大きい有理数）により離散ポイント数をＥ・Ｎに増加した離散伝達関数Ｓ０（ｋ）を取得するステップと、
前記離散伝達関数Ｓ０（ｋ）にＥ・Ｎポイントで逆離散フーリエ変換を施し、Ｅ・Ｎポイントの複素数列ｓ０（ｎ）を取得するステップと、
各補正ステップｉ（０≦ｉ≦Ｊ−１，ｉ，Ｊともに整数）について、複素数列ｓｉ（ｎ）（ｎ：離散ポイント（整数），０≦ｎ≦Ｅ・Ｎ−１）の振幅値の最大値｜Ｐｉ｜、および該最大値｜Ｐｉ｜をとる離散ポイントｎとしての最大値ポイントｎｉを検出するステップと（ただし、ｎ≠ｎｉｐｒｅｖ，ｎｉｐｒｅｖ：現補正ステップｉより前の補正ステップｉｐｒｅｖ（０≦ｉｐｒｅｖ＜ｉ）での最大値ポイント）、
各補正ステップｉについて、前記最大値｜Ｐｉ｜および最大値ポイントｎｉに基づいて次の式（２）
【数４】
2004023405.tif 000051
（ここに、ｂ（ｎ）：補正関数）により複素数列ｓｉ＋１（ｎ）を取得するステップと、
を含み、補正ステップｉ（整数、０≦ｉ≦Ｊ−１）をＪ回デクリメントして取得された複素数列ｓＪ（ｎ）を、ＯＦＤＭ信号の遅延プロファイル（複素インパルス応答）として取得する遅延プロファイル測定方法。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
ステップ１：Ｅ・Ｎポイントの離散フーリエ変換の周波数領域において、低域部Ｎポイントに、元の時間分解能Ｔの遅延プロファイル（複素インパルス応答）ｓ（ｎ）の離散フーリエ変換である（離散）伝達関数Ｓ（ｋ）のＥ倍の複素数値を与え、残りの（Ｅ−１）・Ｎポイントの実数部、虚数部ともに零とし、（Ｅ・Ｎ）ポイントの逆離散フーリエ変換を施し、得られた（Ｅ・Ｎ）ポイントの時間領域の複素数列をｓ０（ｎ）とし、ｉ＝０とする。図６にＮ＝６４、遅延波が１波存在（振幅レベル＝（遅延波レベル／希望波レベル）＝−１０［ｄＢ］、遅延時間＝２０．２５Ｔ、位相＝０度）する場合の｜ｓ（ｎ）｜、図６に対応するＥ＝４、Ｎ＝６４の場合の｜ｓ０（ｎ）｜を図７に示す。


【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００２３
【補正方法】変更
【補正の内容】
【００２３】
なお、元の遅延プロファイル（複素インパルス応答）ｓ（ｎ）に対応する伝達関数Ｓ（ｋ）が既知であれば、図１、２のＳ（ｋ）を求める手順（ＮポイントＤＦＴ部１２の処理）を省略し、図３、４に示すように、Ｓ（ｋ）から時間分解能が向上した遅延プロファイル（複素インパルス応答）ｓｌ（ｎ）の離散フーリエ変換Ｓｌ（ｋ）を得ることも可能である。


【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００２９
【補正方法】変更
【補正の内容】
【００２９】
図１５に示すように、逆数化された伝達関数推定部５４において、ＯＦＤＭ時間信号（再送信信号）から逆数化された伝送路伝達関数Ｔ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部５６で、時間分解能がＴ／Ｅに向上したｔｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＴ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｔ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。ここで、逆数化された伝送路伝達関数Ｔ（ｋ）は、元の伝送路伝達関数Ｓ（ｋ）と
【数１６】
2004023405.tif 000052
なる関係にあり、Ａは補正係数である（この補正係数Ａは、例えばＳ（ｋ）の平均値としてもよいし、予め定めた所定値としてもよい）。以上の処理により得られたｔｌ（ｎ）を利用して、更新タップ係数値算出部５８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値（複素タップ係数値）ｈ（ｎ）（０＜ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数１７】
2004023405.tif 000053
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図１５に示すフィルタタップ係数値算出部５２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）６２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）６２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。


【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００３３
【補正方法】変更
【補正の内容】
【００３３】
図１８に示すように、逆数化された伝達関数推定部８４において、Ａ／Ｄ変換後のＯＦＤＭ時間信号（再送信信号または等化処理後の受信信号）から逆数化された伝送路伝達関数Ｔ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部８６で、時間分解能がＴ／Ｅに向上したｔｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＴ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｔ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。ここで、逆数化された伝送路伝達関数Ｔ（ｋ）は、元の伝送路伝達関数Ｓ（ｋ）と
【数１８】
2004023405.tif 000054
なる関係にあり、Ａは補正係数である（この補正係数Ａは、例えばＳ（ｋ）の平均値としてもよいし、予め定めた所定値としてもよい）。以上の処理により得られたｔｌ（ｎ）を利用して、更新タップ係数値算出部８８において更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０＜ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数１９】
2004023405.tif 000055
によって算出する。以上のような方法により、本発明にかかる上記遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図１８に示すフィルタタップ係数値算出部８２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）９２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）９２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。


【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００３７
【補正方法】変更
【補正の内容】
【００３７】
図２１に示すように、伝達関数推定部１１４において、ＯＦＤＭ時間信号（再送信信号または等化処理後の受信信号）から逆数化された伝送路伝達関数Ｓ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、さらに時間分解能向上複素インパルス応答（遅延プロファイル）推定部１１６で、時間分解能がＴ／Ｅに向上したｓｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＳ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｓ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。このｓｌ（ｎ）を利用して、更新タップ係数値算出部１１８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０＜ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、
【数２０】
2004023405.tif 000056
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）の時間分解能向上方法を図２１に示すフィルタタップ係数値算出部１１２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１２２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善される（小さくなる）。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１２２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。


【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００４１
【補正方法】変更
【補正の内容】
【００４１】
図２４に示すように、伝達関数推定部１４４において、Ａ／Ｄ変換された後のＯＦＤＭ時間信号（再送信信号、等化処理後の受信信号）から伝送路伝達関数Ｓ（ｋ）（ただし、０≦ｋ≦Ｎ−１）が得られ、次いで時間分解能向上複素インパルス応答（遅延プロファイル）推定部１４６で、時間分解能がＴ／Ｅに向上したｓｌ（ｎ）（０≦ｎ≦Ｅ・Ｎ−１）が得られる（単純にＳ（ｋ）を逆離散フーリエ変換して得られる複素インパルス応答ｓ（ｎ）（０≦ｎ≦Ｎ−１）の時間分解能はＴである。）。このｓｌ（ｎ）を利用して、更新タップ係数値算出部１４８において、更新係数等を用いた各種更新アルゴリズムを用いてフィルタタップ係数値ｈ（ｎ）（０＜ｎ≦Ｌ、Ｌ≦Ｅ・Ｎ−１）を算出する。ここでは、例えば、更新前のフィルタタップ係数をｈｏｌｄ（ｎ）、更新係数をλ（０＜λ≦１）としたとき、次の式
【数２１】
2004023405.tif 000057
によって算出する。以上のような方法により、本発明にかかる遅延プロファイル（複素インパルス応答）時間分解能向上方法を図２４に示すフィルタタップ係数値算出部１４２に用いると、より精度の高い複素インパルス応答を、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１５２のフィルタタップ係数値として与えることになるため、キャンセル可能な回り込み波、および遅延波等の希望波対非希望波電力比レベルが改善されて小さくなる。ただし、時間分解能をＴ／Ｅに向上させた遅延プロファイル（複素インパルス応答）をフィルタタップ係数として用いるために、複素トランスバーサルフィルタ（ＦＩＲフィルタ）１５２の遅延器一つの遅延時間もＴ／Ｅに変更する必要がある。