(57)【要約】【課題】回路規模を容易に低減できる歪補償装置及び歪補償方法を提供する。【解決手段】歪補償装置は、一つの態様において、送信信号を増幅するパワーアンプの出力に応じて、複数の奇数次歪が重ねられた第１の周波数成分と複数の偶数次歪が重ねられた第２の周波数成分とを含むフィードバック信号を生成するフィードバック部と、前記生成されたフィードバック信号から前記第１の周波数成分及び前記第２の周波数成分をそれぞれ分離して信号処理する信号処理部とを有する。【選択図】図１
2018101861.tif 000002
【特許請求の範囲】
【請求項１】
送信信号を増幅するパワーアンプの出力に応じて、複数の奇数次歪が重ねられた第１の周波数成分と複数の偶数次歪が重ねられた第２の周波数成分とを含むフィードバック信号を生成するフィードバック部と、
前記生成されたフィードバック信号から前記第１の周波数成分及び前記第２の周波数成分をそれぞれ分離して信号処理する信号処理部と、
を有することを特徴とする歪補償装置。
【請求項２】
前記フィードバック部は、前記パワーアンプの出力に応じた信号を送信周波数から中間周波数にダウンコンバートしながら前記信号を前記中間周波数で折り返して前記フィードバック信号を生成する第１の周波数変換部を有する
ことを特徴とする請求項１に記載の歪補償装置。
【請求項３】
前記フィードバック部は、
前記第１の周波数変換部で生成されたフィードバック信号を前記中間周波数からベースバンド周波数へダウンコンバートする第２の周波数変換部と、
前記第２の周波数変換部でベースバンド周波数へダウンコンバートされたフィードバック信号をアナログ形式からデジタル形式に変換するＡＤ変換部と、
をさらに有する
ことを特徴とする請求項２に記載の歪補償装置。
【請求項４】
前記信号処理部は、
前記フィードバック部で生成されたフィードバック信号から前記第１の周波数成分を抽出し、前記第１の周波数成分に応じて第１の信号及び第２の信号を生成する第１の生成部と、
前記フィードバック部で生成されたフィードバック信号から前記第２の周波数成分を抽出し、前記第２の周波数成分に応じて第３の信号及び第４の信号を生成する第２の生成部と、
を有する
ことを特徴とする請求項１から３のいずれか１項に記載の歪補償装置。
【請求項５】
前記第１の生成部は、前記フィードバック信号から前記第１の周波数成分を抽出する第３の周波数変換部を有し、
前記第２の生成部は、前記フィードバック信号から前記第２の周波数成分を抽出する第４の周波数変換部を有する
ことを特徴とする請求項４に記載の歪補償装置。
【請求項６】
前記第１の生成部は、
前記フィードバック信号の中心周波数を前記第１の周波数成分に対応した周波数にシフトさせる第３の周波数変換部と、
前記第３の周波数変換部で中心周波数がシフトされたフィードバック信号から前記第１の周波数成分を抽出する第１のローパスフィルタと、
をさらに有し、
前記第２の生成部は、
前記フィードバック信号の中心周波数を前記第２の周波数成分に対応した周波数にシフトさせる第４の周波数変換部と、
前記第４の周波数変換部で中心周波数がシフトされたフィードバック信号から前記第２の周波数成分を抽出する第２のローパスフィルタと、
をさらに有する
ことを特徴とする請求項４に記載の歪補償装置。
【請求項７】
前記第１の信号は、第１の周波数帯に対応し、
前記第２の信号は、第２の周波数帯に対応し、
前記送信信号は、前記第１の周波数帯に対応した第５の信号と前記第２の周波数帯に対応した第６の信号とを含み、
前記第１の生成部は、
前記フィードバック信号から前記第１の周波数成分を抽出する第３の周波数変換部と、
前記第３の周波数変換部で抽出された前記第１の周波数成分を含む信号から前記第６の信号を減算し前記第１の信号を生成する第１の減算部と、
前記第３の周波数変換部で抽出された前記第１の周波数成分を含む信号から前記第５の信号を減算し前記第２の信号を生成する第２の減算部と、
をさらに有する
ことを特徴とする請求項４に記載の歪補償装置。
【請求項８】
前記第１の生成部は、１回目のサンプリングにおいて、前記第１の周波数成分を含む信号を出力し、２回目以降のサンプリングにおいて、前記第１の信号及び前記第２の信号を出力し、
前記歪補償装置は、
１回目のサンプリングにおいて、前記第６の信号及び前記第５の信号と前記第１の周波数成分を含む信号とを用いて、奇数次歪補償係数の初期値を求め、２回目以降のサンプリングにおいて、前記第６の信号及び前記第５の信号と前記第１の信号及び前記第２の信号とを用いて、前記奇数次歪補償係数を更新する奇数次更新部と、
前記第３の信号及び前記第４の信号を用いて、偶数次歪補償係数を更新する偶数次更新部と、
１回目のサンプリングにおいて、前記求められた奇数次歪補償係数の初期値を用いて、前記送信信号における奇数次歪を補償し、２回目以降のサンプリングにおいて、前記更新された奇数次歪補償係数を用いて、前記送信信号における奇数次歪を補償する奇数次歪補償部と、
前記偶数次歪補償係数を用いて、前記送信信号における偶数次歪を補償する偶数次歪補償部と、
をさらに備えた
ことを特徴とする請求項７に記載の歪補償装置。
【請求項９】
前記第６の信号及び前記第５の信号と前記第１の周波数成分を含む信号とを用いて、奇数次歪補償係数を更新する奇数次更新部と、
前記第３の信号及び前記第４の信号を用いて、偶数次歪補償係数を更新する偶数次更新部と、
前記奇数次歪補償係数を用いて、前記送信信号における奇数次歪を補償する奇数次歪補償部と、
前記偶数次歪補償係数を用いて、前記送信信号における偶数次歪を補償する偶数次歪補償部と、
をさらに備えた
ことを特徴とする請求項７に記載の歪補償装置。
【請求項１０】
送信信号を増幅するパワーアンプの出力に応じて、複数の奇数次歪が重ねられた第１の周波数成分と複数の偶数次歪が重ねられた第２の周波数成分とを含むフィードバック信号を生成することと、
前記生成されたフィードバック信号から前記第１の周波数成分及び前記第２の周波数成分をそれぞれ分離して信号処理することと、
を含むことを特徴とする歪補償方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、歪補償装置及び歪補償方法に関する。
【背景技術】
【０００２】
従来、無線移動体通信の分野では、基地局等に用いられるパワーアンプの非線形歪を補償する技術が提案されている。かかる技術において、複数の周波数帯の信号を共通のパワーアンプにより同時に増幅する場合、複数の歪補償部やデュアルバンドＤＰＤ（Ｄｉｇｉｔａｌ  Ｐｒｅ−Ｄｉｓｔｏｒｔｉｏｎ）を用いて、非線形歪を補償するものがある。更に、デュアルバンドＤＰＤを用いたパワーアンプでは、二次元ＬＵＴ（Ｌｏｏｋ  Ｕｐ  Ｔａｂｌｅ）の前段に、ｄｕａｌ−ｉｎｐｕｔ  ｔｒｕｎｃａｔｅｄ  Ｖｏｌｔｅｒｒａ  ｍｏｄｅｌのＶｏｌｔｅｒｒａフィルタを設け、二次歪も補償するものもある。他にも、他のバンドに生じる高調波や相互変調歪を逆特性となる信号に付加することにより歪補償を行うものや、非線形連立方程式を解くことにより歪補償を行うものがある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−０８４１４６号公報
【特許文献２】特開２０１２−２２７８８１号公報
【特許文献３】特開２０１４−００３５２７号公報
【非特許文献】
【０００４】
【非特許文献１】Ｂａｓｓａｍ，  Ｓ．Ａ．，  Ｈｅｌａｏｕｉ，  Ｍ．，  Ｇｈａｎｎｏｕｃｈｉ，  Ｆ．Ｍ．，  “２−Ｄ  Ｄｉｇｉｔａｌ  Ｐｒｅｄｉｓｔｏｒｔｉｏｎ  （２−Ｄ−ＤＰＤ）  Ａｒｃｈｉｔｅｃｔｕｒｅ  ｆｏｒ  Ｃｏｎｃｕｒｒｅｎｔ  Ｄｕａｌ−Ｂａｎｄ  Ｔｒａｎｓｍｉｔｔｅｒｓ  ，”  ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｍｉｃｒｏｗａｖｅ  Ｔｈｅｏｒｙ  ａｎｄ  Ｔｅｃｈｎｉｑｕｅｓ，  ｖｏｌ．  ５９，  ｎｏ．  １０，  ｐｐ．  ２５４７−２５５３，  Ｏｃｔ．  ２０１１．
【非特許文献２】Ｙｏｕｎｅｓ，  Ｍ．，  Ｇｈａｎｎｏｕｃｈｉ，  Ｆ．Ｍ．，  “Ｏｎ  ｔｈｅ  Ｍｏｄｅｌｉｎｇ  ａｎｄ  Ｌｉｎｅａｒｉｚａｔｉｏｎ  ｏｆ  ａ  Ｃｏｎｃｕｒｒｅｎｔ  Ｄｕａｌ−Ｂａｎｄ  Ｔｒａｎｓｍｉｔｔｅｒ  Ｅｘｈｉｂｉｔｉｎｇ  Ｎｏｎｌｉｎｅａｒ  Ｄｉｓｔｏｒｔｉｏｎ  ａｎｄ  Ｈａｒｄｗａｒｅ  Ｉｍｐａｉｒｍｅｎｔｓ，”  ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｃｉｒｃｕｉｔｓ  ａｎｄ  Ｓｙｓｔｅｍｓ  Ｉ，  ｖｏｌ．  ６０，  ｎｏ．  １１，  ｐｐ．  ３０５５−３０６８，  Ｎｏｖ．  ２０１３．
【非特許文献３】大田，川崎，木村，玉野井，馬庭透，“適応的マルチバンド共通増幅歪補償法の検討”，信学技報，ｖｏｌ．  １１６，ｎｏ．  １２，ＭＷ２０１６−２，ｐｐ．  ７−１０，Ａｐｒ．  ２０１６．
【非特許文献４】Ｙｕｅｌｉｎ  Ｍａ，  Ｙａｓｕｓｈｉ  Ｙａｍａｏ，  “Ｓｐｅｃｔｒａ−Ｆｏｌｄｉｎｇ  Ｆｅｅｄｂａｃｋ  Ａｒｃｈｉｔｅｃｔｕｒｅ  ｆｏｒ  Ｃｏｎｃｕｒｒｅｎｔ  Ｄｕａｌ−Ｂａｎｄ  Ｐｏｗｅｒ  Ａｍｐｌｉｆｉｅｒ  Ｐｒｅｄｉｓｔｏｒｔｉｏｎ，”  ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｍｉｃｒｏｗａｖｅ  Ｔｈｅｏｒｙ  ａｎｄ  Ｔｅｃｈｎｉｑｕｅｓ，  ｖｏｌ．  ６３，  ｎｏ．  １０，  ｐｐ．  ３１６４−３１７４，  Ｏｃｔ．  ２０１５．
【非特許文献５】Ｙｏｕｊｉａｎｇ  Ｌｉｕ，  Ｊｏｎｍｅｉ  Ｊ．  Ｙａｎ，  Ｐｅｔｅｒ  Ｍ．  Ａｓｂｅｃｋ，  “Ｃｏｎｃｕｒｒｅｎｔ  Ｄｕａｌ−Ｂａｎｄ  Ｄｉｇｉｔａｌ  Ｐｒｅｄｉｓｔｏｒｔｉｏｎ  Ｗｉｔｈ  ａ  Ｓｉｎｇｌｅ  Ｆｅｅｄｂａｃｋ  Ｌｏｏｐ，”  ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｍｉｃｒｏｗａｖｅ  Ｔｈｅｏｒｙ  ａｎｄ  Ｔｅｃｈｎｉｑｕｅｓ，  ｖｏｌ．  ６３，  ｎｏ．  ５，  ｐｐ．  １５５６−１５６８，  Ｍａｙ  ２０１５．
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、歪補償装置では、十分な歪補償性能が得られず、無線品質の低下を招いてしまうことがある。例えば、複数の周波数帯の信号を共通のパワーアンプにより同時に増幅する場合を想定する。この場合、２つのバンドの中心周波数が約２倍程度の関係（例えば、４．６ＧＨｚと９ＧＨｚ）にある条件下では、これらの中心周波数の差に相当する周波数や、低い方の中心周波数の２倍に相当する周波数において、偶数次の歪（以下、単に「偶数次歪」と記す。）が発生することがある。また、複数の周波数帯の信号を複数のパワーアンプで同時増幅して合成する場合であっても、合成器のアイソレーションが不十分なケースでは同様に偶数次歪が発生することがある。このような偶数次歪を補償するように歪補償装置を構成すると、歪補償装置の回路規模が増大することがある。
【０００６】
開示の技術は、上記に鑑みてなされたものであって、回路規模を容易に低減できる歪補償装置及び歪補償方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本願の開示する歪補償装置は、一つの態様において、送信信号を増幅するパワーアンプの出力に応じて、複数の奇数次歪が重ねられた第１の周波数成分と複数の偶数次歪が重ねられた第２の周波数成分とを含むフィードバック信号を生成するフィードバック部と、前記生成されたフィードバック信号から前記第１の周波数成分及び前記第２の周波数成分をそれぞれ分離して信号処理する信号処理部とを有する。
【発明の効果】
【０００８】
本願の開示する歪補償装置及び歪補償方法の一つの態様によれば、歪補償装置の回路規模を容易に低減できる。
【図面の簡単な説明】
【０００９】
【図１】図１は、実施例にかかる歪補償装置の構成を示す図である。
【図２】図２は、実施例における奇数次歪及び偶数次歪を含む周波数スペクトルを示す図である。
【図３】図３は、実施例における共通フィードバック信号処理部の構成を示す図である。
【図４】図４は、実施例にかかる歪補償装置の動作を示すフローチャートである。
【図５】図５は、変形例１における共通フィードバック信号処理部の構成を示す図である。
【図６】図６は、変形例２における共通フィードバック信号処理部の構成を示す図である。
【図７】図７は、変形例２及び変形例３における奇数次歪補償係数更新部の構成を示す図である。
【図８】図８は、変形例４における奇数次歪補償部の構成を示す図である。
【図９】図９は、変形例４における奇数次歪補償係数更新部の構成を示す図である。
【図１０】図１０は、変形例４にかかる奇数次歪補償係数更新部及び奇数次歪補償部の動作を示すフローチャートである。
【図１１】図１１は、変形例５における奇数次歪補償係数更新部の構成を示す図である。
【発明を実施するための形態】
【００１０】
以下に、本願の開示する歪補償装置及び歪補償方法の実施例を図面に基づいて詳細に説明する。なお、この実施例により開示技術が限定されるものではない。また、実施例において同一の機能を有する構成には同一の符号を付し、重複する説明は省略される。
【実施例】
【００１１】
実施例にかかる歪補償装置について説明する。基地局において、送信信号をパワーアンプで増幅する際に、パワーアンプの非線形性が時間的に変動し得る。歪補償装置は、送信信号をパワーアンプで増幅する際におけるパワーアンプの非線形性を適応的に補償できる。
【００１２】
例えば、歪補償装置１は、図１で示すように、パワーアンプ４０及びアンテナ５０を有する。図１に示す歪補償装置１は、パワーアンプ４０の出力に応じた非定常データの歪補償を適応的に行うＤＰＤ（デジタルプレディストーション）技術を採用している。すなわち、歪補償装置１は、時間的に変動し得るパワーアンプ４０の非線形性に対して適応的に歪補償を行う。
【００１３】
パワーアンプ４０の非線形性に応じた歪のうち、パワーアンプ４０の非線形性を数学的にべき級数で展開して表した際の奇数次の項に対応した歪が奇数次歪と呼ばれ、偶数次の項に対応した歪が偶数次歪と呼ばれる。周波数スペクトルで見た場合、奇数次歪は、ＲＦバンド成分の周波数帯域及びＲＦバンド成分の３倍の周波数帯域などに現れ、偶数次歪は、ベースバンド近くの周波数帯域及びＲＦバンド成分の２倍の周波数帯域などに現れる。
【００１４】
送信信号がマルチバンド信号である場合、複数の周波数帯域の信号を同時に共通のパワーアンプ４０で増幅することがある。この際、ある周波数帯域の信号の非線形歪は他の周波数帯域の信号の影響を含む。例えば、図２に破線で示す周波数スペクトルにおいて、一点鎖線で囲った領域に歪が現れる。バンドＡとバンドＢとのデュアルバンド信号をパワーアンプ４０で共通増幅する場合、バンドＡの中心周波数をｆＡとし、バンドＢの中心周波数をｆＢとすると、図２に示すように歪が現れる。バンドＡ，Ｂの周波数ｆＡ，ｆＢの近傍領域のそれぞれに奇数次歪が現れ、バンドＡ，Ｂの差周波数ｆＢ−ｆＡ（＝ｆＣとする）の近傍領域とバンドＡの２倍の周波数２ｆＡ（＝ｆＤとする）の近傍領域とに偶数次歪が現れる。すなわち、低バンド帯ＢＷＬにおいて、周波数ｆＢ−ｆＡ＝ｆＣの近傍領域に偶数次歪が現れ、周波数ｆＡの近傍領域に奇数次歪が現れる。高バンド帯ＢＷＨにおいて、周波数ｆＢの近傍領域に奇数次歪が現れ、周波数２ｆＡ＝ｆＤの近傍領域に偶数次歪が現れる。
【００１５】
図２に示すように、奇数次歪が送信信号の各ＲＦバンド成分（周波数ｆＡ，ｆＢ）に周波数的に重複し得るのに対して、偶数次歪が各ＲＦバンドの帯域外に発生することがある。このため、偶数次歪は、各ＲＦバンドの帯域外の成分を選択的に減衰させるフィルタ処理で除去可能とみなされることがある。
【００１６】
しかしながら、パワーアンプ４０でマルチバンド信号を増幅する場合に、２つのバンドの中心周波数が約２倍程度の関係（２×ｆＡ≒ｆＢ、例えば、ｆＡ＝４．６ＧＨｚ、ｆＢ＝９ＧＨｚ）にあるような条件下では、偶数次歪がＲＦバンド成分に周波数的に近接し得る。バンドＡ，Ｂの差周波数ｆＢ−ｆＡの近傍に現れる偶数次歪が周波数ｆＡのＲＦバンド成分（バンドＡ）に周波数的に近接し得る。バンドＡの２倍の周波数２ｆＡの近傍に現れる偶数次歪が周波数ｆＢのＲＦバンド成分（バンドＢ）に周波数的に近接し得る。このような偶数次歪が送信信号のＲＦバンド成分に近い周波数に発生すると、偶数次歪をフィルタによりカットすることは困難になり得る。また、カットするとしても、急峻なフィルタを用いるので、歪補償装置１の回路規模が大きくなってしまう可能性がある。
【００１７】
それに対して、図１に示す歪補償装置１は、奇数次歪と偶数次歪とをそれぞれ補償できるように構成されている。歪補償装置１は、補償回路３０を有する。補償回路３０は、奇数次歪補償係数更新部３１、偶数次歪補償係数更新部３２、奇数次歪補償部３３、偶数次歪補償部３４、アップコンバータ３５，３６，３７，３８，４８，４９を有する。補償回路３０は、発振器３９，４１，４２，４３，５１，５２、信号合成器４４，４５，５３、ＤＡＣ（Ｄｉｇｉｔａｌ  ｔｏ  Ａｎａｌｏｇ  Ｃｏｎｖｅｒｔｅｒ）４６，４７を有する。
【００１８】
奇数次歪補償係数更新部３１は、奇数次歪補償係数を更新し、更新後の奇数次歪補償係数（＾Ａ，＾Ｂ）を奇数次歪補償部３３へ供給する。奇数次歪補償部３３は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）を受ける。奇数次歪補償部３３は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、奇数次歪補償係数（＾Ａ）を用いて、バンドＡ（中心周波数ｆＡ）の奇数次歪補償信号＾ｚＬ（ｎ）を生成する。奇数次歪補償部３３は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、奇数次歪補償係数（＾Ｂ）を用いて、バンドＢ（中心周波数ｆＢ）の奇数次歪補償信号＾ｚＨ（ｎ）を生成する。
【００１９】
同様に、偶数次歪補償係数更新部３２は、偶数次歪補償係数を更新し、更新後の偶数次歪補償係数（＾Ｃ，＾Ｄ）を偶数次歪補償部３４へ供給する。偶数次歪補償部３４は、バンドＡの送信信号ｘＬ（ｎ）を受け、バンドＢの送信信号ｘＨ（ｎ）を受ける。偶数次歪補償部３４は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、偶数次歪補償係数（＾Ｃ）を用いて、中心周波数ｆＣ＝ｆＢ−ｆＡの偶数次歪補償信号＾ｚＨ−Ｌ（ｎ）を生成する。偶数次歪補償部３４は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、偶数次歪補償係数（＾Ｄ）を用いて、中心周波数ｆＤ＝２ｆＡの奇数次歪補償信号＾ｚＬ，２（ｎ）を生成する。
【００２０】
デジタル領域において、アップコンバータ３５は、バンドＡ（中心周波数ｆＡ）の奇数次歪補償信号＾ｚＬ（ｎ）を発振器３９からの周波数ｆ’Ａ＝ｆＡ−ｆＢ／２のローカル信号に乗算させて周波数変換する。アップコンバータ３８は、中心周波数ｆＣ＝ｆＢ−ｆＡの偶数次歪補償信号＾ｚＨ−Ｌ（ｎ）を発振器３９からの周波数ｆ’Ｄ＝ｆＢ／２−ｆＡのローカル信号に乗算させて周波数変換する。周波数変換後の奇数次歪補償信号＾ｚＬ（ｎ）及び偶数次歪補償信号＾ｚＨ−Ｌ（ｎ）は、デジタルの信号合成器４４で合成された後、ＤＡＣ４６においてＤＡ変換される。アップコンバータ４８は、ＤＡＣ４６でＤＡ変換されたアナログ信号を発振器５１からの周波数ｆＬ１＝ｆＢ／２のローカル信号に乗算させて無線周波数（周波数ｆＡ）に周波数変換してＲＦ信号とする。
【００２１】
同様に、アップコンバータ３６は、バンドＢ（中心周波数ｆＢ）の奇数次歪補償信号＾ｚＨ（ｎ）を発振器４１からの周波数ｆ’Ｂ＝ｆＢ／２−ｆＡのローカル信号に乗算させて周波数変換する。アップコンバータ３７は、中心周波数ｆＤ＝２ｆＡの偶数次歪補償信号＾ｚＬ，２（ｎ）を発振器４２からの周波数ｆ’Ｃ＝ｆＡ−ｆＢ／２のローカル信号に乗算させて周波数変換する。周波数変換後の奇数次歪補償信号＾ｚＨ（ｎ）及び偶数次歪補償信号＾ｚＬ，２（ｎ）は、デジタルの信号合成器４５で合成された後、ＤＡＣ４７においてＤＡ変換される。アップコンバータ４９は、ＤＡＣ４７でＤＡ変換されたアナログ信号を発振器５２からの周波数ｆＬ２＝ｆＡ＋ｆＢ／２のローカル信号に乗算させて無線周波数（周波数ｆＢ）に周波数変換してＲＦ信号とする。
【００２２】
アップコンバータ４８からのＲＦ信号とアップコンバータ４９からのＲＦ信号とは、信号合成器５３で合成された後、パワーアンプ４０に入力される。パワーアンプ４０に入力されるＲＦ信号には奇数次歪の逆特性に応じた成分と偶数次歪の逆特性に応じた成分とが含まれているので、ＲＦ信号（送信信号）がパワーアンプ４０で電力増幅される際における奇数次歪と偶数次歪とをそれぞれ抑制できる。すなわち、歪補償装置１は、パワーアンプ４０の非線形性に応じた奇数次歪と偶数次歪とを補償できる。
【００２３】
歪補償装置１は、カップラ６０をさらに有する。カップラ６０は、パワーアンプ４０の出力側に具備されており、パワーアンプ４０の出力の一部（パワーアンプ４０の出力に応じた信号）を取り出すことができる。
【００２４】
時間的に変動し得るパワーアンプ４０の非線形性を適応的に補償するために、歪補償装置１において、信号分配器、４つのダウンコンバータ、４つのＡＤＣ（Ａｎａｌｏｇ  ｔｏ  Ｄｉｇｉｔａｌ  Ｃｏｎｖｅｒｔｅｒ）を追加することが考えられる。すなわち、パワーアンプ４０の出力に応じてカップラ６０で取り出された信号を信号分配器で４つの信号に分配し、４つの信号を４つのダウンコンバータで周波数ｆＡ、ｆＢ、ｆＣ＝ｆＢ−ｆＡ、ｆＤ＝２ｆＡからベースバンドへそれぞれ周波数変換する。４つのＡＤＣは、４つのダウンコンバータのうち対応するダウンコンバータから入力されるアナログ信号の各々をＡＤ変換して補償回路３０へ供給することになる。この場合、補償すべき歪の個数分のダウンコンバータやＡＤＣが設けられるため、アナログ部品点数が多くなりやすく、歪補償装置１の回路規模が増大しやすい。
【００２５】
そこで、実施例では、歪補償装置１において、奇数次歪どうしを重ねるとともに偶数次歪どうしを重ねたフィードバック信号を生成して周波数変換及びＡＤ変換を共通で行うようにフィードバック系を構成することで、歪補償装置１の回路規模の低減を図る。
【００２６】
具体的には、歪補償装置１は、図１に示すように、共通フィードバック部１０及び共通フィードバック信号処理部７０をさらに有する。共通フィードバック部１０は、カップラ６０及び共通フィードバック信号処理部７０の間に配されている。共通フィードバック部１０は、送信信号を増幅するパワーアンプの出力に応じて共通フィードバック信号を生成する。共通フィードバック信号は、周波数スペクトル上で複数の奇数次歪が重ねられた奇数次歪成分（第１の周波数成分）と周波数スペクトル上で複数の偶数次歪が重ねられた偶数次歪成分（第２の周波数成分）とを含む。共通フィードバック部１０は、パワーアンプ４０の出力に応じてカップラ６０で取り出された信号をダウンコンバートする際に、バンドＡ（中心周波数：ｆＡ）の近傍に発生する奇数次歪とバンドＢ（中心周波数：ｆＢ）の近傍に発生する奇数次歪とを重ねる。それとともに、共通フィードバック部１０は、中心周波数ｆＤ＝２ｆＡの近傍に発生する偶数次歪と中心周波数ｆＣ＝ｆＢ−ｆＡの近傍に発生する偶数次歪とを重ねる。これにより、共通フィードバック部１０は、共通フィードバック信号を生成する。共通フィードバック部１０は、生成された共通フィードバック信号を共通フィードバック信号処理部７０へ供給する。
【００２７】
共通フィードバック部１０は、共通ダウンコンバータ（第１の周波数変換部）１１、共通ダウンコンバータ（第２の周波数変換部）１２、ＡＤＣ（ＡＤ変換部）１３、発振器１４，１５を有する。
【００２８】
共通ダウンコンバータ１１は、パワーアンプ４０の出力の一部がカップラ６０で取り出された信号に対して、発振器１４からのローカル周波数ｆＣＯ１＝（ｆＡ＋ｆＢ）／２のローカル信号を乗算する。これにより、共通ダウンコンバータ１１は、カップラ６０で取り出された信号を、ＲＦ周波数（送信周波数）からＩＦ周波数（中間周波数）にダウンコンバートする。このとき、バンドＡの中心周波数ｆＡ近傍の奇数次歪は
ｆＣＯ１−ｆＡ＝（ｆＢ−ｆＡ）／２・・・数式１
にダウンコンバートされる。バンドＢの中心周波数ｆＢ近傍の奇数次歪は
ｆＢ−ｆＣＯ１＝（ｆＢ−ｆＡ）／２・・・数式２
にダウンコンバートされる。数式１、数式２により、奇数次歪どうしは、略同一の中間周波数（ｆＢ−ｆＡ）／２にダウンコンバートされる。
【００２９】
また、中心周波数ｆＣ＝ｆＢ−ｆＡ近傍の偶数次歪が
ｆＣＯ１−ｆＣ＝（３ｆＡ−ｆＢ）／２・・・数式３
にダウンコンバートされる。中心周波数ｆＤ＝２ｆＡ近傍の偶数次歪が
ｆＤ−ｆＣＯ１＝（３ｆＡ−ｆＢ）／２・・・数式４
にダウンコンバートされる。数式３、数式４により、偶数次歪どうしは、略同一の中間周波数（３ｆＡ−ｆＢ）／２にダウンコンバートされる。
【００３０】
すなわち、共通ダウンコンバータ１１は、カップラ６０で取り出された信号をＲＦ周波数からＩＦ周波数にダウンコンバートしながらその信号をＩＦ周波数で折り返して共通フィードバック信号を生成する。共通フィードバック信号は、周波数スペクトル上で複数の奇数次歪が重ねられた奇数次歪成分（第１の周波数成分）と周波数スペクトル上で複数の偶数次歪が重ねられた偶数次歪成分（第２の周波数成分）とを含む。
【００３１】
共通ダウンコンバータ１２は、共通ダウンコンバータ１１で生成された共通フィードバック信号に対して、発振器１５からのローカル周波数ｆＣＯ２＝ｆＡ／２のローカル信号を乗算し、中間（ＩＦ）周波数からベースバンド周波数にダウンコンバートする。このとき、奇数次歪どうしが重なった奇数次歪成分が
（ｆＢ−ｆＡ）／２−ｆＣＯ２＝ｆＢ／２−ｆＡ・・・数式５
にダウンコンバートされる。偶数次歪どうしが重なった偶数次歪成分が
（３ｆＡ−ｆＢ）／２−ｆＣＯ２＝ｆＡ−ｆＢ／２・・・数式６
にダウンコンバートされる。数式５、数式６により、ダウンコンバート後の共通フィードバック信号では、次の数式７に示すように、奇数次歪成分及び偶数次歪成分のちょうど中点の周波数が略０Ｈｚとなる。
｛（ｆＢ／２−ｆＡ）＋（ｆＡ−ｆＢ／２）｝／２＝０・・・数式７
すなわち、共通ダウンコンバータ１２は、奇数次歪成分及び偶数次歪成分のちょうど中点の周波数が略０Ｈｚとなるように、共通フィードバック信号の中心周波数をシフトさせる。
【００３２】
ＡＤＣ１３は、共通ダウンコンバータ１２でベースバンド周波数にダウンコンバートされた共通フィードバック信号をアナログ形式からデジタル形式に変換し、変換後の共通フィードバック信号を共通フィードバック信号処理部７０へ供給する。
【００３３】
共通フィードバック信号処理部７０は、共通フィードバック信号から奇数次歪成分（第１の周波数成分）及び偶数次歪成分（第２の周波数成分）をそれぞれ分離して信号処理する。
【００３４】
例えば、共通フィードバック信号処理部７０は、図３に示すように、ゲイン調整部７１、生成部（第１の生成部）９１及び生成部（第２の生成部）９２を有する。生成部９１は、ミキサ（第３の周波数変換部）７２、ＬＰＦ（第１のローパスフィルタ）７４、遅延調整部７６、奇数次歪アンプモデル係数算出部７８、奇数次歪フィードバック信号レプリカ生成部８１、及び発振器８３を有する。生成部９２は、ミキサ（第４の周波数変換部）７３、ＬＰＦ（第２のローパスフィルタ）７５、遅延調整部７７、偶数次歪アンプモデル係数算出部７９、偶数次歪フィードバック信号レプリカ生成部８２、及び発振器８４を有する。
【００３５】
ゲイン調整部７１は、共通フィードバック信号を共通フィードバック部１０のＡＤＣ１３から受ける。ゲイン調整部７１は、共通フィードバック信号に応じて、ゲインを調整する。例えば、ゲイン調整部７１は、共通フィードバック信号の電力の平均値と所定の目標値とを比較して、ゲインを求める。ゲイン調整部７１は、求められたゲインに応じてレベル調整された共通フィードバック信号をミキサ７２、ミキサ７３へそれぞれ転送する。
【００３６】
ミキサ７２は、共通フィードバック信号をゲイン調整部７１から受ける。ミキサ７２は、共通フィードバック信号に対して、発振器８３からのローカル周波数ｆｏｄｄ＝ｆＡ−ｆＢ／２のローカル信号を乗算し、共通フィードバック信号の中心周波数をｆＡ−ｆＢ／２だけシフトさせる。このとき、数式５により、共通フィードバック信号における奇数次歪成分の周波数が
（ｆＢ／２−ｆＡ）＋（ｆＡ−ｆＢ／２）＝０・・・数式８
に変換される。数式６により、偶数次歪成分の周波数が
（ｆＡ−ｆＢ／２）＋（ｆＡ−ｆＢ／２）＝２ｆＡ−ｆＢ・・・数式９
近傍に変換される。
【００３７】
ＬＰＦ７４は、ミキサ７２で中心周波数がシフトされた共通フィードバック信号を受ける。ＬＰＦ７４は、共通フィードバック信号における奇数次歪成分を選択的に通過させ偶数次歪成分を除去する。すなわち、ＬＰＦ７４は、共通フィードバック信号から奇数次歪成分を抽出して、奇数次歪成分を含む信号を遅延調整部７６へ供給する。
【００３８】
遅延調整部７６は、奇数次歪成分を含む信号を所定の遅延量で遅延させて、遅延後の信号を奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）として奇数次歪アンプモデル係数算出部７８へ供給する。
【００３９】
同様に、ミキサ７３は、共通フィードバック信号をゲイン調整部７１から受ける。ミキサ７３は、共通フィードバック信号に対して、発振器８４からのローカル周波数ｆｅｖｅｎ＝ｆＢ／２−ｆＡのローカル信号を乗算し、共通フィードバック信号の中心周波数をｆＢ／２−ｆＡだけシフトさせる。このとき、数式５により、共通フィードバック信号における奇数次歪成分の周波数が
（ｆＢ／２−ｆＡ）＋（ｆＢ／２−ｆＡ）＝ｆＢ−２ｆＡ・・・数式１０
に変換される。数式６により、偶数次歪成分の周波数が
（ｆＡ−ｆＢ／２）＋（ｆＢ／２−ｆＡ）＝０・・・数式１１
近傍に変換される。
【００４０】
ＬＰＦ７５は、ミキサ７３で中心周波数がシフトされた共通フィードバック信号を受ける。ＬＰＦ７５は、共通フィードバック信号における偶数次歪成分を選択的に通過させ奇数次歪成分を除去する。すなわち、ＬＰＦ７５は、共通フィードバック信号から偶数次歪成分を抽出して、偶数次歪成分を含む信号を遅延調整部７７へ供給する。
【００４１】
遅延調整部７７は、偶数次歪成分を含む信号を所定の遅延量で遅延させて、遅延後の信号を偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）として偶数次歪アンプモデル係数算出部７９へ供給する。
【００４２】
ここで、簡単のために共通フィードバック信号処理部７０がメモリレスの構成であると仮定すると、奇数次歪アンプモデル係数算出部７８へ供給される奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）は次の数式１２で表される。
【数１】
2018101861.tif 000003
【００４３】
数式１２において、ｘＬ（ｎ）は低バンド帯ＢＷＬ（図２参照）側の送信信号であり、ｘＨ（ｎ）は高バンド帯ＢＷＨ（図２参照）側の送信信号である。また、ｙＬ（ｎ）はパワーアンプ４０から出力されるバンドＡ（中心周波数：ｆＡ）の信号を表し、ｙＨ（ｎ）はパワーアンプ４０から出力されるバンドＢ（中心周波数：ｆＢ）の信号を表す。ａｐ，ｑ、ｂｐ，ｑは、パワーアンプ４０の特性に対応したアンプモデルから導かれるアンプモデル係数である。アンプモデルは、奇数次歪アンプモデル係数算出部７８に予め設定されている。このとき、ｙｏｄｄ（ｎ）を行列形式で表現すると、次の数式１３で表される。
【数２】
2018101861.tif 000004
【００４４】
数式１３中のｙｏｄｄは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における奇数次歪共通フィードバック信号を含むベクトルであり、次の数式１４で表現される。
【数３】
2018101861.tif 000005
【００４５】
数式１３中のＵｏｄｄは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における各バンドの送信信号の奇数次項を含む行列であり、次の数式１５で表現される。
【数４】
2018101861.tif 000006
【００４６】
数式１３中のｈｏｄｄは、パワーアンプ４０の伝達関数のベクトルであり、次の数式１６で表現される。
【数５】
2018101861.tif 000007
【００４７】
数式１３の行列方程式をｈｏｄｄについて解くと、数式１７に示すようになる。数式１７では、ｈｏｄｄを特にベクトル＾ｈｏｄｄと表している。
【数６】
2018101861.tif 000008
【００４８】
数式１７において、添え字Ｈはエルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は逆行列を表す。
【００４９】
すなわち、図３に示す奇数次歪アンプモデル係数算出部７８は、奇数次歪共通フィードバック信号を含むベクトルｙｏｄｄと各バンドの送信信号の奇数次項を含む行列Ｕｏｄｄとを用いて、数式１７に示すように係数ベクトル＾ｈｏｄｄを求める。係数ベクトル＾ｈｏｄｄは、数式１６に示すように、アンプモデル係数＾ａｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）とアンプモデル係数＾ｂｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）とを含む。奇数次歪アンプモデル係数算出部７８は、求められた係数ベクトル＾ｈｏｄｄを奇数次歪フィードバック信号レプリカ生成部８１へ供給する。
【００５０】
奇数次歪フィードバック信号レプリカ生成部８１は、アンプモデル係数＾ａｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）を用いて、次の数式１８に示すようにバンドＡ（中心周波数：ｆＡ）の奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）を生成する。アンプモデル係数＾ａｐ，ｑは、係数ベクトル＾ｈｏｄｄに含まれている。
【数７】
2018101861.tif 000009
【００５１】
奇数次歪フィードバック信号レプリカ生成部８１は、アンプモデル係数＾ｂｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）を用いて、次の数式１９に示すようにバンドＢ（中心周波数：ｆＢ）の奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）を生成する。アンプモデル係数＾ｂｐ，ｑは、係数ベクトル＾ｈｏｄｄに含まれている。
【数８】
2018101861.tif 000010
【００５２】
奇数次歪フィードバック信号レプリカ生成部８１は、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）と奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）とをそれぞれ奇数次歪補償係数更新部３１へ供給する。
【００５３】
同様に、偶数次歪アンプモデル係数算出部７９へ供給される偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）は次の数式２０で表される。
【数９】
2018101861.tif 000011
【００５４】
数式２０において、ｘＬ（ｎ）は低バンド帯ＢＷＬ（図２参照）側の送信信号であり、ｘＨ（ｎ）は高バンド帯ＢＷＨ（図２参照）側の送信信号である。また、ｙＨ−Ｌ（ｎ）は中心周波数ｆＣ＝ｆＢ−ｆＡに発生する偶数次歪成分を表し、ｙＬ，２（ｎ）は中心周波数ｆＤ＝２ｆＡに発生する偶数次歪成分を表す。ｃｐ，ｑ、ｄｐ，ｑは、パワーアンプ４０の特性に対応したアンプモデルから導かれるアンプモデル係数である。アンプモデルは、偶数次歪アンプモデル係数算出部７９に予め設定されている。このとき、ｙｅｖｅｎ（ｎ）を行列形式で表現すると、次の数式２１で表される。
【数１０】
2018101861.tif 000012
【００５５】
数式２１中のｙｅｖｅｎは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における偶数次歪共通フィードバック信号を含むベクトルであり、次の数式２２で表現される。
【数１１】
2018101861.tif 000013
【００５６】
数式２１中のＵｅｖｅｎは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における各バンドの送信信号の偶数次項を含む行列であり、次の数式２３で表現される。
【数１２】
2018101861.tif 000014
【００５７】
数式２１中のｈｅｖｅｎは、パワーアンプ４０の伝達関数のベクトルであり、次の数式２４で表現される。
【数１３】
2018101861.tif 000015
【００５８】
数式２１の行列方程式をｈｅｖｅｎについて解くと、数式２５に示すようになる。数式２５では、ｈｅｖｅｎを特にベクトル＾ｈｅｖｅｎと表している。
【数１４】
2018101861.tif 000016
【００５９】
数式２５において、添え字「Ｈ」はエルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字「−１」は逆行列を表す。
【００６０】
すなわち、図３に示す偶数次歪アンプモデル係数算出部７９は、偶数次歪共通フィードバック信号を含むベクトルｙｅｖｅｎと各バンドの送信信号の偶数次項を含む行列Ｕｅｖｅｎとを用いて、数式２５に示すように係数ベクトル＾ｈｅｖｅｎを求める。係数ベクトル＾ｈｅｖｅｎは、数式２４に示されるように、アンプモデル係数＾ｃｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）とアンプモデル係数＾ｄｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）とを含む。偶数次歪アンプモデル係数算出部７９は、求められた係数ベクトル＾ｈｅｖｅｎを偶数次歪フィードバック信号レプリカ生成部８２へ供給する。
【００６１】
偶数次歪フィードバック信号レプリカ生成部８２は、アンプモデル係数＾ｃｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）を用いて、次の数式２６に示すように中心周波数ｆＣ＝ｆＢ−ｆＡの偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ）を生成する。アンプモデル係数＾ｃｐ，ｑは、係数ベクトル＾ｈｅｖｅｎに含まれている。
【数１５】
2018101861.tif 000017
【００６２】
偶数次歪フィードバック信号レプリカ生成部８２は、アンプモデル係数＾ｄｐ，ｑ（０≦ｐ≦Ｐ、０≦ｑ≦Ｑ）を用いて、次の数式２７に示すように中心周波数ｆＤ＝２ｆＡの偶数次歪フィードバック信号レプリカ＾ｙＬ，２（ｎ）を生成する。アンプモデル係数＾ｄｐ，ｑは、係数ベクトル＾ｈｅｖｅｎに含まれている。
【数１６】
2018101861.tif 000018
【００６３】
偶数次歪フィードバック信号レプリカ生成部８２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ）と偶数次歪フィードバック信号レプリカ＾ｙＬ，２（ｎ）とをそれぞれ偶数次歪補償係数更新部３２へ供給する。
【００６４】
奇数次歪補償係数更新部３１は、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）と奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）とを奇数次歪フィードバック信号レプリカ生成部８１から受ける。奇数次歪補償係数更新部３１は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）を受ける。奇数次歪補償係数更新部３１は、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）と奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）とバンドＡの送信信号ｘＬ（ｎ）とバンドＢの送信信号ｘＨ（ｎ）とを用いて奇数次歪補償信号処理を行う。すなわち、奇数次歪補償係数更新部３１は、次の数式２８に示す計算を行い、歪補償係数ベクトル＾Ａ、＾Ｂを更新する。
【数１７】
2018101861.tif 000019
【００６５】
数式２８中の歪補償係数ベクトル＾Ａ、＾Ｂは、更新後の奇数次歪補償係数を示し、それぞれ数式２９のように表現される。
【数１８】
2018101861.tif 000020
【００６６】
数式２８中の行列ＷＬ、ＷＨは、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）と奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）とにより求められる。行列ＷＬ、ＷＨは、各バンドの奇数次歪フィードバック信号の奇数次項を含む行列であり、それぞれ次の数式３０のように表現される。
【数１９】
2018101861.tif 000021
【００６７】
なお、数式２８において、添え字Ｈはエルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は逆行列を表す。
【００６８】
数式２８中のベクトルｚＬ、ｚＨは、バンドＡの送信信号ｘＬ（ｎ）とバンドＢの送信信号ｘＨ（ｎ）とにより求められ、係数更新前の奇数次歪補償処理後の信号を含むベクトルであり、それぞれ次の数式３１で表現される。
【数２０】
2018101861.tif 000022
【００６９】
図１に示す奇数次歪補償係数更新部３１は、更新後の奇数次歪補償係数（歪補償係数ベクトル＾Ａ，＾Ｂ）を奇数次歪補償部３３へ供給する。奇数次歪補償部３３は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）を受ける。奇数次歪補償部３３は、数式３２に示すように、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、奇数次歪補償係数（＾Ａ）を用いてバンドＡ（中心周波数ｆＡ）の奇数次歪補償信号＾ｚＬ（ｎ）を生成する。奇数次歪補償部３３は、次の数式３２に示すように、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、奇数次歪補償係数（＾Ｂ）を用いてバンドＢ（中心周波数ｆＢ）の奇数次歪補償信号＾ｚＨ（ｎ）を生成する。
【数２１】
2018101861.tif 000023
【００７０】
同様に、偶数次歪補償係数更新部３２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ）と偶数次歪フィードバック信号レプリカ＾ｙＬ，２（ｎ）とを偶数次歪フィードバック信号レプリカ生成部８２から受ける。偶数次歪補償係数更新部３２は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）を受ける。偶数次歪補償係数更新部３２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ）と偶数次歪フィードバック信号レプリカ＾ｙＬ，２（ｎ）とバンドＡの送信信号ｘＬ（ｎ）とバンドＢの送信信号ｘＨ（ｎ）とを用いて偶数次歪補償信号処理を行う。すなわち、偶数次歪補償係数更新部３２は、次の数式３３に示す計算を行い、歪補償係数ベクトル＾Ｃ、＾Ｄを更新する。
【数２２】
2018101861.tif 000024
【００７１】
数式３３中の歪補償係数ベクトル＾Ｃ、＾Ｄは、更新後の偶数次歪補償係数を示し、それぞれ次の数式３４のように表現される。
【数２３】
2018101861.tif 000025
【００７２】
数式３３中の行列ＵＨ−Ｌ、ＵＬ，２は、バンドＡの送信信号ｘＬ（ｎ）とバンドＢの送信信号ｘＨ（ｎ）とにより求められ、各バンドの送信信号の偶数次項を含む行列であり、それぞれ次の数式３５のように表現される。
【数２４】
2018101861.tif 000026
【００７３】
なお、数式３３において、添え字Ｈはエルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は逆行列を表す。
【００７４】
数式３３中のベクトル＾ｙＨ−Ｌ、＾ｙＬ，２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ）と偶数次歪フィードバック信号レプリカ＾ｙＬ，２（ｎ）とにより求められる。ベクトル＾ｙＨ−Ｌ、＾ｙＬ，２は、分離後の偶数次歪フィードバック信号を含むベクトルであり、それぞれ次の数式３６で表現される。
【数２５】
2018101861.tif 000027
【００７５】
図１に示す偶数次歪補償係数更新部３２は、更新後の偶数次歪補償係数（歪補償係数ベクトル＾Ｃ，＾Ｄ）を偶数次歪補償部３４へ供給する。偶数次歪補償部３４は、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）を受ける。偶数次歪補償部３４は、次の数式３７に示すように、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、偶数次歪補償係数（＾Ｃ）を用いて中心周波数ｆＣ＝ｆＢ−ｆＡの偶数次歪補償信号＾ｚＨ−Ｌ（ｎ）を生成する。偶数次歪補償部３４は、次の数式３７に示すように、バンドＡの送信信号ｘＬ（ｎ）、バンドＢの送信信号ｘＨ（ｎ）、偶数次歪補償係数（＾Ｄ）を用いて中心周波数ｆＤ＝２ｆＡの偶数次歪補償信号＾ｚＬ，２（ｎ）を生成する。
【数２６】
2018101861.tif 000028
【００７６】
なお、数式１３〜数式１６と数式２２〜数式２６とを比較すると、次の数式３８に示す関係が成り立つ。
【数２７】
2018101861.tif 000029
【００７７】
偶数次歪補償部３４は、数式３８の関係で求めた係数＾Ｃｐ，ｑ、＾Ｄｐ，ｑを用いて、偶数次歪補償処理後の信号＾ｚＨ−Ｌ（ｎ）、＾ｚＬ，２（ｎ）を数式３７により生成してもよい。
【００７８】
次に、歪補償装置１の動作について図４を用いて説明する。図４は、歪補償装置１の動作を示すフローチャートである。
【００７９】
歪補償装置１において、共通フィードバック部１０は、パワーアンプ４０の出力に応じた信号を受けると、その信号から、複数の奇数次歪が重ねられた奇数次歪成分と複数の偶数次歪が重ねられた偶数次歪成分とを含む共通フィードバック信号を生成する（Ｓ１）。共通フィードバック信号処理部７０が共通フィードバック信号を共通フィードバック部１０から受けると、奇数次歪についての処理（Ｓ１１〜Ｓ１４）と偶数次歪についての処理（Ｓ２１〜Ｓ２４）とが並行して行われる。
【００８０】
奇数次歪についての処理では、共通フィードバック信号処理部７０の生成部９１が、共通フィードバック信号から奇数次歪成分を抽出する（Ｓ１１）。生成部９１は、抽出後の信号と各バンドの送信信号とに応じて、数式１７に示す計算により、アンプモデル係数＾ａｐ，ｑ，＾ｂｐ，ｑを含む係数ベクトル＾ｈｏｄｄを求める（Ｓ１２）。生成部９１は、アンプモデル係数＾ａｐ，ｑ，＾ｂｐ，ｑを用いて、数式１８、数式１９に示す計算により、各バンドＡ，Ｂに対応した奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）を生成する（Ｓ１３）。生成部９１は、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）を奇数次歪補償係数更新部３１へ供給する。奇数次歪補償係数更新部３１は、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）を用いて、奇数次歪補償係数を更新し、更新後の奇数次歪補償係数（＾Ａ，＾Ｂ）を奇数次歪補償部３３へ供給する。奇数次歪補償部３３は、更新後の奇数次歪補償係数（＾Ａ，＾Ｂ）を用いて、奇数次歪を補償するＤＰＤ処理を行う（Ｓ１４）。
【００８１】
偶数次歪についての処理では、共通フィードバック信号処理部７０の生成部９２が、共通フィードバック信号から偶数次歪成分を抽出する（Ｓ２１）。生成部９２は、抽出後の信号と各バンドの送信信号とに応じて、数式２５に示す計算により、アンプモデル係数＾ｃｐ，ｑ，＾ｄｐ，ｑを含む係数ベクトル＾ｈｅｖｅｎを求める（Ｓ２２）。生成部９２は、アンプモデル係数＾ｃｐ，ｑ，＾ｄｐ，ｑを用いて、数式２６、数式２７に示す計算により、各バンドＡ，Ｂに対応した偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ），＾ｙＬ，２（ｎ）を生成する（Ｓ２３）。生成部９２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ），＾ｙＬ，２（ｎ）を偶数次歪補償係数更新部３２へ供給する。偶数次歪補償係数更新部３２は、偶数次歪フィードバック信号レプリカ＾ｙＨ−Ｌ（ｎ），＾ｙＬ，２（ｎ）を用いて、偶数次歪補償係数を更新し、更新後の偶数次歪補償係数（＾Ｃ，＾Ｄ）を偶数次歪補償部３４へ供給する。偶数次歪補償部３４は、更新後の偶数次歪補償係数（＾Ｃ，＾Ｄ）を用いて、偶数次歪を補償するＤＰＤ処理を行う（Ｓ２４）。
【００８２】
歪補償装置１は、終了の指示を受けるまで（Ｓ３０でＮｏ）、Ｓ１〜Ｓ２４までの処理を繰り返し、終了の指示を受けると（Ｓ３０でＹｅｓ）、処理を終了する。
【００８３】
以上のように、実施例では、パワーアンプ４０でマルチバンド信号を増幅する際に発生する非線形歪を補償する歪補償装置１において、奇数次歪（主信号帯域）および偶数次歪をそれぞれ重ねて共通フィードバックする共通フィードバック部１０を設ける。すなわち、歪補償装置１において、奇数次歪どうしを重ねるとともに偶数次歪どうしを重ねた共通フィードバック信号を生成して周波数変換及びＡＤ変換を共通で行うようにフィードバック系を構成する。これにより、ダウンコンバータ（周波数変換器）やＡＤＣなどのアナログ部品点数を少なくすることができ、歪補償装置１におけるフィードバック系の回路規模を容易に低減できる。
【００８４】
また、実施例では、歪補償装置１において、共通フィードバック信号処理部７０が、共通フィードバック部１０で生成された共通フィードバック信号から奇数次歪成分と偶数次歪成分とをそれぞれ分離して信号処理する。これにより、歪補償装置１におけるフィードバック系の回路規模を低減できるとともに、奇数次歪と偶数次歪とをそれぞれ補償できる。
【００８５】
（変形例１）
なお、共通フィードバック信号処理部７０ａの生成部９１ａ，９２ａは、図５に示すように、ＬＰＦ７４，７５（図３参照）が省略された構成であってもよい。図５は、変形例１における共通フィードバック信号処理部７０ａの構成を示す図である。例えば、奇数次歪と偶数次歪とが互いに相関性が低いと想定できると、奇数次歪及び偶数次歪において一方が他方の雑音に見える。このため、奇数次歪アンプモデル係数算出部７８、偶数次歪アンプモデル係数算出部７９における係数の算出に用いるフィードバック信号のサンプル数を多く確保すれば、ＬＰＦによるフィルタリング処理を省略しても他方の歪をキャンセルできる。
【００８６】
例えば、生成部９１ａにおいて、奇数次歪アンプモデル係数算出部７８が係数ベクトル＾ｈｏｄｄの計算に用いる奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）のサンプル数を多く（数式１４，１５におけるＮを大きく）確保する。この場合、０Ｈｚ付近にある奇数次歪成分に比べて、ばらついた周波数位置にある偶数次歪成分の振幅が時間平均した場合に小さくなり得る。これにより、生成部９１ａでは、実質的に、ミキサ７２が、共通フィードバック信号から奇数次歪成分を抽出していると見なすことができる。
【００８７】
例えば、生成部９２ａにおいて、偶数次歪アンプモデル係数算出部７９が係数ベクトル＾ｈｅｖｅｎの計算に用いる偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）のサンプル数を多く（数式２２，２３におけるＮを大きく）確保する。この場合、０Ｈｚ付近にある偶数次歪成分に比べて、ばらついた周波数位置にある偶数次歪成分の振幅が時間平均した場合に小さくなり得る。これにより、生成部９２ａでは、実質的に、ミキサ７３が、共通フィードバック信号から偶数次歪成分を抽出していると見なすことができる。
【００８８】
このように、奇数次歪アンプモデル係数算出部７８、偶数次歪アンプモデル係数算出部７９における係数の算出に用いる信号のサンプル数（Ｎ）を多く確保することで実質的にミキサ７２，７３により奇数次歪成分、偶数次歪成分がそれぞれ分離される。これにより、共通フィードバック信号において奇数次歪（主信号帯域）と偶数次歪とが実質的に重なる場合に、共通フィードバック信号から奇数次歪成分、偶数次歪成分をそれぞれ分離することができる。
【００８９】
また、ＬＰＦ７４，７５（図３参照）が省略された構成であるので、共通フィードバック信号処理部７０ａの構成を簡略化でき低コスト化できる。
【００９０】
（変形例２）
あるいは、共通フィードバック信号処理部７０ｂの生成部９１ｂ，９２ｂは、図６に示すように、フィードバック信号レプリカを用いない構成であってもよい。図６は、変形例２における共通フィードバック信号処理部７０ｂの構成を示す図である。例えば、変形例１の「奇数次歪と偶数次歪が互いに相関性が低い」という想定に加えて、「奇数次歪どうしが互いに相関性が低い」かつ「偶数次歪どうしが互いに相関性が低い」という想定が可能な条件であれば、着目する成分以外は雑音とみなせる。このため、サンプル数を多く確保することでキャンセルできる。
【００９１】
具体的には、生成部９１ｂは、生成部９１ａ（図５参照）に対して、遅延調整部７６、奇数次歪アンプモデル係数算出部７８、及び奇数次歪フィードバック信号レプリカ生成部８１が省略されている。生成部９１ｂは、生成部９１ａ（図５参照）に対して、減算器（第１の減算部）８５ｂ、減算器（第２の減算部）８６ｂ、遅延調整部７６１、及び遅延調整部７６２が追加されている。サンプル数を多く確保することで、実質的に、ミキサ７２が、共通フィードバック信号から奇数次歪成分を抽出していると見なすことができる点は、変形例１と同様である。ミキサ７２で周波数変換された奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）は、次の数式３９で表すことができる。
【数２８】
2018101861.tif 000030
【００９２】
数式３９において、低バンド帯ＢＷＬ（図２参照）側の奇数次歪フィードバック信号＾ｙＬ（ｎ）からは、高バンド帯ＢＷＨ（図２参照）側の歪成分ΔｙＨ（ｎ）が雑音とみなせる。このため、数式３９を奇数次歪フィードバック信号＾ｙＬ（ｎ）について解くと、近似的に次の数式４０が成り立つ。
【数２９】
2018101861.tif 000031
【００９３】
そのため、生成部９１ｂにおいて、減算器８５ｂは、奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）から高バンド帯ＢＷＨ側の送信信号ｘＨ（ｎ）を減算することで、奇数次歪フィードバック信号＾ｙＬ（ｎ）を生成することができる。
【００９４】
遅延調整部７６１は、減算器８５ｂから出力された奇数次歪フィードバック信号＾ｙＬ（ｎ）を所定の遅延量で遅延させて奇数次歪補償係数更新部３１へ供給する。
【００９５】
数式３９において、高バンド帯ＢＷＨ（図２参照）側の奇数次歪フィードバック信号＾ｙＨ（ｎ）からは、低バンド帯ＢＷＬ（図２参照）側の歪成分ΔｙＬ（ｎ）が雑音とみなせる。このため、数式３９を奇数次歪フィードバック信号＾ｙＨ（ｎ）について解くと、近似的に次の数式４１が成り立つ。
【数３０】
2018101861.tif 000032
【００９６】
そのため、生成部９１ｂにおいて、減算器８６ｂは、奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）から低バンド帯ＢＷＬ側の送信信号ｘＬ（ｎ）を減算することで、奇数次歪フィードバック信号＾ｙＨ（ｎ）を生成することができる。
【００９７】
遅延調整部７６２は、減算器８６ｂから出力された奇数次歪フィードバック信号＾ｙＨ（ｎ）を所定の遅延量で遅延させて奇数次歪補償係数更新部３１へ供給する。
【００９８】
生成部９２ｂは、生成部９２ａ（図５参照）に対して、遅延調整部７７、偶数次歪アンプモデル係数算出部７９、及び偶数次歪フィードバック信号レプリカ生成部８２が省略され、遅延調整部７７１、及び遅延調整部７７２が追加されている。サンプル数を多く確保することで、実質的に、ミキサ７３が、共通フィードバック信号から偶数次歪成分を抽出していると見なすことができる点は、変形例１と同様である。ミキサ７３で周波数変換された偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）は、遅延調整部７７１及び遅延調整部７７２にそれぞれ供給される。
【００９９】
遅延調整部７７１は、偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）を所定の遅延量で遅延させ、遅延後の信号を偶数次歪フィードバック信号＾ｙＨ−Ｌ（ｎ）として偶数次歪補償係数更新部３２へ供給する。
【０１００】
遅延調整部７７２は、偶数次歪共通フィードバック信号ｙｅｖｅｎ（ｎ）を所定の遅延量で遅延させ、遅延後の信号を偶数次歪フィードバック信号＾ｙＬ，２（ｎ）として偶数次歪補償係数更新部３２へ供給する。
【０１０１】
このように、フィードバック信号レプリカを用いずに、周波数シフトした共通フィードバック信号から他方の送信信号を減算したものを用いて奇数次歪補償信号処理を行い、周波数シフトした共通フィードバック信号を用いて偶数次歪補償信号処理を行う。これにより、奇数次歪アンプモデル係数算出部７８、奇数次歪フィードバック信号レプリカ生成部８１、偶数次歪アンプモデル係数算出部７９、及び偶数次歪フィードバック信号レプリカ生成部８２等の構成を省略できる。この結果、共通フィードバック信号処理部７０ｂの構成をさらに簡略化できさらに低コスト化できる。
【０１０２】
（変形例３）
変形例２では、第１の条件及び第２の条件が両方成り立つ場合に適切に動作する。第１の条件は、減算する送信信号ｘＨ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる高バンド帯ＢＷＨ（図２参照）側の送信信号成分ｘＨ（ｎ）の振幅および位相がほぼ完全に一致することである。第２の条件は、減算する送信信号ｘＬ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる低バンド帯ＢＷＬ（図２参照）側の送信信号成分ｘＬ（ｎ）の振幅および位相がほぼ完全に一致することである。しかしながら、第１の条件及び第２の条件の少なくとも一方が成り立たない場合は適切に動作することが保証されていない。
【０１０３】
そこで、このような問題を回避するために、変形例３では、フィードバック信号レプリカに代えて送信信号を用いて、歪補償係数を計算する。具体的には、数式２８中の行列ＷＬ、ＷＨを、数式３０に代えて、次の数式４２により求める。
【数３１】
2018101861.tif 000033
【０１０４】
数式４２において、数式３０と異なる点は、絶対値の中身を奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）に代えて送信信号ｘＬ（ｎ），ｘＨ（ｎ）とする点である。
【０１０５】
すなわち、変形例２では、図７（ａ）に示すように、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）の振幅情報を用いて逆特性を求めている。それに対して、変形例３では、図７（ｂ）に示すように、送信信号ｘＬ（ｎ），ｘＨ（ｎ）の振幅情報を用いて逆特性を求める。
【０１０６】
図７（ａ）に示す奇数次歪補償係数更新部３１では、振幅算出部３１１が、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ）の振幅｜＾ｙＬ（ｎ）｜を求めて歪補償係数更新部３１３，３１４へ供給する。振幅算出部３１２が、奇数次歪フィードバック信号レプリカ＾ｙＨ（ｎ）の振幅｜＾ｙＨ（ｎ）｜を求めて歪補償係数更新部３１３，３１４へ供給する。
【０１０７】
歪補償係数更新部３１３は、奇数次歪補償信号ｚＬ（ｎ）を奇数次歪補償部３３から受ける。歪補償係数更新部３１３は、数式３０により行列ＷＬを求め、奇数次歪補償信号ｚＬ（ｎ）及び行列ＷＬを用いて、数式２８により、奇数次歪補償係数（＾Ａ）を求める。
【０１０８】
歪補償係数更新部３１４は、奇数次歪補償信号ｚＨ（ｎ）を奇数次歪補償部３３から受ける。歪補償係数更新部３１４は、数式３０により行列ＷＨを求め、奇数次歪補償信号ｚＨ（ｎ）及び行列ＷＨを用いて、数式２８により、奇数次歪補償係数（＾Ｂ）を求める。
【０１０９】
一方、図７（ｂ）に示す奇数次歪補償係数更新部３１ｃでは、振幅算出部３１１ｃが、バンドＡの送信信号ｘＬ（ｎ）の振幅｜ｘＬ（ｎ）｜を求めて歪補償係数更新部３１３，３１４へ供給する。振幅算出部３１２ｃが、バンドＢの送信信号ｘＨ（ｎ）の振幅｜ｘＨ（ｎ）｜を求めて歪補償係数更新部３１３，３１４へ供給する。
【０１１０】
歪補償係数更新部３１３は、奇数次歪補償信号ｚＬ（ｎ）を奇数次歪補償部３３から受ける。歪補償係数更新部３１３は、数式４２により行列ＷＬを求め、奇数次歪補償信号ｚＬ（ｎ）及び行列ＷＬを用いて、数式２８により、奇数次歪補償係数（＾Ａ）を求める。
【０１１１】
歪補償係数更新部３１４は、奇数次歪補償信号ｚＨ（ｎ）を奇数次歪補償部３３から受ける。歪補償係数更新部３１４は、数式４２により行列ＷＨを求め、奇数次歪補償信号ｚＨ（ｎ）及び行列ＷＨを用いて、数式２８により、奇数次歪補償係数（＾Ｂ）を求める。
【０１１２】
このように、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）に代えて送信信号ｘＬ（ｎ），ｘＨ（ｎ）を用いて歪補償係数を計算するので、より送信信号ｘＬ（ｎ），ｘＨ（ｎ）に追従させながら歪補償係数を更新できる。
【０１１３】
（変形例４）
変形例２では、第１の条件及び第２の条件が両方成り立つ場合に適切に動作する。第１の条件は、減算する送信信号ｘＨ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる高バンド帯ＢＷＨ（図２参照）側の送信信号成分ｘＨ（ｎ）の振幅および位相がほぼ完全に一致することである。第２の条件は、減算する送信信号ｘＬ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる低バンド帯ＢＷＬ（図２参照）側の送信信号成分ｘＬ（ｎ）の振幅および位相がほぼ完全に一致することである。しかしながら、第１の条件及び第２の条件の少なくとも一方が成り立たない場合は適切に動作することが保証されていない。
【０１１４】
そこで、このような問題を回避するために、変形例４では、初期値検出のために最初（１回目のサンプリング）は１次項のみで奇数次歪補償信号処理を行い、２回目以降のサンプリングでは変形例２と同様の処理を行う。
【０１１５】
具体的には、図８に示すように、奇数次歪補償部３３ｄは、歪補償部３３１、歪補償部３３２、振幅算出部３３３、振幅算出部３３４、スイッチ３３５ｄ、スイッチ３３６ｄ、スイッチ３３７ｄ、スイッチ３３８ｄ、乗算器３３９ｄ、及び乗算器３４１ｄを有する。図９に示すように、奇数次歪係数更新部３１ｄは、変形例２の歪補償係数更新部３１３、及び歪補償係数更新部３１４（図７（ａ）参照）に代えて、歪補償係数更新部３１３ｄ、及び歪補償係数更新部３１４ｄを有する。共通フィードバック信号処理部７０ｄは、変形例２の共通フィードバック信号処理部７０ｂ（図６参照）に対して、スイッチ８７ｄ、スイッチ８８ｄ、スイッチ９３ｄ、及びスイッチ９４ｄがさらに追加されている。
【０１１６】
なお、図８は、奇数次歪補償部３３ｄの構成を示す図である。図９は、奇数次歪係数更新部３１ｄの構成を示す図である。図８及び図９において、各スイッチ３３５ｄ〜３３８ｄ，８７ｄ，８８ｄ，９３ｄ，９４ｄは、１回目のサンプリングにおいて実線の状態に切り替えられている。これにより、奇数次歪係数更新部３１ｄ及び奇数次歪補償部３３ｄは、初期値を検出する処理（図１０に示すＳ４１，Ｓ４２）を行う。各スイッチ３３５ｄ〜３３８ｄ，８７ｄ，８８ｄ，９３ｄ，９４ｄは、２回目以降のサンプリングにおいて破線の状態に切り替えられている。これにより、奇数次歪係数更新部３１ｄ及び奇数次歪補償部３３ｄは変形例２と同様の処理（図１０に示すＳ４３〜Ｓ４６）を行う。図１０は、奇数次歪係数更新部３１ｄ及び奇数次歪補償部３３ｄの動作を示すフローチャートである。
【０１１７】
例えば、初期値検出時の１回目のサンプリングにおいて、低バンド帯ＢＷＬ（図２参照）側の奇数次歪補償信号ｚＬ（ｎ）と高バンド帯ＢＷＨ（図２参照）側の奇数次歪補償信号ｚＨ（ｎ）とは、それぞれ次の数式４３で表される。
【数３２】
2018101861.tif 000034
【０１１８】
数式４３において、ｙｏｄｄ（ｎ）は奇数次歪共通フィードバック信号である。このとき、ｚＬ（ｎ）、ｚＨ（ｎ）をそれぞれ行列形式で表現すると、次の数式４４で表される。
【数３３】
2018101861.tif 000035
【０１１９】
数式４４中のｚＬ、ｚＨはそれぞれ次の数式４５で表現される。
【数３４】
2018101861.tif 000036
【０１２０】
数式４４中のＵＬ、ＵＨはそれぞれ次の数式４６で表現される。
【数３５】
2018101861.tif 000037
【０１２１】
数式４４中のｈＬ、ｈＨはそれぞれ次の数式４７で表現される。
【数３６】
2018101861.tif 000038
【０１２２】
従って、初期条件として、ｚＬ（ｎ）＝ｘＬ（ｎ）、ｚＨ（ｎ）＝ｘＨ（ｎ）とし、次の数式４８に示すように係数ベクトル＾ｈＬ、＾ｈＨを求める（Ｓ４１）。
【数３７】
2018101861.tif 000039
【０１２３】
そして、数式４８で求めた係数を用いて、次の数式４９に示すように低バンド帯ＢＷＬ（図２参照）側の奇数次歪補償信号＾ｚＬ（ｎ）と高バンド帯ＢＷＨ（図２参照）側の奇数次歪補償信号＾ｚＨ（ｎ）をそれぞれ生成する（Ｓ４２）。
【数３８】
2018101861.tif 000040
【０１２４】
このように、初期値検出として１次項のみで奇数次歪補償信号処理を行うことで、減算する送信信号ｘＨ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる高バンド帯ＢＷＨ（図２参照）側の送信信号成分ｘＨ（ｎ）の振幅および位相をほぼ完全に一致させ、かつ減算する送信信号ｘＬ（ｎ）と奇数次歪共通フィードバック信号ｙｏｄｄ（ｎ）に含まれる低バンド帯ＢＷＬ（図２参照）側の送信信号成分ｘＬ（ｎ）の振幅および位相をほぼ完全に一致させたうえで歪補償係数を更新できる。
【０１２５】
（変形例５）
変形例５では、変形例３の奇数次歪係数更新部３１ｃ（図７（ｂ）参照）と変形例４の奇数次歪係数更新部３１ｄ（図９参照）とをハイブリッドさせて、奇数次歪係数更新部３１ｅを構成する。具体的には、奇数次歪係数更新部３１ｅは、図１１に示すように、振幅算出部３１１ｃ、振幅算出部３１２ｃ、歪補償係数更新部３１３ｅ、及び歪補償係数更新部３１４ｅを有する。
【０１２６】
例えば、歪補償係数更新部３１３ｅ及び歪補償係数更新部３１４ｅは、係数更新前の歪補償処理後の信号を含むベクトルｚＬ、ｚＨと共通フィードバック信号と各送信信号の奇数次項を含む行列ＷＬ、ＷＨを用いて、歪補償係数ベクトル＾Ａ、＾Ｂを更新する。歪補償係数更新部３１３ｅ及び歪補償係数更新部３１４ｅは、次の数式５０に示すように、歪補償係数ベクトル＾Ａ、＾Ｂを更新する。
【数３９】
2018101861.tif 000041
【０１２７】
数式５０中の歪補償係数ベクトル＾Ａ、＾Ｂは、それぞれ次の数式５１で表現される。
【数４０】
2018101861.tif 000042
【０１２８】
数式５０中の行列ＷＬ、ＷＨは、それぞれ次の数式５２で表現される。
【数４１】
2018101861.tif 000043
【０１２９】
数式５０中のベクトルｚＬ、ｚＨは、それぞれ次の数式５３で表現される。
【数４２】
2018101861.tif 000044
【０１３０】
歪補償係数更新部３１３ｅ及び歪補償係数更新部３１４ｅは、数式５０で更新された歪補償係数ベクトル＾Ａ、＾Ｂを奇数次歪補償部３３へ供給する。奇数次歪補償部３３は、数式５０で更新された歪補償係数ベクトル＾Ａ、＾Ｂを用いて、歪補償処理後の信号＾ｚＬ（ｎ），＾ｚＨ（ｎ）を次の数式５４のように生成する。
【数４３】
2018101861.tif 000045
【０１３１】
このとき、歪補償係数ベクトルの成分のうち半分を歪補償処理後の信号生成に用いる。つまり、歪補償係数ベクトル＾Ａのうち＾ＡＬの部分を用い、歪補償係数ベクトル＾Ｂのうち＾ＢＨの部分を用いる。
【０１３２】
このように、奇数次歪フィードバック信号レプリカ＾ｙＬ（ｎ），＾ｙＨ（ｎ）に代えて送信信号ｘＬ（ｎ），ｘＨ（ｎ）を用いて歪補償係数を計算するので、より送信信号ｘＬ（ｎ），ｘＨ（ｎ）に追従させながら歪補償係数を更新できる。また、歪補償処理には直接用いない項及び係数（歪補償係数ベクトル＾ＡのうちｂＬの部分、歪補償係数ベクトル＾Ｂのうち＾ａＨの部分）を含めたべき級数に基づいて歪補償係数を更新することにより、実質的に奇数次共通フィードバック信号から他方の送信信号を、振幅及び位相を一致させたうえでキャンセルすることができる。
【符号の説明】
【０１３３】
１  歪補償装置
１０  共通フィードバック部
１１  共通ダウンコンバータ
１２  共通ダウンコンバータ
１３  ＡＤＣ
７０  共通フィードバック信号処理部
７２  ミキサ
７３  ミキサ
７４  ＬＰＦ
７５  ＬＰＦ
８５ｂ  減算器
８６ｂ  減算器
９１，９１ａ，９１ｂ  生成部
９２，９２ａ，９２ｂ  生成部
【図１】
2018101861.tif 000046
【図２】
2018101861.tif 000047
【図３】
2018101861.tif 000048
【図４】
2018101861.tif 000049
【図５】
2018101861.tif 000050
【図６】
2018101861.tif 000051
【図７】
2018101861.tif 000052
【図８】
2018101861.tif 000053
【図９】
2018101861.tif 000054
【図１０】
2018101861.tif 000055
【図１１】
2018101861.tif 000056