(57)【要約】      （修正有）【課題】歪補償性能を向上する。【解決手段】歪補償装置１は、第１の送信系３０の送信信号を増幅する第１のパワーアンプ３３の出力に応じた第１の信号と、第２の送信系５０の送信信号を増幅する第２のパワーアンプ５３の出力に応じた第２の信号とを合成したフィードバック信号を生成する共通フィードバック部９０と、生成されたフィードバック信号を用いて、第１のパワーアンプ及び第２のパワーアンプの歪を補償する歪補償信号処理を行う歪補償信号処理部１０とを有する。【選択図】図１
2018113626.tif 000002
【特許請求の範囲】
【請求項１】
第１の送信系の送信信号を増幅する第１のパワーアンプの出力に応じた第１の信号と、第２の送信系の送信信号を増幅する第２のパワーアンプの出力に応じた第２の信号とを合成したフィードバック信号を生成するフィードバック部と、
前記生成されたフィードバック信号を用いて、前記第１のパワーアンプ及び前記第２のパワーアンプそれぞれに対して、歪を補償する歪補償信号処理をそれぞれ行う歪補償信号処理部と
を有することを特徴とする歪補償装置。
【請求項２】
前記フィードバック部は、前記第１のパワーアンプの出力の一部が取り出されて得られる信号と、前記第２のパワーアンプの出力の一部が取り出されて得られる信号とを合成した前記フィードバック信号を生成する信号合成部を有する
ことを特徴とする請求項１に記載の歪補償装置。
【請求項３】
前記フィードバック部は、前記第１のパワーアンプの出力が前記第１の送信系のアンテナから無線送信されて得られる電波と、前記第２のパワーアンプの出力が前記第２の送信系のアンテナから無線送信されて得られる電波とを受信して前記フィードバック信号を生成するアンテナを有する
ことを特徴とする請求項１に記載の歪補償装置。
【請求項４】
前記歪補償信号処理部は、
前記フィードバック部で生成されたフィードバック信号から前記第１の信号のレプリカ及び前記第２の信号のレプリカをそれぞれ分離する信号分離部と、
前記信号分離部で分離された前記第１の信号のレプリカを用いて、前記第１のパワーアンプの歪補償に用いられる第１の歪補償係数を更新する第１の歪補償係数更新部と、
前記信号分離部で分離された前記第２の信号のレプリカを用いて、前記第２のパワーアンプの歪補償に用いられる第２の歪補償係数を更新する第２の歪補償係数更新部と、
前記第１の送信系の送信信号と、前記第１の歪補償係数とを用いて、前記第１のパワーアンプの非線形歪を補償する第１の歪補償信号を生成する第１の歪補償部と、
前記第２の送信系の送信信号と、前記第２の歪補償係数とを用いて、前記第２のパワーアンプの非線形歪を補償する第２の歪補償信号を生成する第２の歪補償部と
を有することを特徴とする請求項１〜３のいずれか一つに記載の歪補償装置。
【請求項５】
前記信号分離部は、
前記フィードバック部で生成されたフィードバック信号と、前記第１の送信系の送信信号と、前記第２の送信系の送信信号とを用いて、前記第１のパワーアンプのアンプモデル係数と、前記第２のパワーアンプのアンプモデル係数とを算出する係数算出部と、
前記係数算出部で算出された前記第１のパワーアンプのアンプモデル係数及び前記第２のパワーアンプのアンプモデル係数を用いて、前記第１の信号のレプリカと前記第２の信号のレプリカとを生成するレプリカ生成部と
を有し、
前記第１の歪補償係数更新部は、
前記レプリカ生成部で生成された前記第１の信号のレプリカを用いて、前記第１の歪補償係数を更新し、
前記第２の歪補償係数更新部は、
前記レプリカ生成部で生成された前記第２の信号のレプリカを用いて、前記第２の歪補償係数を更新する
ことを特徴とする請求項４に記載の歪補償装置。
【請求項６】
前記信号分離部は、
前記フィードバック部で生成されたフィードバック信号から前記第２の送信系の送信信号を減算して前記第１の信号のレプリカを生成する第１の減算部と、
前記フィードバック部で生成されたフィードバック信号から前記第１の送信系の送信信号を減算して前記第２の信号のレプリカを生成する第２の減算部と
を有し、
前記第１の歪補償係数更新部は、
前記第１の減算部で生成された前記第１の信号のレプリカを用いて、前記第１の歪補償係数を更新し、
前記第２の減算部で生成された前記第２の信号のレプリカを用いて、前記第２の歪補償係数を更新する
ことを特徴とする請求項４に記載の歪補償装置。
【請求項７】
前記第１の歪補償係数更新部は、
前記第１の歪補償係数を更新する場合に、前記第１の送信系の送信信号の振幅を算出し、前記第１の信号のレプリカの振幅に代えて前記第１の送信系の送信信号の振幅を用いて、更新後の前記第１の歪補償係数を算出し、
前記第２の歪補償係数更新部は、
前記第２の歪補償係数を更新する場合に、前記第２の送信系の送信信号の振幅を算出し、前記第２の信号のレプリカの振幅に代えて前記第２の送信系の送信信号の振幅を用いて、更新後の前記第２の歪補償係数を算出する
ことを特徴とする請求項６に記載の歪補償装置。
【請求項８】
前記第１の歪補償係数更新部は、
前記第１の送信系の送信信号及び前記第２の送信系の送信信号の１回目のサンプリングにおいて、前記フィードバック部で生成されたフィードバック信号と、前記第２の送信系の送信信号とを用いて、前記第１の歪補償係数の初期値を算出し、前記第１の送信系の送信信号及び前記第２の送信系の送信信号の２回目以降のサンプリングにおいて、前記第１の減算部で生成された前記第１の信号のレプリカを用いて、前記第１の歪補償係数を更新し、
前記第２の歪補償係数更新部は、
前記１回目のサンプリングにおいて、前記フィードバック部で生成されたフィードバック信号と、前記第１の送信系の送信信号とを用いて、前記第２の歪補償係数の初期値を算出し、前記２回目以降のサンプリングにおいて、前記第２の減算部で生成された前記第２の信号のレプリカを用いて、前記第２の歪補償係数を更新し、
前記第１の歪補償部は、
前記１回目のサンプリングにおいて、前記第１の歪補償係数の初期値を用いて、前記第１の送信系の送信信号の１次項のみを含む前記第１の歪補償信号を生成し、前記２回目以降のサンプリングにおいて、前記第１の歪補償係数更新部で更新された前記第１の歪補償係数を用いて、前記１次項と前記１次項以外の項とを含む前記第１の歪補償信号を生成し、
前記第２の歪補償部は、
前記１回目のサンプリングにおいて、前記第２の歪補償係数の初期値を用いて、前記第２の送信系の送信信号の１次項のみを含む前記第２の歪補償信号を生成し、前記２回目以降のサンプリングにおいて、前記第２の歪補償係数更新部で更新された前記第２の歪補償係数を用いて、前記１次項と前記１次項以外の項とを含む前記第２の歪補償信号を生成する
ことを特徴とする請求項６に記載の歪補償装置。
【請求項９】
前記歪補償信号処理部は、
前記フィードバック部で生成されたフィードバック信号と、前記第１の送信系の送信信号と、前記第２の送信系の送信信号とを用いて、前記第１のパワーアンプの歪補償に用いられる第１の歪補償係数を更新する第１の歪補償係数更新部と、
前記フィードバック部で生成されたフィードバック信号と、前記第１の送信系の送信信号と、前記第２の送信系の送信信号とを用いて、前記第２のパワーアンプの歪補償に用いられる第２の歪補償係数を更新する第２の歪補償係数更新部と、
前記第１の送信系の送信信号と、前記第１の歪補償係数を用いて、前記第１のパワーアンプの非線形歪を補償する第１の歪補償信号を生成する第１の歪補償部と、
前記第２の送信系の送信信号と、前記第２の歪補償係数を用いて、前記第２のパワーアンプの非線形歪を補償する第２の歪補償信号を生成する第２の歪補償部と
を有することを特徴とする請求項１〜３のいずれか一つに記載の歪補償装置。
【請求項１０】
第１の送信系の送信信号を増幅する第１のパワーアンプの出力に応じた第１の信号と、第２の送信系の送信信号を増幅する第２のパワーアンプの出力に応じた第２の信号とを合成したフィードバック信号を生成し、
前記生成されたフィードバック信号を用いて、前記第１のパワーアンプ及び前記第２のパワーアンプの歪を補償する歪補償信号処理を行う
ことを特徴とする歪補償方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、歪補償装置及び歪補償方法に関する。
【背景技術】
【０００２】
無線移動体通信の分野では、複数のアンテナを利用する通信方式として、ＭＩＭＯ（Multiple-Input  and  Multiple-Output）方式が知られている。ＭＩＭＯ方式が適用された基地局等の通信装置は、複数のアンテナに対応する複数の送信系（複数の送信回路）を有し、各送信系の送信信号を、各送信系に設けられた各パワーアンプにより増幅する。また、無線移動体通信の分野では、基地局等の歪補償装置に用いられるパワーアンプの非線形歪を補償する技術が提案されている。非線形歪を補償する歪補償方式の一つに「プリディストーション（以下では「ＰＤ」と呼ぶことがある）方式」がある。ＰＤ方式においては、パワーアンプの非線形歪の逆特性を有する歪補償係数がパワーアンプに入力される送信信号に予め乗算される。これにより、パワーアンプの出力の線形性が高まり、パワーアンプの出力の歪が抑制される。また、歪補償性能を高めるために、パワーアンプの出力の一部をフィードバックして、歪補償係数を適応的に更新する方式が知られている。
【０００３】
上記のＭＩＭＯ方式のように複数のアンテナが用いられる場合、歪補償装置は、各アンテナに対応する各送信系に設けられた各パワーアンプの出力から一部を取り出してフィードバックした信号を用いて、各パワーアンプに対応する歪補償係数を更新する。このように、各パワーアンプの出力から取り出される信号を用いて歪補償係数が更新される場合、一般に、各パワーアンプの出力から取り出される信号をフィードバックするための回路がパワーアンプの数に応じて複数設けられることになる。ただし、各パワーアンプの出力から取り出される信号をフィードバックするための回路がパワーアンプの数に応じて複数設けられる場合、回路規模が増大してしまう可能性がある。
【０００４】
これに対し、回路規模を削減する観点から、各パワーアンプの出力から取り出される信号をスイッチにより時分割で切り替えて、各パワーアンプに対応する歪補償係数を更新する技術が知られている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許５１５１７８５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、各パワーアンプの出力から取り出される信号をスイッチにより時分割で切り替える技術では、各パワーアンプに関して信号が選択されない期間が存在するため、この期間において、各パワーアンプに対応する歪補償係数の更新が不十分になる。結果として、歪補償性能が低下する恐れがある。
【０００７】
開示の技術は、歪補償性能を向上することができる歪補償装置及び歪補償方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本願の開示する歪補償装置は、一つの態様において、第１の送信系の送信信号を増幅する第１のパワーアンプの出力に応じた第１の信号と、第２の送信系の送信信号を増幅する第２のパワーアンプの出力に応じた第２の信号とを合成したフィードバック信号を生成するフィードバック部と、生成されたフィードバック信号を用いて、第１のパワーアンプ及び第２のパワーアンプそれぞれに対して、歪を補償する歪補償信号処理をそれぞれ行う歪補償信号処理部とを有する。
【発明の効果】
【０００９】
本願の開示する歪補償装置の一つの態様によれば、歪補償性能を向上することができるという効果を奏する。
【図面の簡単な説明】
【００１０】
【図１】図１は、本実施例に係る歪補償装置の構成を示す図である。
【図２】図２は、本実施例における共通フィードバック歪補償信号処理部の構成を示す図である。
【図３】図３は、本実施例に係る歪補償装置の動作を示すフローチャートである。
【図４】図４は、変形例１における共通フィードバック歪補償信号処理部の構成を示す図である。
【図５】図５は、変形例２における共通フィードバック歪補償信号処理部の構成を示す図である。
【図６】図６は、変形例３における共通フィードバック歪補償信号処理部の構成を示す図である。
【図７】図７は、変形例４における共通フィードバック歪補償信号処理部の構成を示す図である。
【図８】図８は、変形例４における歪補償係数更新部及び歪補償部の動作を示すフローチャートである。
【図９】図９は、変形例５における共通フィードバック歪補償信号処理部の構成を示す図である。
【図１０】図１０は、変形例６に係る歪補償装置の構成を示す図である。
【発明を実施するための形態】
【００１１】
以下に、本願の開示する歪補償装置及び歪補償方法の実施例を図面に基づいて詳細に説明する。なお、この実施例により開示技術が限定されるものではない。
【実施例】
【００１２】
図１は、本実施例に係る歪補償装置１の構成を示す図である。歪補償装置１は、ＭＩＭＯ方式が適用された基地局等の歪補償装置である。歪補償装置１は、複数のアンテナに対応する複数の送信系（複数の送信回路）を有し、各送信系の送信信号を、各送信系に設けられた各パワーアンプにより増幅する。図１に示すように、歪補償装置１は、共通フィードバック歪補償信号処理部１０、送信回路３０，５０、発振器７０及び共通フィードバック部９０を有する。なお、以下では、説明を簡潔にするために、歪補償装置１のアンテナの数が２であり、かつ、送信系（送信回路）の数が２であるものとしているが、これらの数はこれに限定されるものではない。
【００１３】
図１では、送信回路３０の送信信号及び送信回路５０の送信信号を、各送信回路３０，５０に設けられた各パワーアンプ３３，５３により増幅する場合を例にとり説明する。共通フィードバック歪補償信号処理部１０は、送信回路３０の送信信号ｘ１（ｎ）及び送信回路５０の送信信号ｘ２（ｎ）の２つの信号を受け、該信号に基づき、送信回路３０の歪補償信号＾ｚ１（ｎ）及び送信回路５０の歪補償信号＾ｚ２（ｎ）を生成する。
【００１４】
送信回路３０は、アンテナＡ１に対応する送信系であり、ＤＡＣ（Digital  to  Analog  Converter）３１、アップコンバータ３２、パワーアンプ３３及びカップラ３４を有する。以下では、送信回路３０を適宜「第１の送信系」と呼ぶことがある。ＤＡＣ３１は、共通フィードバック歪補償信号処理部１０から入力される第１の送信系の歪補償信号＾ｚ１（ｎ）をデジタル信号からアナログ信号に変換する。
【００１５】
アップコンバータ３２は、ＤＡＣ３１から入力されるアナログ信号を発振器７０からの所定周波数のローカル信号に乗算させて周波数変換し、ＲＦ（Radio  Frequency）信号を生成する。
【００１６】
アップコンバータ３２からのＲＦ信号は、パワーアンプ３３に入力される。パワーアンプ３３に入力されるＲＦ信号にはパワーアンプ３３の非線形歪の逆特性に応じた成分が含まれているので、ＲＦ信号（第１の送信系の送信信号）がパワーアンプ３３で電力増幅される際の歪が抑制される。すなわち、歪補償装置１は、パワーアンプ３３の非線形歪を補償することができる。パワーアンプ３３は、第１のパワーアンプの一例である。
【００１７】
カップラ３４は、パワーアンプ３３の出力端に具備されており、パワーアンプ３３の出力の一部を取り出す。パワーアンプ３３の出力の一部がカップラ３４で取り出されて得られる信号は、共通フィードバック部９０へフィードバックされ、該信号を除くパワーアンプ３３の出力は、アンテナＡ１から無線送信される。以下では、パワーアンプ３３の出力の一部がカップラ３４で取り出されて得られる信号を適宜「第１の送信系のアンプ出力信号」と呼ぶことがある。第１の送信系のアンプ出力信号は、パワーアンプ３３の出力に応じた第１の信号の一例である。
【００１８】
送信回路５０は、アンテナＡ２に対応する送信系であり、ＤＡＣ５１、アップコンバータ５２、パワーアンプ５３及びカップラ５４を有する。以下では、送信回路５０を適宜「第２の送信系」と呼ぶことがある。ＤＡＣ５１は、共通フィードバック歪補償信号処理部１０から入力される第２の送信系の歪補償信号＾ｚ２（ｎ）をデジタル信号からアナログ信号に変換する。
【００１９】
アップコンバータ５２は、ＤＡＣ５１から入力されるアナログ信号を発振器７０からの所定周波数のローカル信号に乗算させて周波数変換し、ＲＦ信号を生成する。
【００２０】
アップコンバータ５２からのＲＦ信号は、パワーアンプ５３の非線形歪の逆特性に応じた成分が含まれているので、ＲＦ信号（第２の送信系の送信信号）がパワーアンプ５３で電力増幅される際の歪が抑制される。すなわち、歪補償装置１は、パワーアンプ５３の非線形歪を補償することができる。パワーアンプ５３は、第２のパワーアンプの一例である。
【００２１】
カップラ５４は、パワーアンプ５３の出力端に具備されており、パワーアンプ５３の出力の一部を取り出す。パワーアンプ５３の出力の一部がカップラ５４で取り出されて得られる信号は、共通フィードバック部９０へフィードバックされ、該信号を除くパワーアンプ５３の出力は、アンテナＡ２から無線送信される。以下では、パワーアンプ５３の出力の一部がカップラ５４で取り出されて得られる信号を適宜「第２の送信系のアンプ出力信号」と呼ぶことがある。第２の送信系のアンプ出力信号は、パワーアンプ５３の出力に応じた第２の信号の一例である。
【００２２】
本実施例では、歪補償装置１において、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とが合成したフィードバック信号を生成して、アンプ出力信号が選択されない期間を無くすことで、歪補償性能の向上を図る。
【００２３】
共通フィードバック部９０は、第１の送信系の送信信号を増幅するパワーアンプ３３の出力に応じた第１の信号と、第２の送信系の送信信号を増幅するパワーアンプ５３の出力に応じた第２の信号とを合成した共通フィードバック信号ｙ（ｎ）を生成する。すなわち、共通フィードバック部９０は、パワーアンプ３３の出力がカップラ３４で取り出され、かつ、パワーアンプ５３の出力がカップラ５４で取り出される際に、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とを重ねる。ＭＩＭＯ方式が適用された歪補償装置１において、第１の送信系のアンプ出力信号の周波数ｆｃと、第２の送信系のアンプ出力信号の周波数ｆｃとは同一である。このため、共通フィードバック部９０は、第１の送信系のアンプ出力信号及び第２の送信系のアンプ出力信号に共通する周波数ｆｃにおいて、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とを重ねる。これにより、共通フィードバック部９０は、共通フィードバック信号ｙ（ｎ）を生成する。共通フィードバック部９０は、生成された共通フィードバック信号ｙ（ｎ）を共通フィードバック歪補償信号処理部１０へ供給する。
【００２４】
共通フィードバック部９０は、信号合成器９１、ダウンコンバータ９２及びＡＤＣ（Analog  to  Digital  Converter）９３を有する。
【００２５】
信号合成器９１は、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とを合成した共通フィードバック信号ｙ（ｎ）を生成する。これにより、共通フィードバック信号ｙ（ｎ）には、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とが重ねられる。
【００２６】
ダウンコンバータ９２は、信号合成器９１で生成された共通フィードバック信号ｙ（ｎ）に対して発振器７０からの所定周波数のローカル信号を乗算して、共通フィードバック信号ｙ（ｎ）をＲＦ周波数からベースバンド周波数にダウンコンバートする。
【００２７】
ＡＤＣ９３は、ダウンコンバータ９２でベースバンド周波数にダウンコンバートされた共通フィードバック信号ｙ（ｎ）をアナログ形式からデジタル形式に変換し、変換後の共通フィードバック信号ｙ（ｎ）を共通フィードバック歪補償信号処理部１０へ供給する。
【００２８】
図１に示す共通フィードバック歪補償信号処理部１０は、共通フィードバック信号ｙ（ｎ）を用いて、パワーアンプ３３及びパワーアンプ５３の非線形歪を補償する歪補償信号処理を行う。
【００２９】
例えば、共通フィードバック歪補償信号処理部１０は、図２に示すように、信号分離部１１、歪補償係数更新部１２，１３及び歪補償部１４，１５を有する。図２は、本実施例における共通フィードバック歪補償信号処理部１０の構成を示す図である。
【００３０】
信号分離部１１は、共通フィードバック信号ｙ（ｎ）から第１の送信系のアンプ出力信号のレプリカ及び第２の送信系のアンプ出力信号のレプリカをそれぞれ分離する。
【００３１】
信号分離部１１は、図２に示すように、アンプモデル係数算出部２１及びフィードバック信号レプリカ生成部２２を有する。
【００３２】
アンプモデル係数算出部２１は、共通フィードバック信号ｙ（ｎ）を共通フィードバック部９０のＡＤＣ９３から受ける。アンプモデル係数算出部２１は、第１の送信系の送信信号ｘ１（ｎ）及び第２の送信系の送信信号ｘ２（ｎ）を受ける。アンプモデル係数算出部２１は、共通フィードバック信号ｙ（ｎ）と、第１の送信系の送信信号ｘ１（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とを用いて、パワーアンプ３３のアンプモデル係数と、パワーアンプ５３のアンプモデル係数とを算出する。
【００３３】
ここで、簡単のために共通フィードバック歪補償信号処理部１０がメモリレスの構成であると仮定すると、アンプモデル係数算出部２１が受ける共通フィードバック信号ｙ（ｎ）は、次の数式（１）で表される。なお、本実施例では、共通フィードバック歪補償信号処理部１０がメモリレスの構成であると仮定したが、共通フィードバック歪補償信号処理部１０は、一般的に知られているメモリ多項式の構成を用いることもできる。
【００３４】
【数１】
2018113626.tif 000003
【００３５】
数式（１）において、ｘ１（ｎ）は、第１の送信系の送信信号であり、ｘ２（ｎ）は、第２の送信系の送信信号である。また、ｙ１（ｎ）は、第１の送信系のアンプ出力信号を表し、ｙ２（ｎ）は、第２の送信系のアンプ出力信号を表す。また、ａｐは、パワーアンプ３３の特性に対応したアンプモデルから導かれるアンプモデル係数であり、ｂｐは、パワーアンプ５３の特性に対応したアンプモデルから導かれるアンプモデル係数である。パワーアンプ３３の特性及びパワーアンプ５３の特性にそれぞれ対応した２つのアンプモデルは、アンプモデル係数算出部２１に予め設定されている。このとき、ｙ（ｎ）を行列形式で表現すると、次の数式（２）で表される。
【００３６】
【数２】
2018113626.tif 000004
【００３７】
数式（２）中のｙは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における共通フィードバック信号を含むベクトルであり、次の数式（３）で表現される。
【００３８】
【数３】
2018113626.tif 000005
【００３９】
数式（２）中のＵは、サンプリングタイミングｎ〜ｎ＋Ｎ−１における各送信系の送信信号の奇数次項を含む行列であり、次の数式（４）で表現される。
【００４０】
【数４】
2018113626.tif 000006
【００４１】
数式（２）中のｈは、パワーアンプ３３の伝達関数及びパワーアンプ５３の伝達関数のベクトルであり、次の数式（５）で表現される。
【００４２】
【数５】
2018113626.tif 000007
【００４３】
数式（２）の行列方程式をｈについて解くと、次の数式（６）が得られる。数式（６）では、ｈを特にベクトル＾ｈと表している。
【００４４】
【数６】
2018113626.tif 000008
【００４５】
数式（６）において、添え字Ｈは、エルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は、逆行列を表す。
【００４６】
すなわち、図２に示すアンプモデル係数算出部２１は、共通フィードバック信号を含むベクトルｙと、各送信系の送信信号の奇数次項を含む行列Ｕとを用いて、数式（６）に示すように係数ベクトル＾ｈを算出する。係数ベクトル＾ｈは、数式（５）に示すように、パワーアンプ３３のアンプモデル係数＾ａｐ（０≦ｐ≦Ｐ）と、パワーアンプ５３のアンプモデル係数＾ｂｐ（０≦ｐ≦Ｐ）とを含む。アンプモデル係数算出部２１は、算出された係数ベクトル＾ｈをフィードバック信号レプリカ生成部２２へ供給する。
【００４７】
フィードバック信号レプリカ生成部２２は、パワーアンプ３３のアンプモデル係数＾ａｐ（０≦ｐ≦Ｐ）を用いて、次の数式（７）に示すように第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を生成する。第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）は、共通フィードバック信号ｙ（ｎ）から分離された第１の送信系のアンプ出力信号のレプリカに相当する。
【００４８】
【数７】
2018113626.tif 000009
【００４９】
フィードバック信号レプリカ生成部２２は、パワーアンプ５３のアンプモデル係数＾ｂｐ（０≦ｐ≦Ｐ）を用いて、次の数式（８）に示すように第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を生成する。第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）は、共通フィードバック信号ｙ（ｎ）から分離された第２の送信系のアンプ出力信号のレプリカに相当する。
【００５０】
【数８】
2018113626.tif 000010
【００５１】
フィードバック信号レプリカ生成部２２は、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を歪補償係数更新部１２へ供給し、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を歪補償係数更新部１３へ供給する。
【００５２】
歪補償係数更新部１２は、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）をフィードバック信号レプリカ生成部２２から受ける。歪補償係数更新部１２は、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を用いて、第１の送信系の歪補償信号処理を行う。すなわち、歪補償係数更新部１２は、次の数式（９）に示す計算を行い、パワーアンプ３３の非線形歪の補償に用いられる第１の歪補償係数（歪補償係数ベクトル＾Ａ）を更新する。
【００５３】
【数９】
2018113626.tif 000011
【００５４】
数式（９）中の歪補償係数ベクトル＾Ａは、更新後の第１の歪補償係数を示し、次の数式（１０）のように表現される。
【００５５】
【数１０】
2018113626.tif 000012
【００５６】
数式（９）中の行列Ｗ１は、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）により求められる。行列Ｗ１は、第１の送信系のフィードバック信号レプリカの奇数次項を含む行列であり、次の数式（１１）のように表現される。
【００５７】
【数１１】
2018113626.tif 000013
【００５８】
なお、数式（９）において、添え字Ｈは、エルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は、逆行列を表す。
【００５９】
数式（９）中のベクトルｚ１は、第１の送信系の送信信号ｘ１（ｎ）により求められ、係数更新前の歪補償処理後の信号を含むベクトルであり、次の数式（１２）で表現される。
【００６０】
【数１２】
2018113626.tif 000014
【００６１】
歪補償係数更新部１２は、更新後の第１の歪補償係数（歪補償係数ベクトル＾Ａ）を歪補償部１４へ供給する。
【００６２】
歪補償部１４は、第１の送信系の送信信号ｘ１（ｎ）を受ける。歪補償部１４は、次の数式（１３）に示すように、第１の送信系の送信信号ｘ１（ｎ）と、第１の歪補償係数（＾Ａ）とを用いて、第１の送信系の歪補償信号＾ｚ１（ｎ）を生成する。
【００６３】
【数１３】
2018113626.tif 000015
【００６４】
歪補償係数更新部１３は、第１の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）をフィードバック信号レプリカ生成部２２から受ける。歪補償係数更新部１３は、第１の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を用いて、第２の送信系の歪補償信号処理を行う。すなわち、歪補償係数更新部１３は、次の数式（１４）に示す計算を行い、パワーアンプ５３の非線形歪の補償に用いられる第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を更新する。
【００６５】
【数１４】
2018113626.tif 000016
【００６６】
数式（１４）中の歪補償係数ベクトル＾Ｂは、更新後の第２の歪補償係数を示し、次の数式（１５）のように表現される。
【００６７】
【数１５】
2018113626.tif 000017
【００６８】
数式（１４）中の行列Ｗ２は、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）により求められる。行列Ｗ２は、第２の送信系のフィードバック信号レプリカの奇数次項を含む行列であり、次の数式（１６）のように表現される。
【００６９】
【数１６】
2018113626.tif 000018
【００７０】
なお、数式（１４）において、添え字Ｈは、エルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は、逆行列を表す。
【００７１】
数式（１４）中のベクトルｚ２は、第２の送信系の送信信号ｘ２（ｎ）により求められ、係数更新前の歪補償処理後の信号を含むベクトルであり、次の数式（１７）で表現される。
【００７２】
【数１７】
2018113626.tif 000019
【００７３】
歪補償係数更新部１３は、更新後の第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を歪補償部１５へ供給する。
【００７４】
歪補償部１５は、第２の送信系の送信信号ｘ２（ｎ）を受ける。歪補償部１５は、次の数式（１８）に示すように、第２の送信系の送信信号ｘ２（ｎ）と、第２の歪補償係数（＾Ｂ）とを用いて、第２の送信系の歪補償信号＾ｚ２（ｎ）を生成する。
【００７５】
【数１８】
2018113626.tif 000020
【００７６】
次に、歪補償装置１の動作について図３を用いて説明する。図３は、本実施例に係る歪補償装置１の動作を示すフローチャートである。
【００７７】
歪補償装置１において、共通フィードバック部９０は、パワーアンプ３３の出力及びパワーアンプ５３の出力に応じた第１の送信系のアンプ出力信号及び第２の送信系のアンプ出力信号を受けると、共通フィードバック信号ｙ（ｎ）を生成する（Ｓ１１）。共通フィードバック信号ｙ（ｎ）には、第１の送信系のアンプ出力信号と第２の送信系のアンプ出力信号とが重ねられる。共通フィードバック歪補償信号処理部１０が共通フィードバック信号ｙ（ｎ）を共通フィードバック部９０から受けると、共通フィードバック信号ｙ（ｎ）を用いる歪補償信号処理（Ｓ１２〜Ｓ１４）が行われる。
【００７８】
共通フィードバック歪補償信号処理部１０の信号分離部１１は、共通フィードバック信号ｙ（ｎ）と各送信系の送信信号とに応じて、数式（６）に示す計算により、アンプモデル係数＾ａｐ，＾ｂｐを含む係数ベクトルを算出する（Ｓ１２）。
【００７９】
信号分離部１１は、アンプモデル係数＾ａｐ，＾ｂｐを用いて、数式（７）、数式（８）に示す計算により、各送信系に対応したフィードバック信号レプリカ＾ｙ１（ｎ），＾ｙ２（ｎ）を生成する（Ｓ１３）。信号分離部１１は、フィードバック信号レプリカ＾ｙ１（ｎ），＾ｙ２（ｎ）を歪補償係数更新部１２，１３へそれぞれ供給する。
【００８０】
歪補償係数更新部１２は、フィードバック信号レプリカ＾ｙ１（ｎ）を用いて、第１の歪補償係数を更新し、更新後の第１の歪補償係数（＾Ａ）を歪補償部１４へ供給する。歪補償係数更新部１３は、フィードバック信号レプリカ＾ｙ２（ｎ）を用いて、第２の歪補償係数を更新し、更新後の第２の歪補償係数（＾Ｂ）を歪補償部１５へ供給する。歪補償部１４，１５は、更新後の第１の歪補償係数及び更新後の第２の歪補償係数（＾Ａ，＾Ｂ）を用いて、パワーアンプ３３及びパワーアンプ５３の非線形歪を補償するＤＰＤ（デジタルプリディストーション）処理を行う（Ｓ１４）。
【００８１】
歪補償装置１は、終了の指示を受けるまで（Ｓ１５でＮｏ）、Ｓ１１〜Ｓ１４までの処理を繰り返し、終了の指示を受け取ると（Ｓ１５でＹｅｓ）、処理を終了する。
【００８２】
以上のように、本実施例では、歪補償装置１は、共通フィードバック部９０と、共通フィードバック歪補償信号処理部１０とを有する。共通フィードバック部９０は、第１の送信系の送信信号を増幅するパワーアンプ３３の出力に応じた第１の信号と、第２の送信系の送信信号を増幅するパワーアンプ５３の出力に応じた第２の信号とを合成した共通フィードバック信号を生成する。共通フィードバック歪補償信号処理部１０は、共通フィードバック部９０で生成された共通フィードバック信号を用いて、パワーアンプ３３及びパワーアンプ５３それぞれの非線形歪を同時に補償する歪補償処理を行う。
【００８３】
この歪補償装置１の構成により、各パワーアンプ３３、５３に関してアンプ出力信号が選択されない期間を無くすことができ、各パワーアンプ３３、５３に対応する歪補償係数の更新が継続的に行われるので、各パワーアンプにおける歪成分の変動が多い場合に歪補償性能を向上できる。また、アンプ出力信号を時分割で切り替えるスイッチや、このスイッチを制御するための回路が省略されるので、歪補償装置１におけるフィードバック系の回路規模を低減することができる。
【００８４】
（変形例１）
なお、共通フィードバック歪補償信号処理部１０ａの信号分離部１１ａは、図４に示すように、アンプモデル係数を用いずにフィードバック信号レプリカを生成する構成であってもよい。この場合、共通フィードバック歪補償信号処理部１０ａの歪補償係数更新部１２ａ，１３ａは、図４に示すように、フィードバック信号レプリカを用いて、ダイレクトラーニング方式によって歪補償係数を更新する構成であってもよい。図４は、変形例１における共通フィードバック歪補償信号処理部１０ａの構成を示す図である。例えば、第１の送信系の歪成分と第２の送信系の歪成分とが互いに相関性が低いと想定できると、第１の送信系の歪成分及び第２の送信系の歪成分において一方が他方の雑音に見える。このため、共通フィードバック部９０で生成される共通フィードバック信号ｙ（ｎ）のサンプル数（例えば数式（３）、数式（４）におけるＮ）を多く確保すれば、他方の歪成分をキャンセルできる。
【００８５】
具体的には、信号分離部１１ａは、信号分離部１１（図２参照）におけるアンプモデル係数算出部２１及びフィードバック信号レプリカ生成部２２に代えて、減算器２３及び２４を有する。また、歪補償係数更新部１２ａには、減算器１１６が接続され、歪補償係数更新部１３ａには、減算器１１７が接続されている。共通フィードバック部９０で生成される共通フィードバック信号ｙ（ｎ）は、次の数式（１９）で表すことができる。
【００８６】
【数１９】
2018113626.tif 000021
【００８７】
数式（１９）において、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）からは、第２の送信系の歪成分Δｙ２が雑音とみなせる。このため、数式（１９）を第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）について解くと、近似的に次の数式（２０）が成り立つ。
【００８８】
【数２０】
2018113626.tif 000022
【００８９】
そのため、信号分離部１１ａにおいて、減算器２３は、共通フィードバック信号ｙ（ｎ）から第２の送信系の送信信号ｘ２（ｎ）を減算することで、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を生成することができる。
【００９０】
減算器２３は、生成された第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を減算器１１６及び歪補償係数更新部１２ａへ供給する。
【００９１】
数式（１９）において、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）からは、第１の送信系の歪成分Δｙ１が雑音とみなせる。このため、数式（１９）を第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）について解くと、近似的に次の数式（２１）が成り立つ。
【００９２】
【数２１】
2018113626.tif 000023
【００９３】
そのため、信号分離部１１ａにおいて、減算器２４は、共通フィードバック信号ｙ（ｎ）から第１の送信系の送信信号ｘ１（ｎ）を減算することで、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を生成することができる。
【００９４】
減算器２４は、生成された第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を減算器１１７及び歪補償係数更新部１３ａへ供給する。
【００９５】
減算器１１６は、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）と第１の送信系の送信信号ｘ１（ｎ）との差分を算出して第１の送信系の誤差信号を得る。
【００９６】
歪補償係数更新部１２ａは、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）と、上記第１の送信系の誤差信号とを用いて、第１の歪補償係数（歪補償係数ベクトル＾Ａ）を更新する。第１の歪補償係数（歪補償係数ベクトル＾Ａ）の更新には、例えば数式（９）が用いられる。歪補償係数更新部１２ａは、更新後の第１の歪補償係数（歪補償係数ベクトル＾Ａ）を歪補償部１４へ供給する。
【００９７】
減算器１１７は、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）と第２の送信系の送信信号ｘ２（ｎ）との差分を算出して第２の送信系の誤差信号を得る。
【００９８】
歪補償係数更新部１３ａは、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）と、上記第２の送信系の誤差信号とを用いて、第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を更新する。第２の歪補償係数（歪補償係数ベクトル＾Ｂ）の更新には、例えば数式（１４）が用いられる。歪補償係数更新部１３ａは、更新後の第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を歪補償部１５へ供給する。
【００９９】
このように、変形例１によれば、共通フィードバック信号から他方の送信信号を減算して生成されたフィードバック信号レプリカを用いて、ダイレクトラーニング方式に基づく歪補償信号処理を行う。これにより、アンプモデル係数算出部２１及びフィードバック信号レプリカ生成部２２等の構成を省略することができる。この結果、共通フィードバック歪補償信号処理部１０ａの構成を簡略化することができる。
【０１００】
（変形例２）
変形例１では、共通フィードバック歪補償信号処理部１０ａの歪補償係数更新部１２ａ，１３ａは、フィードバック信号レプリカを用いて、ダイレクトラーニング方式によって歪補償係数を更新する構成を採る。これに対して、変形例２では、共通フィードバック歪補償信号処理部１０ｂの歪補償係数更新部１２ｂ，１３ｂは、図５に示すように、フィードバック信号レプリカを用いて、インダイレクトラーニング方式によって歪補償係数を更新する構成を採る。信号分離部１１ａにおいて、減算器２３が、共通フィードバック信号ｙ（ｎ）から第２の送信系の送信信号ｘ２（ｎ）を減算することで、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を生成する点は、変形例１と同様である。信号分離部１１ａにおいて、減算器２４が、共通フィードバック信号ｙ（ｎ）から第１の送信系の送信信号ｘ１（ｎ）を減算することで、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を生成する点は、変形例１と同様である。図５は、変形例２における共通フィードバック歪補償信号処理部１０ｂの構成を示す図である。
【０１０１】
図５に示すように、歪補償係数更新部１２ｂは、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）と、係数更新前の歪補償処理後の信号ｚ１（ｎ）とを用いて、第１の歪補償係数（歪補償係数ベクトル＾Ａ）を更新する。第１の歪補償係数（歪補償係数ベクトル＾Ａ）の更新には、例えば数式（９）が用いられる。歪補償係数更新部１２ｂは、更新後の第１の歪補償係数（歪補償係数ベクトル＾Ａ）を歪補償部１４にコピーする。
【０１０２】
歪補償係数更新部１３ｂは、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）と、係数更新前の歪補償処理後の信号ｚ２（ｎ）とを用いて、第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を更新する。第２の歪補償係数（歪補償係数ベクトル＾Ｂ）の更新には、例えば数式（１４）が用いられる。歪補償係数更新部１３ｂは、更新後の第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を歪補償部１５にコピーする。
【０１０３】
このように、変形例２によれば、共通フィードバック信号から他方の送信信号を減算して生成されたフィードバック信号レプリカを用いて、インダイレクトラーニング方式に基づく歪補償信号処理を行う。これにより、アンプモデル係数算出部２１及びフィードバック信号レプリカ生成部２２等の構成を省略することができる。この結果、共通フィードバック歪補償信号処理部１０ｂの構成を簡略化することができる。
【０１０４】
（変形例３）
変形例２では、第１の条件及び第２の条件が両方成り立つ場合に適切に動作する。第１の条件は、減算する送信信号ｘ２（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第２の送信系の送信信号成分ｘ２（ｎ）の振幅及び位相がほぼ完全に一致することである。第２の条件は、減算する送信信号ｘ１（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第１の送信系の送信信号成分ｘ１（ｎ）の振幅及び位相がほぼ完全に一致することである。しかしながら、第１の条件及び第２の条件の少なくとも一方が成り立たない場合は適切に動作することが保証されていない。
【０１０５】
そこで、このような問題を回避するために、変形例３では、歪補償係数を更新する場合に、フィードバック信号レプリカの振幅に代えて送信信号の振幅を用いて、更新後の歪補償係数を算出する。具体的には、変形例３では、数式（９）中の行列Ｗ１を、数式（１１）に代えて、次の数式（２２）により算出し、数式（１４）中の行列Ｗ２を、数式（１６）に代えて、次の数式（２３）により算出する。
【０１０６】
【数２２】
2018113626.tif 000024
【数２３】
2018113626.tif 000025
【０１０７】
数式（２２）において、数式（１１）と異なる点は、絶対値の中身を第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）に代えて第１の送信系の送信信号ｘ１（ｎ）とする点である。
【０１０８】
数式（２３）において、数式（１６）と異なる点は、絶対値の中身を第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）に代えて第２の送信系の送信信号ｘ２（ｎ）とする点である。
【０１０９】
すなわち、変形例３では、フィードバック信号レプリカ＾ｙ１（ｎ），＾ｙ２（ｎ）の振幅情報に代えて、送信信号ｘ１（ｎ），ｘ２（ｎ）の振幅情報を用いて逆特性を算出する。図６は、変形例３における共通フィードバック歪補償信号処理部１０ｃの構成を示す図である。
【０１１０】
図６に示す共通フィードバック歪補償信号処理部１０ｃでは、歪補償係数更新部１２ｃが、第１の歪補償係数（＾Ａ）を更新する場合に、第１の送信系の送信信号ｘ１（ｎ）の振幅｜ｘ１（ｎ）｜を算出する。歪補償係数更新部１２ｃは、第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）の振幅｜＾ｙ１（ｎ）｜に代えて第１の送信系の送信信号ｘ１（ｎ）の振幅｜ｘ１（ｎ）｜を用いて、更新後の第１の歪補償係数（＾Ａ）を算出する。すなわち、歪補償係数更新部１２ｃは、数式（２２）により行列Ｗ１を算出し、係数更新前の歪補償処理後の信号ｚ１（ｎ）及び行列Ｗ１を用いて、数式（９）により、更新後の第１の歪補償係数（＾Ａ）を算出する。
【０１１１】
図６に示す共通フィードバック歪補償信号処理部１０ｃでは、歪補償係数更新部１３ｃが、第２の歪補償係数（＾Ｂ）を更新する場合に、第２の送信系の送信信号ｘ２（ｎ）の振幅｜ｘ２（ｎ）｜を算出する。歪補償係数更新部１３ｃは、第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）の振幅｜＾ｙ２（ｎ）｜に代えて第２の送信系の送信信号ｘ２（ｎ）の振幅｜ｘ２（ｎ）｜を用いて、更新後の第２の歪補償係数（＾Ｂ）を算出する。すなわち、歪補償係数更新部１３ｃは、数式（２３）により行列Ｗ２を算出し、係数更新前の歪補償処理後の信号ｚ２（ｎ）及び行列Ｗ２を用いて、数式（１４）により、更新後の第２の歪補償係数（＾Ｂ）を算出する。
【０１１２】
このように、変形例３では、フィードバック信号レプリカ＾ｙ１（ｎ），＾ｙ２（ｎ）の振幅に代えて送信信号ｘ１（ｎ），ｘ２（ｎ）の振幅を用いて歪補償係数を算出するので、送信信号ｘ１（ｎ），ｘ２（ｎ）に追従させながら歪補償係数を適切に更新できる。
【０１１３】
（変形例４）
変形例２では、第１の条件及び第２の条件が両方成り立つ場合に適切に動作する。第１の条件は、減算する送信信号ｘ２（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第２の送信系の送信信号成分ｘ２（ｎ）の振幅及び位相がほぼ完全に一致することである。第２の条件は、減算する送信信号ｘ１（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第１の送信系の送信信号成分ｘ１（ｎ）の振幅及び位相がほぼ完全に一致することである。しかしながら、第１の条件及び第２の条件の少なくとも一方が成り立たない場合は適切に動作することが保証されていない。
【０１１４】
そこで、このような問題を回避するために、変形例４では、初期値検出のために最初（送信信号の１回目のサンプリング）に送信信号の１次項のみで歪補償信号処理を行い、２回目以降のサンプリングでは変形例２と同様の処理を行う。
【０１１５】
具体的には、図７に示す共通フィードバック歪補償信号処理部１０ｄは、変形例２の歪補償係数更新部１２ｂ，１３ｂ（図５参照）に代えて、歪補償係数更新部１２ｄ，１３ｄを有する。図７に示す共通フィードバック歪補償信号処理部１０ｄは、変形例２の歪補償部１４，１５（図５参照）に代えて、歪補償部１４ｄ，１５ｄを有する。歪補償部１４ｄには、スイッチ１４１ｄ，１４３ｄを介して乗算器１４５ｄが並列に接続され、歪補償部１５ｄには、スイッチ１５１ｄ，１５３ｄを介して乗算器１５５ｄが並列に接続されている。信号分離部１１ｄは、変形例２の信号分離部１１ａ（図５参照）に対して、スイッチ１１１ｄ、スイッチ１１３ｄ、スイッチ１１５ｄ及びスイッチ１１７ｄがさらに追加されている。なお、図７は、変形例４における共通フィードバック歪補償信号処理部１０ｄの構成を示す図である。
【０１１６】
図７において、各スイッチ１４１ｄ，１４３ｄ，１５１ｄ，１５３ｄ，１１１ｄ，１１３ｄ，１１５ｄ，１１７ｄは、第１の送信系の送信信号及び第２の送信系の送信信号の１回目のサンプリングにおいて実線の状態に切り替えられている。これにより、歪補償係数更新部１２ｄ，１３ｄ及び歪補償部１４ｄ，１５ｄは、初期値を検出する処理（図８に示すＳ４１，Ｓ４２）を行う。各スイッチ１４１ｄ，１４３ｄ，１５１ｄ，１５３ｄ，１１１ｄ，１１３ｄ，１１５ｄ，１１７ｄは、第１の送信系の送信信号及び第２の送信系の送信信号の２回目以降のサンプリングにおいて破線の状態に切り替えられている。これにより、歪補償係数更新部１２ｄ，１３ｄ及び歪補償部１４ｄ，１５ｄは、変形例２と同様の処理（図８に示すＳ４３〜Ｓ４６）を行う。図８は、変形例４における歪補償係数更新部１２ｄ，１３ｄ及び歪補償部１４ｄ，１５ｄの動作を示すフローチャートである。
【０１１７】
例えば、初期値検出時の１回目のサンプリングにおいて、第１の送信系の歪補償信号ｚ１（ｎ）と第２の送信系の歪補償信号ｚ２（ｎ）とは、それぞれ次の数式（２４）で表される。
【０１１８】
【数２４】
2018113626.tif 000026
【０１１９】
数式（２４）において、ｘ１（ｎ）は、第１の送信系の送信信号であり、ｘ２（ｎ）は、第２の送信系の送信信号である。ｙ（ｎ）は、共通フィードバック信号である。このとき、ｚ１（ｎ），ｚ２（ｎ）をそれぞれ行列形式で表現すると、次の数式（２５）で表される。
【０１２０】
【数２５】
2018113626.tif 000027
【０１２１】
数式（２５）中のｚ１，ｚ２は、それぞれ次の数式（２６）で表現される。
【０１２２】
【数２６】
2018113626.tif 000028
【０１２３】
数式（２５）中のＵ１，Ｕ２は、それぞれ次の数式（２７）で表現される。
【０１２４】
【数２７】
2018113626.tif 000029
【０１２５】
数式（２５）中のｈ１，ｈ２は、それぞれ次の数式（２８）で表現される。
【０１２６】
【数２８】
2018113626.tif 000030
【０１２７】
従って、歪補償係数更新部１２ｄ，１３ｄは、初期条件として、ｚ１（ｎ）＝ｘ１（ｎ）、ｚ２（ｎ）＝ｘ２（ｎ）を設定し、次の数式（２９）に示すように係数ベクトル＾ｈ１，＾ｈ２を求める（Ｓ４１）。
【０１２８】
【数２９】
2018113626.tif 000031
【０１２９】
すなわち、歪補償係数更新部１２ｄは、１回目のサンプリングにおいて、共通フィードバック信号ｙ（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とを用いて、第１の歪補償係数の初期値（係数ベクトル＾ｈ１）を算出する。歪補償係数更新部１３ｄは、１回目のサンプリングにおいて、共通フィードバック信号ｙ（ｎ）と、第１の送信系の送信信号ｘ１（ｎ）とを用いて、第２の歪補償係数の初期値（係数ベクトル＾ｈ２）を算出する。
【０１３０】
そして、乗算器１４５ｄ，１５５ｄは、数式（２９）で求めた係数α１，α２を用いて、次の数式（３０）に示すように第１の送信系の歪補償信号＾ｚ１（ｎ）と第２の送信系の歪補償信号＾ｚ２（ｎ）をそれぞれ生成する（Ｓ４２）。
【０１３１】
【数３０】
2018113626.tif 000032
【０１３２】
すなわち、乗算器１４５ｄは、１回目のサンプリングにおいて、歪補償部１４ｄとして機能し、第１の歪補償係数の初期値（係数ベクトル＾ｈ１のうちα１）を用いて、第１の送信系の送信信号ｘ１（ｎ）の１次項のみを含む第１の送信系の歪補償信号＾ｚ１（ｎ）を生成する。乗算器１５５ｄは、１回目のサンプリングにおいて、歪補償部１５ｄとして機能し、第２の歪補償係数の初期値（係数ベクトル＾ｈ２のうちα２）を用いて、第２の送信系の送信信号ｘ２（ｎ）の１次項のみを含む第２の送信系の歪補償信号＾ｚ２（ｎ）を生成する。
【０１３３】
続いて、２回目以降のサンプリングにおいて、歪補償係数更新部１２ｄ，１３ｄは、数式（２０），数式（２１），数式（９），数式（１４）に示すように歪補償係数（＾Ａ，＾Ｂ）を更新する（Ｓ４３，Ｓ４４）。
【０１３４】
すなわち、歪補償係数更新部１２ｄは、２回目以降のサンプリングにおいて、減算器２３で生成された第１の送信系のフィードバック信号レプリカ＾ｙ１（ｎ）を用いて、第１の歪補償係数（歪補償係数ベクトル＾Ａ）を更新する。歪補償係数更新部１３ｄは、２回目以降のサンプリングにおいて、減算器２４で生成された第２の送信系のフィードバック信号レプリカ＾ｙ２（ｎ）を用いて、第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を更新する。
【０１３５】
２回目以降のサンプリングにおいて、歪補償部１４ｄ，１５ｄは、数式（１３），数式（１８）に示すように第１の送信系の歪補償信号＾ｚ１（ｎ）と第２の送信系の歪補償信号＾ｚ２（ｎ）をそれぞれ生成する（Ｓ４５）。
【０１３６】
すなわち、歪補償部１４ｄは、２回目以降のサンプリングにおいて、歪補償係数更新部１２ｄで更新された第１の歪補償係数（歪補償係数ベクトル＾Ａ）を用いて、ｘ１（ｎ）の１次項とｘ１（ｎ）の１次項以外の項とを含む第１の送信系の歪補償信号＾ｚ１（ｎ）を生成する。歪補償部１５ｄは、２回目以降のサンプリングにおいて、歪補償係数更新部１３ｄで更新された第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を用いて、ｘ２（ｎ）の１次項とｘ２（ｎ）の１次項以外の項とを含む第２の送信系の歪補償信号＾ｚ２（ｎ）を生成する。
【０１３７】
このように、変形例４では、初期値検出時に１次項のみで歪補償信号処理を行うので、減算する送信信号ｘ２（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第２の送信系の送信信号成分ｘ２（ｎ）の振幅及び位相をほぼ完全に一致させることができる。さらに、変形例４では、減算する送信信号ｘ１（ｎ）と共通フィードバック信号ｙ（ｎ）に含まれる第１の送信系の送信信号成分ｘ１（ｎ）の振幅及び位相をほぼ完全に一致させることができる。結果として、送信信号ｘ１（ｎ），ｘ２（ｎ）に追従させながら歪補償係数を適切に更新できる。
【０１３８】
（変形例５）
変形例５では、変形例３の歪補償係数更新部１２ｃ（図６参照）と変形例４の歪補償係数更新部１２ｄ（図７参照）とをハイブリッドさせて、図９に示すように、歪補償係数更新部１２ｅを構成する。さらに、変形例５では、変形例３の歪補償係数更新部１３ｃ（図６参照）と変形例４の歪補償係数更新部１３ｄ（図７参照）とをハイブリッドさせて、図９に示すように、歪補償係数更新部１３ｅを構成する。図９は、変形例５における共通フィードバック歪補償信号処理部１０ｅの構成を示す図である。
【０１３９】
歪補償係数更新部１２ｅは、共通フィードバック信号ｙ（ｎ）と、第１の送信系の送信信号ｘ１（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とを用いて、第１の歪補償係数（歪補償係数ベクトル＾Ａ）を更新する。歪補償係数更新部１３ｅは、共通フィードバック信号ｙ（ｎ）と、第１の送信系の送信信号ｘ１（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とを用いて、第２の歪補償係数（歪補償係数ベクトル＾Ｂ）を更新する。すなわち、歪補償係数更新部１２ｅ，１３ｅは、次の数式（３１）に示す計算を行い、歪補償係数ベクトル＾Ａ，＾Ｂを更新する。
【０１４０】
【数３１】
2018113626.tif 000033
【０１４１】
数式（３１）中の歪補償係数ベクトル＾Ａ，＾Ｂは、更新後の第１の歪補償係数及び更新後の第２の歪補償係数を示し、それぞれ次の数式（３２）で表現される。
【０１４２】
【数３２】
2018113626.tif 000034
【０１４３】
数式（３１）中の行列Ｗ１，Ｗ２は、共通フィードバック信号ｙ（ｎ）と、第１の送信系の送信信号ｘ１（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とにより求められる。行列Ｗ１，Ｗ２は、共通フィードバック信号と各送信系の送信信号の奇数次項とを含む行列であり、それぞれ次の数式（３３）で表現される。
【０１４４】
【数３３】
2018113626.tif 000035
【０１４５】
なお、数式（３１）において、添え字Ｈは、エルミート転置（行列を転置しその成分の複素共役をとること）を表し、添え字−１は、逆行列を表す。
【０１４６】
数式（３１）中のベクトルｚ１，ｚ２は、第１の送信系の送信信号ｘ１（ｎ）と、第２の送信系の送信信号ｘ２（ｎ）とにより求められ、係数更新前の歪補償処理後の信号を含むベクトルであり、それぞれ次の数式（３４）で表現される。
【０１４７】
【数３４】
2018113626.tif 000036
【０１４８】
歪補償係数更新部１２ｅ，１３ｅは、数式（３１）で更新された歪補償係数ベクトル＾Ａ，＾Ｂをそれぞれ歪補償部１４，１５へコピーする。
【０１４９】
歪補償部１４は、第１の送信系の送信信号ｘ１（ｎ）と、第１の歪補償係数（＾Ａ）とを用いて、第１の送信系の歪補償信号＾ｚ１（ｎ）を生成する。歪補償部１５は、第２の送信系の送信信号ｘ２（ｎ）と、第２の歪補償係数（＾Ｂ）とを用いて、第２の送信系の歪補償信号＾ｚ２（ｎ）を生成する。すなわち、歪補償部１４，１５は、次の数式（３５）に示す計算を行い、歪補償信号＾ｚ１（ｎ），＾ｚ２（ｎ）を生成する。
【０１５０】
【数３５】
2018113626.tif 000037
【０１５１】
このとき、歪補償部１４，１５は、歪補償係数ベクトルの成分のうち半分を歪補償信号の生成に用いる。つまり、歪補償部１４，１５は、歪補償係数ベクトル＾Ａのうち＾Ａ１の部分を用い、歪補償係数ベクトル＾Ｂのうち＾Ｂ２の部分を用いる。
【０１５２】
このように、変形例５では、フィードバック信号レプリカ＾ｙ１（ｎ），＾ｙ２（ｎ）の振幅に代えて送信信号ｘ１（ｎ），ｘ２（ｎ）の振幅を用いて歪補償係数を算出するので、送信信号ｘ１（ｎ），ｘ２（ｎ）に追従させながら歪補償係数を適切に更新できる。また、変形例５では、歪補償処理には用いられない項及び係数を含むべき級数に基づいて歪補償係数を更新することにより、実質的に共通フィードバック信号から他方の送信信号を、振幅及び位相を一致させたうえでキャンセルすることができる。なお、歪補償処理には用いられない項及び係数とは、歪補償係数ベクトル＾Ａのうちｂ１の部分と、歪補償係数ベクトル＾Ｂのうちａ２の部分とを指す。
【０１５３】
（変形例６）
なお、歪補償装置１ａの共通フィードバック部９０ａは、図１０に示すように、信号合成器を用いずに共通フィードバック信号を生成する構成であってもよい。図１０は、変形例６に係る歪補償装置１ａの構成を示す図である。
【０１５４】
図１０に示すように、共通フィードバック部９０ａは、共通フィードバック部９０（図１参照）における信号合成器９１に代えて、アンテナ９１ａ及び低雑音パワーアンプ９１ｂを有する。また、図１０に示す送信回路３０，５０においてカップラ３４，５４（図１参照）がそれぞれ省略されている。
【０１５５】
アンテナ９１ａは、パワーアンプ３３の出力が第１の送信系のアンテナＡ１から無線送信されて得られる電波と、パワーアンプ５３の出力が第２の送信系のアンテナＡ２から無線送信されて得られる電波とを受信して共通フィードバック信号ｙ（ｎ）を生成する。これにより、共通フィードバック信号ｙ（ｎ）には、第１の送信系のアンプ出力信号と、第２の送信系のアンプ出力信号とが重ねられる。パワーアンプ３３の出力が第１の送信系のアンテナＡ１から無線送信されて得られる電波は、上記の「第１の送信系のアンプ出力信号」に相当する。パワーアンプ５３の出力が第２の送信系のアンテナＡ２から無線送信されて得られる電波は、上記の「第２の送信系のアンプ出力信号」に相当する。アンテナ９１ａで生成された共通フィードバック信号ｙ（ｎ）は、低雑音パワーアンプ９１ｂで電力増幅され、ダウンコンバータ９２へ供給される。
【０１５６】
このように、変形例６によれば、信号合成器を用いずに共通フィードバック信号を生成する。これにより、送信回路３０，５０からカップラ３４，５４を省略することができ、カップラ３４，５４における各送信系の送信信号の減衰を回避することができる。
【符号の説明】
【０１５７】
１、１ａ  歪補償装置
１０，１０ａ，１０ｂ，１０ｃ，１０ｄ，１０ｅ  共通フィードバック歪補償信号処理部
１１，１１ａ，１１ｄ  信号分離部
１２，１２ａ，１２ｂ，１２ｃ，１２ｄ，１２ｅ  歪補償係数更新部
１３，１３ａ，１３ｂ，１３ｃ，１３ｄ，１３ｅ  歪補償係数更新部
１４，１４ｄ  歪補償部
１５，１５ｄ  歪補償部
２１  アンプモデル係数算出部
２２  フィードバック信号レプリカ生成部
２３  減算器
２４  減算器
３０  送信回路
３３  パワーアンプ
５０  送信回路
５３  パワーアンプ
９０，９０ａ  共通フィードバック部
９１  信号合成器
９１ａ  アンテナ
【図１】
2018113626.tif 000038
【図２】
2018113626.tif 000039
【図３】
2018113626.tif 000040
【図４】
2018113626.tif 000041
【図５】
2018113626.tif 000042
【図６】
2018113626.tif 000043
【図７】
2018113626.tif 000044
【図８】
2018113626.tif 000045
【図９】
2018113626.tif 000046
【図１０】
2018113626.tif 000047