(57)【要約】【課題】合成後の出力高周波信号の劣化を抑制する歪補償装置を提供する。【解決手段】入力信号に対し、前記入力信号に依存する歪補償係数に基づいて、歪補償を行う歪補償部と、前記歪補償部の出力を、前記歪補償部の出力の振幅に基づく位相差を有する定振幅の第１信号及び第２信号に分離する信号生成部と、前記第１信号に対し、第１係数を乗算する第１調整部と、前記第１調整部の出力を増幅する第１増幅部と、前記第２信号に対し、第２係数を乗算する第２調整部と、前記第２調整部の出力を増幅する第２増幅部と、前記第１増幅部の出力と前記第２増幅部の出力とを合成する合成部と、前記歪補償部の出力と、前記第１調整部の出力と、前記第２調整部の出力と、前記合成部の出力とに基づいて、前記歪補償係数、前記第１係数、前記第２係数を算出する演算部と、を備える歪補償装置とする。【選択図】図７
2014011653.tif 000002
【特許請求の範囲】
【請求項１】
入力信号に対し、前記入力信号に依存する歪補償係数に基づいて、歪補償を行う歪補償部と、
前記歪補償部の出力を、前記歪補償部の出力の振幅に基づく位相差を有する定振幅の第１信号及び第２信号に分離する信号生成部と、
前記第１信号に対し、第１係数を乗算する第１調整部と、
前記第１調整部の出力を増幅する第１増幅部と、
前記第２信号に対し、第２係数を乗算する第２調整部と、
前記第２調整部の出力を増幅する第２増幅部と、
前記第１増幅部の出力と前記第２増幅部の出力とを合成する合成部と、
前記歪補償部の出力と、前記第１調整部の出力と、前記第２調整部の出力と、前記合成部の出力とに基づいて、前記歪補償係数、前記第１係数、前記第２係数を算出する演算部と、
を備える歪補償装置。
【請求項２】
前記演算部は、前記合成部の出力のうち前記合成部の出力の大きさが所定値以上であるものに基づいて前記歪補償係数を算出し、前記合成部の出力のうち前記合成部の出力の大きさが所定値未満であるものに基づいて、前記第１係数及び前記第２係数を算出する請求項１に記載の歪補償装置。
【請求項３】
入力信号に対し、前記入力信号に依存する第１歪補償係数に基づいて、歪補償を行う第１歪補償部と、
前記第１歪補償部の出力を、前記歪補償部の出力の振幅に基づく位相差を有する定振幅の第１信号及び第２信号に分離する信号生成部と、
前記第１信号の高周波成分を除去する第１フィルタと、
前記第１フィルタの出力に対し、前記第１フィルタの出力の振幅に依存する第２歪補償係数に基づいて、歪補償を行う第２歪補償部と、
前記第２歪補償部の出力を増幅する第１増幅部と、
前記第２信号の高周波成分を除去する第２フィルタと、
前記第２フィルタの出力に対し、前記第２フィルタの出力の振幅に依存する第３歪補償係数に基づいて、歪補償を行う第３歪補償部と、
前記第３歪補償部の出力を増幅する第２増幅部と、
前記第１増幅部の出力と前記第２増幅部の出力とを合成する合成部と、
前記合成部の出力のうち前記合成部の出力の大きさが所定値以上であるものに基づいて前記歪補償係数を算出し、前記合成部の出力のうち前記合成部の出力の大きさが所定値未満であるものに基づいて、前記第１係数及び前記第２係数を算出する演算部と、
を備える歪補償装置。
【請求項４】
入力信号に対し、前記入力信号に依存する歪補償係数に基づいて、歪補償を行い、
前記歪補償を行った信号を、前記歪補償を行った信号の振幅に基づく位相差を有する定振幅の第１信号及び第２信号に分離し、
前記第１信号に対し、第１係数を乗算し、
前記第１係数を乗算された信号を増幅し、
前記第２信号に対し、第２係数を乗算し、
前記第２係数を乗算された信号を増幅し、
前記第１係数を乗算された信号を増幅した信号と、前記第２係数を乗算された信号を増幅した信号とを合成し、
前記歪補償を行った信号と、前記第１係数を乗算された信号と、前記第２係数を乗算された信号と、前記合成された信号とに基づいて、前記歪補償係数、前記第１係数、前記第
２係数を算出する、
歪補償方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、歪補償装置および歪補償方法に関する。
【背景技術】
【０００２】
高効率な線形増幅器を実現する手段の１つとして、Linear Amplification with Nonlinear Components（ＬＩＮＣ）による増幅器を用いた高周波増幅回路が知られている。
【０００３】
図１は、ＬＩＮＣによる増幅器の例を示す図である。ＬＩＮＣによる増幅器では、ＬＩＮＣ信号生成部が、入力変調信号Ｓｉｎ（ｔ）をその振幅に応じた位相差を有する位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）に分離して出力する。例えば、入力変調信号Ｓｉｎ（ｔ）は、振幅変調及び位相変調（角度変調）を伴う変調信号であり、位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）は、定包絡線となる定振幅位相変調信号である。ここでの入力変調信号Ｓｉｎと位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）とは、いずれもベースバンド信号であってもよいし、Intermediate Frequency（ＩＦ）信号であってもよい。ＬＩＮＣ信号生成部は、位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）をディジタル信号として出力する。
【０００４】
ここで、信号Ｓｉｎ（ｔ）、信号Ｓｃ１（ｔ）、信号Ｓｃ２（ｔ）は、例えば、次のように表される。
【０００５】
【数１】
2014011653.tif 000003
【０００６】
ａ（ｔ）は、入力変調信号Ｓｉｎ（ｔ）の振幅成分である。θ（ｔ）は、入力変調信号Ｓｉｎ（ｔ）の位相成分である。振幅ａ（ｔ）に応じた位相差２×ψ（ｔ）が生じるように位相変調が与えられる。ａｍａｘは、振幅ａ（ｔ）の最大値であり、定数である。信号Ｓｃ１（ｔ）、信号Ｓｃ２（ｔ）は、定包絡線の信号である。即ち、信号Ｓｃ１（ｔ）、信号Ｓｃ２（ｔ）の振幅は、一定である。
【０００７】
ＬＩＮＣ信号生成部から出力された位相変調信号対の一方Ｓｃ１（ｔ）は、ＤＡＣ（Digital to Analog Converter）でディジタル信号からアナログ信号に変換される。さらに、変換されたアナログ信号が、ローパスフィルタを通過することで、位相変調信号対の一方Ｓｃ１（ｔ）の周波数帯域に対応する成分が抽出されるとともに、それ以外の周波数成分が抑圧される。同様に、位相変調信号対の他方Ｓｃ２（ｔ）は、ＤＡＣでディジタル信号からアナログ信号に変換される。さらに、変換されたアナログ信号が、ローパスフィルタを通過することで、位相変調信号対の他方Ｓｃ２（ｔ）の周波数帯域に対応する成分が
抽出されるとともにそれ以外の周波数成分が抑圧される。
【０００８】
ＬＩＮＣによる増幅器では、直交変調器は、ローパスフィルタを通過した位相変調信号対の一方Ｓｃ１（ｔ）を直交変調する。周波数変換器は、発振器から出力された高周波信号（発振信号）を用いて、ＲＦ（Radio Frequency）信号である高周波信号対の一方Ｓ１（ｔ）を生成して出力する。同様に、直交変調器は、ローパスフィルタを通過した位相変調信号対の他方Ｓｃ２（ｔ）を直交変調する。周波数変換器は、発振器から出力された高周波信号を用いて、ＲＦ信号である高周波信号対の他方Ｓ２（ｔ）を生成して出力する。
【０００９】
高周波信号Ｓ１（ｔ）及び高周波信号Ｓ２（ｔ）は、次のように表される。ここで、無線周波数（発振器の周波数）は、ｆｃとする。
【００１０】
【数２】
2014011653.tif 000004
【００１１】
増幅器対は、並列に設けられた２つの増幅器を含む。２つの増幅器の利得および位相特性は、ほぼ同一である。増幅器は、周波数変換器から出力された高周波信号を増幅する。ＤＡＣから増幅器（増幅器Ａ又は増幅器Ｂ）までを１つのブランチともいう。
【００１２】
合成器は、増幅器対により増幅された高周波信号対を合成し、合成後の信号を高周波信号Ｓｏｕｔ（ｔ）として、出力する。合成器から出力される信号Ｓｏｕｔ（ｔ）は、増幅器の利得をＧとすると、次のように表される。
【００１３】
【数３】
2014011653.tif 000005
【００１４】
ここで、φは、高周波信号対Ｓ１（ｔ）、Ｓ２（ｔ）の通過位相である。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】特表２００９−５３３９４７号公報
【特許文献２】特開２００３−２９８３６１号公報
【特許文献３】特開２００３−１５２４６４号公報
【特許文献４】特開平５−３７２６３号公報
【特許文献５】特開平９−７４３２０号公報
【特許文献６】特開２００６−３３９８８号公報
【発明の概要】
【発明が解決しようとする課題】
【００１６】
ＬＩＮＣによる増幅器を用いた高周波増幅回路に、入力変調信号としてＰＳＫ信号等のキャリア極性が反転する信号を入力すると、ＬＩＮＣ信号生成部で生成される位相変調信号対Ｓｃ１（ｔ）、Ｓｃ２（ｔ）には、位相が１８０度反転する点が存在し、信号の周波数帯域幅が広くなる。
【００１７】
図２は、入力変調信号が２トーンの場合のディジタル信号系列コンステレーションの例を示す図である。図２は、例えば、信号Ｓｃ１を表す。図２の例では、信号は定包絡線の信号である。図２の例では、（Ｉ，Ｑ）＝（０，−１）の点と（Ｉ，Ｑ）＝（０，１）の点との間で、位相が１８０度反転している。
【００１８】
しかし、ディジタル信号の位相変調信号対Ｓｃ１（ｔ）、Ｓｃ２（ｔ）はナイキストの定理により、サンプリング周波数の半分の周波数までしか表現できない。このため、ＤＡＣでアナログ信号に変換され、ローパスフィルタで折り返し成分を除去された信号は、大きなリンギングを生じ、定包絡線の信号とは異なる信号となる。
【００１９】
図３は、図２の信号がＤＡＣでアナログ信号に変換され、ローパスフィルタで高周波成分が除去された後のコンステレーションの例を示す図である。図３は、例えば、信号Ｓｃ１がＤＡＣでアナログ変換された後の信号（アナログ信号）である。
【００２０】
このアナログ信号の振幅成分は、リンギングにより変化する。即ち、このアナログ信号の振幅成分は、一定ではない。このアナログ信号を増幅器対で増幅する際に、このアナログ信号は増幅器対のＡＭ／ＡＭ特性やＡＭ／ＰＭ特性（ＡＭ／ＡＭ歪やＡＭ／ＰＭ歪）の影響を受ける。これらの歪により、合成後の出力高周波信号Ｓｏｕｔ（ｔ）が劣化する。即ち、出力高周波信号Ｓｏｕｔ（ｔ）に、歪が発生する。ＡＭ／ＡＭ特性は、入力信号の振幅に対する出力信号の振幅を表す。ＡＭ／ＡＭ歪は、ＡＭ／ＡＭ特性による歪である。ＡＭ／ＰＭ特性は、入力信号の振幅に対する出力信号の位相回転を表す。ＡＭ／ＰＭ歪は、ＡＭ／ＰＭ特性による歪である。
【００２１】
このアナログ信号の振幅成分は、時間変化するが、ＤＡＣに入力されるディジタル信号は定包絡線信号である。即ち、ＤＡＣに入力されるディジタル信号の振幅成分は時間変化しない。このため、増幅器の非線形特性を補償するために、ＤＡＣに入力されるディジタル信号にディジタルプリディストーション処理を行うことは、困難である。
【００２２】
図４は、図１の増幅器Ａの出力のコンステレーションの例である。図４は、図３のような信号が、増幅器Ａに入力された時の出力を表す。増幅器Ａの出力は、ＡＭ／ＡＭ歪及びＡＭ／ＰＭ歪の影響を受けている。
【００２３】
図５は、図１の増幅器Ｂの出力のコンステレーションの例である。図５は、図３と対になる信号が、増幅器Ｂに入力された時の出力を表す。増幅器Ｂの出力は、ＡＭ／ＡＭ歪及びＡＭ／ＰＭ歪の影響を受けている。
【００２４】
図６は、図４の信号と図５の信号とを合成した合成器の出力のコンステレーションの例を示す図である。増幅器における歪の影響がなければ、合成器の出力は、（Ｉ，Ｑ）＝（−２，０）の点と（Ｉ，Ｑ）＝（２，０）の点とを結ぶ直線となる。しかし、図６のように、合成器の出力は、ＡＭ／ＡＭ歪及びＡＭ／ＰＭ歪により、原点を中心として回転して
いる。
【００２５】
増幅器で使用される代表的な合成器は、線形性の優れたウィルキンソン型電力合成器である。また、電力増幅器である場合、効率を高めるためにChireix型電力合成器が使用されることがある。Chireix型電力合成器が使用されると、線形性は低下するが、効率は高くなる。
【００２６】
合成器として、ＬＩＮＣ型増幅器においてChireix型電力合成器のような線形性が低い合成器が使用されると、ブランチ間の位相差と合成電力の線形性が低下し、高周波出力に歪が生じる。
【００２７】
また、ブランチ間アンバランスの補償を合成器出力信号を使用して行う場合、合成器の非線形性が影響して、ブランチ間アンバランスの補償誤差が発生し、高周波出力が劣化することがある。
【００２８】
本件開示の技術は、合成後の出力高周波信号の劣化を抑制する歪補償装置を提供することを課題とする。
【課題を解決するための手段】
【００２９】
開示の技術は、上記課題を解決するために、以下の手段を採用する。
【００３０】
即ち、第１の態様は、
入力信号に対し、前記入力信号に依存する歪補償係数に基づいて、歪補償を行う歪補償部と、
前記歪補償部の出力を、前記歪補償部の出力の振幅に基づく位相差を有する定振幅の第１信号及び第２信号に分離する信号生成部と、
前記第１信号に対し、第１係数を乗算する第１調整部と、
前記第１調整部の出力を増幅する第１増幅部と、
前記第２信号に対し、第２係数を乗算する第２調整部と、
前記第２調整部の出力を増幅する第２増幅部と、
前記第１増幅部の出力と前記第２増幅部の出力とを合成する合成部と、
前記歪補償部の出力と、前記第１調整部の出力と、前記第２調整部の出力と、前記合成部の出力とに基づいて、前記歪補償係数、前記第１係数、前記第２係数を算出する演算部と、
を備える歪補償装置である。
【００３１】
開示の態様は、プログラムが情報処理装置によって実行されることによって実現されてもよい。即ち、開示の構成は、上記した態様における各手段が実行する処理を、情報処理装置に対して実行させるためのプログラム、或いは当該プログラムを記録した記録媒体として特定することができる。また、開示の構成は、上記した各手段が実行する処理を情報処理装置が実行する方法をもって特定されてもよい。
【発明の効果】
【００３２】
開示の技術によれば、合成後の出力高周波信号の劣化を抑制する歪補償装置を提供することができる。
【図面の簡単な説明】
【００３３】
【図１】図１は、ＬＩＮＣによる増幅器の例を示す図である。
【図２】図２は、入力変調信号が２トーンの場合のディジタル信号系列コンステレーションの例を示す図である。
【図３】図３は、図２の信号がＤＡＣでアナログ信号に変換された後のコンステレーションの例を示す図である。
【図４】図４は、図１の増幅器Ａの出力のコンステレーションの例である。
【図５】図５は、図１の増幅器Ｂの出力のコンステレーションの例である。
【図６】図６は、図４の信号と図５の信号とを合成した合成器の出力のコンステレーションの例を示す図である。
【図７】図７は、実施形態１の歪補償装置の構成例を示す図である。
【図８】図８は、歪補償装置１００の動作フローの例を示す図である。
【図９】図９は、実施形態１の歪補償装置のハードウェア構成例を示す図である。
【図１０】図１０は、実施形態２の歪補償装置の構成例を示す図である。
【図１１】図１１は、歪補償装置２００の動作フローの例を示す図である。
【発明を実施するための形態】
【００３４】
以下、図面を参照して実施形態について説明する。実施形態の構成は例示であり、開示の構成は、開示の実施形態の具体的構成に限定されない。開示の構成の実施にあたって、実施形態に応じた具体的構成が適宜採用されてもよい。各実施形態は、可能な限り組み合わせて実施され得る。
【００３５】
〔実施形態１〕
（構成例）
図７は、実施形態１の歪補償装置の構成例を示す図である。歪補償装置１００は、ＤＰＤ１０２、ＬＩＮＣ信号生成部１０４、位相振幅調整部１１４、ＤＡＣ１１６、ローパスフィルタ１１８、直交変調器１２０、周波数変換器１２２、増幅器１２４、逆数演算部１２６を含む。また、歪補償装置１００は、位相振幅調整部ＤＰＤ１３４、ＤＡＣ１３６、ローパスフィルタ１３８、直交変調器１４０、周波数変換器１４２、増幅器１４４、逆数演算部１４６を含む。さらに、歪補償装置１００は、合成器１５２、周波数変換器１５４、直交復調器１５６、ＡＤＣ１５８、モデリング部１８０、逆特性演算部１８２を含む。
【００３６】
ＤＰＤ（Digital PreDistortion）１０２は、入力変調信号Ｓｉｎ（ｔ）に対し、歪補償処理を行う。ここでの歪補償処理は、ＬＩＮＣによる増幅器で増幅する前の信号に対して、プリディストーション係数を乗算してあらかじめ歪を発生させ、ＬＩＮＣによる増幅器で生じる歪を打ち消す処理である。プリディストーション係数は、増幅器の歪を補償する歪補償係数である。入力変調信号に乗算されるプリディストーション係数は、入力変調信号Ｓｉｎ（ｔ）に依存する。プリディストーション係数は、増幅器における出力信号の振幅が入力信号の振幅に比例するように、かつ、入力信号と出力信号との位相差が０になるように補償する係数である。入力信号、プリディストーション係数は、例えば、Ｉ（In-phase）成分及びＱ（Quadrature）成分で表される。プリディストーション係数は、逆特性演算部１８２により求められる。ＤＰＤ１０２は、処理した信号を、ＬＩＮＣ信号生成部１０４及びモデリング部１８０に出力する。入力変調信号Ｓｉｎ（ｔ）は、例えば、ベースバンド信号である。入力変調信号Ｓｉｎ（ｔ）は、中間周波数（ＩＦ： Intermediate Frequency）信号であってもよい。ＤＰＤ１０２は、歪補償部の１例である。プリディストーション係数は、入力信号の大きさに依存する。
【００３７】
ＬＩＮＣ信号生成部１０４は、ＤＰＤ１０２の出力をその振幅に応じた位相差を有する位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）に分離して出力する。信号Ｓｃ１（ｔ）は、位相振幅調整部１１４に出力される。信号Ｓｃ２（ｔ）は、位相振幅調整部１３４に出力される。
【００３８】
位相振幅調整部１１４は、逆数演算部１２６で演算される値を、ＬＩＮＣ信号生成部１０４が出力した信号Ｓｃ１（ｔ）に乗算する。ここでは、当該係数は、信号Ｓｃ１（ｔ）
の大きさに依存しない。また、当該係数は、時間に依存してもよい。位相振幅調整部１１４は、調整部の１例である。
【００３９】
ＤＡＣ（Digital to Analog Converter）１１６は、位相振幅調整部１１４の出力であるディジタル信号をアナログ信号に変換する。
【００４０】
ローパスフィルタ（ＬＰＦ）１１８は、ＤＡＣ１１６から出力された信号の高周波成分をカットする。
【００４１】
直交変調器１２０は、ローパスフィルタ１１８が出力した信号を直交変調して、出力する。
【００４２】
周波数変換器１２２は、発振器を含む。周波数変換器１２２は、直交変調器１２０で直交変調された信号を、無線（ＲＦ： Radio Frequency）周波数にアップコンバートする。周波数変換器１２２は、発振器から出力される高周波信号（発振信号）を用いて、高周波信号を生成して出力する。
【００４３】
増幅器１２４は、周波数変換器１２２から出力された信号を電力増幅する。増幅器１２４は、増幅した信号を、合成器１５２に出力する。
【００４４】
位相振幅調整部１３４、ＤＡＣ１３６は、それぞれ、位相振幅調整部１１４、ＤＡＣ１１６と同様である。ローパスフィルタ１３８、直交変調器１４０、周波数変換器１４２、増幅器１４４は、それぞれ、ローパスフィルタ１１８、直交変調器１２０、周波数変換器１２２、増幅器１２４と同様である。
【００４５】
合成器１５２は、増幅器１２４の出力と増幅器１４４の出力とを合成し、Ｓｏｕt（ｔ）として出力する。合成器１５２から出力された信号は、アンテナ等から送信される。
【００４６】
周波数変換器１５４は、合成器１５２の出力信号を、無線周波数からベースバンド周波数にダウンコンバートする。直交復調器１５６は、周波数変換器１５４の出力を、同相信号および直交信号に復調する。ＡＤＣ１５８は、直交復調器１５６が出力する信号を、アナログ信号からディジタル信号に変換する。
【００４７】
モデリング部１８０は、ＤＰＤ１０２が出力したディジタル信号、位相振幅調整部１１４が出力したディジタル信号、位相振幅調整部１３４が出力したディジタル信号、ＡＤＣ１５８が出力したディジタル信号を受信する。モデリング部１８０は、受信したディジタル信号が所定のサンプル数に達するまで、ディジタル信号を受信する。モデリング部１８０は、所定のサンプル数の当該ディジタル信号を受信する毎に、合成器１５２の出力を模擬する。
【００４８】
ここで、ＤＰＤ１０２が出力した信号をｘｉ、位相振幅調整部１１４が出力したディジタル信号をｐｉ、位相振幅調整部１３４が出力したディジタル信号をｑｉ、ＡＤＣ１５８が出力したディジタル信号ｒｉとする。添字のｉは、モデリング部１８０がｉ番目に受信したディジタル信号であることを示す。ここでは、所定のサンプル数をＮとする。Ｎとして、例えば、Ｎ＝１００、Ｎ＝１０００を取りうる。Ｎの値はこれらに限定されるものではない。
【００４９】
モデリング部１８０は、次の式で示されるε１が０になるように（ε１の絶対値が最小になるように）、係数ｃｎの値を決定する。係数ｃｎを決定する方法として、どのような方法が用いられてもよい。モデリング部１８０は、所定のサンプル数Ｎ毎に、係数ｃｎを
算出する。ｘｉは、ｒｉに対応する入力信号である。
【００５０】
【数４】
2014011653.tif 000006
【００５１】
ここで、Ｍは、定数であり、自然数である。Ｍが大きいほど、合成器１５２の出力を精度よく模擬する。Ｍとして、例えば、Ｍ＝１、Ｍ＝２を取りうる。Ｍの値はこれらに限定されるものではない。関数ｆは、１次以上の項を表す。ｆ（ｘi）は、モデリング部１８０が合成器１５２のｉ番目の出力を模擬したものである。
【００５２】
さらに、モデリング部１８０は、次の式で示されるε２が０になるように（ε２の絶対値が最小になるように）、係数ａ、係数ｂの値を決定する。係数ａ、係数ｂを決定する方法として、どのような方法が用いられてもよい。モデリング部１８０は、所定のサンプル数Ｎ毎に、係数ａ、係数ｂを算出する。
【００５３】
【数５】
2014011653.tif 000007
【００５４】
モデリング部１８０は、算出した係数ａを逆数演算部１２６に出力する。モデリング部１８０は、算出した係数ｂを逆数演算部１４６に出力する。モデリング部１８０は、算出した係数ｃｎを逆特性演算部１８２に出力する。
【００５５】
逆特性演算部１８２は、プリディストーション係数を算出する。逆特性演算部１８２は、モデリング部１８０が出力する関数ｆに基づいて、プリディストーション係数を算出する。逆特性演算部１８２は、算出したプリディストーション係数をＤＰＤ１０２に出力する。逆特性演算部１８２は、入力信号とプリディストーション係数との対応関係を示すテーブルを、所定時間毎にＤＰＤ１０２に送信してもよい。ＤＰＤ１０２は、当該テーブルが逆特性演算部１８２から送信された場合、当該テーブルを格納する。
【００５６】
逆数演算部１２６は、係数ａの逆数（１／ａ）を算出する。逆数演算部１２６は、算出した値を位相振幅調整部１１４に出力する。
【００５７】
逆数演算部１４６は、係数ｂの逆数（１／ｂ）を算出する。逆数演算部１４６は、算出した値を位相振幅調整部１３４に出力する。
【００５８】
係数ａ及び係数ｂによって、ＬＩＮＣによる増幅器のブランチ間のアンバランスが補正される。
【００５９】
モデリング部１８０は、係数ａ及び係数ｂを逆数演算部１４６に出力してもよい。このとき、逆数演算部１４６は、係数ｂの逆数（１／ｂ）を算出する。さらに、逆数演算部１４６は、算出した逆数を、係数ａの逆数（１／ａ）で正規化した値を算出する。即ち、逆数演算部１４６が算出する値は、ａ／ｂである。逆数演算部１４６は、算出した値を位相振幅調整部１３４に出力する。このとき、逆数演算部１２６は、値１を位相振幅調整部１１４に出力する。この値１は、係数ａの逆数（１／ａ）を、係数ａの逆数（１／ａ）で正規化した値である。
【００６０】
また、モデリング部１８０は、係数ａ及び係数ｂを逆数演算部１２６に出力してもよい。このとき、逆数演算部１２６は、係数ａの逆数（１／ａ）を算出する。さらに、逆数演算部１２６は、算出した逆数を、係数ｂの逆数（１／ｂ）で正規化した値を算出する。即ち、逆数演算部１２６が算出する値は、ｂ／ａである。逆数演算部１２６は、算出した値を位相振幅調整部１１４に出力する。このとき、逆数演算部１４６は、値１を位相振幅調整部１３４に出力する。この値１は、係数ｂの逆数（１／ｂ）を、係数ｂの逆数（１／ｂ）で正規化した値である。
【００６１】
逆数演算部１２６、逆数演算部１４６、逆特性演算部１８２は、モデリング部１８０に含まれてもよい。即ち、モデリング部１８０が、逆数演算部１２６、逆数演算部１４６、逆特性演算部２８２として動作してもよい。
【００６２】
（歪補償装置の動作）
歪補償装置１００の動作について説明する。
【００６３】
図８は、歪補償装置１００の動作フローの例を示す図である。
【００６４】
歪補償装置１００のＤＰＤ１０２は、送信対象のディジタル信号Ｓｉｎ（ｔ）を受信する。ここで使用される信号は、複素信号（complex信号）である。ＤＰＤ１０２は、ディジタル信号Ｓｉｎ（ｔ）に、プリディストーション係数を乗算して出力する（Ｓ１０１）。プリディストーション係数は、ＬＩＮＣ信号生成部１０４から合成器１５２までで発生する歪を補償する係数である。プリディストーション係数は、逆特性演算部１８２によって算出される。ＤＰＤ１０２では、プリディストーション係数によって、入力信号に対し、ＬＩＮＣによる増幅器で与えられる歪特性と逆の特性が与えられる。プリディストーション係数は、入力信号に依存する。ＤＰＤ１０２は、逆特性演算部１８２が算出した、信号とプリディストーション係数との対応関係を、テーブルとして有してもよい。ＤＰＤ１０２は、入力信号に、プリディストーション係数を乗算した信号を、ＬＩＮＣ信号生成部１０４に出力する。
【００６５】
ＬＩＮＣ信号生成部１０４は、ＤＰＤ１０２が出力した信号を受信する。ＬＩＮＣ信号生成部１０４は、ディジタル信号を、当該ディジタル信号の振幅に応じた位相差を有する位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）に分離する（Ｓ１０２）。ＬＩＮＣ信号生成部１０２は、信号Ｓｃ１（ｔ）を、位相振幅調整部１１４に出力する。ＬＩＮＣ信号生成部１０２は、信号Ｓｃ２（ｔ）を、位相振幅調整部１３４に出力する。ＬＩＮＣ信号生成部１０２から出力される信号は、例えば、振幅および位相により表現される。また、ＬＩＮＣ信号生成部１０４から出力される信号は、Ｉ（In-phase）成分およびＱ（Quadrature phase）成分により表現されてもよい。
【００６６】
位相振幅調整部１１４は、ＬＩＮＣ信号生成部１０４が出力したディジタル信号に、所定の係数を乗算して出力する（Ｓ１０３）。所定の係数は、逆数演算部１２６によって、算出される。位相振幅調整部１１４では、所定の係数によって、信号の位相及び振幅が調整される。所定の係数によって、増幅器１２４の出力と増幅器１４４の出力とがアンバランスにならないように調整される。
【００６７】
ＤＡＣ１１６は、位相振幅調整部１１４が出力したディジタル信号を、アナログ信号に変換する（Ｓ１０４）。
【００６８】
ローパスフィルタ１１８は、ＤＡＣ１１６で変換されたアナログ信号の高周波成分をカットする（Ｓ１０５）。高周波成分がカットされた信号は、直交変調器１２０に出力される。
【００６９】
直交変調器１２０は、ローパスフィルタ１１８が出力した信号を、直交変調する。周波数変換器１２２は、発振器から出力される高周波信号を用いて、直交変調器から出力される信号の周波数を無線周波数に変換して出力する（Ｓ１０６）。
【００７０】
増幅器１２４は、周波数変換器１２２が出力した信号を、増幅する（Ｓ１０７）。増幅器１２４で増幅される信号は、あらかじめ、ＤＰＤ１０２等で、歪補償処理がされている。
【００７１】
ＬＩＮＣ信号生成部１０２が出力した信号Ｓｃ２（ｔ）は、信号Ｓｃ１（ｔ）と同様に、位相振幅調整部１３４、ＤＡＣ１３６、ローパスフィルタ１３８、直交変調器１４０、周波数変換器１４２、増幅器１４４で処理される。
【００７２】
合成器１５２は、増幅器１２４が出力した信号と、増幅器１４４が出力した信号とを、合成して、出力する（Ｓ１０８）。出力された信号は、アンテナ等により、他の装置に向けて送信される。また、出力された信号の一部は、周波数変換器１５４に入力される。
【００７３】
周波数変換器１５４は、合成器１５２の出力信号を、無線周波数からベースバンド周波数にダウンコンバートする。直交復調器１５６は、周波数変換器１５４の出力を、同相信号および直交信号に復調する。ＡＤＣ１５８は、直交復調器１５６が出力する信号を、ディジタル信号に変換する。
【００７４】
モデリング部１８０は、ＤＰＤ１０２が出力したディジタル信号、位相振幅調整部１１４が出力したディジタル信号、位相振幅調整部１３４が出力したディジタル信号、ＡＤＣ１５８が出力したディジタル信号を受信する。モデリング部１８０は、所定のサンプル数の当該ディジタル信号を受信する毎に、合成器１５２の出力を模擬する関数ｆ（ｘ）を算出する。さらに、モデリング部１８０は、ＬＩＮＣによる増幅器のブランチ間のアンバランスを調整する係数ａ及び係数ｂを算出する（Ｓ１０９）。モデリング部１８０は、算出した係数ａを逆数演算部１２６に出力する。モデリング部１８０は、算出した係数ｂを逆数演算部１４６に出力する。モデリング部１８０は、算出した係数ｃｎ（関数ｆ（ｘ））を逆特性演算部１８２に出力する。
【００７５】
逆数演算部１２６は、モデリング部１８０が出力した係数ａの逆数を算出し、位相振幅調整部１１４に出力する。逆数演算部１４６は、モデリング部１８０が出力した係数ｂの逆数を算出し、位相振幅調整部１３４に出力する。
【００７６】
逆特性演算部１８２は、プリディストーション係数を算出する（Ｓ１０９）。逆特性演算部１８２は、モデリング部１８０が出力した係数Ｃｎ（関数ｆ（ｘ））を受信する。逆
特性演算部１８２は、係数Ｃｎに基づいて、ＡＭ／ＡＭ歪及びＡＭ／ＰＭ歪を補償するプリディストーション係数を算出する。逆特性演算部１８２は、算出したプリディストーション係数をＤＰＤ１０２に出力する。
【００７７】
逆特性演算部１８２は、例えば、次のようにプリディストーション係数を算出する。逆特性演算部１８２は、関数ｆ（ｘ）を係数Ｃｎの和で、正規化する。関数ｆ（ｘ）を係数Ｃｎの和で正規化した関数を、関数ｆＮ（ｘ）とする。
【００７８】
【数６】
2014011653.tif 000008
【００７９】
逆特性演算部１８２は、関数ｆを正規化した関数ｆＮを用いて、ＬＩＮＣによる増幅器の逆特性を求め、プリディストーション係数を算出する。ＤＰＤ１０２は、逆特性演算部１８２が算出したプリディストーション係数を使用して、入力信号に対して、歪補償を行う。
【００８０】
また、プリディストーション係数（歪補償係数）の算出方法として、例えば、ＬＭＳ（Least Mean Square最小二乗法）アルゴリズムによるもの、指数重み付きＲＬＳ（Recursive Least Mean）アルゴリズムによるもの等が採用され得る。プリディストーション係数の算出方法として採用されるアルゴリズムは、これらに限定されるものではない。プリディストーション係数の更新は、例えば、所定時間毎に行われる。
【００８１】
合成器、増幅器の特性は、経年変化、動作温度、環境温度、入力信号等によって、変化し得る。よって、プリディストーション係数の更新を所定時間毎に行うことで、合成器、増幅器の特性の変化にプリディストーション係数が追随するようにできる。プリディストーション係数の更新（算出）は、この動作フローから独立して、個別に実行されてもよい。
【００８２】
一連の処理は、ハードウェアにより実行させることもできるが、ソフトウェアにより実行させることもできる。
【００８３】
プログラムを記述するステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくても、並列的または個別に実行される処理を含む。
【００８４】
（実施形態１の歪補償装置のハードウェア構成例）
図９は、実施形態１の歪補償装置のハードウェア構成例を示す図である。歪補償装置１０００は、入力されたディジタル信号を無線周波数にアップコンバートし、増幅し、出力する。歪補償装置１０００は、プロセッサ１００２、記憶装置１００４、ＤＡＣ１１０２、ＬＰＦ１１０４、アップコンバータ１１０６、増幅器１１０８を含む。歪補償装置１０００は、さらに、ＤＡＣ１２０２、ＬＰＦ１２０４、アップコンバータ１２０６、増幅器１２０８、合成器１０１２、アンテナ１０１４を含む。歪補償装置１００は、例えば、歪補償装置１０００のようなハードウェア構成によって実現される。
【００８５】
プロセッサ１００２は、例えば、Central Processing Unit（ＣＰＵ）やDigital Signal Processor（ＤＳＰ）である。プロセッサ１００２は、歪補償装置１０００の全体を制御する。プロセッサ１００２として、特定用途向け集積回路（ＡＳＩＣ: Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）が使用されてもよい。
【００８６】
記憶装置１００４は、例えば、Random Access Memory（ＲＡＭ）やRead Only Memory（ＲＯＭ）である。また、記憶装置１００４は、例えば、Erasable Programmable Read Only Memory（ＥＰＲＯＭ）、ハードディスクドライブ（ＨＤＤ：Hard Disk Drive）である。また、二次記憶装置は、リムーバブルメディア、即ち可搬記録媒体を含むことができる。リムーバブルメディアは、例えば、Universal Serial Bus（ＵＳＢ）メモリ、あるいは、Compact Disk（ＣＤ）やDigital Versatile Disk（ＤＶＤ）のようなディスク記録媒体である。記憶装置１００４には、入力信号とプリディストーション係数との対応関係を示す対応テーブル等が格納されうる。
【００８７】
歪補償装置１０００は、プロセッサ１００２が記憶装置１００４に記憶されるプログラムを実行することにより、ＤＰＤ１０２、ＬＩＮＣ信号生成部１０４、位相振幅調整部１１４、逆数演算部１２６、モデリング部１８０、逆特性演算部１８２等の機能を実現する。
【００８８】
ＤＡＣ１１０２は、プロセッサ１００２から出力されるディジタル信号を、アナログ信号に変換する。ＤＡＣ１１０２は、ＤＡＣ１１６の機能を実現する。
【００８９】
ＬＰＦ１１０４は、ＤＡＣ１１０２が出力するアナログ信号から、高周波成分を除去する。ＬＰＦ１１０４は、ローパスフィルタ１１８の機能を実現する。
【００９０】
アップコンバータ１１０６は、直交変調器１２０、周波数変換器１２２の機能を実現する。
【００９１】
増幅器１１０８は、アップコンバータ１１０６が出力するアナログ信号を、増幅する。増幅器１１０８として、様々な増幅器が使用されうるが、増幅器１１０８と増幅器１２０８とは同一の特性を有する増幅器が使用されることが好ましい。増幅器１１０８及び増幅器１２０８は、それぞれ、増幅器１２４および増幅器１４４の機能を実現する。
【００９２】
ＤＡＣ１２０２、ＬＰＦ１２０４は、アップコンバータ１２０６、増幅器１２０８は、それぞれ、ＤＡＣ１１０２、ＬＰＦ１１０４は、アップコンバータ１１０６、増幅器１１０８と同様の機能を有する。
【００９３】
合成器１０１２は、増幅器１１０８の出力と、増幅器１２０８の出力と合成する。合成器１０１２は、合成器１５２の機能を実現する。合成器１０１２として、例えば、Chireix型電力合成器が使用される。合成器１０１２は、Chireix型電力合成器に限定されるものではない。
【００９４】
アンテナ１０１４は、合成器１０１２で合成された信号を、他の装置に向けて送信する。
【００９５】
（変形例）
ここでは、モデリング部１８０の変形例について説明する。ここでは、係数Ｃｎ（関数ｆ（ｘ））を求める際に使用するｒｉは、所定の閾値以上であるとする。また、係数ａ及
び係数ｂを求める際に使用するｒｉは、所定の閾値未満であるとする。モデリング部１８０以外の部分については、上記の構成と同様である。
【００９６】
モデリング部１８０は、ＤＰＤ１０２が出力したディジタル信号、位相振幅調整部１１４が出力したディジタル信号、位相振幅調整部１３４が出力したディジタル信号、ＡＤＣ１５８が出力したディジタル信号を受信する。モデリング部１８０は、受信したディジタル信号が所定のサンプル数に達するまで、ディジタル信号を受信する。モデリング部１８０は、所定のサンプル数の当該ディジタル信号を受信する毎に、合成器１５２の出力を模擬する。
【００９７】
ここで、ＤＰＤ１０２が出力した信号をｘｉ、位相振幅調整部１１４が出力したディジタル信号をｐｉ、位相振幅調整部１３４が出力したディジタル信号をｑｉ、ＡＤＣ１５８が出力したディジタル信号ｒｉとする。添字のｉは、モデリング部１８０がｉ番目に受信したディジタル信号であることを示す。ここでは、所定のサンプル数をＮとする。Ｎとして、例えば、Ｎ＝１００、Ｎ＝１０００を取りうる。Ｎの値はこれらに限定されるものではない。ここで、ディジタル信号ｒｉ（ｉ＝１，・・・，Ｎ）のうち、ディジタル信号ｒｉの大きさが所定値Ｓ以上のものがＪ個、ディジタル信号ｒｉの大きさが所定値Ｓ未満のものがＫ個であるとする。このとき、Ｎ、Ｊ、Ｋの関係は、Ｎ＝Ｊ＋Ｋである。また、ディジタル信号ｒｉの大きさが所定値Ｓ以上のものをｒｊ（ｊ＝１，・・・，Ｊ）とする。ディジタル信号ｒｊに対応する、ＤＰＤ１０２が出力した信号をｘｊ、位相振幅調整部１１４が出力したディジタル信号をｐｊ、位相振幅調整部１３４が出力したディジタル信号をｑｊとする。さらに、ディジタル信号ｒｉの大きさが所定値Ｓ未満のものをｒｋ（ｋ＝１，・・・，Ｋ）とする。ディジタル信号ｒｋに対応する、ＤＰＤ１０２が出力した信号をｘｋ、位相振幅調整部１１４が出力したディジタル信号をｐｋ、位相振幅調整部１３４が出力したディジタル信号をｑｋとする。ディジタル信号ｒｉの大きさが所定値Ｓ未満であれば歪補償装置１００の出力信号の大きさが入力信号の大きさに比例するとみなされることを基準として、所定値Ｓがあらかじめ決定される。
【００９８】
モデリング部１８０は、次の式で示されるε１が０になるように（ε１の絶対値が最小になるように）、係数ｃｎの値を決定する。係数ｃｎを決定する方法として、どのような方法が用いられてもよい。モデリング部１８０は、所定のサンプル数Ｎ毎に、係数ｃｎを算出する。ｘｊは、ｒｊに対応する入力信号である。
【００９９】
【数７】
2014011653.tif 000009
【０１００】
ここで、Ｍは、定数であり、自然数である。Ｍが大きいほど、合成器１５２の出力を精度よく模擬する。Ｍとして、例えば、Ｍ＝１、Ｍ＝２を取りうる。Ｍの値はこれらに限定
されるものではない。関数ｆは、１次以上の項を表す。ｆ（ｘｊ）は、モデリング部１８０が合成器１５２の出力ｒｊを模擬したものである。ｆ（ｘ）は、所定値Ｓ以上であるディジタル信号ｒｊを使用して、算出される。歪補償装置１００の出力信号の大きさが入力信号の大きさに比例するとみなされないディジタル信号ｒｊが使用されて、係数ｃｎが算出される。
【０１０１】
さらに、モデリング部１８０は、次の式で示されるε２が０になるように（ε２の絶対値が最小になるように）、係数ａ、係数ｂの値を決定する。係数ａ、係数ｂを決定する方法として、どのような方法が用いられてもよい。モデリング部１８０は、所定のサンプル数Ｎ毎に、係数ａ、係数ｂを算出する。
【０１０２】
【数８】
2014011653.tif 000010
【０１０３】
モデリング部１８０は、算出した係数ａを逆数演算部１２６に出力する。モデリング部１８０は、算出した係数ｂを逆数演算部１４６に出力する。モデリング部１８０は、算出した係数ｃｎを逆特性演算部１８２に出力する。係数ａ及び係数ｂは、所定値Ｓ未満であるディジタル信号ｒｋを使用して、算出される。歪補償装置１００の出力信号の大きさが入力信号の大きさに比例するとみなされるディジタル信号ｒｋが使用されて、係数ａおよび係数ｂが算出される。
【０１０４】
出力信号は、ＬＩＮＣによる増幅器のブランチ間のアンバランス、および、合成器の非線形性によって劣化しうる。出力信号の大きさが小さい場合には、ＬＩＮＣによる増幅器のブランチ間のアンバランスによる影響を受けて、出力信号は劣化する。また、出力信号の大きさが大きい場合には、合成器の非線形性による影響を受けて、出力信号は劣化する。従って、変形例のモデリング部１８０は、出力信号の大きさが所定値以上のディジタル信号を用いて、合成器の非線形性を調整し、出力信号の大きさが所定値未満のディジタル信号を用いて、ＬＩＮＣによる増幅器のブランチ間のアンバランスを調整する。
【０１０５】
（実施形態１の作用、効果）
歪補償装置１００は、合成器１５２の出力を用いて、増幅器、合成器等の特性を、多項式を用いて、模擬する。歪補償装置１００は、ＬＩＮＣ信号に分離する前の入力信号に対して歪補償を行うことで、ＬＩＮＣによる増幅器における歪を補償する。ＬＩＮＣによる増幅器の歪を補償するプリディストーション係数（歪補償係数）は、出力信号に基づいて算出される。また、歪補償装置１００は、ＬＩＮＣ信号に分離した後の信号に対して、所定の係数を乗算することで、ＬＩＮＣによる増幅器のブランチ間のバランスを調整する。所定の係数は、出力信号と、ＬＩＮＣによる増幅器を模擬した関数ｆの３次以上の項との差に基づいて、算出される。歪補償装置１００は、合成器の特性と、増幅器の歪特性と、信号Ｓｃ１を処理する増幅器１２４等と信号Ｓｃ２を処理する増幅器１４４等とのアンバランスとを補償することができる。また、歪補償装置１００は、合成器１５２の出力の劣化を抑制する。
【０１０６】
歪補償装置２００は、ＬＩＮＣによる増幅器に対する歪補償と、ＬＩＮＣによる増幅器のブランチ間のアンバランスの調整とにより、高周波出力の劣化を抑制することができる
。
【０１０７】
歪補償装置１００は、合成器１５２による出力に基づいて歪を補償するので、合成器１５２として線形性の低い合成器を使用することができる。即ち、歪補償装置１００において、線形性が低いが、効率がよい合成器を使用することができる。
【０１０８】
合成器１５２、増幅器１２４、増幅器１４４の特性は、経年劣化、動作環境（温度等）等によって、変化し得る。歪補償装置１００は、ＤＰＤの出力と、位相振幅調整部の出力、合成器の出力とを使用して、プリディストーション係数を所定時間毎などに算出することにより、増幅器等の特性が変化した場合でも、より適切な歪補償を行うことができる。
【０１０９】
〔実施形態２〕
次に実施形態２について説明する。実施形態２は、実施形態１との共通点を有する。従って、主として相違点について説明し、共通点については、説明を省略する。
【０１１０】
（構成例）
図１０は、実施形態２の歪補償装置の構成例を示す図である。歪補償装置２００は、ＤＰＤ２０２、ＬＩＮＣ信号生成部２０４、第１ローパスフィルタ２１２、ＤＰＤ２１４、ＤＡＣ２１６、第２ローパスフィルタ２１８、直交変調器２２０、周波数変換器２２２、増幅器２２４、逆特性演算部２２６を含む。また、歪補償装置２００は、第１ローパスフィルタ２３２、ＤＰＤ２３４、ＤＡＣ２３６、第２ローパスフィルタ２３８、直交変調器２４０、周波数変換器２４２、増幅器２４４、逆特性演算部２４６を含む。歪補償装置３００は、さらに、合成器２５２、周波数変換器２５４、直交復調器２５６、ＡＤＣ２５８、モデリング部２８０、逆特性演算部２８２を含む。
【０１１１】
ＤＰＤ（Digital PreDistortion）２０２は、入力変調信号Ｓｉｎ（ｔ）に対し、歪補償処理を行う。ここでの歪補償処理は、ＬＩＮＣによる増幅器で増幅する前の信号に対して、プリディストーション係数を乗算してあらかじめ歪を発生させ、ＬＩＮＣによる増幅器で生じる歪を打ち消す処理である。プリディストーション係数は、増幅器の歪を補償する歪補償係数である。入力変調信号に乗算されるプリディストーション係数は、入力変調信号Ｓｉｎ（ｔ）に依存する。プリディストーション係数は、増幅器における出力信号の振幅が入力信号の振幅に比例するように、かつ、入力信号と出力信号との位相差が０になるように補償する係数である。プリディストーション係数は、逆特性演算部２８２により求められる。ＤＰＤ２０２は、処理した信号を、ＬＩＮＣ信号生成部２０４及びモデリング部２８０に出力する。入力変調信号Ｓｉｎ（ｔ）は、例えば、ベースバンド信号である。入力変調信号Ｓｉｎ（ｔ）は、中間周波数（ＩＦ： Intermediate Frequency）信号であってもよい。ＤＰＤ２０２は、歪補償部の１例である。
【０１１２】
ＬＩＮＣ信号生成部２０４は、ＤＰＤ２０２の出力をその振幅に応じた位相差を有する位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）に分離して出力する。信号Ｓｃ１（ｔ）は、第１ローパスフィルタ２１２に出力される。信号Ｓｃ２（ｔ）は、第１ローパスフィルタ２３２に出力される。
【０１１３】
第１ローパスフィルタ（ＬＰＦ： Low Pass Filter）２１２は、ＬＩＮＣ信号生成部２０４から出力された信号の高周波成分をカットする。
【０１１４】
ＤＰＤ（Digital PreDistortion）２１４は、第１ローパスフィルタ２１２の出力に対し、歪補償処理を行う。ここでの歪補償処理は、増幅器２２４で増幅する前の信号に対して、プリディストーション係数を乗算してあらかじめ歪を発生させ、増幅器２２４で生じる歪を打ち消す処理である。プリディストーション係数は、増幅器の歪を補償する歪補償
係数である。第１ローパスフィルタ２１２の出力に乗算されるプリディストーション係数は、第１ローパスフィルタ２１２の出力に依存する。プリディストーション係数は、出力信号の振幅が入力信号の振幅に比例するように、かつ、入力信号と出力信号との位相差が０になるように補償する係数である。プリディストーション係数は、逆特性演算部２２６により算出される。ＤＰＤ２１４は、歪補償部の１例である。
【０１１５】
ＤＡＣ（Digital to Analog Converter）２１６は、ＤＰＤ部２１４の出力であるディジタル信号をアナログ信号に変換する。
【０１１６】
第２ローパスフィルタ（ＬＰＦ）２１８は、ＤＡＣ２１６から出力された信号の高周波成分をカットする。
【０１１７】
直交変調器２２０は、ローパスフィルタ２１８が出力した信号を直交変調して、出力する。
【０１１８】
周波数変換器２２２は、発振器を含む。周波数変換器２２２は、直交変調器２２０で直交変調された信号を、無線（ＲＦ： Radio Frequency）周波数にアップコンバートする。周波数変換器２２２は、発振器から出力される高周波信号（発振信号）を用いて、高周波信号を生成して出力する。
【０１１９】
増幅器２２４は、周波数変換器２２２から出力された信号を電力増幅する。増幅器２２４は、増幅した信号を、合成器２５２に出力する。
【０１２０】
逆特性演算部２２６は、プリディストーション係数を算出する。逆特性演算部２２６は、モデリング部２８０が出力する関数ｇに基づいて、プリディストーション係数を算出する。関数ｇによって、ＤＰＤ２１４の出力と増幅器２２４の出力との関係が分かる。これを利用することにより、逆特性演算部２２６は、プリディストーション係数を算出する。逆特性演算部２２６は、算出したプリディストーション係数をＤＰＤ２１４に出力する。逆特性演算部２２６は、入力信号とプリディストーション係数との対応関係を示すテーブルを、所定時間毎にＤＰＤ２１４に送信してもよい。ＤＰＤ２１４は、当該テーブルが逆特性演算部２２６から送信された場合、当該テーブルを格納する。
【０１２１】
第１ローパスフィルタ２３２、ＤＰＤ２３４、ＤＡＣ２３６は、それぞれ、第１ローパスフィルタ２１２、ＤＰＤ２１４、ＤＡＣ２１６と同様である。第２ローパスフィルタ２３８、直交変調器２４０、周波数変換器２４２、増幅器２４４、逆特性演算部２４６は、それぞれ、第２ローパスフィルタ２１８、直交変調器２２０、周波数変換器２２２、増幅器２２４、逆特性演算部２２６と同様である。
【０１２２】
合成器２５２は、増幅器２２４の出力と増幅器２４４の出力とを合成し、Ｓｏｕt（ｔ）として出力する。合成器２５２から出力された信号は、アンテナ等から送信される。
【０１２３】
周波数変換器２５４は、合成器２５２の出力信号を、無線周波数からベースバンド周波数にダウンコンバートする。直交復調器２５６は、周波数変換器２５４の出力を、同相信号および直交信号に復調する。ＡＤＣ２５８は、直交復調器２５６が出力する信号を、アナログ信号からディジタル信号に変換する。
【０１２４】
モデリング部２８０は、ＤＰＤ２０２が出力したディジタル信号、ＤＰＤ２１４が出力したディジタル信号、ＤＰＤ２３４が出力したディジタル信号、ＡＤＣ２５８が出力したディジタル信号を受信する。モデリング部２８０は、受信したディジタル信号が所定のサンプル数に達するまで、ディジタル信号を受信する。モデリング部２８０は、所定のサン
プル数の当該ディジタル信号を受信する毎に、合成器２５２の出力を模擬する。
【０１２５】
ここで、ＤＰＤ２０２が出力した信号をｘｉ、位相振幅調整部２１４が出力したディジタル信号をｐｉ、位相振幅調整部２３４が出力したディジタル信号をｑｉ、ＡＤＣ２５８が出力したディジタル信号ｒｉとする。添字のｉは、モデリング部２８０がｉ番目に受信したディジタル信号であることを示す。ここでは、所定のサンプル数をＮとする。Ｎとして、例えば、Ｎ＝１００、Ｎ＝１０００を取りうる。Ｎの値はこれらに限定されるものではない。ここで、ディジタル信号ｒｉ（ｉ＝１，・・・，Ｎ）のうち、ディジタル信号ｒｉの大きさが所定値Ｓ以上のものがＪ個、ディジタル信号ｒｉの大きさが所定値Ｓ未満のものがＫ個であるとする。このとき、Ｎ、Ｊ、Ｋの関係は、Ｎ＝Ｊ＋Ｋである。また、ディジタル信号ｒｉの大きさが所定値Ｓ以上のものをｒｊ（ｊ＝１，・・・，Ｊ）とする。ディジタル信号ｒｊに対応する、ＤＰＤ２０２が出力した信号をｘｊ、位相振幅調整部２１４が出力したディジタル信号をｐｊ、位相振幅調整部２３４が出力したディジタル信号をｑｊとする。さらに、ディジタル信号ｒｉの大きさが所定値Ｓ未満のものをｒｋ（ｋ＝１，・・・，Ｋ）とする。ディジタル信号ｒｋに対応する、ＤＰＤ２０２が出力した信号をｘｋ、位相振幅調整部２１４が出力したディジタル信号をｐｋ、位相振幅調整部２３４が出力したディジタル信号をｑｋとする。ディジタル信号ｒｉの大きさが所定値Ｓ未満であれば歪補償装置２００の出力信号の大きさが入力信号の大きさに比例するとみなされることを基準として、所定値Ｓがあらかじめ決定される。
【０１２６】
モデリング部２８０は、次の式で示されるε１が０になるように（ε１の絶対値が最小になるように）、係数ｃｎの値を決定する。係数ｃｎを決定する方法として、どのような方法が用いられてもよい。モデリング部２８０は、所定のサンプル数Ｎ毎に、係数ｃｎを算出する。ｘｊは、ｒｊに対応する入力信号である。
【０１２７】
【数９】
2014011653.tif 000011
【０１２８】
ここで、Ｍは、定数であり、自然数である。Ｍが大きいほど、合成器１５２の出力を精度よく模擬する。Ｍとして、例えば、Ｍ＝１、Ｍ＝２を取りうる。Ｍの値はこれらに限定されるものではない。関数ｆは、１次以上の項を表す。ｆ（ｘｊ）は、モデリング部２８０が合成器２５２の出力ｒｊを模擬したものである。ｆ（ｘ）は、所定値Ｓ以上であるディジタル信号ｒｊを使用して、算出される。歪補償装置２００の出力信号の大きさが入力信号の大きさに比例するとみなされないディジタル信号ｒｊが使用されて、係数ｃｎが算出される。
【０１２９】
さらに、モデリング部２８０は、次の式で示されるε２が０になるように（ε２の絶対値が最小になるように）、関数ｇ（係数ｕｎ）、関数ｈ（係数ｖｎ）の値を決定する。関
数ｇ、関数ｈを決定する方法として、どのような方法が用いられてもよい。モデリング部２８０は、所定のサンプル数Ｎ毎に、関数ｇ、関数ｈを算出する。ここでは、関数ｇ、関数ｈは、増幅器の振幅、位相及び非線形特性を模擬する。
【０１３０】
【数１０】
2014011653.tif 000012
【０１３１】
ここで、Ｍ１は、定数であり、自然数である。Ｍ１が大きいほど、関数ｇ、関数ｈは、増幅器の振幅、位相および非線形特性を精度よく模擬する。Ｍ１として、例えば、Ｍ１＝１、Ｍ１＝２を取りうる。
【０１３２】
モデリング部２８０は、算出した係数ｕｎを逆数演算部２２６に出力する。モデリング部２８０は、算出した係数ｖｎを逆数演算部２４６に出力する。モデリング部２８０は、算出した係数ｃｎを逆特性演算部２８２に出力する。係数ｕｎおよび係数ｖｎは、所定値Ｓ未満であるディジタル信号ｒｋを使用して、算出される。歪補償装置２００の出力信号の大きさが入力信号の大きさに比例するとみなされるディジタル信号ｒｋが使用されて、係数ｕｎおよび係数ｖｎが算出される。
【０１３３】
逆特性演算部２８２は、プリディストーション係数を算出する。逆特性演算部２８２は、モデリング部２８０が出力する関数ｆに基づいて、プリディストーション係数を算出する。逆特性演算部２８２は、算出したプリディストーション係数をＤＰＤ２０２に出力する。逆特性演算部２８２は、入力信号とプリディストーション係数との対応関係を示すテーブルを、所定時間毎にＤＰＤ２０２に送信してもよい。ＤＰＤ２０２は、当該テーブルが逆特性演算部２８２から送信された場合、当該テーブルを格納する。
【０１３４】
逆特性演算部２８２は、例えば、次のようにプリディストーション係数を算出する。逆特性演算部２８２は、関数ｆ（ｘ）を係数Ｃｎの和で、正規化する。関数ｆ（ｘ）を係数Ｃｎの和で正規化した関数を、関数ｆＮ（ｘ）とする。
【０１３５】
【数１１】
2014011653.tif 000013
【０１３６】
逆特性演算部２８２は、関数ｆを正規化した関数ｆＮを用いて、ＬＩＮＣによる増幅器の逆特性を求め、プリディストーション係数を算出する。ＤＰＤ２０２は、逆特性演算部２８２が算出したプリディストーション係数を使用して、入力信号に対して、歪補償を行う。
【０１３７】
逆特性演算部２２６は、は、プリディストーション係数を算出する。逆特性演算部２２６は、モデリング部２８０が出力する関数ｇに基づいて、プリディストーション係数を算出する。逆特性演算部２２６は、算出したプリディストーション係数をＤＰＤ２１４に出力する。逆特性演算部２２６は、入力信号とプリディストーション係数との対応関係を示すテーブルを、所定時間毎にＤＰＤ２１４に送信してもよい。ＤＰＤ２１４は、当該テーブルが逆特性演算部２２６から送信された場合、当該テーブルを格納する。
【０１３８】
逆特性演算部２４６は、は、プリディストーション係数を算出する。逆特性演算部２４６は、モデリング部２８０が出力する関数ｈに基づいて、プリディストーション係数を算出する。逆特性演算部２４６は、算出したプリディストーション係数をＤＰＤ２３４に出力する。逆特性演算部２４６は、入力信号とプリディストーション係数との対応関係を示すテーブルを、所定時間毎にＤＰＤ２３４に送信してもよい。ＤＰＤ２３４は、当該テーブルが逆特性演算部２４６から送信された場合、当該テーブルを格納する。
【０１３９】
逆特性演算部２２６及び逆特性演算部２４６は、関数ｇ又は関数ｈの逆特性のゲインを基準として、各関数を正規化して、プリディストーション係数を算出してもよい。
【０１４０】
逆特性演算部２２６、逆特性演算部２４６、逆特性演算部２８２は、モデリング部２８０に含まれてもよい。即ち、モデリング部２８０が、逆特性演算部２２６、逆特性演算部２４６、逆特性演算部２８２として動作してもよい。
【０１４１】
（歪補償装置の動作）
歪補償装置２００の動作について説明する。
【０１４２】
図１１は、歪補償装置２００の動作フローの例を示す図である。
【０１４３】
歪補償装置２００のＤＰＤ２０２は、送信対象のディジタル信号Ｓｉｎ（ｔ）を受信する。ここで使用される信号は、複素信号（complex信号）である。ＤＰＤ２０２は、ディジタル信号Ｓｉｎ（ｔ）に、プリディストーション係数を乗算して出力する（Ｓ２０１）。プリディストーション係数は、ＬＩＮＣ信号生成部２０４から合成器２５２までで発生する歪を補償する係数である。プリディストーション係数は、逆特性演算部２８２によって算出される。ＤＰＤ２０２では、プリディストーション係数によって、入力信号に対し、ＬＩＮＣによる増幅器で与えられる歪特性と逆の特性が与えられる。プリディストーション係数は、入力信号に依存する。ＤＰＤ２０２は、逆特性演算部２８２が算出した、信
号とプリディストーション係数との対応関係を、テーブルとして有してもよい。ＤＰＤ２０２は、入力信号に、プリディストーション係数を乗算した信号を、ＬＩＮＣ信号生成部２０４に出力する。
【０１４４】
ＬＩＮＣ信号生成部２０４は、ＤＰＤ２０２が出力した信号を受信する。ＬＩＮＣ信号生成部２０４は、ディジタル信号を、当該ディジタル信号の振幅に応じた位相差を有する位相変調信号対Ｓｃ１（ｔ）及びＳｃ２（ｔ）に分離する（Ｓ２０２）。ＬＩＮＣ信号生成部２０２は、信号Ｓｃ１（ｔ）を、第１ローパスフィルタ２１２に出力する。ＬＩＮＣ信号生成部２０２は、信号Ｓｃ２（ｔ）を、第１ローパスフィルタ２３２に出力する。ＬＩＮＣ信号生成部１０２から出力される信号は、例えば、振幅および位相により表現される。また、ＬＩＮＣ信号生成部１０４から出力される信号は、Ｉ（In-phase）成分およびＱ（Quadrature phase）成分により表現されてもよい。
【０１４５】
第１ローパスフィルタ２１２は、信号Ｓｃ１（ｔ）の高周波成分をカットする（Ｓ２０３）。高周波成分をカットされた信号は、ＤＰＤ２１４に入力される。第１ローパスフィルタ２１２による帯域制限の周波数帯域は、第２ローパスフィルタ２１８によるアナログ信号の周波数帯域と同様もしくは、第２ローパスフィルタ２１８によるアナログ信号の周波数帯域より狭くする。第１ローパスフィルタ２１２を通過した信号には、リンギングが生じる。即ち、第１ローパスフィルタ２１２に入力される信号が定包絡線の信号であっても、第１ローパスフィルタ２１２から出力される信号の振幅成分は一定ではない。
【０１４６】
ＤＰＤ２１４は、第１ローパスフィルタ２１４が出力したディジタル信号に、所定の係数を乗算して出力する（Ｓ２０４）。所定の係数は、逆特性演算部２２６によって、算出される。ＤＰＤ２１４では、所定の係数によって、信号の位相及び振幅が調整される。所定の係数によって、増幅器２２４の出力と増幅器２４４の出力とがアンバランスにならないように調整される。
【０１４７】
ＤＡＣ２１６は、ＤＰＤ２１４が出力したディジタル信号を、アナログ信号に変換する（Ｓ２０５）。
【０１４８】
第２ローパスフィルタ２１８は、ＤＡＣ２１６で変換されたアナログ信号の高周波成分をカットする（Ｓ２０６）。高周波成分がカットされた信号は、直交変調器２２０に出力される。
【０１４９】
直交変調器２２０は、第２ローパスフィルタ２１８が出力した信号を、直交変調する。周波数変換器２２２は、発振器から出力される高周波信号を用いて、直交変調器から出力される信号の周波数を無線周波数に変換して出力する（Ｓ２０７）。
【０１５０】
増幅器２２４は、周波数変換器２２２が出力した信号を、増幅する（Ｓ２０８）。増幅器２２４で増幅される信号は、あらかじめ、ＤＰＤ２０２等で、歪補償処理がされている。
【０１５１】
ＬＩＮＣ信号生成部２０２が出力した信号Ｓｃ２（ｔ）は、信号Ｓｃ１（ｔ）と同様に、第１ローパスフィルタ２３２、ＤＰＤ２３４、ＤＡＣ２３６、第２ローパスフィルタ２３８、直交変調器２４０、周波数変換器２４２、増幅器２４４で処理される。
【０１５２】
合成器２５２は、増幅器２２４が出力した信号と、増幅器２４４が出力した信号とを、合成して、出力する（Ｓ２０８）。出力された信号は、アンテナ等により、他の装置に向けて送信される。また、出力された信号の一部は、周波数変換器２５４に入力される。
【０１５３】
周波数変換器２５４は、合成器２５２の出力信号を、無線周波数からベースバンド周波数にダウンコンバートする。直交復調器２５６は、周波数変換器２５４の出力を、同相信号および直交信号に復調する。ＡＤＣ２５８は、直交復調器２５６が出力する信号を、ディジタル信号に変換する。
【０１５４】
モデリング部２８０は、ＤＰＤ２０２が出力したディジタル信号、ＤＰＤ２１４が出力したディジタル信号、ＤＰＤ２３４が出力したディジタル信号、ＡＤＣ２５８が出力したディジタル信号を受信する。モデリング部２８０は、受信したディジタル信号が所定のサンプル数に達するまで、ディジタル信号を受信する。モデリング部２８０は、所定のサンプル数の当該ディジタル信号を受信する毎に、合成器２５２の出力を模擬する関数ｆ（ｘ）を算出する。さらに、モデリング部２８０は、ＬＩＮＣによる増幅器のブランチ間のアンバランス及び増幅器２２４及び増幅器２４４の非線形性を調整する関数ｇ（p）、ｈ（q）を算出する（Ｓ２１０）。モデリング部２８０は、算出した係数ｕｎ（関数ｇ（p））を逆特性演算部２２６に出力する。モデリング部２８０は、算出した係数ｖｎ（関数ｈ（q））を逆特性演算部２２６に出力する。モデリング部２８０は、算出した係数ｃｎ（関数ｆ（ｘ））を逆特性演算部２８２に出力する。
【０１５５】
逆特性演算部２８２は、プリディストーション係数を算出する。逆特性演算部２８２は、モデリング部２８０が出力した係数Ｃｎ（関数ｆ（ｘ））を受信する。逆特性演算部２８２は、係数Ｃｎに基づいて、ＡＭ／ＡＭ歪及びＡＭ／ＰＭ歪を補償するプリディストーション係数を算出する。逆特性演算部２８２は、算出したプリディストーション係数をＤＰＤ２０２に出力する。ＤＰＤ２０２は、逆特性演算部２８２が算出したプリディストーション係数を使用して、入力信号に対して、歪補償を行う。
【０１５６】
逆特性演算部２２６は、モデリング部２８０が出力した関数ｇの逆特性を算出し、ＤＰＤ２１４に出力する。逆特性演算部２４６は、モデリング部２８０が出力した関数ｈの逆特性を算出し、ＤＰＤ１３４に出力する。
【０１５７】
一連の処理は、ハードウェアにより実行させることもできるが、ソフトウェアにより実行させることもできる。
【０１５８】
プログラムを記述するステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくても、並列的または個別に実行される処理を含む。
【０１５９】
（実施形態２の歪補償装置のハードウェア構成例）
実施形態２の歪補償装置２００は、実施形態１の歪補償装置１０００と同様なハードウェア構成で実現される。ここでは、歪補償装置１０００を用いて、歪補償装置２００のハードウェア構成例を説明する。
【０１６０】
プロセッサ１００２が記憶装置１００４に記憶されるプログラムを実行することにより、ＤＰＤ２０２、ＬＩＮＣ信号生成部２０４、第１ローパスフィルタ２１２、位相振幅調整部２１４等の機能を実現する。さらに、プロセッサ１００２が記憶装置１００４に記憶されるプログラムを実行することにより、逆数演算部２２６、モデリング部２８０、逆特性演算部２８２等の機能を実現する。
【０１６１】
記憶装置１００２には、プロセッサ１００２で使用されるプリディストーション係数等の情報が格納される。
【０１６２】
ＤＡＣ１１０２は、プロセッサ１００２から出力されるディジタル信号を、アナログ信
号に変換する。ＤＡＣ１１０２は、ＤＡＣ２１６の機能を実現する。
【０１６３】
ＬＰＦ１１０４は、ＤＡＣ１１０２が出力するアナログ信号から、高周波成分を除去する。ＬＰＦ１１０４は、第２ローパスフィルタ２１８の機能を実現する。
【０１６４】
アップコンバータ１１０６は、直交変調器２２０、周波数変換器２２２の機能を実現する。
【０１６５】
増幅器１１０８は、アップコンバータ１１０６が出力するアナログ信号を、増幅する。増幅器１１０８として、様々な増幅器が使用されうるが、増幅器１１０８と増幅器１２０８とは同一の特性を有する増幅器が使用されることが好ましい。増幅器１１０８及び増幅器１２０８は、それぞれ、増幅器１２４および増幅器１４４の機能を実現する。
【０１６６】
ＤＡＣ１２０２、ＬＰＦ１２０４は、アップコンバータ１２０６、増幅器１２０８は、それぞれ、ＤＡＣ１１０２、ＬＰＦ１１０４は、アップコンバータ１１０６、増幅器１１０８と同様の機能を有する。
【０１６７】
合成器１０１２は、増幅器１１０８の出力と、増幅器１２０８の出力と合成する。合成器１０１２は、合成器１５２の機能を実現する。
【０１６８】
アンテナ１０１４は、合成器１０１２で合成された信号を、他の装置に向けて送信する。
【０１６９】
（実施形態２の作用、効果）
歪補償装置２００は、ＬＩＮＣ信号生成部２０４が生成した定包絡線のディジタル信号の高周波成分をカットする。高周波成分がカットされた信号には、リンギングが生じる。即ち、定包絡線の信号から高周波成分がカットされた信号の振幅成分は、一定ではない。歪補償装置２００は、高周波成分がカットされた信号に対し、増幅器の非線形歪を補償する歪補償処理を行う。歪補償装置２００は、リンギングが生じた信号に対して歪補償処理をすることで、リンギング及び増幅器の歪特性による出力信号の劣化を抑制することができる。また、歪補償装置２００は、ディジタル信号に対して周波数帯域制限を行う（高周波成分をカットする）ことで、ディジタル信号に対して、ディジタルプリディストーション処理（歪補償処理）を行うことができる。
【０１７０】
ＬＩＮＣによる増幅器の歪を補償するプリディストーション係数（歪補償係数）は、出力信号に基づいて算出される。増幅器２２４および増幅器２４４に対するプリディストーション係数（歪補償係数）は、出力信号とＬＩＮＣによる増幅器を模擬した関数ｆの３次以上の項との差に基づいて算出される。
【０１７１】
出力信号は、ＬＩＮＣによる増幅器のブランチ間のアンバランス、および、合成器の非線形性によって劣化しうる。出力信号の大きさが小さい場合には、ＬＩＮＣによる増幅器のブランチ間のアンバランスによる影響を受けて、出力信号は劣化する。また、出力信号の大きさが大きい場合には、合成器の非線形性による影響を受けて、出力信号は劣化する。従って、モデリング部２８０は、出力信号の大きさが所定値以上のディジタル信号を用いて、合成器の非線形性を調整し、出力信号の大きさが所定値未満のディジタル信号を用いて、ＬＩＮＣによる増幅器のブランチ間のアンバランス等を調整する。
【０１７２】
また、歪補償装置２００は、多項式を用いて増幅器２２４及び増幅器２４４の同定を行うことにより、増幅器２２４及び増幅器２４４の非線形特性を補償することができる。
【０１７３】
歪補償装置２００は、ＬＩＮＣによる増幅器に対する歪補償と、ＬＩＮＣによる増幅器の各ブランチの増幅器に対する歪補償とにより、高周波出力の劣化を抑制することができる。
【符号の説明】
【０１７４】
１００    歪補償装置
１０２      ＤＰＤ
１０４      ＬＩＮＣ信号生成部
１１４      位相振幅調整部
１１６      ＤＡＣ
１１８      ローパスフィルタ
１２０      直交変調器
１２２      周波数変換器
１２４      増幅器
１３４      位相振幅調整部
１３６      ＤＡＣ
１３８      ローパスフィルタ
１４０      直交変調器
１４２      周波数変換器
１４４      増幅器
１５２      合成器
１８０      モデリング部
１８２      逆特性演算部
１０００    歪補償装置
１００２      プロセッサ
１００４      記憶装置
１０１２      合成器
１０１４      アンテナ
１０２２      ダウンコンバータ
１０２４      ＡＤＣ
１１０２      ＤＡＣ
１１０４      ＬＰＦ
１１０６      アップコンバータ
１１０８      増幅器
１２０２      ＤＡＣ
１２０４      ＬＰＦ
１２０６      アップコンバータ
１２０８      増幅器
２００    歪補償装置
２０２      ＤＰＤ
２０４      ＬＩＮＣ信号生成部
２１２      第１ローパスフィルタ
２１４      ＤＰＤ
２１６      ＤＡＣ
２１８      第２ローパスフィルタ
２２０      直交変調器
２２２      周波数変換器
２２４      増幅器
２２６      逆特性演算部
２３２      第１ローパスフィルタ
２３４      ＤＰＤ
２３６      ＤＡＣ
２３８      第２ローパスフィルタ
２４０      直交変調器
２４２      周波数変換器
２４４      増幅器
２４６      逆特性演算部
２５２      合成器
２８０      モデリング部
２８２      逆特性演算部
【図１】
2014011653.tif 000014
【図２】
2014011653.tif 000015
【図３】
2014011653.tif 000016
【図４】
2014011653.tif 000017
【図５】
2014011653.tif 000018
【図６】
2014011653.tif 000019
【図７】
2014011653.tif 000020
【図８】
2014011653.tif 000021
【図９】
2014011653.tif 000022
【図１０】
2014011653.tif 000023
【図１１】
2014011653.tif 000024