(57)【特許請求の範囲】
【請求項１】
時間領域信号である入力信号を周波数領域信号に変換して周波数領域で歪を検出する歪検出装置であって、
所定の時間長のｎ（ｎは２以上の整数）倍の周期を有する時間窓の少なくとも一部の波形を記憶する窓関数メモリを有し、前記記憶された波形を有する前記時間窓を生成し、前記所定の時間長を一単位とした前記入力信号に前記時間窓を乗算して、前記所定の時間長のｎ倍の時間長を有する時間領域信号を出力する時間窓処理部と、
前記時間窓処理部の出力信号に含まれるｎ個の前記所定の時間長を有する時間領域信号を前記一単位毎に全て加算して、前記所定の時間長を有する平均化された一単位の時間領域信号を出力する平均化処理部と、
前記平均化処理部から出力される前記一単位の時間領域信号を、所定数のポイントにおいて１度に時間領域信号から周波数領域信号に変換する時間／周波数変換処理部と、
前記時間／周波数変換処理部の出力から歪成分を抽出する歪抽出部と、
を備える歪検出装置。
【請求項２】
前記窓関数メモリは、複数の時間窓の少なくとも一部の波形を記憶し、前記時間窓処理部は、前記複数の時間窓を切り替えて前記乗算する時間窓に用いることを特徴とする
請求項１記載の歪検出装置。
【請求項３】
前記窓関数メモリは、複数の時間窓の少なくとも一部の波形を記憶し、前記時間窓処理部は、前記複数の時間窓を合成した時間窓を前記乗算する時間窓に用いることを特徴とする
請求項１記載の歪検出装置。
【請求項４】
前記窓関数メモリは、前記時間窓の一部の波形を記憶し、前記時間窓処理部は、前記時間窓の一部の波形に基づいて、所定時刻に対して対称な波形を有する前記時間窓を生成することを特徴とする
請求項１記載の歪検出装置。
【請求項５】
前記時間／周波数変換処理部は、ＦＦＴ（Fast Fourier Transform）演算により前記平均化処理部の出力を時間領域から周波数領域に変換し、
前記所定の時間長は、前記時間／周波数変換処理部におけるＦＦＴ演算で変換処理を行う単位である前記所定数のポイントに対応する時間長であることを特徴とする請求項１乃至４のいずれかに記載の歪検出装置。
【請求項６】
前記平均化処理部は、前記所定の時間長の時間領域信号を記憶するメモリを有し、前記時間窓が乗算された前記所定の時間長のｎ倍の長さを有する時間領域信号について、前記所定の長さの最初の一単位の時間領域信号を前記メモリに記憶し、さらに前記メモリに記憶された時間領域信号に、次の一単位の時間領域信号を加算して前記メモリに記憶することをｎ−１回繰り返し、前記メモリから前記所定の時間長を有する一単位の時間領域信号を出力することを特徴とする請求項１乃至４のいずれかに記載の歪検出装置。
【請求項７】
外部から入力された信号を増幅し、歪を補償する歪補償増幅装置であって、
第１の信号の増幅を行う増幅部と、
前記増幅部の出力を分配する分配部と、
前記分配部の一方の出力に含まれる歪成分を抽出して出力する歪検出装置と、
前記外部から入力された信号に対して、前記歪検出装置の出力を用いて適応アルゴリズムにより更新された逆特性に基づいて予歪を与えたものを前記第１の信号とする歪補償部とを備え、
前記歪検出装置は請求項１乃至６のいずれかに記載の歪検出装置であることを特徴とする歪補償増幅装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、増幅時の歪を検出し、補償するためのスペクトル解析方法、歪検出装置、歪補償増幅装置に関するものである。
【背景技術】
【０００２】
無線通信において、送信信号の増幅のために増幅器が必要である。増幅器は、入力レベルに対して非線形な出力レベル、出力位相となる非線形歪を発生する。非線形歪を表す入力対出力特性は、入力レベルに対する出力レベルを示すＡＭ−ＡＭ特性と、入力レベルに対する出力位相を示すＡＭ−ＰＭ特性を持つ。しかし、電力効率の向上のために、敢えて非線形歪が発生する領域で増幅器が使用されている。非線形歪は一般に、送信信号の帯域幅の３倍から５倍（悪いときには７倍）の帯域にわたって不要な電力放出を起こす。図１１は、歪を含む送信信号の一例を示すスペクトルである。横軸はキャリア周波数に対する相対周波数を表し、縦軸は希望波に対する相対レベルを表す。希望波のほかに、−ＩＭ７（負側の７次高調波），−ＩＭ５（負側の５次高調波），−ＩＭ３（負側の５次高調波帯域），＋ＩＭ３（正側の３次高調波），＋ＩＭ５（正側の５次高調波），＋ＩＭ７（正側の７次高調波）を示す。【０００３】
従って、非線形歪を補償する歪補償増幅装置が必要となる。非線形歪を補償する信号処理の方式として、ＦＦ（Feed Forward）方式や、ＡＰＤ（Adaptive Pre-Distortion）方式があり、いずれも商用レベルや試作レベルで実施され、有効とされている。ＦＦ方式は、増幅器出力の歪だけを抜き出し、増幅器出力から差し引くものである。ＡＰＤ方式は、増幅器入力前に予め歪特性を与え、増幅器と相殺するものである。ＡＰＤ方式は、デジタル信号処理で全てを実現可能であるため、小型化、低コスト化において有望である。また、ＡＰＤ方式においては、温度特性や経年変化等、増幅器の特性変化に追随するための適応アルゴリズムが使用されている。【０００４】
図１２は、従来の歪補償増幅装置の構成の一例を示すブロック図である。この歪補償増幅装置は、増幅部１、分配部２、歪検出装置１０３、補償係数演算部１０４、包絡線検出部１１、逆特性記憶部１２、遅延部１３、ＰＤ（予歪）部１４を備える。【０００５】
適応アルゴリズムには複数の種類があるが、例えば歪信号をモニタしながら、歪補償係数を変更する方法がある（例えば特許文献１参照）。逆特性記憶部１２は、増幅部１が持つＡＭ−ＡＭ特性とＡＭ−ＰＭ特性の逆特性を予め格納している。包絡線検出部１１は、送信信号のエンベロープを逆特性記憶部１２のアドレスとして出力する。逆特性記憶部１２は、このアドレスに格納された逆特性をＰＤ部１４に出力する。ＰＤ部１４は、遅延部１３で遅延された送信信号に対して逆特性を与える。増幅部１は、ＰＤ部１４の出力を増幅する。分配部２は、増幅部２の出力を分配し、一方を外部へ出力し、もう一方を歪検出装置１０３へ出力する。補償係数演算部１０４は、歪検出装置１０３により検出される歪が最小となるように、適応アルゴリズムを用いて逆特性記憶部１２の内容を更新する。【０００６】
なお、本発明の関連ある従来技術として、例えば、下記に示す特許文献２が知られている。この脈波検出装置は、計算能力が低い演算手段を用いた場合でも、正確な脈拍数の測定を可能にするものである。
【特許文献１】特開２００４−１２８９２１号公報
【特許文献２】特開２００４−１２１６２５号公報  （第４−７頁、第１図）
【特許文献３】特開平９−１３８２５１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、上述した歪補償装置は、ハードウェアのコストが掛かる、または適応アルゴリズムの収束時間があまりにも遅くなる、といった問題がある。【０００８】
ハードウェアのコストが掛かる方法の例として、同時に複数の歪帯域をモニタする方法がある。図１１のような周波数特性においては、少なくとも−ＩＭ５，−ＩＭ３，＋ＩＭ３，＋ＩＭ５の４箇所の歪帯域をモニタする必要がある。また、このように４箇所の歪帯域をアナログ信号においてモニタするには、４系統のアナログ受信系が必要となる。しかも、４系統のアナログ受信系の間にレベル偏差があってはならないため、コスト面で困難となる。【０００９】
一方、ハードウェアのコストを抑えた方法の例として、ＦＦＴ（Fast Fourier Transform）を用いた周波数領域での歪検出の方法がある。この方法において、受信系は１系統で良い。高速Ａ／Ｄ変換器で増幅器出力をデジタル信号に変換し、ＦＦＴで時間／周波数変換することにより、スペクトルが得られる。【００１０】
しかし、近年の送信信号は、２０ＭＨｚ程度の帯域が必要となっている。これは、例えばＷ−ＣＤＭＡ（Wideband-Code Division Multiple Access）において、１キャリアあたり５ＭＨｚの帯域で４キャリアを想定した場合である。図１１における送信信号はＷ−ＣＤＭＡの２キャリアのものであるが、この図から想像されるように、４キャリアで非線形歪が発生すると、１００ＭＨｚの帯域幅が必要となる。サンプリング定理により、２００ＭＨｚ以上のサンプリング動作を行わなければ、正確な歪検出を行うことができない。【００１１】
このような条件において、ＦＦＴをハードウェアで実現しようとすると回路規模が膨大であることから、コスト面で困難となる。また、ソフトウェアで実現しようとすると、処理に時間が掛かる。例えば、３００ＭＨｚで動作するＤＳＰ（Digital Signal Processor）を用いて、１０２４ポイント（サンプル）のＦＦＴを実施すると、約４ｍｓｅｃ必要となる。また、この場合、ＦＦＴを用いて歪検出できるのは約４ｍｓｅｃの時間領域信号のうち１０２４サンプルとわずかな時間だけであり、平均的なスペクトルを得るためには長時間の平均が必要となる。また、この歪補償増幅装置は、補償係数が１度更新される度に、その結果が正しいか否かを判断するために長時間平均を行うため、収束時間が２０分程度必要となることもある。一般に、収束時間は１０秒程度とすることが要求されていることから、１２０倍以上の高速化を図る必要がある。このように、ソフトウェアで実現する場合、時間／周波数変換処理中にリアルタイムで入力される信号は時間／周波数変換処理できないため、補償係数の収束時間が多大となってしまう。【００１２】
本発明は上述した問題点を解決するためになされたものであり、時間／周波数変換処理を高速化し、収束時間を低減するスペクトル解析方法、歪検出装置、歪補償増幅装置を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上述した課題を解決するため、本発明は、スペクトルの時間平均を計算するスペクトル解析方法であって、入力信号にＦＦＴポイント数の整数倍の周期の窓関数を乗算する第１ステップと、前記第１ステップの出力信号をＦＦＴポイント数間隔で平均して、ＦＦＴポイント数の信号を出力する第２ステップと、前記第２ステップの出力信号をＦＦＴする第３ステップとを有するものである。【００１４】
また、本発明は、入力信号について周波数領域で歪を検出する歪検出装置であって、前記入力信号を平均化する平均化処理部と、前記平均化処理部の出力を時間領域から周波数領域に変換する時間／周波数変換処理部と、前記時間／周波数変換処理部の出力から歪成分を抽出する歪抽出部とを備えたものである。【００１５】
また、本発明に係る歪検出装置において、さらに、前記時間領域信号に時間窓を乗算し、前記平均化処理部へ出力する時間窓処理部を備えることを特徴とするものである。【００１６】
また、本発明は、入力信号について周波数領域で歪を検出する歪検出装置であって、前記入力信号に時間窓を乗算する時間窓処理部と、前記時間窓処理部の出力を平均化する平均化処理部と、前記平均化処理部の出力を時間領域から周波数領域に変換する時間／周波数変換処理部と、前記時間／周波数変換処理部の出力から複数の歪成分を抽出する歪抽出部と、を備え、前記平均化処理部は、前記時間／周波数変換処理部が変換を行う単位毎に、平均を行い、前記時間窓の幅は、前記単位の整数倍であることを特徴とするものである。【００１７】
また、本発明に係る歪検出装置において、前記時間／周波数変換処理部はＦＦＴを用いることを特徴とするものである。【００１８】
また、本発明に係る歪検出装置において、前記平均化処理部は前記ＦＦＴの長さ毎に前記時間領域信号の加算を行うことを特徴とするものである。【００１９】
また、本発明は、外部から入力された信号を増幅し、歪を補償する歪補償増幅装置であって、第１の信号の増幅を行う増幅部と、前記増幅部の出力を分配する分配部と、前記分配部の一方の出力を平均化する平均化処理部と、前記平均化処理部の出力を時間領域から周波数領域に変換する時間／周波数変換処理部と、前記時間／周波数変換処理部の出力から歪成分を抽出する歪抽出部と、前記外部から入力された信号に対して、前記歪抽出部の出力に基づく予歪を与えたものを前記第１の信号とする歪補償部とを備えたものである。
【発明の効果】
【００２０】
本発明によれば、ハードウェアによる平均化処理を行った後に、ソフトウェアによる時間／周波数変換処理を行うことにより、収束時間を大幅に低減することができる。従って、広帯域信号であっても低コストで歪検出を行うことができる。
【図面の簡単な説明】
【００２１】
【図１】本発明に係る歪補償増幅装置の構成の一例を示すブロック図である。
【図２】本発明に係る歪検出装置の構成の一例を示すブロック図である。
【図３】本発明に係る時間／周波数変換を表す式１の一例を示す式である。
【図４】本発明に係る時間／周波数変換を表す式２の一例を示す式である。
【図５】本発明に係る時間／周波数変換を表す式３の一例を示す式である。
【図６】本発明に係るローゼンフィールド窓を表す式４の一例を示す式である。
【図７】本発明に係る第１の時間窓の一例を示す波形である。
【図８】本発明に係る第２の時間窓の一例を示す波形である。
【図９】本発明に係る第３の時間窓の一例を示す波形である。
【図１０】本発明に係る時間窓の周波数特性の一例を示すスペクトルである。
【図１１】歪を含む送信信号の一例を示すスペクトルである。
【図１２】従来の歪補償増幅装置の構成の一例を示すブロック図である。
【符号の説明】
【００２２】
１  増幅部
２  分配部
３  歪検出装置
４  補償係数演算部
１１  包絡線検出部
１２  逆特性記憶部
１３  遅延部
１４  ＰＤ部
２１  周波数変換部
２２  ＢＰＦ
２３  Ａ／Ｄ
３１  アドレス生成部
３２  時間窓生成部
３３  乗算部
３４  加算部
４１  カウンタ
４２  切り替え部
５１  デュアルポートメモリ
５２  ＦＦＴ処理部
５３  歪抽出部
【発明を実施するための最良の形態】
【００２３】
以下、本発明の実施の形態について図面を参照しつつ説明する。【００２４】
まず、本発明に係る歪補償増幅装置の構成について説明する。図１は、本発明に係る歪補償増幅装置の構成の一例を示すブロック図である。図１において、図１２と同一符号は図１２に示された対象と同一又は相当物を示しており、ここでの説明を省略する。また、図１２の歪補償増幅装置と比較すると図１の歪補償増幅装置は、歪検出装置１０３の代わりに歪検出装置３を備える。【００２５】
次に、本発明に係る歪検出装置３の構成について説明する。図２は、本発明に係る歪検出装置の構成の一例を示すブロック図である。この歪検出装置は、周波数変換部２１、ＢＰＦ（Band-Pass Filter）２２、Ａ／Ｄ変換部（Ａ／Ｄ）２３、アドレス生成部３１、時間窓生成部３２、乗算部３３、加算部３４、カウンタ４１、切り替え部４２、デュアルポートメモリ５１、ＦＦＴ処理部５２、歪抽出部５３を備える。本実施の形態では、歪検出装置３が、１０２４ポイントＦＦＴを用いて、４倍の時間である４０９６ポイントの時間領域信号から１度に歪検出を行う例について説明する。【００２６】
分配部２の出力は、周波数変換部２１でダウンコンバートされ、ＢＰＦ２２で帯域制限され、Ａ／Ｄ変換部２３でアナログからデジタルの時間領域信号に変換され、Ａ／Ｄ変換後にデジタル直交検波やフィルタリング等、必要な処理が適宜施される。アドレス生成部３１は、ＦＦＴのポイント数（時間／周波数変換を行う単位）の整数倍（例えば４倍）に対応するアドレスを出力する。時間窓生成部３２は、アドレス生成部３１からのアドレスに従って時間窓を出力する。乗算部３３は、Ａ／Ｄ変換部２３の出力と時間窓生成部３２の出力を乗算することにより、時間領域信号に対して時間窓処理を行う。【００２７】
加算部３４は、乗算部３３の出力とデュアルポートメモリ５１のDATA OUTの出力を加算する。カウンタ４１は、ＦＦＴポイント数毎に平均化回数をカウントし、平均化回数に１回切り替え部４２を下側に接続する。本実施の形態ではＦＦＴポイント数を１０２４、平均化回数を４とする。切り替え部４２は、カウンタ４１に従って、乗算部３３または加算部３４の出力をデュアルポートメモリ５１のDATA INへ入力する。デュアルポートメモリ５１は、アドレス生成部３１からのアドレスに従って時間領域信号を読み出し、DATA OUTから出力すると共に、DATA INに入力された時間領域信号をアドレス生成部３１からADDRESSへ入力されたアドレスに従って書き込む。結果として、１０２４ポイントの時間領域信号が４つ加算されたことにより平均化処理された平均化時間領域信号がデュアルポートメモリ５１に格納される。このように、時間窓処理、平均化処理はハードウェアで実現される。【００２８】
さらに、デュアルポートメモリ５１は、ＦＦＴ処理部５２からADDRESS2へ入力されたアドレスに従って平均化時間領域信号を読み出し、DATA OUT2から出力する。ＦＦＴ処理部５２は、平均化時間領域信号の時間／周波数変換処理を行い、得られた周波数領域信号を歪抽出部５３へ出力する。ＦＦＴ処理部５２による時間／周波数変換処理は、複素ＦＦＴであり、ＤＳＰを用いたバタフライ演算で実現される。歪抽出部５３は、周波数領域信号から−ＩＭ５，−ＩＭ３，＋ＩＭ３，＋ＩＭ５について対キャリア相対レベルを補償係数演算部４へ出力する。【００２９】
補償係数演算部４は、歪検出装置３により検出される歪が最小となるように、適応アルゴリズムを用いて逆特性記憶部１２の内容を更新する。具体的な構成は例えば特許文献１と同じでよい。また、メモリ効果を補償するための制御量をＰＤ部１４が有しているときは、−ＩＭ３と＋ＩＭ３のアンバランス等に基づき、その制御量を更新しても良い。【００３０】
次に、本発明に係る歪検出装置３の動作原理として、平均化処理、時間／周波数変換処理、時間窓処理について説明する。【００３１】
まず、加算部３４における４０９６ポイントから１０２４ポイントへの平均化処理と、ＦＦＴ処理部５２における１０２４ポイントの時間／周波数変換処理について説明する。ＦＦＴ処理部５２は、等価的には周波数に応じた回転因子を乗算し、積分することにより周波数領域信号を算出する。ここでは、規格化周波数が４Ｈｚの周波数領域信号Ａ４の例について示す。周波数領域信号Ａ４は、図３の式１で表される。ここで、ｓ（ｔ）は時間領域信号、ｔは時間、ｊは虚数単位である。式１における回転因子は、４０９６ポイントを４周期とする正弦波と余弦波であり、１０２４ポイントで１周期となる。【００３２】
次に、Ａ４を展開し、簡単のためにＡ４の第１項であるＡ４ｃだけを表した式を、図４の式２に示す。回転因子の周期は１０２４ポイントであり、１０２４ポイント離れた時間領域信号に乗算される回転因子は等しくなることから、式２における等しい回転因子の項をまとめた式を、図５の式３に示す。式３は、回転因子の周期だけ離れた時間領域信号を加算してから、回転因子を乗算し、積分しても式１と同じ結果が得られることを表している。従って本実施例の１０２４ポイントＦＦＴの結果は、通常の４０９６ポイントＦＦＴの結果を４ポイントおきに取り出したものと等しい。【００３３】
実際には、まず時間領域信号の平均化処理を加算部３４やデュアルポートメモリ５１等のハードウェアで行い、次に平均化時間領域信号の時間／周波数変換処理をＦＦＴ処理部５２のソフトウェアで行うことにより、時間／周波数変換処理に掛かる時間が低減され、適応アルゴリズムの更新周期を短くできる。また、リアルタイム信号の取りこぼしが減り、より長い時間信号に対して平均化処理を行うことができる。従って、送信信号の電力の変動の影響が軽減された、より平均的な歪を補償係数演算部４にフィードバックすることができ、適応アルゴリズムの収束が高速化する。この相乗効果により、従来のソフトウェアによるＦＦＴでは１０分程度掛かっていた収束時間が、本実施の形態によれば１０秒程度まで低減できる。また、本実施の形態では４０９６ポイントの時間領域信号毎に処理を行うが、周波数分解能は１０２４ポイントのＦＦＴに対応する。ただし、各周波数間の直交度は改善されている。【００３４】
なお、本実施の形態における平均化処理は、平均化回数だけ加算を行ってリセットするようにしたが、移動平均や、忘却係数を使った重み付き平均を行っても良い。【００３５】
次に、時間窓処理について説明する。ＦＦＴは周期信号を前提とする変換なので、リアルタイム信号を処理する際は、ＦＦＴ長（ＦＦＴポイント数）分の信号の始めと終わりを０に減衰させる窓関数を用いるのが通常である。ここでは、時間窓の例を挙げる。第１の時間窓は、平均化長の窓関数である。平均化長は、ＦＦＴ長×平均化回数である。窓関数であるローゼンフィールド窓ｗ（ｔ）は、図６の式４で表される。ここで、Ｍは窓関数長である。図７は、本発明に係る第１の時間窓の一例を示す波形である。横軸は時間をポイント数で表し、縦軸は振幅を表す。本実施の形態において第１の時間窓に用いるローゼンフィールド窓は、Ｍ＝４０９６とする。【００３６】
第２の時間窓は、ＦＦＴ長の窓関数を平均化回数だけ繰り返した時間窓である。図８は、本発明に係る第２の時間窓の一例を示す波形である。本実施の形態において第２の時間窓に用いるローゼンフィールド窓は、Ｍ＝１０２４とする。【００３７】
第３の時間窓は、立ち上がりにＦＦＴ長の窓関数の前半を用い、立ち下がりにＦＦＴ長の窓関数の後半を用い、途中が最大値１でホールドされた時間窓である。図９は、本発明に係る第３の時間窓の一例を示す波形である。本実施の形態において第３の時間窓の一部に用いるローゼンフィールド窓は、Ｍ＝１０２４とする。【００３８】
上述した３種類の時間窓の周波数特性は異なるため、信号のスペクトルに応じて選択する。図１０は、本発明に係る時間窓の周波数特性の一例を示すスペクトルである。また、時間窓生成部３２は、窓関数の波形を記憶させたメモリで実現することができ、窓関数は左右対称であることからメモリ容量は窓関数長の半分にすることもできる。特に、第２の時間窓、第３の時間窓は、５１２ポイント分のメモリ容量で済む。なお、本実施の形態においては、窓関数としてローゼンフィールド窓を用いたが、ブラックマン窓等、他の窓関数を用いても良い。【００３９】
窓関数は、帯域が広く平坦であるほど良いとは限らない。Ｗ−ＣＤＭＡのスペクトルは、ある時間で平均すれば図１１のようになるが、瞬時スペクトルはキャリア数のトーンに近いものであり、キャリアトーンの離調間隔とＦＦＴの周波致分解能との関係によっては、検出レベルが均一でなくなったり、特定のキャリアや歪が検出できなくなる可能性がある。従って、Ｗ−ＣＤＭＡのスペクトル測定には３０ｋＨｚ程度の分解能が推奨されている。一方、図１１のＩＭ３とＩＭ５の一部を取得するためにサンプルレートを１２３ＭＨｚとしても、１０２４ポイントＦＦＴでは周波数分解能は１２０ｋＨｚとなる。しかし、変調された実際のキャリアは少なからず帯域を持つことや、周波数成分を持つ窓関数と乗算されスペクトルが適度に拡散されることで、分解できなかった周波数成分に対しても問題にならない感度で検出はできる。従って４０９６ポイント未満のＦＦＴが使用できる。【００４０】
最適な窓関数はシステム構成に依存するが、複数の窓関数を切り替えて使う、或いは複数の窓関数を合成した窓関数を使う、などが有効な場合がある。例えば、元の窓関数（４０９６ポイント）に対し窓幅正規化周波数１Ｈｚ，２Ｈｚ，…（実周波数３０ｋＨｚ，６０ｋＨｚ…）などをそれぞれ乗算し、連結、或いは合成しても良い。また、ＦＦＴで得られる周波数領域信号は、該当する周波数を中心に、ナイキスト周波数の範囲でｓｉｎｃ関数状に感度を持つので、（特に位相の）周波数特性の異なる複数の窓関数を用いると、所望の分解周波数以外の周波数に対する感度を打ち消すことができる。これは窓関数にランダムなノイズを重畳することでも容易に達成でき、これによりスペクトル検出のダイナミックレンジが増大するので、本実施の形態のようにキャリアに比べ−６０ｄＢ以下のＩＭ５を精度良く検出する用途では極めて有効である。【００４１】
またどうしても分解能が不足するのであれば、回転因子を異ならせたＦＦＴ処理部５３、或いは予め１，２，３Ｈｚといった回転因子を乗算させて入カするＦＦＴ処理部５３を並列に設ければ良い。ＦＦＴの一般的な演算量のオーダ（ｎ×ｌｏｇ（ｎ））に比べ少ない演算量で処理できる可能性がある。【００４２】
次に歪抽出部５３について説明する。歪抽出部５３には、ポイント数に対応するスペクトル強度情報が入カされる。それらを、−ＩＭ５、−ＩＭ３、キャリア、＋ＩＭ３、＋ＩＭ５の領域に区分する。各領域は４キャリアを全て使っている時とそうでない時とでは異なるので、必要に応じてキャリアの使用状況を上位装置から取得するものとする。区分されたスペクトル強度情報は、２乗により電カ化し、各領域内で加算し、それぞれの領域の総電力となる。その後、各歪成分の総電力をキャリアの総電力で除算し、対キャリア相対レベルを出カする。各総電力或いは対キャリア相対レベルは複数のＦＦＴ回数分、更に平均化しても良い。【００４３】
或いは歪抽出部５３は、以下の方法で独自にキャリアの使用状況を把握して各ＩＭを求めても良い。即ち、キャリアは一定の周波数間隔Δｆで配置されるとして、予めキャリアとして使用可能な複数の周波数ｆCN（Ｎ＝１…ｎ）、及びキャリアの帯域幅ｆBWを記憶しているものとする。【００４４】
ステップ１として、スペクトル強度情報をスキャンし、最大の強度を検索する。ステップ２として、最大強度に予め定める１以下の数を乗じた値を闘値とする。ステップ３として、スペクトル強度情報をスキャンし、闘値を超えたスペクトルの周波数を含むキャリアを特定する。含むか否かは、スペクトルの周波数ｆに対し、ｆ＜ｆCN＋ｆBW／２且つｆ＞ｆCN−ｆBW／２を満たすｆCNが存在するか否かを判断して行う。これにより各キャリアの使用情況がわかる。【００４５】
ステップ４として、使用キャリアから、目的の各ＩＭの帯域を決定する。例えば使用キャリア数が２以上であれば、各使用キャリアの周波数にΔｆを加算もしくは減算した周波数が別のキャリアに使用されていなければ、その周波数が１つのＩＭ３の中心周波数である。ＩＭの帯域幅は例えばΔｆとする。複数のＩＭ３の内、周波数の（最も）高いものを＋ＩＭ３、低いものを−ＩＭ３とする。同様に２Δｆを加算もしくは減算した周波数がＩＭ５となり得る。使用キャリア数が１のとき、当該キャリアの中心周波数を基準として、予め記憶する周波数範囲のものを各ＩＭとする。【００４６】
なお、本実施の形態においては、時間窓処理、平均化処理をハードウェアで行い、時間／周波数変換処理をソフトウェアで行うとしたが、全てをソフトウェアで行っても良いし、全てをハードウェアで行っても良い。またＦＦＴのみならず同様のたたみこみを行う各種変換（ウェーブレット変換等）に応用してもよい。【００４７】
また、本発明に係る歪補償増幅装置は、無線通信装置に容易に適用することができ、無線通信装置の性能をより高めることができる。ここで、無線通信装置には、例えば携帯電話等の無線通信システムにおける基地局、端末、中継局、また無線ＬＡＮ（Local Area Network）等が含まれ得る。【００４８】
なお、時間窓処理部とは、実施の形態における時間窓生成部３２、乗算部３３に対応する。また、平均化処理部とは、実施の形態における加算部３４、カウンタ４１、切り替え部４２、デュアルポートメモリ５１に対応する。また、時間／周波数変換処理部とは、実施の形態におけるＦＦＴ処理部５２に対応する。
【図１】
0004619402.tif 000002
【図２】
0004619402.tif 000003
【図３】
0004619402.tif 000004
【図４】
0004619402.tif 000005
【図５】
0004619402.tif 000006
【図６】
0004619402.tif 000007
【図７】
0004619402.tif 000008
【図８】
0004619402.tif 000009
【図９】
0004619402.tif 000010
【図１０】
0004619402.tif 000011
【図１１】
0004619402.tif 000012
【図１２】
0004619402.tif 000013