(57)【要約】【課題】歪補償係数の収束時間を短縮する技術を提供する。【解決手段】アドレス毎に歪補償係数を有するルックアップテーブルを格納する記憶部と、入力信号に基づいてルックアップテーブルからアドレスを選択し、選択されたアドレスに格納された歪補償係数を取得し、取得された歪補償係数に基づいて入力信号のプリディストーションを行うプリディストーション部と、電力増幅器の出力を示すフィードバック信号から、入力信号に対する誤差を算出する誤差算出部と、誤差と取得された歪補償係数とから複数の新たな歪補償係数を適応アルゴリズムにより算出する係数算出部と、選択されたアドレス毎に、誤差に基づいて複数の新たな歪補償係数の中から複数の適合係数を選択する係数選択部と、アドレス毎に、複数の適合係数の平均値を算出し、歪補償係数を平均値に書き換える係数平均化部とを備える歪補償装置である。【選択図】図１
2011160043.tif 000002
【特許請求の範囲】
【請求項１】
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、
アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、
入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、
前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部と
を備える歪補償装置。
【請求項２】
前記係数選択部は、前記選択されたアドレス毎に、前記誤差から評価値を算出すると共に前記評価値の範囲を決定し、前記算出された評価値が前記範囲内にある場合に前記誤差から算出された歪補償係数を前記複数の適合係数の一つとして選択する、
請求項１に記載の歪補償装置。
【請求項３】
前記係数選択部は、前記誤差の絶対値を算出し、前記選択されたアドレス毎に、前記絶対値の平均値である基準値を算出し、前記基準値に対する前記絶対値の偏差を算出し、
前記評価値は、前記偏差である、
請求項２に記載の歪補償装置。
【請求項４】
前記係数選択部は、前記範囲を定義するしきい値を、前記基準値に基づいて決定する、
請求項３に記載の歪補償装置。
【請求項５】
前記係数選択部は、前記誤差の絶対値の標準偏差を算出し、前記標準偏差に基づいて前記しきい値を決定する、
請求項４に記載の歪補償装置。
【請求項６】
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、
アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、
入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、
前記プリディストーション結果の増幅を行う電力増幅器と、
前記電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部と
を備える送信装置。
【請求項７】
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償方法であって、
記憶部に格納されアドレス毎に歪補償係数を有するルックアップテーブルから、入力信号に基づいてアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行い、
前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出し、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出し、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択し、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える、
歪補償方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プリディストーションに用いられる歪補償係数を更新する技術に関する。
【背景技術】
【０００２】
無線通信の高速化に伴い、送信信号の広帯域化と高ダイナミックレンジ化が進んでいる。信号品質の劣化を最小限に抑えるために、送信装置には高い線形性が要求される。更に、装置の小型化、運用コスト削減、環境問題等の観点から、高効率で動作する電力増幅器が要求される。
【０００３】
しかしながら、一般的な電力増幅器において、線形性と電力変換効率は、相反関係にある。飽和電力から十分にバックオフした線形領域で電力増幅器を動作させることにより、帯域外歪の発生を小さくすることができる。しかし、電力変換効率が著しく低下し、装置の消費電力が増えてしまう。そのため、電力増幅器を効率の高い非線形領域で動作させる際に発生する非線形歪を除去する歪補償回路が使用される。
【０００４】
歪補償方式の一つであるプリディストーション方式は、電力増幅器の非線形歪の逆特性を予め送信信号に乗算することにより、電力増幅器出力の線形性を高める技術である。
【０００５】
ディジタル信号処理を用いるプリディストーション方式は、ディジタル・プリディストーション（ＤＰＤ）方式と呼ばれる。ＤＰＤ方式の一つに、メモリ内のルック・アップ・テーブル（ＬＵＴ）に歪補償係数を保持するＬＵＴ型ＤＰＤが知られている。各ＬＵＴアドレスには、送信信号の振幅に応じた歪補償係数が格納される。ＬＵＴ型ＤＰＤは、送信信号の振幅に基づくＬＵＴアドレスを決定し、そのＬＵＴアドレスに格納された歪補償係数を読み出し、その歪補償係数を送信信号に適用してプリディストーション信号を生成する。更に、ＬＵＴ型ＤＰＤは、電力増幅器の出力の一部をフィードバックした信号（フィードバック信号）に基づいて歪補償係数を更新することにより、電力増幅器の特性変動や経年変化に応じて歪補償係数を変化させることができる。
【０００６】
関連する技術として、歪補償係数の収束時間を短縮する非線形歪補償送信装置や、送信信号がＤＡ変換器のダイナミックレンジを越えないように補正する歪補償装置が知られている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００２−２２３１７１号公報
【特許文献２】特開２００１−２５１１４８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、ＬＵＴ型ＤＰＤにおいて、フィードバックループの雑音成分が大きい場合や、電力増幅器を低バックオフ動作する際に生じる瞬時的な特性変動が大きい場合、更新後の歪補償係数の精度が低くなり、歪補償係数が最適値に収束するまでの時間が長くなる場合がある。
【０００９】
本発明は上述した問題点を解決するためになされたものであり、歪補償係数の収束時間を短縮する技術を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上述した課題を解決するため、本発明の一態様は、ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部とを備える歪補償装置である。
【００１１】
また、本発明の別の一態様は、ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、前記プリディストーション結果の増幅を行う増幅器と、前記電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部とを備える送信装置である。
【００１２】
また、本発明の別の一態様は、ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、記憶部に格納されアドレス毎に歪補償係数を有するルックアップテーブルから、入力信号に基づいてアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行い、前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出し、前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出し、前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択し、前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える、歪補償方法である。
【発明の効果】
【００１３】
この出願に開示された技術によれば、歪補償係数の収束時間を短縮することができる。
【図面の簡単な説明】
【００１４】
【図１】送信装置１ａの構成を示すブロック図である。
【図２】歪補償処理を示すフローチャートである。
【図３】係数選択部３７ａと係数平均化部３３の構成を示すブロック図である。
【図４】ＤＰＤ部２ｘの構成を示すブロック図である。
【図５】実施の形態２における係数選択部３７ａの構成を示すブロック図である。
【図６】実施の形態２における第１特定アドレス演算を示すフローチャートである。
【図７】誤差信号の振幅頻度分布を示す図である。
【図８】複素平面における誤差信号の分布を示す図である。
【図９】実施の形態３における係数選択部３７ａの構成を示すブロック図である。
【図１０】実施の形態３における第１特定アドレス演算を示すフローチャートである。
【図１１】実施の形態４における係数選択部３７ａの構成を示すブロック図である。
【図１２】実施の形態５における係数選択部３７ａの構成を示すブロック図である。
【図１３】実施の形態５における第１特定アドレス演算を示すフローチャートである。
【図１４】偏差しきい値算出関数ｆ（ｍａk）の特性を示すグラフである。
【図１５】実施の形態６における係数選択部３７ａの構成を示すブロック図である。
【図１６】偏差しきい値算出関数ｇ（ｎ）の特性を示すグラフである。
【発明を実施するための形態】
【００１５】
以下、本発明の実施の形態について図面を参照しつつ説明する。
【００１６】
＜実施の形態１＞
以下、本発明の送信装置の適用例である送信装置１ａについて説明する。
【００１７】
図１は、送信装置１ａの構成を示すブロック図である。送信装置１ａは、ＤＰＤ部２ａ、ＤＡＣ（digital analog converter）３、変調部４、電力増幅器５、復調部６、ＡＤＣ（analog digital converter）７を有する。ＤＰＤ部２ａは、本発明の歪補償装置の適用例である。
【００１８】
送信装置１ａにより送信されるためのディジタル信号である送信信号は、ＤＰＤ部２ａへ入力される。送信信号は例えば、ベースバンド信号や中間周波（intermediate frequency：ＩＦ）信号の波形が所定のサンプリング時間間隔でサンプリングされたサンプル列であり、実部（Ｉ成分）および虚部（Ｑ成分）を有する複素サンプル列である。
【００１９】
ＤＰＤ部２ａは、送信信号に歪補償のためのプリディストーションを施し、プリディストーションを施された送信信号である補償後送信信号を出力する。
【００２０】
ＤＡＣ３は、ＤＰＤ部２ａからの補償後送信信号を、アナログ信号へ変換する。変調部４は、キャリア発生器１１ａ、乗算器１２ａを有する。キャリア発生器１１ａは、所定の送信キャリア周波数を有する正弦波を発生する。乗算器１２ａは、ＤＡＣ３の出力とキャリア発生器１１ａの出力とを乗算することによりアップコンバートを行う。
【００２１】
電力増幅器５は、変調部４の出力を増幅し、出力信号として出力する。出力信号は、電力増幅器５により与えられた非線形歪を含む。送信装置１ａが無線送信を行う場合、電力増幅器５の出力は、アンテナへ出力される。送信装置１ａが有線送信を行う場合、電力増幅器５の出力は、伝送路へ出力される。
【００２２】
復調部６は、キャリア発生器１１ｂ、乗算器１２ｂを有する。キャリア発生器１１ｂは、所定の受信キャリア周波数を有する正弦波を発生する。乗算器１２ｂは、電力増幅器５の出力の一部とキャリア発生器１１ｂの出力とを乗算することによりダウンコンバートを行う。ＡＤＣ７は、復調部６の出力をディジタル信号へ変換することにより、フィードバック信号を生成する。
【００２３】
以下、ＤＰＤ部２ａについて説明する。
【００２４】
ＤＰＤ部２ａは、ＲＥＦ（reference：参照）信号バッファ２１、ＦＢ（feedback：フィードバック）信号バッファ２２、減算器２３、アドレス生成部３１、アドレスバッファ３２、係数平均化部３３、ＬＵＴ３４、係数バッファ３５、乗算器３６、係数選択部３７ａ、係数算出部３８を有する。係数算出部３８は、複素共役演算器２４、乗算器２５、乗算器２６、加算器２７を有する。
【００２５】
アドレスバッファ３２の長さとＲＥＦ信号バッファ２１の長さとＦＢ信号バッファ２２の長さと係数バッファ３５の長さとは、バッファ長Ｎである。ここで、Ｎは正の整数である。アドレスバッファ３２は、アドレス生成部３１からのＬＵＴアドレスのＮサンプルを格納する。ＲＥＦ信号バッファ２１は、送信装置１ａへの入力である送信信号のＮサンプルを格納する。ＦＢ信号バッファは、ＡＤＣ７の出力であるフィードバック信号のＮサンプルを格納する。係数バッファ３５は、ＬＵＴ３４の出力である歪補償係数のＮサンプルを格納する。以後、Ｎサンプルの時間は１バッファ時間と定義される。
【００２６】
ＬＵＴ３４は、Ｍ個のＬＵＴアドレスにそれぞれＭ個の歪補償係数を格納する。あるＬＵＴアドレスに格納された歪補償係数は、そのＬＵＴアドレスに対応する送信信号に乗算される。ここで、ＭはＬＵＴアドレスの数であり、正の整数である。ＬＵＴアドレスはｋと定義される。ここで、ｋは０からＭ−１までの整数である。
【００２７】
以下、ＤＰＤ部２ａがＬＵＴ３４内の歪補償係数から新たな歪補償係数を算出するための係数算出方法について説明する。
【００２８】
この係数算出方法は、適応アルゴリズムを用いる。ここでは、説明を簡単にするため、ループ遅延や処理遅延等は、無視される。１バッファ時間を単位とする時刻ｎにおいて、ＲＥＦ信号バッファ２１内の１バッファ時間分の送信信号はｘ［ｎ］、ＦＢ信号バッファ２２内の１バッファ時間分のフィードバック信号はｙ［ｎ］、係数バッファ３５内の１バッファ時間分の歪補償係数はｈ［ｎ］、と定義される。参照信号である送信信号ｘ［ｎ］とフィードバック信号との差分である誤差信号はｅ［ｎ］と定義される。ＤＰＤ部２ａは、ＲＥＦ信号バッファ２１に蓄えられた参照信号、ＦＢ信号バッファ２２に蓄えられたフィードバック信号、および歪補償係数を使用して、次の式１により新たな歪補償係数である更新後係数ｈ［ｎ＋１］を算出する。
【００２９】
ｈ［ｎ＋１］＝ｈ［ｎ］＋μ×ｅ［ｎ］×ｈ［ｎ］×ｙ［ｎ］*
・・・（式１）
【００３０】
ここでμはステップサイズパラメータ、ｙ［ｎ］*はｙ［ｎ］の複素共役を表す。式１に示された歪補償係数の更新により、１バッファ時間分の更新後係数ｈ［ｎ＋１］が生成される。
【００３１】
アドレスバッファ３２内の１バッファ時間分のＬＵＴアドレスはｋ［ｎ］と定義される。ｋ［ｎ］のそれぞれは０からＭ−１までの整数のいずれか一つである。アドレスバッファ３２とＲＥＦ信号バッファ２１とＦＢ信号バッファ２２とにおいて参照されるサンプル番号はｉと定義される。ここで、ｉは０からＮ−１までの整数である。
【００３２】
アドレスバッファ３２内のＮサンプルのＬＵＴアドレスｋ［ｎ］のうちサンプル番号ｉに従って読み出されるＬＵＴアドレスは添え字ｉを用いてｋ＿ｉと定義される。同様に、ＲＥＦ信号バッファ２１内のＮサンプルの送信信号ｘ［ｎ］のうちサンプル番号ｉに従って読み出される送信信号はｘ＿ｉ、ＦＢ信号バッファ２２内のＮサンプルのフィードバック信号ｙ［ｎ］のうちサンプル番号ｉに従って読み出されるフィードバック信号はｙ＿ｉ、係数バッファ３５内のＮサンプルの歪補償係数ｈ［ｎ］のうちサンプル番号ｉに従って読み出される歪補償係数はｈ＿ｉと、それぞれ定義される。同様に、更新後係数ｈ［ｎ＋１］のうちサンプル番号ｉに従って算出される更新後係数はｈｕ＿ｉ、誤差信号ｅ［ｎ］のうちサンプル番号ｉに従って算出される誤差信号はｅ＿ｉと、それぞれ定義される。
【００３３】
送信信号ｘ＿ｉとＬＵＴアドレスｋ＿ｉとフィードバック信号ｙ＿ｉと歪補償係数ｈ＿ｉは、サンプル番号ｉに従って互いに同期するように、ＲＥＦ信号バッファ２１、ＦＢ信号バッファ２２、アドレスバッファ３２、係数バッファ３５からそれぞれ読み出される。アドレス生成部３１により送信信号ｘ＿ｉから決定されたＬＵＴアドレスはｋ＿ｉである。乗算器３６により送信信号ｘ＿ｉに乗算された歪補償係数はｈ＿ｉである。送信信号ｘ＿ｉと歪補償係数ｈ＿ｉの乗算結果が電力増幅器５により増幅された結果を示すフィードバック信号はｙ＿ｉである。
【００３４】
以下、ＤＰＤ部２ａによる歪補償処理について説明する。
【００３５】
図２は、歪補償処理を示すフローチャートである。この歪補償処理は、１バッファ時間毎に行われる。
【００３６】
アドレス生成部３１は、送信信号ｘ＿ｉの所定パラメータを示すＬＵＴアドレスｋ＿ｉを選択する（Ｓ２１１）。この例における所定パラメータは、送信信号ｘ＿ｉの振幅（絶対値）である。この場合、アドレス生成部３１は、予め複数の振幅の範囲と複数のＬＵＴアドレスとの対応関係を有する。乗算器３６は、ＬＵＴ３４のＬＵＴアドレスｋ＿ｉに格納された歪補償係数ｈ＿ｉを読み出し（Ｓ２１２）、歪補償係数ｈ＿ｉと送信信号ｘ＿ｉとを乗算することにより、送信信号ｘ＿ｉのプリディストーションを行い、プリディストーション結果である補償後送信信号を生成し、ＤＡＣ３へ出力する（Ｓ２１３）。
【００３７】
減算器２３は、ＲＥＦ信号バッファ２１内の送信信号ｘ＿ｉからＦＢ信号バッファ２２内のフィードバック信号ｙ＿ｉを減算することにより、誤差信号ｅ＿ｉを算出する（Ｓ２２１）。係数算出部３８は、前述の係数算出方法により、歪補償係数ｈ＿ｉと誤差信号ｅ＿ｉとフィードバック信号ｙ＿ｉとから更新後係数ｈｕ＿ｉを算出する（Ｓ２２２）。係数選択部３７ａは、誤差信号ｅ＿ｉに基づいて、更新後係数ｈｕ＿ｉの中から適合係数ｈｓ＿ｉを選択する（Ｓ２３１）。係数平均化部３３は、ＬＵＴアドレスｋに対応して算出された適合係数ｈｓ＿ｉの平均である平均係数ｈａkを算出する（Ｓ２４２）。係数平均化部３３は、ＬＵＴ３４内のＬＵＴアドレスｋに格納された歪補償係数を平均係数ｈａkに書き換える（Ｓ２４３）。
【００３８】
以上でこの歪補償処理は終了する。以後、１バッファ時間毎にこの歪補償処理が繰り返される。この歪補償処理のうちＳ２２１〜Ｓ２４３は、ＬＵＴ３４内の歪補償係数を更新する係数更新処理と定義される。
【００３９】
アドレス生成部３１により測定される所定パラメータは、送信信号ｘ＿ｉの位相であっても良いし、振幅および位相であっても良い。
【００４０】
以下、係数算出部３８について説明する。
【００４１】
複素共役演算器２４は、フィードバック信号ｙ＿ｉの複素共役信号ｙ＿ｉ*を算出する。乗算器２５は、係数バッファ３５内のｈ＿ｉと複素共役信号ｙ＿ｉ*とを乗算する。乗算器２６は、減算器２３により算出された誤差信号ｅ＿ｉと所定のステップサイズパラメータμと乗算器２５の出力とを乗算する。加算器２７は、係数バッファ３５内のｈ＿ｉと乗算器２６の出力とを加算することにより、更新後係数ｈｕ＿ｉを算出する。誤差信号ｅ＿ｉと、アドレスバッファ３２に格納されたＬＵＴアドレスｋ＿ｉと、更新後係数ｈｕ＿ｉとは、係数選択部３７ａへ入力される。
【００４２】
以下、係数選択部３７ａと係数平均化部３３について説明する。
【００４３】
図３は、係数選択部３７ａと係数平均化部３３の構成を示すブロック図である。係数選択部３７ａは、１入力Ｍ出力のセレクタ（デマルチプレクサ）５１，５２、Ｍ個の第１特定アドレス演算部５３ａを有する。係数平均化部３３は、Ｍ個の第２特定アドレス演算部６１、Ｍ入力１出力のセレクタ（マルチプレクサ）６２を有する。
【００４４】
以下、係数選択部３７ａについて説明する。
【００４５】
サンプル番号ｉの増加につれて、アドレスバッファ３２からのＬＵＴアドレスｋ＿ｉは、セレクタ５１およびセレクタ５２の制御端子へ順次入力される。同様に、誤差信号ｅ＿ｉは、セレクタ５１の入力端子へ順次入力される。同様に、更新後係数ｈｕ＿ｉは、セレクタ５２の入力端子へ順次入力される。
【００４６】
セレクタ５１は、１個の制御端子、１個の入力端子、Ｍ個のＬＵＴアドレスにそれぞれ対応するＭ個の出力端子ｅｏ0，ｅｏ1，…，ｅｏk，…，ｅｏM-1を有する。セレクタ５１において、制御端子にＬＵＴアドレスｋが入力された場合、ＬＵＴアドレスｋに対応する出力端子ｅｏkが選択され、入力端子へ入力された信号は選択された出力端子ｅｏkへ出力される。同様に、セレクタ５２は、１個の制御端子、１個の入力端子、Ｍ個のＬＵＴアドレスにそれぞれ対応するＭ個の出力端子ｈｕｏ0，ｈｕｏ1，…，ｈｕｏk，…，ｈｕｏM-1を有する。セレクタ５２において、制御端子にＬＵＴアドレスｋが入力された場合、ＬＵＴアドレスｋに対応する出力端子ｈｕｏkが選択され、入力端子へ入力された信号は選択された出力端子ｈｕｏkへ出力される。
【００４７】
Ｍ個の第１特定アドレス演算部５３ａは、Ｍ個のＬＵＴアドレスにそれぞれ対応して設けられる。Ｍ個の第１特定アドレス演算部５３ａのそれぞれは、条件決定部７１ａ、判定部７２、２入力１出力のセレクタ（マルチプレクサ）７３を有する。
【００４８】
以下、ＬＵＴアドレスｋに対する第１特定アドレス演算について説明する。
【００４９】
係数選択部３７ａへ入力されるＬＵＴアドレスｋ＿ｉがｋである場合、セレクタ５１の入力端子へ入力された誤差信号ｅ＿ｉは、ｋに対応して設けられた出力端子ｅｏkから出力される。出力端子ｅｏkからの誤差信号ｅ＿ｉは、ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ａへ入力される。同様に、係数選択部３７ａへ入力されるＬＵＴアドレスｋ＿ｉがｋである場合、セレクタ５２の入力端子へ入力された更新後係数ｈｕ＿ｉは、ｋに対応して設けられた出力端子ｈｕｏkから出力される。セレクタ５２の出力端子ｈｕｏkからの更新後係数ｈｕ＿ｉは、ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ａへ入力される。
【００５０】
セレクタ７３の２入力のうち、第１入力には固定値０が入力され、第２入力にはセレクタ５２の出力端子ｈｕｏkからの更新後係数ｈｕ＿ｉが入力される。条件決定部７１ａは、１バッファ時間内に入力される誤差信号ｅ＿ｉに基づいて選択条件を決定する。判定部７２は、誤差信号ｅ＿ｉおよび更新後係数ｈｕ＿ｉが入力される度に、誤差信号ｅ＿ｉが選択条件を満たすか否かを判定する。判定部７２は、誤差信号ｅ＿ｉが選択条件を満たす場合、セレクタ７３に第１入力を選択させ、誤差信号ｅ＿ｉが選択条件を満たさない場合、セレクタ７３に第２入力を選択させる。これにより、セレクタ７３は、誤差信号ｅ＿ｉが選択条件を満たす場合に更新後係数ｈｕ＿ｉを出力端子ｈｓｏkへ出力し、誤差信号ｅ＿ｉが選択条件を満たさない場合に更新後係数ｈｕ＿ｉを出力端子ｈｓｏkへ出力しない。
【００５１】
以下、係数平均化部３３について説明する。
【００５２】
Ｍ個の第２特定アドレス演算部６１は、Ｍ個のＬＵＴアドレスにそれぞれ対応して設けられる。
【００５３】
ここで、ＬＵＴアドレスｋに対する第２特定アドレス演算について説明する。更新後係数ｈｕ＿ｉのうち係数選択部３７ａから係数平均化部３３へ出力されたものが、適合係数ｈｓ＿ｉと定義される。ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ａの出力端子ｈｓｏkから適合係数ｈｓ＿ｉが出力された場合、その適合係数ｈｓ＿ｉはＬＵＴアドレスｋに対応する第２特定アドレス演算部６１へ入力される。ＬＵＴアドレスｋに対応する第２特定アドレス演算部６１は、１バッファ時間内に入力される適合係数ｈｓ＿ｉを平均化することにより、平均係数ｈａkを算出する。係数平均化部３３から出力されるｋが０からＭ−１までの全ての平均係数ｈａkは、ｈａで表される。
【００５４】
第２特定アドレス演算部６１による１バッファ時間の平均化が終了すると、０からＭ−１までのＬＵＴアドレスｋがセレクタ６２の制御端子とＬＵＴ３４の書き込みアドレスとへ順次入力される。セレクタ６２は、１個の制御端子、Ｍ個の入力端子ｈａｉ0，ｈａｉ1，…，ｈａｉk，…，ｈａｉM-1、１個の出力端子を有する。セレクタ６２において、制御端子にＬＵＴアドレスｋが入力された場合、ＬＵＴアドレスｋに対応する入力端子へ入力された信号が選択されて出力端子へ出力される。
【００５５】
次に、ＬＵＴアドレスｋに対応する第２特定アドレス演算部６１により算出された平均係数ｈａkは、セレクタ６２のＬＵＴアドレスｋに対応する入力端子ｈａｉkへ入力される。セレクタ６２は、制御端子へ入力されたＬＵＴアドレスｋに従って入力端子ｈａｉkを選択し、選択された入力端子ｈａｉkに入力された平均係数ｈａkを出力端子へ出力する。これに同期してＬＵＴアドレスｋがＬＵＴ３４の書き込みアドレスへ入力される。これにより、０からＭ−１までのＬＵＴアドレスｋに対応する平均係数ｈａkがＬＵＴ３４のＬＵＴアドレスｋへ順次書き込まれることにより、ＬＵＴ３４の歪補償係数が更新される。以上で、１回の係数更新処理が終了する。
【００５６】
このように、係数選択部３７ａは、更新後係数ｈｕ＿ｉのうち、誤差信号ｅ＿ｉが選択条件を満たす更新後係数ｈｕ＿ｉである適合係数ｈｓ＿ｉだけを、係数平均化部３３へ出力する。係数平均化部３３は、ＬＵＴアドレスｋ毎に、１バッファ時間内の適合係数ｈｓ＿ｉをＬＵＴアドレス毎に平均化する。
【００５７】
なお、ＤＰＤ部２ａは、回路等により実現されても良いし、ＤＳＰ（digital signal processor）やコンピュータのようなプロセッサとメモリにより実現されても良い。ＤＰＤ部２ａがプロセッサとメモリにより実現される場合、メモリはＤＰＤ部２ａの機能をプロセッサに実行させるソフトウェアを格納し、プロセッサはこのソフトウェアに従ってＤＰＤ部２ａの機能を実行する。
【００５８】
以下、ＤＰＤ部２ａの比較例であるＤＰＤ部２ｘについて説明する。
【００５９】
図４は、ＤＰＤ部２ｘの構成を示すブロック図である。ＤＰＤ部２ｘにおいて、ＤＰＤ部２ａと同一符号はＤＰＤ部２ａに示された対象と同一又は相当物を示しており、ここでの説明を省略する。ＤＰＤ部２ａと比較すると、ＤＰＤ部２ｘは、係数選択部３７ａを持たない。係数平均化部３３から出力されるｋが０からＭ−１までの全ての平均係数ｈａkは、ｈａｘで表される。
【００６０】
つまり、比較例のＤＰＤ２ｘにおける係数平均化部３３は、ＬＵＴアドレス毎に全ての更新後係数ｈｕ＿ｉを平均化することにより平均係数ｈａｘを算出する。そのため、雑音や電力増幅器の瞬時的な特性変動を含むフィードバック信号から生成された更新後係数の影響により、平均係数ｈａｘの精度が下がり、歪補償係数が最適値に収束するまでの時間が長くなる場合がある。
【００６１】
一方、実施の形態１のＤＰＤ部２ａは、係数平均化部３３の前段に係数選択部３７ａを有する。これにより、係数平均化部３３は、ＬＵＴアドレス毎に適合係数ｈｓ＿ｉのみを平均化することにより平均係数ｈａを算出する。つまり、ＤＰＤ部２ａにおける係数平均化部３３は、雑音や電力増幅器５の瞬時的な特性変動を含むフィードバック信号から生成された更新後係数を平均化に使用しない。従って、平均係数ｈａの精度が向上し、歪補償係数が最適値に収束するまでの時間を短縮することができる。これにより、電力増幅器５の特性変動が発生する動作条件においても、帯域外歪を良好に抑圧することが可能となる。
【００６２】
以下、この実施の形態の歪補償装置の他の形態について述べる。この歪補償装置（例えば、ＤＰＤ部２ａ）は、歪補償装置へ連続して入力される入力信号値（例えば、送信信号）に乗算される歪補償係数を所定期間毎に更新する。更にこの歪補償装置は、入力信号値の所定パラメータ（例えば、振幅）の値を示す複数のアドレスの夫々に複数の歪補償係数の夫々を記憶部に格納する記憶部（例えば、ＬＵＴ３４）を有する。更にこの歪補償装置は、入力信号値毎に、複数のアドレスの中から、入力信号値の所定パラメータの値を示すアドレスである適用アドレスを選択し、記憶部内の適用アドレスに格納された歪補償係数である適用係数と前記入力信号値とを乗算するプリディストーション部（例えば、乗算器３６）を有する。更にこの歪補償装置は、乗算の結果に基づく信号が増幅器（例えば、電力増幅器５）により増幅されると、入力信号値毎に、増幅の結果に基づく増幅信号値を取得し、入力信号値に対する増幅信号値の誤差を算出する誤差算出部（例えば、減算器２３）を有する。更にこの歪補償装置は、適用係数と誤差と増幅信号値とに基づいて新たな歪補償係数である推定係数を算出する係数算出部（例えば、係数算出部３８）を有する。更にこの歪補償装置は、誤差の算出により複数の入力信号値から算出された複数の誤差と推定係数の算出により複数の入力信号値から算出された複数の推定係数とを取得し、複数の誤差に基づいて複数の推定係数の中の複数の歪補償係数である複数の適合係数を選択する係数選択部（例えば、係数選択部３７ａ）を有する。更にこの歪補償装置は、複数のアドレスの中の対象アドレス毎に、複数の適合係数の中から、記憶部内の対象アドレスに格納された歪補償係数から算出された適合係数である複数の対象係数の夫々を選択し、複数の対象係数の平均値である平均係数を算出し、前記記憶部内の前記対象アドレスへ前記平均係数を書き込む係数平均化部（例えば、係数平均化部３３）を有する。
【００６３】
＜実施の形態２＞
以下、係数選択部３７ａの別の形態について説明する。
【００６４】
図５は、実施の形態２における係数選択部３７ａの構成を示すブロック図である。この図において、実施の形態１で示された係数選択部３７ａ内の要素と同一符号は実施の形態１で示された対象と同一又は相当物を示しており、ここでの説明を省略する。実施の形態１における係数選択部３７ａと比較すると、実施の形態２における係数選択部３７ａは、複数の第１特定アドレス演算部５３ａのそれぞれの代わりに複数の第１特定アドレス演算部５３ｂのそれぞれを有する。第１特定アドレス演算部５３ａと比較すると、第１特定アドレス演算部５３ｂは条件決定部７１ａの代わりに条件決定部７１ｂを有する。
【００６５】
以下、第１特定アドレス演算部５３ｂによる第１特定アドレス演算について説明する。
【００６６】
図６は、実施の形態２における第１特定アドレス演算を示すフローチャートである。ここでは、ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ｂによる第１特定アドレス演算について説明する。
【００６７】
まず、条件決定部７１ｂは、サンプル番号ｉ毎に、セレクタ５１の出力端子ｅｏkからの誤差信号ｅ＿ｉの振幅（絶対値）である瞬時誤差振幅ｍ＿ｉを算出する（Ｓ１１）。次に、条件決定部７１ｂは、瞬時誤差振幅ｍ＿ｉを１バッファ時間にわたって平均化することにより、平均誤差振幅ｍａkを算出する（Ｓ１２）。
【００６８】
次に、条件決定部７１ｂは、瞬時誤差振幅ｍ＿ｉと平均誤差振幅ｍａkに基づいて、瞬時誤差振幅ｍ＿ｉの標準偏差σkを算出する（Ｓ２１）。次に、条件決定部７１ｂは、標準偏差σkに基づいて選択条件を決定する（Ｓ２２）。ここでの選択条件は、偏差ｄ＿ｉが偏差しきい値Ｗkより小さいことである。偏差ｄ＿ｉは、瞬時誤差振幅ｍ＿ｉと平均誤差振幅ｍａkの差の絶対値である。偏差しきい値Ｗkは、Ｋ×σkと定義される。ここで、Ｋは、所定の正数である。即ち、条件決定部７１ｂは、標準偏差σkに基づいて偏差しきい値Ｗkを決定する
【００６９】
次に、判定部７２は、サンプル番号ｉを０として初期化する（Ｓ１１０）。
【００７０】
次に、判定部７２は、サンプル番号ｉがバッファ長Ｎ以上であるか否かを判定する（Ｓ１１１）。即ち、判定部７２は、１バッファ時間内の全てのサンプルに対する処理を終了したか否かを判定する。
【００７１】
サンプル番号ｉがＮ以上である場合（Ｓ１１１，Ｙｅｓ）、このフローは終了する。
【００７２】
サンプル番号ｉがＮより小さい場合（Ｓ１１１，Ｎｏ）、判定部７２は、セレクタ５１の出力端子ｅｏkから第１特定アドレス演算部５３ａへ誤差信号ｅ＿ｉが入力されたか否かを判定する（Ｓ１１２）。
【００７３】
誤差信号ｅ＿ｉが入力されていない場合（Ｓ１１２，Ｎｏ）、このフローはＳ１３０へ移行する。
【００７４】
誤差信号ｅ＿ｉが入力された場合（Ｓ１１２，Ｙｅｓ）、判定部７２は、瞬時誤差振幅ｍ＿ｉと平均誤差振幅ｍａkの差の絶対値である偏差ｄ＿ｉを算出する（Ｓ１２１）。次に、判定部７２は、偏差ｄ＿ｉが偏差しきい値Ｗkより小さいか否かを判定する（Ｓ１２２）。
【００７５】
偏差ｄ＿ｉが偏差しきい値Ｗkより小さい場合（Ｓ１２２，Ｙｅｓ）、判定部７２は、加算器２７からの更新後係数ｈｕ＿ｉを適合係数ｈｓ＿ｉとして係数平均化部３３へ出力し（Ｓ１２３）、このフローはＳ１３０へ移行する。
【００７６】
偏差ｄ＿ｉが偏差しきい値Ｗkより小さくない場合（Ｓ１２２，Ｎｏ）、判定部７２は、加算器２７からの更新後係数ｈｕ＿ｉを係数平均化部３３へ出力せず、このフローはＳ１３０へ移行する。
【００７７】
Ｓ１１２において誤差信号ｅ＿ｉが入力されていない場合（Ｓ１１２，Ｎｏ）、またはＳ１２３において判定部７２が適合係数ｈｓ＿ｉを係数平均化部３３へ出力した場合、またはＳ１２２において偏差ｄ＿ｉが偏差しきい値Ｗkより小さくない場合（Ｓ１２２，Ｎｏ）、判定部７２は、サンプル番号ｉを１増加させ（Ｓ１３０）、このフローはＳ１１１へ移行し、次のサンプル番号に対する処理が行われる。
【００７８】
前述のように、係数選択部３７ａは、偏差ｄ＿ｉが偏差しきい値Ｗkより小さい場合の更新後係数ｈｕ＿ｉのみを、適合係数ｈｓ＿ｉとして係数平均化部３３へ渡す。これにより、係数平均化部３３は、係数選択部３７ａから出力された適合係数ｈｓ＿ｉのみを用いて、更新後係数の平均化と歪補償係数の更新とを行うことができる。
【００７９】
以下、選択条件について説明する。
【００８０】
図７は、誤差信号の振幅頻度分布を示す図である。この図において、横軸は瞬時誤差振幅ｍ＿ｉを示し、縦軸は頻度（サンプル数）を示す。更にこの図において、横軸上のｍａkは平均誤差振幅を示し、斜線で示された範囲は選択条件を満たす範囲を示す。選択条件を満たす範囲は、偏差しきい値Ｗkにより定義される。この選択条件を満たす範囲は、瞬時誤差振幅ｍ＿ｉが（平均誤差振幅ｍａk−偏差しきい値Ｗk）より大きく且つ（平均誤差振幅ｍａk＋偏差しきい値Ｗk）より小さくなる範囲である。
【００８１】
前述のように送信信号ｘ＿ｉは複素信号であるため、誤差信号ｅ＿ｉは複素信号である。図８は、複素平面における誤差信号の分布を示す図である。この図において、横軸は誤差信号ｅ＿ｉの実部（Ｉ成分）を示し、縦軸は誤差信号ｅ＿ｉの虚部（Ｑ成分）を示す。更にこの図において、「×」印で示された複数の点は１バッファ時間内の誤差信号ｅ＿ｉを示し、斜線で示された範囲は選択条件を満たす範囲を示す。更にこの図において、半径Ｗkの円は選択条件を示す。この円の内部の領域は、選択条件を満たす誤差信号ｅ＿ｉの範囲を示す。この円の中心は誤差信号ｅ＿ｉの平均値を示す。
【００８２】
この実施の形態に示された選択条件によれば、係数選択部３７ａは、ばらつきに基づいて選択条件を決定し、選択条件を満たす誤差信号ｅ＿ｉから生成された更新後係数ｈｕ＿ｉだけを、適合係数ｈｓ＿ｉとして係数平均化部３３へ渡すことができる。
【００８３】
なお、選択条件は、誤差信号ｅ＿ｉの範囲、瞬時誤差振幅ｍ＿ｉの範囲、等により定義されても良い。
【００８４】
なお、係数選択部３７ａは、ＬＵＴアドレス毎の選択条件を決定するのではなく、ＬＵＴアドレスに関わらず１バッファ時間内の全ての誤差信号ｅ＿ｉから一つの選択条件を決定しても良い。この場合、係数選択部３７ａは例えば、ＬＵＴアドレスに関わらず１バッファ時間内の全ての誤差信号ｅ＿ｉから一つの偏差しきい値を決定する。
【００８５】
＜実施の形態３＞
以下、係数選択部３７ａの別の形態について説明する。
【００８６】
図９は、実施の形態３における係数選択部３７ａの構成を示すブロック図である。この図において、実施の形態２で示された係数選択部３７ａ内の要素と同一符号は実施の形態２で示された対象と同一又は相当物を示しており、ここでの説明を省略する。実施の形態２における係数選択部３７ａと比較すると、実施の形態３における係数選択部３７ａは、複数の第１特定アドレス演算部５３ｂのそれぞれの代わりに複数の第１特定アドレス演算部５３ｃのそれぞれを有し、新たに選択条件記憶部７４ｃを有する。第１特定アドレス演算部５３ｂと比較すると、第１特定アドレス演算部５３ｃは条件決定部７１ｂの代わりに条件決定部７１ｃを有する。
【００８７】
以下、第１特定アドレス演算部５３ｃによる第１特定アドレス演算について説明する。
【００８８】
選択条件記憶部７４ｃは、例えばメモリであり、選択条件テーブルを格納する。選択条件テーブルは、Ｍ個のＬＵＴアドレスにそれぞれ対応して予め定義されたＭ個の選択条件に関する情報を有する。この例において、ＬＵＴアドレスｋに対応する選択条件は、偏差ｄ＿ｉが偏差しきい値Ｗｔkより小さいことである。この例において、選択条件テーブルは、Ｍ個のＬＵＴアドレスにそれぞれ対応して予め定義されたＭ個の偏差しきい値Ｗｔkを有する。
【００８９】
図１０は、実施の形態３における第１特定アドレス演算を示すフローチャートである。ここでは、ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ｃによる第１特定アドレス演算について説明する。まず、条件決定部７１ｃは、実施の形態２と同様の処理Ｓ１１、Ｓ１２を行う。次に、条件決定部７１ｃは、選択条件記憶部７４ｃ内の選択条件テーブルからＬＵＴアドレスｋ＿ｉに対応して格納されたＷｔkを読み出し、読み出したＷｔkを偏差しきい値Ｗkとして決定する（Ｓ３１）。次に、判定部７２は、実施の形態２と同様の処理Ｓ１１０〜Ｓ１３０を行う。
【００９０】
この実施の形態に示された選択条件によれば、第１特定アドレス演算は、ＬＵＴアドレスに応じて予め設定された選択条件を用いることにより、条件決定部７１ｃの処理が簡単になる。
【００９１】
＜実施の形態４＞
以下、係数選択部３７ａの別の形態について説明する。
【００９２】
図１１は、実施の形態４における係数選択部３７ａの構成を示すブロック図である。この図において、実施の形態３で示された係数選択部３７ａ内の要素と同一符号は実施の形態３で示された対象と同一又は相当物を示しており、ここでの説明を省略する。実施の形態３における係数選択部３７ａと比較すると、実施の形態４における係数選択部３７ａは、複数の第１特定アドレス演算部５３ｃのそれぞれの代わりに複数の第１特定アドレス演算部５３ｄのそれぞれを有し、選択条件記憶部７４ｃの代わりに選択条件記憶部７４ｄを有し、新たに特性値算出部７５ｄを有する。第１特定アドレス演算部５３ｃと比較すると、第１特定アドレス演算部５３ｄは条件決定部７１ｃの代わりに条件決定部７１ｄを有する。
【００９３】
特性値算出部７５ｄは、ＲＥＦ信号バッファ２１内の送信信号を読み出し、送信信号の特性値を算出する。特性値は例えば、送信信号の平均電力、送信信号の帯域幅、送信信号のＰＡＰＲ（peak to average power ratio：ピーク電力対平均電力比）等である。更に、この特性値の範囲として、異なる複数の特性値範囲が定義される。次に、特性値算出部７５ｄは、特性値算出部７５ｄにより算出された特性値が含まれる特性値範囲を決定する。
【００９４】
選択条件記憶部７４ｄは、前述の複数の特性値範囲にそれぞれ対応して予め定義された複数の選択条件テーブルを格納する。複数の選択条件テーブルのそれぞれは例えば、実施の形態３の選択条件テーブルと同様、Ｍ個のＬＵＴアドレスにそれぞれ対応するＭ個の偏差しきい値を格納する。
【００９５】
前述の第１特定アドレス演算部５３ｃと比較すると、第１特定アドレス演算部５３ｄによる第１特定アドレス演算は、Ｓ３１における偏差しきい値Ｗkの決定の処理内容が異なる。Ｓ３１において、第１特定アドレス演算部５３ｄにおける条件決定部７１ｄは、選択条件記憶部７４ｄ内の複数の選択条件テーブルの中から、特性値算出部７５ｄにより決定された特性値範囲に対応して定義された一つの選択条件テーブルを選択する。次に、条件決定部７１ｄは、条件決定部７１ｃと同様、選択された選択条件テーブルに基づいて選択条件を決定する。
【００９６】
例えば、送信信号が広帯域である場合、一般的にはアンプのメモリ効果により電力増幅部５のアンプ特性の変動は大きくなる。この場合、条件決定部７１ｄは例えば、ＬＵＴアドレスｋの全範囲にわたって偏差しきい値Ｗｔkが大きいセットを使用する。あるいは例えば、電力増幅部５が飽和電力からバックオフされた平均電力で動作する場合、非線形歪量は小さくなる。この場合、条件決定部７１ｄは例えば、偏差しきい値Ｗkが全体的に小さい選択条件テーブルや、送信信号の振幅が小さい場合を示すＬＵＴアドレスｋに対応して定義された選択条件テーブルであって偏差しきい値Ｗkが小さいもの等を使用する。
【００９７】
この実施の形態に示された選択条件によれば、係数選択部３７ａは、複数の選択条件テーブルの中から送信信号の特性値に適した選択条件テーブルを選択することにより、送信信号の特性に応じて選択条件を変化させることができる。
【００９８】
＜実施の形態５＞
以下、係数選択部３７ａの別の形態について説明する。
【００９９】
図１２は、実施の形態５における係数選択部３７ａの構成を示すブロック図である。この図において、実施の形態２で示された係数選択部３７ａ内の要素と同一符号は実施の形態２で示された対象と同一又は相当物を示しており、ここでの説明を省略する。実施の形態２における係数選択部３７ａと比較すると、実施の形態５における係数選択部３７ａは、複数の第１特定アドレス演算部５３ｂのそれぞれの代わりに複数の第１特定アドレス演算部５３ｅのそれぞれを有する。第１特定アドレス演算部５３ｂと比較すると、第１特定アドレス演算部５３ｅは条件決定部７１ｂの代わりに条件決定部７１ｅを有する。
【０１００】
条件決定部７１ｅは、平均誤差振幅ｍａkから偏差しきい値Ｗkを決定する偏差しきい値算出関数Ｗk＝ｆ（ｍａk）を有する。
【０１０１】
以下、第１特定アドレス演算部５３ｅによる第１特定アドレス演算について説明する。
【０１０２】
図１３は、実施の形態５における第１特定アドレス演算を示すフローチャートである。ここでは、ＬＵＴアドレスｋに対応する第１特定アドレス演算部５３ｅによる第１特定アドレス演算について説明する。まず、条件決定部７１ｅは、実施の形態２と同様の処理Ｓ１１、Ｓ１２を行う。次に、条件決定部７１ｅは、算出された平均誤差振幅ｍａkと偏差しきい値算出関数ｆ（ｍａk）から偏差しきい値算出Ｗkを算出する（Ｓ５１）。次に、判定部７２は、実施の形態２と同様の処理Ｓ１１０〜Ｓ１３０を行う。
【０１０３】
図１４は、偏差しきい値算出関数ｆ（ｍａk）の特性を示すグラフである。この図において、横軸は平均誤差振幅ｍａkを示し、縦軸は偏差しきい値Ｗkを示す。この偏差しきい値算出関数ｆ（ｍａk）の例においては、平均誤差振幅ｍａkが小さくなるにつれて、偏差しきい値Ｗkは小さくなり、平均誤差振幅ｍａkに対する偏差しきい値Ｗkの傾きは大きくなる。
【０１０４】
係数更新処理が繰り返されることにより歪補償係数の収束が進むと、フィードバック信号が参照信号（送信信号）に近づく。そのため、誤差信号ｅ［ｎ］について、平均誤差振幅ｍａkおよび標準偏差σkは小さくなる。そこで前述のような特性を有する偏差しきい値算出関数ｆ（ｍａk）を用いることにより、歪補償係数の収束が進むにつれて、偏差しきい値Ｗkは小さくなり、平均化に用いられる更新後係数の精度は高くなり、収束時間を短縮することができる。
【０１０５】
＜実施の形態６＞
以下、係数選択部３７ａの別の形態について説明する。
【０１０６】
図１５は、実施の形態６における係数選択部３７ａの構成を示すブロック図である。この図において、実施の形態２で示された係数選択部３７ａ内の要素と同一符号は実施の形態２で示された対象と同一又は相当物を示しており、ここでの説明を省略する。実施の形態２における係数選択部３７ａと比較すると、実施の形態６における係数選択部３７ａは、複数の第１特定アドレス演算部５３ｂのそれぞれの代わりに複数の第１特定アドレス演算部５３ｆのそれぞれを有し、新たに更新回数カウンタ７６ｆを有する。第１特定アドレス演算部５３ｂと比較すると、第１特定アドレス演算部５３ｆは条件決定部７１ｂの代わりに条件決定部７１ｆを有する。
【０１０７】
更新回数カウンタ７６ｆは、所定時点から係数更新処理が実行された回数である更新回数ｎをカウントする。条件決定部７１ｆは、更新回数ｎから偏差しきい値Ｗkを決定する偏差しきい値算出関数Ｗk＝ｇ（ｎ）を有する。なお、条件決定部７１ｆは、更新回数ｎの代わりに、所定時点からの時間を示すパラメータを用いて偏差しきい値Ｗkを決定しても良い。
【０１０８】
前述の第１特定アドレス演算部５３ｅと比較すると、第１特定アドレス演算部５３ｆによる第１特定アドレス演算は、Ｓ５１における偏差しきい値算出Ｗkの算出の処理内容が異なる。Ｓ５１において、第１特定アドレス演算部５３ｆにおける条件決定部７１ｆは、更新回数カウンタ７６ｆによりカウントされた更新回数ｎと偏差しきい値算出関数ｇ（ｎ）から偏差しきい値算出Ｗkを算出する。
【０１０９】
図１６は、偏差しきい値算出関数ｇ（ｎ）の特性を示すグラフである。この図において、横軸は更新回数ｎを示し、縦軸は偏差しきい値Ｗkを示す。この偏差しきい値算出関数ｇ（ｎ）の例においては、更新回数ｎが大きくなるにつれて、偏差しきい値Ｗkは小さくなり、更新回数ｎに対する偏差しきい値Ｗkの傾きは小さくなる。
【０１１０】
係数更新処理が繰り返されることにより歪補償係数の収束が進むと、フィードバック信号が参照信号に近づく。そのため、平均誤差振幅ｍａkおよび標準偏差σkは小さくなる。そこで前述のような特性を有する偏差しきい値算出関数ｇ（ｎ）を用いることにより、歪補償係数の収束が進むにつれて、偏差しきい値Ｗkは小さくなり、平均化に用いられる更新後係数の精度は高くなり、収束時間を短縮することができる。
【０１１１】
以上の全ての実施の形態に関し、更に以下の付記を開示する。
（付記１）
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、
アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、
入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、
前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部と
を備える歪補償装置。
（付記２）
前記係数選択部は、前記選択されたアドレス毎に、前記誤差から評価値を算出すると共に前記評価値の範囲を決定し、前記算出された評価値が前記範囲内にある場合に前記誤差から算出された歪補償係数を前記複数の適合係数の一つとして選択する、
付記１に記載の歪補償装置。
（付記３）
前記係数選択部は、前記誤差の絶対値を算出し、前記選択されたアドレス毎に、前記絶対値の平均値である基準値を算出し、前記基準値に対する前記絶対値の偏差を算出し、
前記評価値は、前記偏差である、
付記２に記載の歪補償装置。
（付記４）
前記係数選択部は、前記範囲を定義するしきい値を、前記基準値に基づいて決定する、
付記３に記載の歪補償装置。
（付記５）
前記係数選択部は、前記誤差の絶対値の標準偏差を算出し、前記標準偏差に基づいて前記しきい値を決定する、
付記４に記載の歪補償装置。
（付記６）
前記係数選択部は、前記基準値と前記しきい値の関係を有し、前記関係に基づいて前記基準値から前記しきい値を決定する、
付記４に記載の歪補償装置。
（付記７）
前記係数選択部は、所定時点からの時間を示すパラメータを取得し、前記時間と前記範囲を定義するしきい値との関係を有し、前記関係に基づいて前記時間から前記しきい値を決定する、
付記３に記載の歪補償装置。
（付記８）
前記係数選択部は、前記ルックアップテーブル内のアドレスと前記範囲を定義するしきい値との関係を有し、前記関係に基づいて前記選択されたアドレスから前記しきい値を決定する、
付記２に記載の歪補償装置。
（付記９）
前記係数選択部は、前記入力信号の特性を示す特性値を算出し、前記特定値と前記ルックアップテーブル内のアドレスと前記範囲を定義するしきい値との関係を有し、前記関係に基づいて前記特性値と前記選択されたアドレスから前記しきい値を決定する、
付記２に記載の歪補償装置。
（付記１０）
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償装置であって、
アドレス毎に歪補償係数を有する前記ルックアップテーブルを格納する記憶部と、
入力信号に基づいて前記ルックアップテーブルからアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行うプリディストーション部と、
前記プリディストーション結果の増幅を行う電力増幅器と、
前記電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出する誤差算出部と、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出する係数算出部と、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択する係数選択部と、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える係数平均化部と
を備える送信装置。
（付記１１）
ルックアップテーブルを用いてプリディストーションを行うと共に、適応アルゴリズムを用いて前記ルックアップテーブルを最適化する歪補償方法であって、
記憶部に格納されアドレス毎に歪補償係数を有するルックアップテーブルから、入力信号に基づいてアドレスを選択し、前記選択されたアドレスに格納された歪補償係数を取得し、前記取得された歪補償係数に基づいて前記入力信号のプリディストーションを行い、
前記プリディストーション結果が入力された電力増幅器の出力を示すフィードバック信号から、前記入力信号に対する誤差を算出し、
前記誤差と前記取得された歪補償係数とから複数の新たな歪補償係数を前記適応アルゴリズムにより算出し、
前記選択されたアドレス毎に、前記誤差に基づいて前記複数の新たな歪補償係数の中から複数の歪補償係数である複数の適合係数を選択し、
前記ルックアップテーブル内のアドレス毎に、前記複数の適合係数の平均値を算出し、前記ルックアップテーブル内の歪補償係数を前記平均値に書き換える、
歪補償方法。
【０１１２】
記憶部は例えば、ＬＵＴ３４を含む。
電力増幅器は例えば、電力増幅器５を含む。
プリディストーションは例えば、アドレス生成部３１、乗算器３６を含む。
誤差算出部は例えば、減算器２３を含む。
係数算出部は例えば、係数算出部３８を含む。
係数選択部は例えば、係数選択部３７ａを含む。
係数平均化部は例えば、係数平均化部３３を含む。
アドレスは例えば、ＬＵＴアドレスを含む。
入力信号は例えば、送信信号を含む。
誤差は例えば、誤差信号を含む。
【符号の説明】
【０１１３】
１ａ  送信装置
２ａ  ＤＰＤ部
３  ＤＡＣ
４  変調部
５  電力増幅器
６  復調部
７  ＡＤＣ
１１ａ，１１ｂ  キャリア発生器
１２ａ，１２ｂ  乗算器
２１  ＲＥＦ信号バッファ
２２  ＦＢ信号バッファ
２３  減算器
２４  複素共役演算器
２５  乗算器
２６  乗算器
２７  加算器
３１  アドレス生成部
３２  アドレスバッファ
３３  係数平均化部
３４  ＬＵＴ
３５  係数バッファ
３６  乗算器
３７ａ  係数選択部
３８  係数算出部
５１，５２  セレクタ
５３ａ，５３ｂ，５３ｃ，５３ｄ，５３ｅ，５３ｆ  第１特定アドレス演算部
６１  第２特定アドレス演算部
６２  セレクタ
７１ａ，７１ｂ，７１ｃ，７１ｄ，７１ｅ，７１ｆ  条件決定部
７２  判定部
７３  セレクタ
７４ｃ，７４ｄ  選択条件記憶部
７５ｄ  特性値算出部
７５ｆ  更新回数カウンタ
【図１】
2011160043.tif 000003
【図２】
2011160043.tif 000004
【図３】
2011160043.tif 000005
【図４】
2011160043.tif 000006
【図５】
2011160043.tif 000007
【図６】
2011160043.tif 000008
【図７】
2011160043.tif 000009
【図８】
2011160043.tif 000010
【図９】
2011160043.tif 000011
【図１０】
2011160043.tif 000012
【図１１】
2011160043.tif 000013
【図１２】
2011160043.tif 000014
【図１３】
2011160043.tif 000015
【図１４】
2011160043.tif 000016
【図１５】
2011160043.tif 000017
【図１６】
2011160043.tif 000018