JP5115979B2

JP5115979B2 - プリディストータ

Info

Publication number: JP5115979B2
Application number: JP2008191239A
Authority: JP
Inventors: 孝基柴田; 拓也舩山; 孝義佐々木
Original assignee: Japan Radio Co Ltd
Current assignee: Japan Radio Co Ltd
Priority date: 2008-07-24
Filing date: 2008-07-24
Publication date: 2013-01-09
Anticipated expiration: 2028-07-24
Also published as: JP2010028766A

Description

本発明は、信号増幅器などの被補償回路からの出力信号の歪を補償するため、入力信号に予め加えるプリディストーション信号を生成するプリディストータに関するものである。

従来のプリディストータは、歪補償用のプリディストーション信号を生成するためのルックアップテーブル（以下、「ルックアップテーブル」を「ＬＵＴ」と略記する。）を備えており、ＬＵＴのテーブル値をアドレス毎に更新している（例えば、特許文献１を参照。）。特許文献１のプリディストータは、適応信号処理を用いてアドレス毎にＬＵＴのテーブル値の更新をするため、アドレス数（入力信号の振幅又は電力を量子化した値の数）に比例してメモリ量及び計算量が大幅に増えることになる。

このため、適応信号処理によるテーブル値の計算量及びメモリ量を減らす目的で、更新する際にアドレス数を減らしてテーブル値を算出し、後にアドレスの間のテーブル値を補間するプリディストータが知られている（例えば、特許文献２を参照。）。図１は、特許文献２のプリディストータにおけるＬＵＴ作成方法のフローチャートを示した図である。また、図２は、特許文献２のプリディストータが作成する歪補償値の例を示した図である。

ステップＳ１０１で、入力信号及び出力信号をサンプリングする。ステップＳ１０２で、入力信号と出力信号との差分を計算して誤差を算出する。ステップＳ１０３で、ＬＵＴの歪補償値と出力信号と誤差と更新係数を乗じて歪補償値更新量を算出する。ステップＳ１０３において入力信号と出力信号を量子化した各アドレスで求められた歪補償値の例を図２（ａ）に示す。このとき、図２（ｂ）のように、各アドレスに適応信号処理を適用して収束させておく。ステップＳ１０４及びステップＳ１０５でステップＳ１０３で計算した歪補償値更新量を用いてＬＵＴを更新し、新たなＬＵＴを作成する。ステップＳ１０６で、図２（ｃ）のように、歪補償値が存在しないアドレスの歪補償値を隣接するアドレスの歪補償値で直線補間し、図２（ｄ）のように、当該アドレスの歪補償値を直線補間した値から取得する。更に、ステップＳ１０６では図２（ｄ）の処理の後に図３のように、図２（ｄ）よりも入力信号の振幅値（または電力値）の量子化ステップを細かく設定したアドレスに対する歪補償値を直線補間した値から取得する。ステップＳ１０７で、当該アドレスを補間した歪補償値を用いたＬＵＴに差し替える。

一方、被補償回路の歪補償値そのものを推定するプリディストータも知られている（例えば、特許文献３を参照。）。図４は、特許文献３のプリディストータにおけるＬＵＴ作成方法のフローチャートを示した図である。また、図５は、特許文献３のプリディストータが作成する歪補償値の例を示した図である。

ステップＳ３０１及びステップＳ３０２は図１のステップＳ１０１及びステップＳ１０２と同じである。このときの各アドレスの歪補償値は図５（ａ）及び図５（ｂ）である。ステップＳ３０３で歪補償値を算出し、ステップＳ３０４で取り込んだ全ての入出力信号に対する歪補償値の総和を算出する。さらにステップＳ３０５で歪補償値が存在しないアドレスの歪補償値を補間し、ステップＳ３０６で平均化した結果で歪補償特性を更新する（ステップＳ３０７）。ステップＳ３０８で行列要素を算出し、ステップＳ３０９で歪補償多項式の係数を算出する。図５（ｃ）は最小二乗法で歪補償多項式を推定した例である。ステップＳ３１０で、ステップＳ３０９で求めた歪補償多項式により全アドレスの歪補償値を求め、ＬＵＴを作成する（図５（ｄ））。ステップＳ３１１で、作成したＬＵＴに差し替える。
特許３５６０３９８号特開２００２−２２３１７１号公報特開２００６−０９３９４７号公報

特許文献２のようにテーブル値を補間する方法は、メモリ量及び計算量を少なくすることはできるが、補間したテーブル値の精度が悪くなり歪補償が不十分であるという課題があった。一方、特許文献３のように歪補償値を推定する方法は、計算量が多くなるという課題があった。

そこで、本発明は、上記課題を解決するためになされたもので、歪補償するためのメモリ量及び計算量を低減でき、且つ歪補償の精度の高いプリディストータを提供することを目的とする。

前記目的を達成するために、本発明に係るプリディストータは、歪補償多項式を保有し、その係数を定期的に更新することとした。

具体的には、本発明に係るプリディストータは、歪補償値で入力信号を歪補償して被補償回路へ出力する歪補償部と、前記入力信号の強度を測定する強度算出部と、前記入力信号及び前記被補償回路の出力信号が入力され、歪補償多項式の係数を算出する多項式係数算出部と、前記歪補償多項式の係数を保管する歪補償多項式係数テーブルを有しており、前記強度算出部が測定した前記入力信号の強度に基づく前記歪補償多項式係数テーブルが保管する係数を適用した前記歪補償多項式で前記歪補償値を算出し、前記歪補償部へ前記歪補償値を出力する歪補償値算出部と、を備えるプリディストータであって、前記多項式係数算出部は、前記歪補償値算出部の前記歪補償多項式係数テーブルが保管する係数を適用した前記歪補償多項式に入力信号を代入して得た値と前記入力信号とを乗算して入力レプリカ信号を生成し、前記歪補償値算出部の前記歪補償多項式係数テーブルが保管する係数を適用した前記歪補償多項式に出力信号を代入して得た値と前記出力信号とを乗算して出力レプリカ信号を生成し、前記入力レプリカ信号と前記出力レプリカ信号との誤差が最小となる前記歪補償多項式の係数を算出し、前記歪補償値算出部は、前記多項式係数算出部が算出した前記歪補償多項式の係数で、前記歪補償多項式係数テーブルが保管する係数を更新することを特徴とする。

本発明に係るプリディストータは、更新及び保管するのは歪補償多項式の係数のみでありメモリ量及び計算量を大幅に減らせる。また、本発明に係るプリディストータでは歪補償多項式係数テーブルのテーブル値を歪補償多項式を用いて更新するので、任意の分解能で正確な歪補償値が求まり、補間処理も必要としない。結果として追従性能が向上し、収束も速くなる。

従って、本発明に係るプリディストータは、歪補償するためのメモリ量及び計算量を低減でき、歪補償の精度を高くすることができる。

本発明に係るプリディストータの前記歪補償値算出部は、前記歪補償部が歪補償を開始する際に、前記歪補償多項式係数テーブルが保管する係数を適用した前記歪補償多項式で前記歪補償値を算出する。

従来のプリディストータのように予め測定しておいたＬＵＴのテーブル値そのものではなく、歪補償多項式係数テーブルに保持しておいた歪補償多項式の係数を次回のプリディストーション開始の初期値とするので、少ないメモリ量で収束を速くできる。

本発明に係るプリディストータの前記歪補償多項式係数テーブルが保管する係数の数は、前記歪補償多項式の項数以下となる。

歪補償多項式の係数のうち、振幅（または電力）が大きい係数のみを歪補償多項式係数テーブルに保管することで、メモリ量を低減でき、歪補償の精度も維持することができる。

本発明は、歪補償するためのメモリ量及び計算量を低減でき、歪補償値の補間処理が不要のため歪補償の精度の高いプリディストータを提供することができる。また、プリディストーションを行うための初期値を保持しておくメモリ量を低減することができる。

添付の図面を参照して本発明の実施の形態を説明する。以下に説明する実施の形態は本発明の構成の例であり、本発明は、以下の実施の形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。

（実施の形態１）
図６は本実施形態のプリディストータ３０１の構成を説明するブロック図である。プリディストータ３０１は、歪補償値で入力信号Ｘを歪補償して被補償回路４０１へ出力する歪補償部１１と、入力信号Ｘの強度を測定する強度算出部１２と、入力信号Ｘ及び被補償回路４０１の出力信号Ｙが入力され、歪補償多項式の係数を算出する多項式係数算出部１３と、歪補償多項式の係数を保管する歪補償多項式係数テーブル１５を有しており、強度算出部１２が測定した入力信号Ｘの強度に基づく歪補償多項式係数テーブル１５が保管する係数を適用した歪補償多項式で歪補償値を算出し、歪補償部１１へ歪補償値を出力する歪補償値算出部１４と、を備える。例えば、被補償回路４０１は信号増幅器である。

また、歪補償多項式は予めプリディストータ３０１に設定されている。例えば、プリディストータ３０１は図示しない記憶回路を備えており、記憶回路が歪補償多項式を記憶していてもよい。多項式係数算出部１３及び歪補償値算出部１４は、必要に応じて記憶回路にアクセスし、記憶回路が記憶する歪補償多項式を取得することができる。

多項式係数算出部１３は、歪補償値算出部１４の歪補償多項式係数テーブル１５が保管する係数を適用した歪補償多項式に入力信号を代入して得た値と入力信号とを乗算して入力レプリカ信号を生成する。また、多項式係数算出部１３は、歪補償値算出部１４の歪補償多項式係数テーブルが保管する係数を適用した歪補償多項式に出力信号を代入して得た値と出力信号とを乗算して出力レプリカ信号を生成する。さらに、多項式係数算出部１３は、入力レプリカ信号と出力レプリカ信号との差分を誤差とし、この誤差が最小となる歪補償多項式の係数を算出する。

多項式係数算出部１３が生成する入力レプリカ信号とは、入力信号Ｘ及び歪補償多項式に入力信号を代入して得た値を乗算して被補償回路４０１に入力する予歪補償信号Ａについて算出したものである。一方、多項式係数算出部１３が生成する出力レプリカ信号とは、出力信号Ｙ及び歪補償多項式に出力信号を代入して得た値を乗算して被補償回路４０１に入力する予歪補償信号Ａについて算出したものである。入力レプリカ信号と出力レプリカ信号との誤差が最小となる歪補償多項式の係数を算出する歪補償アルゴリズムについては後述する。

強度算出部１２は、入力信号Ｘの強度を測定する。強度とは、例えば、入力信号Ｘの振幅や電力である。

歪補償値算出部１４は、多項式係数算出部１３が算出した歪補償多項式の係数で、歪補償多項式係数テーブル１５が保管する係数を更新する。また、歪補償値算出部１４は、歪補償多項式係数テーブル１５が保管する係数を適用した歪補償多項式を用いて、強度算出部１２が測定した入力信号Ｘの強度に基づく歪補償値を歪補償部１１に出力する。なお、歪補償多項式係数テーブル１５が保管する係数の数は、歪補償多項式の項数以下でもよい。

歪補償部１１は、歪補償値算出部１４からの歪補償値で入力信号Ｘを歪補償し、予歪補償信号Ａを生成して被補償回路４０１へ出力する。このように、プリディストータ３０１は、入力信号Ｘの強度に応じた歪補償値で歪補償できるため、出力信号Ｙの歪を小さくすることができる。

図７は、プリディストータ３０１の動作を説明する図である。取り込み時間ｔ_ｓ１に、多項式係数算出部１３は入力信号Ｘ及び出力信号Ｙを取り込む。Ｎは取り込み時間ｔ_ｓｉ（ｉは自然数）における入力信号Ｘ及び出力信号Ｙの取り込み点数である。その後、時間ｔ_ｃ１に、多項式係数算出部１３は歪補償多項式の係数を算出し、歪補償値算出部１４は歪補償多項式係数テーブル１５の係数を更新し、歪補償部１１は更新された係数に基づく歪補償値で歪補償を行う。Ｍは時間ｔ_ｃｉにおけるサンプリング間隔数である。更新された係数に基づく歪補償値の歪補償は、時間ｔ_ｃ１で歪補償多項式係数テーブル１５の係数が更新され、次の取り込み時間ｔ_ｓ２が終了し、時間ｔ_ｃ２で歪補償多項式係数テーブル１５の係数が更新されるまで行われる。このように、プリディストータ３０１は、取り込み時間ｔ_ｓｉと時間ｔ_ｃｉを繰り返すことで、温度変化、湿度変化又は経年変化する被補償回路４０１の歪特性に追従して適正な歪補償をすることができる。

なお、歪補償部１１が歪補償を開始する際には、歪補償値算出部１４は、歪補償多項式係数テーブル１５が予め保管する係数を適用した歪補償多項式で歪補償値を算出するように動作する。

（歪補償アルゴリズム）
次に、プリディストータ３０１が歪補償多項式係数を算出する歪補償アルゴリズムについて説明する。ここでは、被補償回路４０１がゲインＧの増幅装置であるとして説明する。図８は、歪補償アルゴリズムについてのフローチャートを示した図である。図９は、プリディストータ３０１が作成する歪補償値の例を示した図である。

ここで、増幅装置に入力する信号をｘ（ｎ）とすると、出力信号ｙ（ｎ）は次式で表される。増幅装置の入出力信号の関係は次式のように線形であるのが理想的である。

Ｇは実数定数である。なお、ここでは、以降の議論を簡単にする目的で、Ｇ＝１とおくこととする。

しかし、実際の増幅装置では、入力信号Ｘの振幅や電力が大きくなると入出力信号の関係は数式１のような線形ではなく、例えば、数式２のように非線形となることが知られている。

ここでのｇ_ｊは複素数であり、０≦ｊ≦Lである。

そこで、増幅装置の入出力関係を線形にする目的で、増幅装置の非線形な歪特性（歪値）に対する逆の歪特性（歪補償値）を生成するプリディストータ３０１を増幅装置の前段に設けている。プリディストータ３０１は、入力信号Ｘに歪補償値を予め与えた予歪補償信号Ａを生成して増幅装置に入力する。すなわち、入力信号Ｘを周期Ｔ_ｓでサンプリングした離散時間信号をｘ（ｎＴ_ｓ）、出力信号Ｙを周期Ｔ_ｓでサンプリングした離散時間信号をｙ（ｎＴ_ｓ）とし、表記を簡単にするためにそれぞれｘ（ｎ）、ｙ（ｎ）で表すこととする。このとき、増幅装置の出力ｙ（ｎ）をｙ（ｎ）＝ｘ（ｎ）とする増幅装置への入力信号を得るために、プリディストータ３０１は、入力信号Ｘに歪補償値を予め与えた次式の予歪補償信号Ａの関数ａ（ｎ）を生成する。

但し、ここでのｗ_ｊは複素数であり、０≦ｊ≦Ｋである。

数式２で与えられる非線形な歪を補償する逆特性（歪補償値）を与えるａ（ｎ）を得るには、数式２において、ｘ（ｎ）とｙ（ｎ）を入れ替えた数式４の関係を満足する係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを得ればよい。

数式４における係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを求めるには、Ｋ＋１個以上のｘ（ｎ）とｙ（ｎ）の組が必要である。ここで、数式４を行列表現すると数式５となる。

歪補償値を推定する際に使用する入出力信号の取り込み点数はＮ（Ｎ≧Ｋ＋１）なので、数式４をＮ個まとめると数式５は数式６の連立方程式となる。

である。また、Ｔは行列の転置を表す。

数式６又は数式７の連立方程式を解くことで、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋが求まる。数式６又は数式７の連立方程式を解くには掃き出し法を用いても良いし、最小二乗法を用いて数式８としても良い。

但し、Hは行列の複素共役転置を表す。このようにして得られた係数ｗ_０、ｗ_１、・・・、ｗ_Ｋをａ（ｎ）に適用し、ｘ（ｎ）に応じた歪補償値を算出すれば歪補償ができる。

増幅装置の歪特性が変化せず、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを求める際に用いる取り込み点数Ｎが十分に大きく確保できる場合は数式６の連立方程式を一回だけ解けばよい

しかし、一方で実際には増幅装置の歪特性は温度変化、湿度変化又は経年変化などにより変動し、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを求める際に用いる取り込み点数Ｎが十分に確保できない場合が多く、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを更新していく必要がある。そこで、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを更新する場合について述べる。

歪補償部１１が出力する予歪補償信号Ａについて、入力信号Ｘに歪補償値を乗算して入力レプリカ信号Ａ’（関数ａ’（ｎ））を作成する。

また、予歪補償信号Ａについて、出力信号Ｙに歪補償値を乗算した出力レプリカ信号Ａ”（関数ａ”（ｎ））を作成する。

ｙ（ｎ）＝ｘ（ｎ）となるとき、

が成り立つので、

が満たされる。すなわち、プリディストータ３０１と増幅装置とを含む全体の入出力関係を線形とする条件ｙ（ｎ）＝ｘ（ｎ）を満たす係数ｗ_０、ｗ_１、・・・、ｗ_Ｋは数式１２を満たす。従って、ａ’（ｎ）とａ”（ｎ）とが一致するように、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを求めればよい。

ここで、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋの更新について説明する。まず、ｉ回目（ｉ≧０）の更新時における係数の値をそれぞれｗ_０（ｉ）、ｗ_１（ｉ）、・・・、ｗ_Ｋ（ｉ）とおく。ｉ回目の更新時における係数ｗ（ｉ）の行列は数式１３となる。

なお、数式１３の係数ｗ（ｉ）の行列の初期値としては、ｗ（０）＝（１，０，・・・，０）^Ｔを用いるか、同じ被補償部に対して十分な取り込み点数Ｎのときに数式６、数式７又は数式８を用いて予め求めておいた係数の行列ｗをｗ（０）＝ｗとして用いるか、同じ被補償部に対して十分に更新回数を確保して予め求めておいた係数の行列ｗ（∞）をｗ（０）＝ｗ（∞）として用いればよい。

図７の取り込み時間ｔ_ｓｉにおいてステップＳ７０１を行う。図７の時間ｔ_ｃｉでステップＳ７０２〜ステップＳ７０８を行う。ステップＳ７０１でｉ回目の更新時における係数ｗ_０（ｉ）、ｗ_１（ｉ）、・・・、ｗ_Ｋ（ｉ）に対応した入力信号Ｘ及び出力信号Ｙを取り込み、ステップＳ７０２で入力信号Ｘ及び出力信号Ｙの行列を生成する。入力信号Ｘ及び出力信号Ｙの行列は数式１４である。

ステップＳ７０３で入力レプリカ信号Ａ’と出力レプリカ信号Ａ”を生成する。数式９及び数式１０の入力レプリカ信号Ａ’と出力レプリカ信号Ａ”を行列で表すと数式１５のようになる。

次に、ステップＳ７０４で、入力レプリカ信号Ａ’と出力レプリカ信号Ａ”の差分である誤差ｅ（ｉ）を算出する。誤差ｅ（ｉ）は数式１６で表される。

ここで、係数ｗ_０、ｗ_１、・・・、ｗ_Ｋを更新する量を更新量Δｗ（ｉ）とすると、数式６の連立方程式は数式１７となる。同様に数式７の連立方程式は数式１８となる。

さらに、ステップＳ７０５で行列要素を計算し、ステップＳ７０６で多項式係数を更新する量Δｗ（ｉ）を算出する。数式１７より行列Δｗ（ｉ）は数式１９となる。

ステップＳ７０７で数式２０となるように係数ｗ（ｉ）を更新する。

但し、μは０＜μ≦１．０を満たす。最後にステップＳ７０８で歪補償多項式係数テーブル１５が保管する係数をステップＳ７０７で更新した係数に差し替える。図９（ｂ）は、更新した係数に基づく歪補償多項式による、入力信号電力を量子化した全てのアドレスの歪補償値の曲線の例である。

また、入力信号に対応した歪補償値を得るのにＬＵＴを用いることもできる。図１０は、入力信号に対応した歪補償値を得るのにＬＵＴを用いた場合のフローチャートである。ステップＳ７０１〜ステップＳ７０７までは図８の説明と同じである。歪補償値を得る際にＬＵＴを用いる場合は、ステップＳ９０１でｉ回目の更新時における係数の値ｗ（ｉ）に、数式２１のように想定される入力信号の全ての振幅値又は電力値を量子化した値｜ｘ｜を代入してＬＵＴを作成する。

ここでの｜ｘ｜は、０≦｜ｘ｜≦｜ｘ｜_ｍａｘを満たす全ての量子化された値である。例えば、０≦｜ｘ｜≦８１９１を満たす全ての整数値とする。

その後、ステップＳ９０２でＬＵＴを差し替える。図９（ｃ）は、更新した係数に基づく歪補償多項式を利用して求めた各アドレスの歪補償値（ＬＵＴ）の例である。

（メモリ量の比較）
次に、従来の方法と本発明において使用するメモリ量について比較する。従来の方法では、ＬＵＴにおいて入力信号の振幅又は電力を量子化したアドレス数をＬビットで表現し、歪補償値をＢビットで量子化して表現している場合、歪補償値としては実数部と虚数部が必要なので、２×２^ＬＢ＝２^Ｌ＋１Ｂビット必要となる。さらに、従来の方法では、歪補償部に用いるＬＵＴに加えて、ＬＵＴの更新に用いるＬＵＴを別に保持しておく必要があるので、必要なビット数は２倍となり、２×２^Ｌ＋１Ｂ＝２^Ｌ＋２Ｂビットのメモリが必要である。

一方、プリディストータ３０１はＬＵＴではなく、歪補償多項式の係数を保持するのみであるので、歪補償多項式の次数をＫ（すなわち項数を（Ｋ＋１））とし、係数をＣビットで表現している場合、係数としては実数部と虚数部が必要なので、メモリ量は２（Ｋ＋１）Ｃビット必要である。プリディストータ３０１においても歪補償多項式の係数の更新に用いる歪補償多項式の係数を別に保持しておく必要があるので、必要なビット数は２倍となり、４（Ｋ＋１）Ｃビットのメモリが必要である。

例えば、Ｌ＝１３、Ｂ＝１６、Ｋ＝１００およびＣ＝３２の場合について考える。このとき、従来の方法では
２^１３＋２×１６＝５２４２８８
ビットのメモリが必要なのに対し、プリディストータ３０１では
４×（１００＋１）×３２＝１２９２８
ビットのメモリが必要である。従って、プリディストータ３０１を用いた場合、必要なメモリ量は従来の
１２９２８／５２４２８８＝０．０２４６６
程度で済む。

本発明に係るプリディストータは、移動体通信基地局などに用いられる無線送信機の電力増幅器に適用することができる。

従来のプリディストータにおけるＬＵＴ作成方法のフローチャートを示した図である。従来のプリディストータが作成する歪補償値の例を示した図である。（ａ）は入力信号と出力信号を量子化した各アドレスで求められた歪補償値の例である。（ｂ）は各アドレスに適応信号処理を適用して収束させた歪補償値である。（ｃ）は歪補償値が存在しないアドレスの歪補償値を隣接するアドレスの歪補償値で直線補間したものである。（ｄ）は歪補償値が存在しないアドレスの歪補償値を直線補間した値から取得したものである。従来のプリディストータが作成する歪補償値の例を示した図である。従来のプリディストータにおけるＬＵＴ作成方法のフローチャートを示した図である。従来のプリディストータが作成する歪補償値の例を示した図である。（ａ）は入力信号と出力信号を量子化した各アドレスで求められた歪補償値の例である。（ｂ）は各アドレスの歪補償値の平均値である。（ｃ）は最小二乗法で歪補償多項式を推定した例である。（ｄ）は（ｃ）で求めた歪補償多項式から全アドレスの歪補償値を求めた例である。本発明に係るプリディストータの構成を説明するブロック図である。本発明に係るプリディストータの動作を説明する図である。本発明に係るプリディストータの歪補償アルゴリズムについてのフローチャートを示した図である。本発明に係るプリディストータが作成する歪補償値の例を示した図である。（ａ）は入力信号を量子化した各アドレスで求められた歪補償値の例である。（ｂ）は（ａ）のサンプルを用いて入力信号電力を量子化した全てのアドレスの歪補償値に適応信号処理を適用して求めた歪補償値を表す歪補償多項式による曲線の例である。（ｃ）は（ｂ）の歪補償多項式を利用して求めた各アドレスの歪補償値（ＬＵＴ）の例である。本発明に係るプリディストータの歪補償アルゴリズムについてのフローチャートを示した図である。

符号の説明

３０１：プリディストータ
１１：歪補償部
１２：強度算出部
１３：多項式係数算出部
１４：歪補償値算出部
１５：歪補償多項式係数テーブル
４０１：被補償回路
Ｘ：入力信号
Ｙ：出力信号
Ａ：予歪補償信号
Ｎ：取り込み点数
Ｍ：サンプリング間隔数
ｔ_ｓｉ：取り込み時間
ｔ_ｃｉ：時間

Claims

歪補償値で入力信号を歪補償して被補償回路へ出力する歪補償部と、
前記入力信号の強度を測定する強度算出部と、
前記入力信号及び前記被補償回路の出力信号が入力され、歪補償多項式の係数を算出す
る多項式係数算出部と、
前記歪補償多項式の係数を保管する歪補償多項式係数テーブルを有しており、前記強度
算出部が測定した前記入力信号の強度に基づく前記歪補償多項式係数テーブルが保管する
係数を数Ｃ１で表される前記歪補償多項式に適用して予歪補償信号を算出し、前記予歪補償信号を前記歪補償値として前記歪補償部へ出力する歪補償値算出部と、
を備えるプリディストータであって、
前記多項式係数算出部は、
現在の入力信号及び出力信号から数Ｃ２で表されるｉ回目の更新における入力信号行列及び出力信号行列を生成し、前記歪補償値算出部の前記歪補償多項式係数テーブルが保管する（ｉ−１）回目に更新した係数の行列と前記入力信号行列及び出力信号行列を用いて数Ｃ３で表される入力レプリカ信号及び出力レプリカ信号を生成し、前記入力レプリカ信号と前記出力レプリカ信号との誤差行列が最小となる前記歪補償多項式の係数の行列を算出し、
前記歪補償値算出部は、前記多項式係数算出部が算出した前記歪補償多項式の係数の行列で、前記歪補償多項式係数テーブルが保管する係数を更新することを特徴とするプリディストータ。

ここで、ａ（ｎ）は予歪補償信号、ｘ（ｎ）は入力信号サンプリングした離散時間信号、ｗ _ｊ（ｊ＝１、２、３、・・・、Ｋ）は複素数であり歪補償多項式の係数である。

Ｎはサンプリング番号、Ｍはサンプリング間隔数である。

ｗ（ｉ−１）は歪補償多項式の係数（ｉ−１回目の更新）である。
前記歪補償値算出部は、前記歪補償部が歪補償を開始する際に、前記歪補償多項式係数テーブルが保管する係数を適用した前記歪補償多項式で前記歪補償値を算出することを特徴とする請求項１に記載のプリディストータ。
前記歪補償多項式係数テーブルが保管する係数の数は、前記歪補償多項式の項数以下となることを特徴とする請求項１又は２に記載のプリディストータ。