JP2001268150A - リニアライザ - Google Patents
リニアライザInfo
- Publication number
- JP2001268150A JP2001268150A JP2000077353A JP2000077353A JP2001268150A JP 2001268150 A JP2001268150 A JP 2001268150A JP 2000077353 A JP2000077353 A JP 2000077353A JP 2000077353 A JP2000077353 A JP 2000077353A JP 2001268150 A JP2001268150 A JP 2001268150A
- Authority
- JP
- Japan
- Prior art keywords
- data
- input
- value
- linearizer
- input signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Amplifiers (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
(57)【要約】
【課題】局部復調器を有するアダプティブプリディスト
ーション方式の非線形補償回路において、回路規模が小
さく、かつ高速動作可能で高精度に補償が可能な回路構
成を提供する。 【解決手段】入力信号レベルの変動範囲を複数個の信号
レベル代表点で表わして、代表点でのみ補償係数を計算
し、補間によって他のレベルでの補償係数を求めること
で、回路規模の縮小化、動作の高速化を図る。また、補
償係数の計算には、累積乗算器を用いたフィードバック
ループにより、漸化的に計算を行うことで、高精度な補
償が可能となる。
ーション方式の非線形補償回路において、回路規模が小
さく、かつ高速動作可能で高精度に補償が可能な回路構
成を提供する。 【解決手段】入力信号レベルの変動範囲を複数個の信号
レベル代表点で表わして、代表点でのみ補償係数を計算
し、補間によって他のレベルでの補償係数を求めること
で、回路規模の縮小化、動作の高速化を図る。また、補
償係数の計算には、累積乗算器を用いたフィードバック
ループにより、漸化的に計算を行うことで、高精度な補
償が可能となる。
Description
【0001】
【発明の属する技術分野】本発明は非線形補償回路に係
わり、特にディジタル無線機に用いる電力増幅器の非線
形特性を補償する回路(リニアライザ)に関するもので
ある。
わり、特にディジタル無線機に用いる電力増幅器の非線
形特性を補償する回路(リニアライザ)に関するもので
ある。
【0002】
【従来の技術】移動無線等のディジタル無線の急速な普
及に伴って送信機の電力増幅器の大電力化、高効率化の
要求が高まっており、デバイス、回路両面で開発が盛ん
に行われている。ディジタル無線機では、電力増幅器の
非線形な入出力振幅特性に起因する高次歪みによる隣接
チャンネル妨害を防ぐため、電力増幅器の入出力振幅特
性の線形性が厳しく要求されている。しかし、良好な線
形性を保ちながら、大電力化、高効率化することは困難
であることから、リニアライザによる非線形補償技術の
適用が重要となっている。電力増幅器の非線形特性を補
償するリニアライザには大きく分けて、フィードフォワ
ード方式、フィードバック方式、プリディストーション
方式がある。
及に伴って送信機の電力増幅器の大電力化、高効率化の
要求が高まっており、デバイス、回路両面で開発が盛ん
に行われている。ディジタル無線機では、電力増幅器の
非線形な入出力振幅特性に起因する高次歪みによる隣接
チャンネル妨害を防ぐため、電力増幅器の入出力振幅特
性の線形性が厳しく要求されている。しかし、良好な線
形性を保ちながら、大電力化、高効率化することは困難
であることから、リニアライザによる非線形補償技術の
適用が重要となっている。電力増幅器の非線形特性を補
償するリニアライザには大きく分けて、フィードフォワ
ード方式、フィードバック方式、プリディストーション
方式がある。
【0003】フィードフォワード方式は、電力増幅器の
出力の一部を取り出し、別途発生した無歪みの信号成分
を減算して歪み成分を作成し、これを電力増幅器の出力
から差し引いて歪み補償する方式で、高周波(Radio Fr
equency、以降RFと称する)帯での回路技術で構成され
る。フィードバック方式は、電力増幅器に負帰還をかけ
て歪を補償する方式で、RF帯、IF(Intermediate Frequ
ency、中間周波数)帯、BB(Base Band:基底帯域)の
様々な箇所で負帰還をかける方式である。プリディスト
ーション方式は、電力増幅器で発生する歪み特性と逆の
非線形特性を予め入力信号に与えてから電力増幅器に入
力する方法で、電力増幅器の前段で信号処理する方法で
ある。プリディストーション方式もRF帯、IF帯、BB帯で
処理が可能である。
出力の一部を取り出し、別途発生した無歪みの信号成分
を減算して歪み成分を作成し、これを電力増幅器の出力
から差し引いて歪み補償する方式で、高周波(Radio Fr
equency、以降RFと称する)帯での回路技術で構成され
る。フィードバック方式は、電力増幅器に負帰還をかけ
て歪を補償する方式で、RF帯、IF(Intermediate Frequ
ency、中間周波数)帯、BB(Base Band:基底帯域)の
様々な箇所で負帰還をかける方式である。プリディスト
ーション方式は、電力増幅器で発生する歪み特性と逆の
非線形特性を予め入力信号に与えてから電力増幅器に入
力する方法で、電力増幅器の前段で信号処理する方法で
ある。プリディストーション方式もRF帯、IF帯、BB帯で
処理が可能である。
【0004】上記3方式のリニアライザの中で、フィー
ドバック方式は、自動的に電力増幅器の特性変動に追随
するが、他の方式はアダプティブに特性変動を検出する
適応補償制御を行う必要がある。適応信号制御のように
複雑な信号処理は、アナログ回路技術よりもディジタル
信号処理の方が容易である。しかしながら、ディジタル
信号処理は高周波信号を扱うのが不得手であるので、デ
ィジタル信号処理によるアダプティブ制御は主にプリデ
ィストーション方式に適用されている。
ドバック方式は、自動的に電力増幅器の特性変動に追随
するが、他の方式はアダプティブに特性変動を検出する
適応補償制御を行う必要がある。適応信号制御のように
複雑な信号処理は、アナログ回路技術よりもディジタル
信号処理の方が容易である。しかしながら、ディジタル
信号処理は高周波信号を扱うのが不得手であるので、デ
ィジタル信号処理によるアダプティブ制御は主にプリデ
ィストーション方式に適用されている。
【0005】プリディストーション非線形補償の原理を
図2を用いて説明する。図2(a)は、電力増幅器の入出
力特性の一例を示す図である。また、図2(b)は、従来
のプリディストーション方式を用いた非線形補償電力増
幅器(リニアライザ)の構成を示すブロック図である。
217はI成分入力端子、218はQ成分入力端子、20は複素掛
算器、21は逆特性計算部、22は直交変調部、23は直交復
調部、24は電力増幅器(HPA)、219は出力端子である。
図2は原理説明のための概略図であるので、BB信号をRF
信号に変換する周波数変換部は省略している。
図2を用いて説明する。図2(a)は、電力増幅器の入出
力特性の一例を示す図である。また、図2(b)は、従来
のプリディストーション方式を用いた非線形補償電力増
幅器(リニアライザ)の構成を示すブロック図である。
217はI成分入力端子、218はQ成分入力端子、20は複素掛
算器、21は逆特性計算部、22は直交変調部、23は直交復
調部、24は電力増幅器(HPA)、219は出力端子である。
図2は原理説明のための概略図であるので、BB信号をRF
信号に変換する周波数変換部は省略している。
【0006】図2において、電力増幅器の非線形特性は
入力信号レベルp(p2 =xI 2 + xQ 2、ただし、xIは信号の
同相成分振幅、xQは直交成分振幅である)の関数とし
て、 A(p)ejφ(p) = A(p)cosφ(p) + jA(p)sinφ(p) ‥‥‥式(1) と与えられる。式(1)でA(p)は振幅の非線形成分を表
し、φ(p)は位相の非線形成分を表している。この非線
形特性により、HPA24の出力信号は、 y = yI+jyQ = A(p)ejφ(p)(xI + jxQ) = A(p)(xIcosφ(p) - xQsinφ(p)) + jA(p)(xIsinφ(p) + xQcosφ(p)) ‥‥‥式(2) となる。
入力信号レベルp(p2 =xI 2 + xQ 2、ただし、xIは信号の
同相成分振幅、xQは直交成分振幅である)の関数とし
て、 A(p)ejφ(p) = A(p)cosφ(p) + jA(p)sinφ(p) ‥‥‥式(1) と与えられる。式(1)でA(p)は振幅の非線形成分を表
し、φ(p)は位相の非線形成分を表している。この非線
形特性により、HPA24の出力信号は、 y = yI+jyQ = A(p)ejφ(p)(xI + jxQ) = A(p)(xIcosφ(p) - xQsinφ(p)) + jA(p)(xIsinφ(p) + xQcosφ(p)) ‥‥‥式(2) となる。
【0007】式(2)の出力信号yの振幅部(yI 2 + yQ 2)
1/2の入力信号レベルpに対する変化の様子を図2(a)に
示す。ただし出力信号レベルは小信号利得で正規化して
おく。一般的に、電力増幅器の振幅非線形特性は、入力
信号レベルpの増加に対して出力振幅が低下し、図2(a)
の45°線より小さくなる傾向を示す。そのため、出力振
幅yを得るには、入力信号レベルとしてx = yでなく、
x′のレベルの信号を入力しなければならない。x′を取
得するには非線形特性(式(1))の逆特性である次の式
(3)を求め、これに入力信号xの信号レベルpを代入する
ことで得られる。 A-1(p)e-jφ(p) = aI + jaQ = A-1(p)cosφ(p) − jA-1(p)sinφ(p) ‥‥‥式(3)
1/2の入力信号レベルpに対する変化の様子を図2(a)に
示す。ただし出力信号レベルは小信号利得で正規化して
おく。一般的に、電力増幅器の振幅非線形特性は、入力
信号レベルpの増加に対して出力振幅が低下し、図2(a)
の45°線より小さくなる傾向を示す。そのため、出力振
幅yを得るには、入力信号レベルとしてx = yでなく、
x′のレベルの信号を入力しなければならない。x′を取
得するには非線形特性(式(1))の逆特性である次の式
(3)を求め、これに入力信号xの信号レベルpを代入する
ことで得られる。 A-1(p)e-jφ(p) = aI + jaQ = A-1(p)cosφ(p) − jA-1(p)sinφ(p) ‥‥‥式(3)
【0008】図2(b)において、入力端子217を介して信
号の同相成分振幅xIを入力し、入力端子218を介して信
号の直交成分振幅yIを入力する。これらの信号振幅は複
素乗算記20と直交変調器22を通って、HPA24に与えられ
増幅される、この増幅されたHPA24の出力信号yは出力端
子219を介して出力されるが、また一方直交復調器23に
与えられる。HPA24の出力信号yを直交復調器23により復
調し、yI + jyQとして、逆特性計算部21に入力する。同
様に入力信号x = xI + jxQも入力端子217と218から逆特
性計算部21に入力する。逆特性計算部21は、これらxとy
のデータから逆特性x/yを計算する。求めた逆特性x/y
に、入力信号レベルpを代入して、式(3)の非線形補正値
A-1(p)e-jφ(p) = aI + jaQを求め、複素掛算器20に与
える。これを複素掛算器20では、入力信号xに複素乗算
し、非線形補正信号x′を求め、HPA24に与える。
号の同相成分振幅xIを入力し、入力端子218を介して信
号の直交成分振幅yIを入力する。これらの信号振幅は複
素乗算記20と直交変調器22を通って、HPA24に与えられ
増幅される、この増幅されたHPA24の出力信号yは出力端
子219を介して出力されるが、また一方直交復調器23に
与えられる。HPA24の出力信号yを直交復調器23により復
調し、yI + jyQとして、逆特性計算部21に入力する。同
様に入力信号x = xI + jxQも入力端子217と218から逆特
性計算部21に入力する。逆特性計算部21は、これらxとy
のデータから逆特性x/yを計算する。求めた逆特性x/y
に、入力信号レベルpを代入して、式(3)の非線形補正値
A-1(p)e-jφ(p) = aI + jaQを求め、複素掛算器20に与
える。これを複素掛算器20では、入力信号xに複素乗算
し、非線形補正信号x′を求め、HPA24に与える。
【0009】上述のプリディストーションの原理によれ
ば、HPA24の逆特性を直交復調器23により常に監視して
いるので、HPA24が温度変動等により特性変動を生じて
も、アダプティブに特性追従が可能であり、高精度な非
線形補償が実現できる。ディジタル信号処理によるアダ
プティブプリディストーション方式リニアライザの従来
例として、高林、折橋、松岡、森井、“ディジタル直交
変復調器を適用した送信系線形補償の検討”、電子情報
通信学会1998年ソサイエティ大会、B-5-4で報告された
補償回路がある。図3にこの従来例の概略構成を示す。
ば、HPA24の逆特性を直交復調器23により常に監視して
いるので、HPA24が温度変動等により特性変動を生じて
も、アダプティブに特性追従が可能であり、高精度な非
線形補償が実現できる。ディジタル信号処理によるアダ
プティブプリディストーション方式リニアライザの従来
例として、高林、折橋、松岡、森井、“ディジタル直交
変復調器を適用した送信系線形補償の検討”、電子情報
通信学会1998年ソサイエティ大会、B-5-4で報告された
補償回路がある。図3にこの従来例の概略構成を示す。
【0010】図3は、従来のアダプティブプリディスト
ーション方式非線形補償回路の構成を示すブロック図で
ある。201はI成分入力端子、202はQ成分入力端子、30は
参照テーブル(LUT)、31はパワー計算部(PER)、32と
33は複素掛算器、34は係数計算部、35はデータ更新部、
10は直交変調器、11は直交復調器、12はD/A変換器、13
はA/D変換器、14と17は周波数混合器、15と16は局部発
振器、18はHPA、19は方向性結合器である。図中、信号
線が2本平行に描かれているのは、複素信号(同相成分
を実数部、直交成分を虚数部とする)で表わした信号で
ある。
ーション方式非線形補償回路の構成を示すブロック図で
ある。201はI成分入力端子、202はQ成分入力端子、30は
参照テーブル(LUT)、31はパワー計算部(PER)、32と
33は複素掛算器、34は係数計算部、35はデータ更新部、
10は直交変調器、11は直交復調器、12はD/A変換器、13
はA/D変換器、14と17は周波数混合器、15と16は局部発
振器、18はHPA、19は方向性結合器である。図中、信号
線が2本平行に描かれているのは、複素信号(同相成分
を実数部、直交成分を虚数部とする)で表わした信号で
ある。
【0011】図3において、変調入力信号のI成分xIが
入力端子201を介して、またQ成分xQが入力端子202を介
して、それぞれ複素掛算器32に与えられる。複素掛算器
32では、入力した変調入力信号に乗算器33から与えられ
る非線形補償係数を乗算し、プリディストーション補償
した信号として直交変調器10に与える。直交変調器10で
は、変調を行いD/A変換器12に与える。更にD/A変換器12
では、変調された信号をアナログ信号に変換して、周波
数混合器14に与える。アナログの変調信号は周波数混合
器14により局部発振器15からの局発信号fIFと混合され
て、RF帯に周波数変換され、HPA18に与えられる。HPA18
は信号を所定のレベルに増幅して出力端子203を介して
出力する。
入力端子201を介して、またQ成分xQが入力端子202を介
して、それぞれ複素掛算器32に与えられる。複素掛算器
32では、入力した変調入力信号に乗算器33から与えられ
る非線形補償係数を乗算し、プリディストーション補償
した信号として直交変調器10に与える。直交変調器10で
は、変調を行いD/A変換器12に与える。更にD/A変換器12
では、変調された信号をアナログ信号に変換して、周波
数混合器14に与える。アナログの変調信号は周波数混合
器14により局部発振器15からの局発信号fIFと混合され
て、RF帯に周波数変換され、HPA18に与えられる。HPA18
は信号を所定のレベルに増幅して出力端子203を介して
出力する。
【0012】HPA18で増幅された信号は、方向性結合器1
9により一部が取り出され、周波数混合器17に与えられ
る。周波数混合器17では、HPA18の出力の一部である信
号を局部発振器16の局発信号fIF′と混合して、IF帯信
号としてA/D変換器13に与える。A/D変換器13では、IF帯
信号をディジタル信号に変換して直交復調器11に与え
る。直交復調器11では、入力した信号を同相成分信号と
直交成分信号に復調して、ディジタル電力増幅器の出力
データ信号としてデータ更新部35に与えられる。
9により一部が取り出され、周波数混合器17に与えられ
る。周波数混合器17では、HPA18の出力の一部である信
号を局部発振器16の局発信号fIF′と混合して、IF帯信
号としてA/D変換器13に与える。A/D変換器13では、IF帯
信号をディジタル信号に変換して直交復調器11に与え
る。直交復調器11では、入力した信号を同相成分信号と
直交成分信号に復調して、ディジタル電力増幅器の出力
データ信号としてデータ更新部35に与えられる。
【0013】データ更新部35には、別に、変調入力信号
のI成分xIとQ成分xQとが、入力端子201と202とを介し
て、それぞれ与えられている。データ更新部35では、こ
の与えられた変調入力信号と出力データ信号とによっ
て、それぞれ、入力データと出力データの更新を行い、
更新したデータを係数計算部34に与える。
のI成分xIとQ成分xQとが、入力端子201と202とを介し
て、それぞれ与えられている。データ更新部35では、こ
の与えられた変調入力信号と出力データ信号とによっ
て、それぞれ、入力データと出力データの更新を行い、
更新したデータを係数計算部34に与える。
【0014】また、変調入力信号のI成分xIとQ成分xQと
は、更にパワー計算部31にも与えられる。パワー計算部
31では、与えられた変調入力信号の入力信号レベルを計
算し、これを固定分非線形補償係数を格納したROM(Rea
d Only Memory)等の記憶デバイスで構成された参照テ
ーブル30と変動分非線形補償係数を計算する係数計算部
34に与える。
は、更にパワー計算部31にも与えられる。パワー計算部
31では、与えられた変調入力信号の入力信号レベルを計
算し、これを固定分非線形補償係数を格納したROM(Rea
d Only Memory)等の記憶デバイスで構成された参照テ
ーブル30と変動分非線形補償係数を計算する係数計算部
34に与える。
【0015】係数計算部34では、パワー計算部31から与
えられた入力信号レベルと、データ更新部35から与えら
れた更新データとによって、変動分の非線形補償係数の
計算を行う。また、参照テーブル30では、パワー計算部
31から与えられた入力レベル値を、前もって作成した参
照テーブルで補償した値に変換した固定分非線形補償係
数として、複素掛算器33に与える。
えられた入力信号レベルと、データ更新部35から与えら
れた更新データとによって、変動分の非線形補償係数の
計算を行う。また、参照テーブル30では、パワー計算部
31から与えられた入力レベル値を、前もって作成した参
照テーブルで補償した値に変換した固定分非線形補償係
数として、複素掛算器33に与える。
【0016】複素掛算器33は、入力した固定分非線形補
償係数と変動分非線形補償係数を、複素乗算して、変調
入力信号に対する非線形補償係数を求め、これを複素掛
算器32に与える。複素掛算器32は、入力端子201と202と
を介して与えられた変調入力信号のI成分xIとQ成分xQと
に、複素掛算器33から与えられた非線形補償係数によっ
て複素乗算し、プリディストーション非線形補償動作を
行う。
償係数と変動分非線形補償係数を、複素乗算して、変調
入力信号に対する非線形補償係数を求め、これを複素掛
算器32に与える。複素掛算器32は、入力端子201と202と
を介して与えられた変調入力信号のI成分xIとQ成分xQと
に、複素掛算器33から与えられた非線形補償係数によっ
て複素乗算し、プリディストーション非線形補償動作を
行う。
【0017】即ち、上述の従来例では、HPA18の非線形
特性を固定分と変動分とに分け、固定分は前もって作成
した参照テーブル30で補償し、固定分と現状の特性との
ずれを変動分として、アダプティブにプリディストーシ
ョン補償動作を行う方法である。また、係数計算部34で
変動分非線形補償係数を計算する際には、非線形特性の
変動分を2次の近似多項式により近似し、近似係数をLM
S(最小二乗誤差)アルゴリズムにより求めている。電
力増幅器の非線形特性は、電源や温度の変動及び経年変
化などの要因によって変動するため、それら全ての非線
形特性の逆特性に対応したテーブル参照を行うことがで
きる記憶デバイスには、膨大な記憶容量を必要とする。
そこで、上述の従来では、逆特性を固定分と変動分とに
分け、固定分のみ参照テーブルを用いている。
特性を固定分と変動分とに分け、固定分は前もって作成
した参照テーブル30で補償し、固定分と現状の特性との
ずれを変動分として、アダプティブにプリディストーシ
ョン補償動作を行う方法である。また、係数計算部34で
変動分非線形補償係数を計算する際には、非線形特性の
変動分を2次の近似多項式により近似し、近似係数をLM
S(最小二乗誤差)アルゴリズムにより求めている。電
力増幅器の非線形特性は、電源や温度の変動及び経年変
化などの要因によって変動するため、それら全ての非線
形特性の逆特性に対応したテーブル参照を行うことがで
きる記憶デバイスには、膨大な記憶容量を必要とする。
そこで、上述の従来では、逆特性を固定分と変動分とに
分け、固定分のみ参照テーブルを用いている。
【0018】
【発明が解決しようとする課題】前述の従来技術には、
電力増幅器の非線形特性が入力信号レベルの関数である
ことから、非線形特性の逆特性を精度よく実現するため
に、入力信号レベルの分解能を上げなければならない
が、固定分の逆特性データだけであっても、参照テーブ
ル用にかなりの記憶容量を必要とする欠点があった。
電力増幅器の非線形特性が入力信号レベルの関数である
ことから、非線形特性の逆特性を精度よく実現するため
に、入力信号レベルの分解能を上げなければならない
が、固定分の逆特性データだけであっても、参照テーブ
ル用にかなりの記憶容量を必要とする欠点があった。
【0019】また非線形特性の変動分については、前述
の従来例では近似式を用い、この近似式の係数をLMSア
ルゴリズムにより求めている。しかし、このアルゴリズ
ムは収束が遅いので、データ更新を何度も(数百回程
度)繰り返さなければならず処理速度が遅い欠点があっ
た。また、ディジタル回路で実現される変調器では、信
号が標本化されているので、入力信号レベルの変動する
全ての範囲にわたって細かくデータを取得するために
は、長い取得時間を必要とする。従って変動分の逆特性
データの更新処理が遅い欠点があった。
の従来例では近似式を用い、この近似式の係数をLMSア
ルゴリズムにより求めている。しかし、このアルゴリズ
ムは収束が遅いので、データ更新を何度も(数百回程
度)繰り返さなければならず処理速度が遅い欠点があっ
た。また、ディジタル回路で実現される変調器では、信
号が標本化されているので、入力信号レベルの変動する
全ての範囲にわたって細かくデータを取得するために
は、長い取得時間を必要とする。従って変動分の逆特性
データの更新処理が遅い欠点があった。
【0020】更に、変動分の逆特性を求める場合には、
電力増幅器の非線形特性を多項式で近似して求めている
が、収束を速くするためには、近似多項式の次数を低く
しなければならず、このため、非線形特性の近似精度が
下がってしまう欠点があった。
電力増幅器の非線形特性を多項式で近似して求めている
が、収束を速くするためには、近似多項式の次数を低く
しなければならず、このため、非線形特性の近似精度が
下がってしまう欠点があった。
【0021】また、固定分の逆特性は予め電力増幅器の
非線形特性を測定して、参照テーブルに書き込んでおか
なければならない。これを行うには、前述のリニアライ
ザとは別の測定回路が必要であり、個々の電力増幅器毎
にこの参照テーブルを作成しなければならない欠点があ
った。
非線形特性を測定して、参照テーブルに書き込んでおか
なければならない。これを行うには、前述のリニアライ
ザとは別の測定回路が必要であり、個々の電力増幅器毎
にこの参照テーブルを作成しなければならない欠点があ
った。
【0022】また、上記変動分の逆特性を取得している
時間は無線機を動作させることができないので、無線機
の通常動作を行う前に準備動作が必要となる欠点があっ
た。
時間は無線機を動作させることができないので、無線機
の通常動作を行う前に準備動作が必要となる欠点があっ
た。
【0023】以上述べたように、従来技術には、以下に
示す欠点が有る。即ち、 (1)大容量のメモリ及びLMSアルゴリズムを用いるた
め回路規模が大きい。 (2)変動分の逆特性計算の収束が遅い。 (3)非線形補償精度が低い。 (4)製造時の初期測定(固定分データの取得)、無線
機の動作モード切替が必要である。
示す欠点が有る。即ち、 (1)大容量のメモリ及びLMSアルゴリズムを用いるた
め回路規模が大きい。 (2)変動分の逆特性計算の収束が遅い。 (3)非線形補償精度が低い。 (4)製造時の初期測定(固定分データの取得)、無線
機の動作モード切替が必要である。
【0024】本発明の目的は、上記のような欠点を除去
し、非線形補償精度の高いリニアライザを提供すること
にある。また、本発明の第2の目的は、上記のような欠
点を除去し、回路規模が小さく、高速動作可能なリニア
ライザを提供することにある。更に、本発明の第3の目
的は、上記のような欠点を除去し、無調整でかつ保守不
要で、自動的に動作可能なリニアライザを提供すること
にある。また更に、本発明の他の目的は、ディジタル信
号処理では一般的に得意としない割算や高精度の関数計
算を用いずに、ディジタル回路で実現しやすい、リニア
ライザを提供することにある。
し、非線形補償精度の高いリニアライザを提供すること
にある。また、本発明の第2の目的は、上記のような欠
点を除去し、回路規模が小さく、高速動作可能なリニア
ライザを提供することにある。更に、本発明の第3の目
的は、上記のような欠点を除去し、無調整でかつ保守不
要で、自動的に動作可能なリニアライザを提供すること
にある。また更に、本発明の他の目的は、ディジタル信
号処理では一般的に得意としない割算や高精度の関数計
算を用いずに、ディジタル回路で実現しやすい、リニア
ライザを提供することにある。
【0025】
【課題を解決するための手段】上記の目的を達成するた
め、本発明のリニアライザは、逆特性の計算に近似式を
用いず、取得した入出力データを直接演算することによ
り、高精度な逆特性を求めるようにしたものである。
め、本発明のリニアライザは、逆特性の計算に近似式を
用いず、取得した入出力データを直接演算することによ
り、高精度な逆特性を求めるようにしたものである。
【0026】また、上記第2の目的を達成するため、本
発明の非線形補償回路は、 全ての信号レベルにわたっ
て逆特性の計算を行わず、予め選定した信号レベル代表
点でのみ逆特性計算を行い、他のレベルについては、レ
ベル代表点での逆特性を補間して計算を行う。また、デ
ータを取得する際にも取得した信号レベルから補間し
て、レベル代表点での逆特性を求める。従って、全信号
レベルに亘る逆特性計算は不要で、レベル代表点でのデ
ータのみ計算、記憶しておけば良いので、回路規模の縮
減と、動作の高速化を図ることができる。
発明の非線形補償回路は、 全ての信号レベルにわたっ
て逆特性の計算を行わず、予め選定した信号レベル代表
点でのみ逆特性計算を行い、他のレベルについては、レ
ベル代表点での逆特性を補間して計算を行う。また、デ
ータを取得する際にも取得した信号レベルから補間し
て、レベル代表点での逆特性を求める。従って、全信号
レベルに亘る逆特性計算は不要で、レベル代表点でのデ
ータのみ計算、記憶しておけば良いので、回路規模の縮
減と、動作の高速化を図ることができる。
【0027】更に、上記第3の目的を達成するため、本
発明の非線形補償回路は、逆特性データの計算におい
て、振幅非線形特性に対しては乗算器と遅延レジスタと
から構成した累積乗算器を用い、位相非線形特性に対し
ては乗算器と加算器と遅延レジスタとから構成した累積
加算器を用いる。累積乗算器では入出力信号レベル比が
入力となり、積の形で累積されていく。また位相の累積
加算器では入出力信号位相差の余弦、正弦値が入力とな
り、位相差の総和の余弦、正弦値が遅延レジスタに累積
されていく。逆特性補償用複素掛算器と、非線形電力増
幅器とで帰還ループを構成すると、累積乗算器の遅延レ
ジスタのデータ値は、HPAの入出力信号レベル比の逆数
に収束して行き、収束した時点では入出力信号レベル比
(すなわち、累積乗算器の入力値)は1に収束する(累
積加算器の入力値は0に収束)。非線形特性が変動する
と、累積乗算器の入力値が1でない(累積加算器の入力
値は0でない)値となるので、ループは変動した値に再
び収束して行くように動作する。こうして、自動的に非
線形補償動作が可能な補償回路が実現できる。
発明の非線形補償回路は、逆特性データの計算におい
て、振幅非線形特性に対しては乗算器と遅延レジスタと
から構成した累積乗算器を用い、位相非線形特性に対し
ては乗算器と加算器と遅延レジスタとから構成した累積
加算器を用いる。累積乗算器では入出力信号レベル比が
入力となり、積の形で累積されていく。また位相の累積
加算器では入出力信号位相差の余弦、正弦値が入力とな
り、位相差の総和の余弦、正弦値が遅延レジスタに累積
されていく。逆特性補償用複素掛算器と、非線形電力増
幅器とで帰還ループを構成すると、累積乗算器の遅延レ
ジスタのデータ値は、HPAの入出力信号レベル比の逆数
に収束して行き、収束した時点では入出力信号レベル比
(すなわち、累積乗算器の入力値)は1に収束する(累
積加算器の入力値は0に収束)。非線形特性が変動する
と、累積乗算器の入力値が1でない(累積加算器の入力
値は0でない)値となるので、ループは変動した値に再
び収束して行くように動作する。こうして、自動的に非
線形補償動作が可能な補償回路が実現できる。
【0028】更に、本発明では、割算演算を用いずに非
線形の逆特性を計算するために、ニュートン法による漸
化式で逆数平方根や逆数を計算する。これにより、割算
命令のないディジタル信号処理プロセッサでも、非線形
補償回路が容易に実現できる。
線形の逆特性を計算するために、ニュートン法による漸
化式で逆数平方根や逆数を計算する。これにより、割算
命令のないディジタル信号処理プロセッサでも、非線形
補償回路が容易に実現できる。
【0029】
【発明の実施の形態】本発明による非線形補償回路の実
施の形態を説明する前に、本発明に関わる基本的事項を
説明する。課題を解決する手段の項で述べたように、電
力増幅器の非線形特性は入力信号レベルpの関数である
が、Pの変動範囲(0〜最大レベルpmax)の全てにわたっ
て計算を行うことは、演算処理量が多すぎて実現困難で
ある。そこで、pの変動範囲を、例えばN個の等間隔(Δ
p = pmax /N、Nは整数値)に分割し、N個のレベル代表
値pi(pi = iΔp 、i = 0〜N)についてのみ逆特性デー
タを計算する。もちろん、レベル分割は一定間隔に限る
ことはなく、任意の予め定められた間隔で分割してもよ
い。
施の形態を説明する前に、本発明に関わる基本的事項を
説明する。課題を解決する手段の項で述べたように、電
力増幅器の非線形特性は入力信号レベルpの関数である
が、Pの変動範囲(0〜最大レベルpmax)の全てにわたっ
て計算を行うことは、演算処理量が多すぎて実現困難で
ある。そこで、pの変動範囲を、例えばN個の等間隔(Δ
p = pmax /N、Nは整数値)に分割し、N個のレベル代表
値pi(pi = iΔp 、i = 0〜N)についてのみ逆特性デー
タを計算する。もちろん、レベル分割は一定間隔に限る
ことはなく、任意の予め定められた間隔で分割してもよ
い。
【0030】一方、取得されるデータ、およびプリディ
ストーションに用いるデータは、pの変動範囲の任意の
値を取り得るので、補間により計算を行う。データ補間
はラグランジのM次補間多項式 f(x) = Σj=0 Mfj[Πk≠j(x - pk)/(pj - pk)] = Σj=0 MfjCj(x) ‥‥‥式(4) によって行う。式(4)、M + 1個の点(p0,f0),(p1,f1),
‥‥‥,(pM,fM)を通るM次多項式である。
ストーションに用いるデータは、pの変動範囲の任意の
値を取り得るので、補間により計算を行う。データ補間
はラグランジのM次補間多項式 f(x) = Σj=0 Mfj[Πk≠j(x - pk)/(pj - pk)] = Σj=0 MfjCj(x) ‥‥‥式(4) によって行う。式(4)、M + 1個の点(p0,f0),(p1,f1),
‥‥‥,(pM,fM)を通るM次多項式である。
【0031】x = pでの値f(p)を精度良く求めるため
に、M + 1個の点p0,p1,‥‥‥,pMを、N + 1個のデー
タp′i( = iΔp、i = 0〜N、M<N)の中から、pが p0
〜pM の中央に位置するように、式(5)または式(6)とな
るように選ぶ。 p0 = p′i―M/2,‥‥‥,pM/2 = p′i,‥‥‥,pM = p′i+M/2 ‥‥‥式(5) ただし、p′i - Δp/2 < p < p′i + Δp/2 (Mは偶数) p0 = p′i―(M-1)/2,‥‥‥,p(M-1)/2 = p′i,‥‥‥,pM = p′i+(M+1)/2 ‥‥‥式(6) ただし、p′i < p < p′i+1 (Mは奇数)補間の
精度は多項式の次数を上げるほど、またレベル代表点の
刻みを小さくするほど良くなる。また、当然のことであ
るが、補間ができるためには、fはxのみの関数で、xに
対して連続な関数でなければならない。
に、M + 1個の点p0,p1,‥‥‥,pMを、N + 1個のデー
タp′i( = iΔp、i = 0〜N、M<N)の中から、pが p0
〜pM の中央に位置するように、式(5)または式(6)とな
るように選ぶ。 p0 = p′i―M/2,‥‥‥,pM/2 = p′i,‥‥‥,pM = p′i+M/2 ‥‥‥式(5) ただし、p′i - Δp/2 < p < p′i + Δp/2 (Mは偶数) p0 = p′i―(M-1)/2,‥‥‥,p(M-1)/2 = p′i,‥‥‥,pM = p′i+(M+1)/2 ‥‥‥式(6) ただし、p′i < p < p′i+1 (Mは奇数)補間の
精度は多項式の次数を上げるほど、またレベル代表点の
刻みを小さくするほど良くなる。また、当然のことであ
るが、補間ができるためには、fはxのみの関数で、xに
対して連続な関数でなければならない。
【0032】次に本発明の電力増幅器の非線形の逆特性
計算法について説明する。従来技術の説明で述べたよう
に、非線形性は、その振幅および位相が、入力信号レベ
ルpの関数として式(1)に示すように表わされる。一方、
取得できるデータは入出力信号の同相成分と直交成分
(xI、xQ、yI、yQ)であり、非線形補償を行う逆特
性データも同相成分と直交成分の形(aI、aQ )であ
る。従って、振幅位相成分と同相直交成分の変換を三角
関数テーブルによる変換でなく演算によって行なうこと
により、精度が低下することを防ぐことができる。
計算法について説明する。従来技術の説明で述べたよう
に、非線形性は、その振幅および位相が、入力信号レベ
ルpの関数として式(1)に示すように表わされる。一方、
取得できるデータは入出力信号の同相成分と直交成分
(xI、xQ、yI、yQ)であり、非線形補償を行う逆特
性データも同相成分と直交成分の形(aI、aQ )であ
る。従って、振幅位相成分と同相直交成分の変換を三角
関数テーブルによる変換でなく演算によって行なうこと
により、精度が低下することを防ぐことができる。
【0033】即ち、電力増幅器の非線形特性A(p)e-jφ
(p)は、電力増幅器の入力信号(x =xI + jxQ )と出力
信号(y = yI + jyQ)の比として与えられる。従って、
逆特性A-1(p)e-jφ(p)は式(7)となる。 A-1(p)e-jφ(p) = (xI + jxQ)/(yI + jyQ) =((xIyI + xQyQ) - j(xIyQ - xQyI))/(yI 2 + yQ 2) ‥‥‥式(7)
(p)は、電力増幅器の入力信号(x =xI + jxQ )と出力
信号(y = yI + jyQ)の比として与えられる。従って、
逆特性A-1(p)e-jφ(p)は式(7)となる。 A-1(p)e-jφ(p) = (xI + jxQ)/(yI + jyQ) =((xIyI + xQyQ) - j(xIyQ - xQyI))/(yI 2 + yQ 2) ‥‥‥式(7)
【0034】式(7)と式(3)とを比較して、逆特性データ
として式(8)を得る。
として式(8)を得る。
【0035】更に、x2 = xI 2 + xQ 2、 y2=yI 2 + yQ 2、
xy = ((xI 2 + xQ 2)(yI 2+yQ 2))1/2と置くと、A-1(p) =
xy/y2となるので、次の式(9)となる。 aI = xyCOSφ(P)/y2 AQ = −xysinφ(p)/y2 ‥‥‥式(9)
xy = ((xI 2 + xQ 2)(yI 2+yQ 2))1/2と置くと、A-1(p) =
xy/y2となるので、次の式(9)となる。 aI = xyCOSφ(P)/y2 AQ = −xysinφ(p)/y2 ‥‥‥式(9)
【0036】式(8)の変数( xI、xQ、yI、yQ )はレベ
ルpだけの関数ではなく、位相にも依存するので、振
幅、位相が様々な値を取り得るデータを補間してpの関
数として求めることはできない。そこで、補間に用いる
変数としては、変数( x2、y2、xy cosφ、−xy sinφ
)を用いて、式(9)によって逆特性を計算する。
ルpだけの関数ではなく、位相にも依存するので、振
幅、位相が様々な値を取り得るデータを補間してpの関
数として求めることはできない。そこで、補間に用いる
変数としては、変数( x2、y2、xy cosφ、−xy sinφ
)を用いて、式(9)によって逆特性を計算する。
【0037】式(8)は1回のデータ取得で計算できる
が、割算演算を備えた高速DSP(ディジタル信号処理プ
ロセッサ)が必要である。また、実際の無線機では、電
源電圧や周囲温度の変動により、電力増幅器の非線形特
性も時間とともに変動している。このような場合の非線
形補償としては、複数回のデータ取得で、目的の特性に
収束するような補償制御で十分であり、DSPも高速動作
が要求されないので、安価で低消費電力のプロセッサで
済む。
が、割算演算を備えた高速DSP(ディジタル信号処理プ
ロセッサ)が必要である。また、実際の無線機では、電
源電圧や周囲温度の変動により、電力増幅器の非線形特
性も時間とともに変動している。このような場合の非線
形補償としては、複数回のデータ取得で、目的の特性に
収束するような補償制御で十分であり、DSPも高速動作
が要求されないので、安価で低消費電力のプロセッサで
済む。
【0038】図4によって、フィードバックループを用
いて非線形補償を行うアダプティブプリディストーショ
ン部の動作例を説明する。図4は、フィードバックルー
プを用いた非線形補償ループ(アダプティブプリディス
トーション)部の構成を示すブロック図である。204は
入力端子、40は逆特性補償器、41と45はは遅延レジス
タ、42は掛算器、43は累積乗算器、44は加算器、46は累
積加算器、47は複素数比演算回路、48は電力増幅器(HP
A)、205は出力端子である。
いて非線形補償を行うアダプティブプリディストーショ
ン部の動作例を説明する。図4は、フィードバックルー
プを用いた非線形補償ループ(アダプティブプリディス
トーション)部の構成を示すブロック図である。204は
入力端子、40は逆特性補償器、41と45はは遅延レジス
タ、42は掛算器、43は累積乗算器、44は加算器、46は累
積加算器、47は複素数比演算回路、48は電力増幅器(HP
A)、205は出力端子である。
【0039】図4において、入力信号xは入力端子204を
介して逆特性補償器40と複素数比演算回路47とに与えら
れる。逆特性補償器40では、逆特性係数An -1e-jφnが
掛けられ、An -1e-jφn・xとなってHPA48に与える。HPA
48では、入力信号An -1e-jφ n・xを非線形利得Ae-jΦ
倍し、出力信号y = Ae-jΦAn -1e-jφn・xを出力す
る。
介して逆特性補償器40と複素数比演算回路47とに与えら
れる。逆特性補償器40では、逆特性係数An -1e-jφnが
掛けられ、An -1e-jφn・xとなってHPA48に与える。HPA
48では、入力信号An -1e-jφ n・xを非線形利得Ae-jΦ
倍し、出力信号y = Ae-jΦAn -1e-jφn・xを出力す
る。
【0040】出力信号yは出力端子205から出力されると
共に、複素数比演算回路47に与えられる。従って複素数
比演算回路48には、入力信号xと出力信号yとが入力し、
複素数比演算回路48はn時点での入出力信号の逆比an -1
e-jφ = x/yを計算する。求めた振幅逆比an -1は、遅延
レジスタ41と掛算器42とで構成した累積乗算器43に与え
られ、掛算器42によって、遅延レジスタ41に記憶されて
いる1時点前の逆特性係数An-1 -1と乗算して、n時点で
の逆特性係数An -1を計算し逆特性補償器40に与える。即
ち、An -1 = an -1An-1 -1 = Πn(an -1)となる。また、
入出力位相差φは遅延レジスタ45と加算器44とで構成し
た累積加算器46に与えられ、加算器44によって、遅延レ
ジスタ45に記憶されている1時点前の位相累算値φn-1
と加算して、n時点での位相累算値φnを計算し、同じ
く逆特性補償器40に与える。即ち、-φn = -φ - φn-1
= Σn(φ)となる。
共に、複素数比演算回路47に与えられる。従って複素数
比演算回路48には、入力信号xと出力信号yとが入力し、
複素数比演算回路48はn時点での入出力信号の逆比an -1
e-jφ = x/yを計算する。求めた振幅逆比an -1は、遅延
レジスタ41と掛算器42とで構成した累積乗算器43に与え
られ、掛算器42によって、遅延レジスタ41に記憶されて
いる1時点前の逆特性係数An-1 -1と乗算して、n時点で
の逆特性係数An -1を計算し逆特性補償器40に与える。即
ち、An -1 = an -1An-1 -1 = Πn(an -1)となる。また、
入出力位相差φは遅延レジスタ45と加算器44とで構成し
た累積加算器46に与えられ、加算器44によって、遅延レ
ジスタ45に記憶されている1時点前の位相累算値φn-1
と加算して、n時点での位相累算値φnを計算し、同じ
く逆特性補償器40に与える。即ち、-φn = -φ - φn-1
= Σn(φ)となる。
【0041】図4の構成で、遅延レジスタ41の初期値A0
-1は1、遅延レジスタ45の初期値φ0は0としておく。非
線形補償ループが収束した時点では、y = AejΦAn -1e
-jφn・x = xとなるから、振幅An -1はA-1に収束し、an
-1は1に収束する。また、位相φnはHPA48の位相特性Φ
に収束し、φは0に収束する。An -1 =A-1=((xI 2 + x
Q 2)/(yI 2 + yQ 2))1/2は、瞬時計算で求めたx/yに
等しいので、これを式(8)に代入すると、式(10)のaIとa
Qとなって、非線形補償ループにより、HPA48の逆特性が
求められる。 aI = (xIyI + xQyQ)An -1/(xy) aQ = (xQyI - xIyQ)An -1/(xy) ただし、xy = ((xI 2 + xQ 2)(yI 2 + yQ 2))1/2 ‥‥‥式(10)
-1は1、遅延レジスタ45の初期値φ0は0としておく。非
線形補償ループが収束した時点では、y = AejΦAn -1e
-jφn・x = xとなるから、振幅An -1はA-1に収束し、an
-1は1に収束する。また、位相φnはHPA48の位相特性Φ
に収束し、φは0に収束する。An -1 =A-1=((xI 2 + x
Q 2)/(yI 2 + yQ 2))1/2は、瞬時計算で求めたx/yに
等しいので、これを式(8)に代入すると、式(10)のaIとa
Qとなって、非線形補償ループにより、HPA48の逆特性が
求められる。 aI = (xIyI + xQyQ)An -1/(xy) aQ = (xQyI - xIyQ)An -1/(xy) ただし、xy = ((xI 2 + xQ 2)(yI 2 + yQ 2))1/2 ‥‥‥式(10)
【0042】累積乗算器43と累積加算器46とを用いた非
線形補償ループでは、各時点での逆比a-1e-jφは正確
にx/yと等しくなくても良い。そこで、逆特性を計算す
るのに必要となる逆数平方根や、逆数などの演算をニュ
ートン法を用いた漸化式計算で行うことができる。ニュ
ートン法は関数f(x)を近似値xnでの接線で近似し、x
軸との交点を新しい近似値xn+1として漸近的にf(x)の
解を求める解法である。xnでの接線の傾きはyn′= f′
(xn) (xnでの微係数)で与えられ、かつ、yn′= yn
/(xn - xn+ 1)が成り立つので、漸化式は、式(11)と
なる。 xn+1 = xn - yn/yn′ ‥‥‥式(11)
線形補償ループでは、各時点での逆比a-1e-jφは正確
にx/yと等しくなくても良い。そこで、逆特性を計算す
るのに必要となる逆数平方根や、逆数などの演算をニュ
ートン法を用いた漸化式計算で行うことができる。ニュ
ートン法は関数f(x)を近似値xnでの接線で近似し、x
軸との交点を新しい近似値xn+1として漸近的にf(x)の
解を求める解法である。xnでの接線の傾きはyn′= f′
(xn) (xnでの微係数)で与えられ、かつ、yn′= yn
/(xn - xn+ 1)が成り立つので、漸化式は、式(11)と
なる。 xn+1 = xn - yn/yn′ ‥‥‥式(11)
【0043】ニュートン法でf(x) = aを解くには、変
形してy = f(x) - a = 0を解く。これを式(11)に代入
すると、式(12)となる。 xn+1 = xn -(f(xn)−a)/f′(xn) ‥‥‥式(12) ニュートン法による漸化式を用いると、割算演算命令の
無いDSPでも、逆数を求めたり、また、平方根演算をテ
ーブル参照などを用いずに精度良く求めることができ
る。
形してy = f(x) - a = 0を解く。これを式(11)に代入
すると、式(12)となる。 xn+1 = xn -(f(xn)−a)/f′(xn) ‥‥‥式(12) ニュートン法による漸化式を用いると、割算演算命令の
無いDSPでも、逆数を求めたり、また、平方根演算をテ
ーブル参照などを用いずに精度良く求めることができ
る。
【0044】以上説明した基本的事項に続き、図11を
用いて、本発明のリニアライザについて更に説明する。
図11は本発明のリニアライザの一実施例の構成を示す
ブロック図である。201は入力信号xの同相成分xI入力端
子、202は入力信号xの直交成分xQ入力端子、2はレベル
検出器(PWR)、3は複素掛算器、4′は第2の補間部、
6′は逆特性計算部、7′は第2の補間部、9′はデータ
入力部、10′は直交変調器、11′は局部直交復調器、18
は電力増幅器(HPA)、19は方向性結合器、245と246は
遅延レジスタ、247は累積乗算器、248は累積加算器、20
3は出力端子である。レベル検出器2、複素掛算器3、第
2の補間部4′、逆特性計算部6′、第2の補間部7′、
データ入力部9′、直交変調器10′、局部直交復調器1
1′、遅延レジスタ245と246、累積乗算器247、累積加算
器248は概アダプティブプリディストーション部を構成
する。
用いて、本発明のリニアライザについて更に説明する。
図11は本発明のリニアライザの一実施例の構成を示す
ブロック図である。201は入力信号xの同相成分xI入力端
子、202は入力信号xの直交成分xQ入力端子、2はレベル
検出器(PWR)、3は複素掛算器、4′は第2の補間部、
6′は逆特性計算部、7′は第2の補間部、9′はデータ
入力部、10′は直交変調器、11′は局部直交復調器、18
は電力増幅器(HPA)、19は方向性結合器、245と246は
遅延レジスタ、247は累積乗算器、248は累積加算器、20
3は出力端子である。レベル検出器2、複素掛算器3、第
2の補間部4′、逆特性計算部6′、第2の補間部7′、
データ入力部9′、直交変調器10′、局部直交復調器1
1′、遅延レジスタ245と246、累積乗算器247、累積加算
器248は概アダプティブプリディストーション部を構成
する。
【0045】図11では、HPA18の非線形特性を取得す
るため、局部直交復調器11′を備えた構成としている。
出力信号yと入力信号xの比( y/x = anejΔφ )を計
算し、これを逆特性計算部6′で累積乗算( A = an・A
n-1 )及び累積加算(φn =Δφ + φn―1 )して非線
形特性( A exp(jφn ))を求める。この逆数を入力
信号xに複素乗算して非線形補償信号x′を得る。逆特性
計算の構成を簡単化するため、HPA入力信号変化範囲をN
分割し、分割点でのみ逆特性を求め、他の点では補間、
及び逆補間によって計算を行う。以上の構成により本方
式のアダプティブプリディストーションでは、累積乗算
器247、及び累積加算器248のレジスタに初期値1、及び0
を与え漸近的に非線形特性を求めていくので、HPA18の
非線形特性を測定するための特別なテスト信号を必要と
しない。
るため、局部直交復調器11′を備えた構成としている。
出力信号yと入力信号xの比( y/x = anejΔφ )を計
算し、これを逆特性計算部6′で累積乗算( A = an・A
n-1 )及び累積加算(φn =Δφ + φn―1 )して非線
形特性( A exp(jφn ))を求める。この逆数を入力
信号xに複素乗算して非線形補償信号x′を得る。逆特性
計算の構成を簡単化するため、HPA入力信号変化範囲をN
分割し、分割点でのみ逆特性を求め、他の点では補間、
及び逆補間によって計算を行う。以上の構成により本方
式のアダプティブプリディストーションでは、累積乗算
器247、及び累積加算器248のレジスタに初期値1、及び0
を与え漸近的に非線形特性を求めていくので、HPA18の
非線形特性を測定するための特別なテスト信号を必要と
しない。
【0046】図11に示した本発明によるリニアライザ
の一実施例を、更に詳しく、図1以下を使用して説明す
る。図1は、本発明のリニアライザの構成を示すブロッ
ク図である。201はI成分入力端子、202はQ成分入力端
子、1は非線形補償部、2はレベル検出器(PWR)、3は複
素掛算器、4は第2の補間部、5は第2のメモリ(逆特性
メモリ)、6は逆特性計算部、7は第1の補間部、8は第
1のメモリ、9はデータ入力部、10は直交変調器、11は
直交復調器、12はD/A変換器、13はA/D変換器、14と17は
周波数混合器、15と16は局部発振器、18はHPA、19は方
向性結合器、203は出力端子である。また、レベル検出
器2、複素掛算器3、第2の補間部4、第2のメモリ5、逆
特性計算部6、第1の補間部7、第1のメモリ8、データ
入力部9までの部分が非線形補償部1を構成しており、全
てディジタル的に信号処理を行っている。図中、2本平
行に描いた信号線は、同相成分を実数部、直交成分を虚
数部とする複素信号である。
の一実施例を、更に詳しく、図1以下を使用して説明す
る。図1は、本発明のリニアライザの構成を示すブロッ
ク図である。201はI成分入力端子、202はQ成分入力端
子、1は非線形補償部、2はレベル検出器(PWR)、3は複
素掛算器、4は第2の補間部、5は第2のメモリ(逆特性
メモリ)、6は逆特性計算部、7は第1の補間部、8は第
1のメモリ、9はデータ入力部、10は直交変調器、11は
直交復調器、12はD/A変換器、13はA/D変換器、14と17は
周波数混合器、15と16は局部発振器、18はHPA、19は方
向性結合器、203は出力端子である。また、レベル検出
器2、複素掛算器3、第2の補間部4、第2のメモリ5、逆
特性計算部6、第1の補間部7、第1のメモリ8、データ
入力部9までの部分が非線形補償部1を構成しており、全
てディジタル的に信号処理を行っている。図中、2本平
行に描いた信号線は、同相成分を実数部、直交成分を虚
数部とする複素信号である。
【0047】図1において、I成分入力端子201を介して
I成分が、そして、Q成分入力端子202を介してQ成分が、
レベル検出器2と複素掛算器3とにそれぞれ与えられる。
これによって変調入力信号x( = xI + jxQ)が、レベル
検出器2と複素掛算器2にそれぞれ与えられることにな
る。
I成分が、そして、Q成分入力端子202を介してQ成分が、
レベル検出器2と複素掛算器3とにそれぞれ与えられる。
これによって変調入力信号x( = xI + jxQ)が、レベル
検出器2と複素掛算器2にそれぞれ与えられることにな
る。
【0048】レベル検出器2では、入力信号レベルp(p2
= xI 2 + xQ 2)を計算し、複素掛算器3では、入力信号
レベルpに応じた非線形逆特性A-1(p)e-jφ(p)を掛け
て、プリディストーション補償する。このプリディスト
ーション補償信号x′(x′ =A-1(p)e-jφ(p)・x)は
直交変調器10に与えられる。直交変調器10は、入力した
プリディストーション補償信号x′を変調してD/A変換器
12に与える。D/A変換器12では、被変調信号をアナログ
値に変換し、このアナログ変調信号は周波数混合器14に
与える。
= xI 2 + xQ 2)を計算し、複素掛算器3では、入力信号
レベルpに応じた非線形逆特性A-1(p)e-jφ(p)を掛け
て、プリディストーション補償する。このプリディスト
ーション補償信号x′(x′ =A-1(p)e-jφ(p)・x)は
直交変調器10に与えられる。直交変調器10は、入力した
プリディストーション補償信号x′を変調してD/A変換器
12に与える。D/A変換器12では、被変調信号をアナログ
値に変換し、このアナログ変調信号は周波数混合器14に
与える。
【0049】周波数混合器14では、アナログ変調信号を
局部発振器15の局発信号fIFと混合して、RF帯信号に周
波数変換し、HPA18与える。HPA18では、入力した信号を
増幅して出力する。HPA18の出力は、出力端子203を介し
て出力されると共に、方向性結合器19により一部が取り
出され、周波数混合器17に与えられる。
局部発振器15の局発信号fIFと混合して、RF帯信号に周
波数変換し、HPA18与える。HPA18では、入力した信号を
増幅して出力する。HPA18の出力は、出力端子203を介し
て出力されると共に、方向性結合器19により一部が取り
出され、周波数混合器17に与えられる。
【0050】周波数混合器17では、入力した出力信号の
一部を、局部発振器16の局発信号fI F′と混合してIF帯
信号へ周波数変換する。このIF帯信号をA/D変換器13で
ディジタル信号に変換し、変換したディジタル信号を直
交復調器11に与える。
一部を、局部発振器16の局発信号fI F′と混合してIF帯
信号へ周波数変換する。このIF帯信号をA/D変換器13で
ディジタル信号に変換し、変換したディジタル信号を直
交復調器11に与える。
【0051】直交復調器11では、入力した信号を復調し
て、出力信号y(y = yI + jyQ)を出力する。この出力
信号yの同相成分yIと直交成分yQと、入力信号レベルp
と、変調入力信号xの同相成分xI及び直交成分xQとは、
データ入力部9にそれぞれ与えられる。
て、出力信号y(y = yI + jyQ)を出力する。この出力
信号yの同相成分yIと直交成分yQと、入力信号レベルp
と、変調入力信号xの同相成分xI及び直交成分xQとは、
データ入力部9にそれぞれ与えられる。
【0052】データ入力部9では、以下の変数変換を行
う。 x2 = xI 2 + xQ 2 、 y2 = yI 2 + yQ 2 、 xy cosφ = xIyI + xQyQ 、 -xy sinφ = xQyI - xIyQ 、 (ただし xy =((xI 2 + xQ 2)(yI 2 + yQ 2))1/2 ) ‥‥‥式(13) この式(13)によって取得された出力データ( x2、y2、x
y cosφ、-xy sinφ)は第1のメモリ(データメモリ)
8に与えられる。
う。 x2 = xI 2 + xQ 2 、 y2 = yI 2 + yQ 2 、 xy cosφ = xIyI + xQyQ 、 -xy sinφ = xQyI - xIyQ 、 (ただし xy =((xI 2 + xQ 2)(yI 2 + yQ 2))1/2 ) ‥‥‥式(13) この式(13)によって取得された出力データ( x2、y2、x
y cosφ、-xy sinφ)は第1のメモリ(データメモリ)
8に与えられる。
【0053】第1のメモリ8では、これらのデータを、
第1のメモリ8の最も近い入力レベル代表値pi(i = 0〜
N、p0 = 0、pN = pmax )の記憶領域に書き込んでお
く。第1のメモリ8の全メモリ領域が書き込まれた時点
で、第1の補間部(補間部1)7により、第1のメモリ8
のデータを補間し、入力レベル代表値でのデータ値( x
2、y2、xy cosφ、-xy sinφ )を計算する。
第1のメモリ8の最も近い入力レベル代表値pi(i = 0〜
N、p0 = 0、pN = pmax )の記憶領域に書き込んでお
く。第1のメモリ8の全メモリ領域が書き込まれた時点
で、第1の補間部(補間部1)7により、第1のメモリ8
のデータを補間し、入力レベル代表値でのデータ値( x
2、y2、xy cosφ、-xy sinφ )を計算する。
【0054】更に、逆特性計算部6において、入力レベ
ル代表値piでの非線形逆特性A-1(pi)e-jφ(pi) = aI
+ jaQを、式(9)によって求め、第2のメモリ(逆特性
メモリ)5に書き込む。第2の補間部(補間部2)4によ
り第2のメモリ5に書き込まれた入力レベル代表値piで
の非線形逆特性データを補間して、レベル検出器2で検
出された入力信号レベルpに対応する非線形逆特性を計
算する。こうして求めた変調入力信号xに対する非線形
逆特性を複素掛算器3に入力して、プリディストーショ
ン非線形補償動作を行う。
ル代表値piでの非線形逆特性A-1(pi)e-jφ(pi) = aI
+ jaQを、式(9)によって求め、第2のメモリ(逆特性
メモリ)5に書き込む。第2の補間部(補間部2)4によ
り第2のメモリ5に書き込まれた入力レベル代表値piで
の非線形逆特性データを補間して、レベル検出器2で検
出された入力信号レベルpに対応する非線形逆特性を計
算する。こうして求めた変調入力信号xに対する非線形
逆特性を複素掛算器3に入力して、プリディストーショ
ン非線形補償動作を行う。
【0055】以上説明したように、図1に示した本発明
の実施例によれば、逆特性の計算に近似式を用いないで
直接演算を行うので、高精度な非線形補償を行うことが
できる。代表レベル点以外では補間によって計算するの
で、若干精度は落ちるが、代表レベル点の数を増やした
り、補間式の次数を上げることで容易に精度を上げるこ
とが可能である。また、逆特性演算を代表レベル点での
み行うことにより、回路規模の縮減と、高速動作可能な
回路が実現できる。また、参照テーブルを用いず、局部
復調器により常に電力増幅器の特性を監視、追従してい
るので、無調整、自動制御が可能である。
の実施例によれば、逆特性の計算に近似式を用いないで
直接演算を行うので、高精度な非線形補償を行うことが
できる。代表レベル点以外では補間によって計算するの
で、若干精度は落ちるが、代表レベル点の数を増やした
り、補間式の次数を上げることで容易に精度を上げるこ
とが可能である。また、逆特性演算を代表レベル点での
み行うことにより、回路規模の縮減と、高速動作可能な
回路が実現できる。また、参照テーブルを用いず、局部
復調器により常に電力増幅器の特性を監視、追従してい
るので、無調整、自動制御が可能である。
【0056】図1の本発明の実施例におけるデータ入力
部9と第1のメモリ(データメモリ)8の詳細を、図5に
よって説明する。図5は本発明のデータ入力部とデータ
メモリ部の一実施例の構成を示すブロック図である。20
6はxI成分入力端子、207はxQ成分入力端子、208はHPA出
力入力端子、50と51は遅延回路、52は直交復調器、53と
54,56と57,59,510,512と513は掛算器、55,58,51
1,514は加算器、515はメモリ制御回路、516はメモリで
ある。
部9と第1のメモリ(データメモリ)8の詳細を、図5に
よって説明する。図5は本発明のデータ入力部とデータ
メモリ部の一実施例の構成を示すブロック図である。20
6はxI成分入力端子、207はxQ成分入力端子、208はHPA出
力入力端子、50と51は遅延回路、52は直交復調器、53と
54,56と57,59,510,512と513は掛算器、55,58,51
1,514は加算器、515はメモリ制御回路、516はメモリで
ある。
【0057】図5において、電力増幅器(HPA)の出力
信号yが入力端子208を介して直行復調器52に与えられ、
また変調入力信号の同相成分xI及び直交成分xQは遅延回
路50と51にそれぞれ与えられる。直交復調器52では復調
を行い、同相成分yIと直交成分yQを出力する。同相成分
yIは、掛算器59,512に与えられ、また、掛算器56に2
乗演算するために同時に2つ与えられる。更に直交成分
yQは、掛算器510,513に与えられ、また、掛算器57に2
乗演算するために同時に2つ与えられる。
信号yが入力端子208を介して直行復調器52に与えられ、
また変調入力信号の同相成分xI及び直交成分xQは遅延回
路50と51にそれぞれ与えられる。直交復調器52では復調
を行い、同相成分yIと直交成分yQを出力する。同相成分
yIは、掛算器59,512に与えられ、また、掛算器56に2
乗演算するために同時に2つ与えられる。更に直交成分
yQは、掛算器510,513に与えられ、また、掛算器57に2
乗演算するために同時に2つ与えられる。
【0058】また変調入力信号は、遅延回路50と51とに
よってそれぞれ遅延され、直交復調器52を介してデータ
入力部に戻ってくるHPA出力信号y = yI + jyQと時間合
わせを行う。遅延された同相成分xIは、掛算器53に2乗
演算するために同時に2つ与えられ、また遅延された同
相成分xQは、掛算器54に2乗演算するために同時に2つ
与えられる。
よってそれぞれ遅延され、直交復調器52を介してデータ
入力部に戻ってくるHPA出力信号y = yI + jyQと時間合
わせを行う。遅延された同相成分xIは、掛算器53に2乗
演算するために同時に2つ与えられ、また遅延された同
相成分xQは、掛算器54に2乗演算するために同時に2つ
与えられる。
【0059】掛算器53と掛算器54とでそれぞれ自乗され
た信号は、加算器55に与えられ、加算器55で加算される
ことにより、式(13)のx2( x2 = xI 2 + xQ 2 )の計算結
果が加算器55から出力される。加算器55の出力x2は、メ
モリ制御回路515とメモリ516とに与えられる。同様に、
掛算器56と掛算器57でそれぞれ自乗された信号は、加算
器58に与えられ、加算器58で加算されることにより、式
(13)のy2( y2 = yI 2 + yQ 2 )の計算結果が加算器58か
ら出力される。加算器58の出力y2はメモリ516に与えら
れる。
た信号は、加算器55に与えられ、加算器55で加算される
ことにより、式(13)のx2( x2 = xI 2 + xQ 2 )の計算結
果が加算器55から出力される。加算器55の出力x2は、メ
モリ制御回路515とメモリ516とに与えられる。同様に、
掛算器56と掛算器57でそれぞれ自乗された信号は、加算
器58に与えられ、加算器58で加算されることにより、式
(13)のy2( y2 = yI 2 + yQ 2 )の計算結果が加算器58か
ら出力される。加算器58の出力y2はメモリ516に与えら
れる。
【0060】また、遅延回路50の出力信号xIは掛算器59
と掛算器513とも与えられ、遅延回路51の出力信号xQは
掛算器510と掛算器512とも与えられる。掛算器59の出力
と掛算器510の出力とはそれぞれ加算器511に与えられ、
加算器511で加算されることにより、式(13)のxy cosφ
( xy cosφ = xIyI + xQyQ )の計算結果cが加算器511
から出力される。加算器511の出力c(c = xy cosφ)は
メモリ516に与えられる。また同様に、掛算器513の出力
と掛算器512の出力とはそれぞれ加算器514に与えられ、
加算器514で加算(ただし、掛算器513の出力は減算)さ
れることにより、式(13)の-xy sinφ( -xy sinφ = xQ
yI - xIyQ )の計算結果sが加算器514から出力される。
加算器514の出力s(s = -xy sinφ)はメモリ516に与え
られる。
と掛算器513とも与えられ、遅延回路51の出力信号xQは
掛算器510と掛算器512とも与えられる。掛算器59の出力
と掛算器510の出力とはそれぞれ加算器511に与えられ、
加算器511で加算されることにより、式(13)のxy cosφ
( xy cosφ = xIyI + xQyQ )の計算結果cが加算器511
から出力される。加算器511の出力c(c = xy cosφ)は
メモリ516に与えられる。また同様に、掛算器513の出力
と掛算器512の出力とはそれぞれ加算器514に与えられ、
加算器514で加算(ただし、掛算器513の出力は減算)さ
れることにより、式(13)の-xy sinφ( -xy sinφ = xQ
yI - xIyQ )の計算結果sが加算器514から出力される。
加算器514の出力s(s = -xy sinφ)はメモリ516に与え
られる。
【0061】次に、メモリ制御回路515では、入力するx
2のデータから、 p =(xI 2 + xQ 2)1/2 =(x2 )1/2 を算出することによって入力信号レベルpを求める。
2のデータから、 p =(xI 2 + xQ 2)1/2 =(x2 )1/2 を算出することによって入力信号レベルpを求める。
【0062】入力信号レベルpは、0から最大入力レベル
pmax間での任意の値を取り得る。そこで、入力信号レベ
ルpの変動範囲を予め定めた間隔でN分割(Nは自然数)
して、メモリ516の記憶領域を、例えば図5(b)のよう
に、0からNのに(N+1個に)分割しておく。( pi-1 + p
i)/2 < p <( pi + pi+1 )/2 となるデータ(
x2,y2,c,s )をi番目の記憶領域に格納し、データが
格納されたことを示すフラグFを1にセットする。新たな
データが既にデータが格納されている記憶領域である場
合には、新旧データの入力信号レベル値pを比較し、レ
ベル代表値piに近い方のデータを残す。こうして、メモ
リ516の全記憶領域のフラグFがセットされるまで、以上
の操作を繰り返す。
pmax間での任意の値を取り得る。そこで、入力信号レベ
ルpの変動範囲を予め定めた間隔でN分割(Nは自然数)
して、メモリ516の記憶領域を、例えば図5(b)のよう
に、0からNのに(N+1個に)分割しておく。( pi-1 + p
i)/2 < p <( pi + pi+1 )/2 となるデータ(
x2,y2,c,s )をi番目の記憶領域に格納し、データが
格納されたことを示すフラグFを1にセットする。新たな
データが既にデータが格納されている記憶領域である場
合には、新旧データの入力信号レベル値pを比較し、レ
ベル代表値piに近い方のデータを残す。こうして、メモ
リ516の全記憶領域のフラグFがセットされるまで、以上
の操作を繰り返す。
【0063】図5の実施例に拠れば、入力信号レベルの
全変動範囲に亘って細かくデータを取得する必要が無い
ので、データ取得に要する時間を短くでき、かつメモリ
の記憶容量を小さくできる。またレベル代表値に近いデ
ータを取得することにより、逆特性計算の精度を上げる
ことが可能となる。また、メモリに書き込むデータを入
出力信号そのもの( xI,xQ,yI,yQ )でなく、(
x2,y2,c = xy cosφ,s = -xy sinφ)にデータ変換
しているので、補間する際の変数が入力信号レベルpの
連続関数となる条件を満たすことができる。
全変動範囲に亘って細かくデータを取得する必要が無い
ので、データ取得に要する時間を短くでき、かつメモリ
の記憶容量を小さくできる。またレベル代表値に近いデ
ータを取得することにより、逆特性計算の精度を上げる
ことが可能となる。また、メモリに書き込むデータを入
出力信号そのもの( xI,xQ,yI,yQ )でなく、(
x2,y2,c = xy cosφ,s = -xy sinφ)にデータ変換
しているので、補間する際の変数が入力信号レベルpの
連続関数となる条件を満たすことができる。
【0064】図6によって図1の第1の補間部の具体的
実施例を説明する。図6は、本発明の第1の補間部の一
実施例の構成を示すブロック図である。209は入力信号
レベルpの入力端子、210はx0 2の入力端子、211はx1 2の
入力端子、212はx2 2の入力端子、213はf0の入力端子、2
14はf1の入力端子、215はf2の入力端子、60〜65,619,
620は加算器、66〜68は割算器、69,610,611,615〜61
8は掛算器、612〜614はデータレジスタ、621〜623は平
方根演算回路、242は値-1の入力端子、216は補間出力f
(p)の出力端子である。
実施例を説明する。図6は、本発明の第1の補間部の一
実施例の構成を示すブロック図である。209は入力信号
レベルpの入力端子、210はx0 2の入力端子、211はx1 2の
入力端子、212はx2 2の入力端子、213はf0の入力端子、2
14はf1の入力端子、215はf2の入力端子、60〜65,619,
620は加算器、66〜68は割算器、69,610,611,615〜61
8は掛算器、612〜614はデータレジスタ、621〜623は平
方根演算回路、242は値-1の入力端子、216は補間出力f
(p)の出力端子である。
【0065】補間には式(4)に示したラグランジの補間
多項式を用いる。図6では補間多項式の次数を2次とし
た例を示している。式(4)で次数M = 2とすると、補間多
項式は次の式(14)となる。 f(p) = f0( p - p1 )( p - p2 )/(( p - p1 )( p - p2 )) + f1( p - p0 )( p - p2 )/(( p1 - p0 )( p1 - p2 )) + f2( p - p0 )( p - p1 )/(( p2 - p0 )( p2 - p1 )) (ただし、pi =(x2)i 1/2 ) ‥‥‥式(14)
多項式を用いる。図6では補間多項式の次数を2次とし
た例を示している。式(4)で次数M = 2とすると、補間多
項式は次の式(14)となる。 f(p) = f0( p - p1 )( p - p2 )/(( p - p1 )( p - p2 )) + f1( p - p0 )( p - p2 )/(( p1 - p0 )( p1 - p2 )) + f2( p - p0 )( p - p1 )/(( p2 - p0 )( p2 - p1 )) (ただし、pi =(x2)i 1/2 ) ‥‥‥式(14)
【0066】式(14)は、(p0,f0),(p1,f1),
(p2,f2)の3点を通る2次多項式である。fiは4個の
データ( x2,y2,xy cosφ,-xy sinφ)を代表して表
わしたものである。また、入力信号レベルpiは,(x2)
i 1/2によって求めることができる。補間を精度良く行う
ために、補間するレベル点pがp1に近くなるようにp0,p
1,p2を選ぶ。
(p2,f2)の3点を通る2次多項式である。fiは4個の
データ( x2,y2,xy cosφ,-xy sinφ)を代表して表
わしたものである。また、入力信号レベルpiは,(x2)
i 1/2によって求めることができる。補間を精度良く行う
ために、補間するレベル点pがp1に近くなるようにp0,p
1,p2を選ぶ。
【0067】式(14)の補間係数の計算において、次の式
(15)とおくと、 補間多項式は f(p)= -f0c0 - f1c1 - f2c2 ‥‥‥式(16) となり、補間係数は4個のデータの補間に共通に使える
ことがわかる。
(15)とおくと、 補間多項式は f(p)= -f0c0 - f1c1 - f2c2 ‥‥‥式(16) となり、補間係数は4個のデータの補間に共通に使える
ことがわかる。
【0068】式(15)と式(16)の演算を図6の構成によっ
て実行する。まず、補間を行うレベルpに対し、( pi-1
+ pi )/2< p <( pi + pi+1)/2となるiを求
め、第1のメモリのi-1、i、i+1番目の記憶領域から、x
2データを読み出し、入力端子210,211,212からそれぞ
れ、平方根演算回路621,621,623に与え、平方根演算
回路621,621,623によって、信号レベルデータp0,
p1,p2を求める。
て実行する。まず、補間を行うレベルpに対し、( pi-1
+ pi )/2< p <( pi + pi+1)/2となるiを求
め、第1のメモリのi-1、i、i+1番目の記憶領域から、x
2データを読み出し、入力端子210,211,212からそれぞ
れ、平方根演算回路621,621,623に与え、平方根演算
回路621,621,623によって、信号レベルデータp0,
p1,p2を求める。
【0069】入力信号レベルpは入力端子209を介して、
加算器60,61,62にそれぞれ与えられ、また、平方根演
算回路621の出力p0が加算器60、平方根演算回路622の出
力p1が加算器61、平方根演算回路623の出力p2が加算器6
2にそれぞれ与えられる。、この結果、加算器60,61,6
2は、それぞれ、p-p0,p-p1,p-p2を計算し、被除数値
としてそれぞれ割算器67,68,66に与える。
加算器60,61,62にそれぞれ与えられ、また、平方根演
算回路621の出力p0が加算器60、平方根演算回路622の出
力p1が加算器61、平方根演算回路623の出力p2が加算器6
2にそれぞれ与えられる。、この結果、加算器60,61,6
2は、それぞれ、p-p0,p-p1,p-p2を計算し、被除数値
としてそれぞれ割算器67,68,66に与える。
【0070】また、加算器65,63,64によってそれぞ
れ、p1-p2,p2-p0,p0-p1を計算し、それぞれ割算器6
7,68,66に与える。これらの入力データによって、割
算器67,68,66では、式(15)のd0,d1,d2を計算する。
割算器67の出力d0は掛算器610と611とに与えられ、割算
器68の出力d1は掛算器69と610とに与えられ、割算器66
の出力d2は掛算器69と611とに与えられる。
れ、p1-p2,p2-p0,p0-p1を計算し、それぞれ割算器6
7,68,66に与える。これらの入力データによって、割
算器67,68,66では、式(15)のd0,d1,d2を計算する。
割算器67の出力d0は掛算器610と611とに与えられ、割算
器68の出力d1は掛算器69と610とに与えられ、割算器66
の出力d2は掛算器69と611とに与えられる。
【0071】掛算器69によって補間係数c0を計算し、求
めた補間係数c0はデータレジスタ612に与えられ記憶さ
れる。同様に、掛算器611によって補間係数c1を計算
し、求めた補間係数c1はデータレジスタ613に与えられ
記憶され、掛算器610によって補間係数c2を計算し、求
めた補間係数c2はデータレジスタ614に与えられ記憶さ
れる。
めた補間係数c0はデータレジスタ612に与えられ記憶さ
れる。同様に、掛算器611によって補間係数c1を計算
し、求めた補間係数c1はデータレジスタ613に与えられ
記憶され、掛算器610によって補間係数c2を計算し、求
めた補間係数c2はデータレジスタ614に与えられ記憶さ
れる。
【0072】次に、第1のメモリ8のi-1,i,i+1番目の
記憶領域から、入力信号の振幅2乗値x2 i-1,x2 i,x2
i+1を読み出し、式(16)のf0,f1,f2とし、入力端子21
3,214,215を介して、掛算器615,616,617にそれぞれ
与える。掛算器615では、入力したf0に補間係数c0を乗
算し、c0f0を求め掛算器618に与える。また、掛算器616
では、入力したf1に補間係数c1を乗算し、c1f1を求め加
算器619に与える。更に掛算器617では、入力したf2に補
間係数c2を乗算し、c2f2を求め加算器620に与える。
記憶領域から、入力信号の振幅2乗値x2 i-1,x2 i,x2
i+1を読み出し、式(16)のf0,f1,f2とし、入力端子21
3,214,215を介して、掛算器615,616,617にそれぞれ
与える。掛算器615では、入力したf0に補間係数c0を乗
算し、c0f0を求め掛算器618に与える。また、掛算器616
では、入力したf1に補間係数c1を乗算し、c1f1を求め加
算器619に与える。更に掛算器617では、入力したf2に補
間係数c2を乗算し、c2f2を求め加算器620に与える。
【0073】掛算器618では、入力端子242から与えられ
た値-1を乗算し、-c0f0を加算器619に与える。そして加
算器619では、-c0f0から乗算器616から与えられたc1f1
を減算し、-c0f0-c1f1として加算器620に与える。更
に、加算器620では、-c0f0-c1f1から乗算器617から与え
られたc2f2を減算し、-c0f0-c1f1-c2f2( = f(p))と
して出力端子216を介して出力する。このように、入力
信号レベルpにおける入力信号振幅2乗値x2(p)として
出力する。同様の演算を出力信号振幅2乗値y2 、及びx
y cosφ、-xy sinφについて行う。
た値-1を乗算し、-c0f0を加算器619に与える。そして加
算器619では、-c0f0から乗算器616から与えられたc1f1
を減算し、-c0f0-c1f1として加算器620に与える。更
に、加算器620では、-c0f0-c1f1から乗算器617から与え
られたc2f2を減算し、-c0f0-c1f1-c2f2( = f(p))と
して出力端子216を介して出力する。このように、入力
信号レベルpにおける入力信号振幅2乗値x2(p)として
出力する。同様の演算を出力信号振幅2乗値y2 、及びx
y cosφ、-xy sinφについて行う。
【0074】図6に示した第1の補間の実施例では、補
間する入力信号レベル点の両側のデータを用いてデータ
補間を行うので、精度良く補間を行うことができる。ま
た補間係数は4個のデータについて共通に用いることが
できるので、演算処理量を大幅に縮減できる。
間する入力信号レベル点の両側のデータを用いてデータ
補間を行うので、精度良く補間を行うことができる。ま
た補間係数は4個のデータについて共通に用いることが
できるので、演算処理量を大幅に縮減できる。
【0075】図6の実施例では、割算器を用いて補間係
数を求めているが、一般に用いられているDSP(ディジ
タル信号処理プロセッサ)では割算器が用意されていな
い場合が多い。その場合にはテーブル参照や、プログラ
ム処理により割算を実行することになるが、演算精度が
落ち、処理速度が遅くなるといった問題点がある。そこ
で、逆数を高速に演算できる逆数演算回路の実施例を図
7によって説明する。
数を求めているが、一般に用いられているDSP(ディジ
タル信号処理プロセッサ)では割算器が用意されていな
い場合が多い。その場合にはテーブル参照や、プログラ
ム処理により割算を実行することになるが、演算精度が
落ち、処理速度が遅くなるといった問題点がある。そこ
で、逆数を高速に演算できる逆数演算回路の実施例を図
7によって説明する。
【0076】図7は、本発明の逆数演算回路の一実施例
の構成を示すブロック図である。図7(a)は逆数演算回
路、図7(b)は割算回路である。220は入力値xを入力す
る入力端子、221は値2の入力端子、70,73,75は掛算
器、71は遅延レジスタ、72は加算器、74は逆数演算回
路、222は入力値bの逆数1/xを出力する出力端子であ
る。
の構成を示すブロック図である。図7(a)は逆数演算回
路、図7(b)は割算回路である。220は入力値xを入力す
る入力端子、221は値2の入力端子、70,73,75は掛算
器、71は遅延レジスタ、72は加算器、74は逆数演算回
路、222は入力値bの逆数1/xを出力する出力端子であ
る。
【0077】図7(a)において逆数を計算するために、
前述したニュートン法による漸化式(式(12))を用い
る。入力値をbとしたときの逆数x = 1/bを求めるため
に、f(x) = 1/x - b = 0 を解く。f(x)を微分し
て、f′(x) = -1/x2とし、式(12)に代入すると、漸
化式(式(17))が得られる。 (ただし、繰り返し回数nは1ずらしてある)。
前述したニュートン法による漸化式(式(12))を用い
る。入力値をbとしたときの逆数x = 1/bを求めるため
に、f(x) = 1/x - b = 0 を解く。f(x)を微分し
て、f′(x) = -1/x2とし、式(12)に代入すると、漸
化式(式(17))が得られる。 (ただし、繰り返し回数nは1ずらしてある)。
【0078】式(17)を回路で実行するには、図7(a)に
示すように、入力値xを掛算器70によって遅延レジスタ7
1に格納したxn-1(前回の計算結果)と乗算し、乗算結果
を加算器72により (入力端子221から与えられる数値)
2から引き、加算器72の出力に再び、掛算器73によってx
n-1を乗算して、新しい計算結果xnを得る。xnは次回の
計算のために遅延レジスタ71に格納する。以上の演算を
繰り返し実行すると、xnは1/x(入力値xの逆数)に収
束する。遅延レジスタに設定する初期値は結果に近い値
にすると、収束が速く3回ぐらいの繰り返しで十分な精
度が得られる。
示すように、入力値xを掛算器70によって遅延レジスタ7
1に格納したxn-1(前回の計算結果)と乗算し、乗算結果
を加算器72により (入力端子221から与えられる数値)
2から引き、加算器72の出力に再び、掛算器73によってx
n-1を乗算して、新しい計算結果xnを得る。xnは次回の
計算のために遅延レジスタ71に格納する。以上の演算を
繰り返し実行すると、xnは1/x(入力値xの逆数)に収
束する。遅延レジスタに設定する初期値は結果に近い値
にすると、収束が速く3回ぐらいの繰り返しで十分な精
度が得られる。
【0079】図7(a)に示した逆数演算回路により、割
算を実行するには、図7(b)に示すように逆数演算回路7
4で除数xの逆数を求め、これを掛算器75で被除数18に乗
算することで、商y/xを計算できる。図7に示した実施
例により、割算命令の無いDSPでも精度良く逆数を求め
ることができる。
算を実行するには、図7(b)に示すように逆数演算回路7
4で除数xの逆数を求め、これを掛算器75で被除数18に乗
算することで、商y/xを計算できる。図7に示した実施
例により、割算命令の無いDSPでも精度良く逆数を求め
ることができる。
【0080】図1の実施例で示した逆特性計算部は、式
(8)を直接計算することにより、構成することができる
が、式(8)から解るように、割算演算が必要であり、図
7の逆数演算回路を用いても高速なDSPが必要となる。
そこで、図4に示すようにフィードバックループを用い
て、漸近的に逆特性の計算を式(10)によって行う逆特性
計算部の一実施例を図8に示す。
(8)を直接計算することにより、構成することができる
が、式(8)から解るように、割算演算が必要であり、図
7の逆数演算回路を用いても高速なDSPが必要となる。
そこで、図4に示すようにフィードバックループを用い
て、漸近的に逆特性の計算を式(10)によって行う逆特性
計算部の一実施例を図8に示す。
【0081】図8は本発明の逆特性計算部の一実施例の
構成を示すブロック図である。226はx2入力端子、227は
y2入力端子、228はxy cosφ入力端子、229は-xy sinφ
入力端子、80,82,83,86〜89,811,812,816,817は
掛算器、810,813は加算器、81は逆数平方根演算回路、
84,814,815は遅延レジスタ、85は累積乗算器、818は
三角関数累算器、230はaI出力端子、231はaQ出力端子、
818は掛算器88,89,811,812と加算器810,813と遅延
レジスタ814,815の回路要素で構成された三角関数累算
器である。
構成を示すブロック図である。226はx2入力端子、227は
y2入力端子、228はxy cosφ入力端子、229は-xy sinφ
入力端子、80,82,83,86〜89,811,812,816,817は
掛算器、810,813は加算器、81は逆数平方根演算回路、
84,814,815は遅延レジスタ、85は累積乗算器、818は
三角関数累算器、230はaI出力端子、231はaQ出力端子、
818は掛算器88,89,811,812と加算器810,813と遅延
レジスタ814,815の回路要素で構成された三角関数累算
器である。
【0082】図8において、第1の補間部7で求めた信
号レベル代表値piでの補間データx2,y2,xy cosφ,-x
y sinφをそれぞれ入力端子226〜229を介して入力す
る。即ち、x2は入力端子226を介して乗算器80と82とに
与えられ、y2は入力端子227を介して乗算器80に与えら
れ、xy cosφは入力端子228を介して乗算器86に与えら
れ、-xy sinφは入力端子229を介して乗算器87に与えら
れる。
号レベル代表値piでの補間データx2,y2,xy cosφ,-x
y sinφをそれぞれ入力端子226〜229を介して入力す
る。即ち、x2は入力端子226を介して乗算器80と82とに
与えられ、y2は入力端子227を介して乗算器80に与えら
れ、xy cosφは入力端子228を介して乗算器86に与えら
れ、-xy sinφは入力端子229を介して乗算器87に与えら
れる。
【0083】掛算器80で乗算し(出力= x2y2)、更に逆
数平方根演算回路81により1/xyを求める。求めた1/xy
は、掛算器82と掛算器86と掛算器87とに与えられる。掛
算器82により、x2及び1/xyを乗算して、n時点での入
出力信号の逆比an -1= x/yを計算する。逆比an -1は、遅
延レジスタ84と掛算器83で構成した累積乗算器85に与え
られ、掛算器83によって、前もって第3のメモリから読
み出して遅延レジスタ84にセットした1時点前の逆特性
係数An-1 -1と乗算して、n時点での逆特性係数An -1を計
算する。
数平方根演算回路81により1/xyを求める。求めた1/xy
は、掛算器82と掛算器86と掛算器87とに与えられる。掛
算器82により、x2及び1/xyを乗算して、n時点での入
出力信号の逆比an -1= x/yを計算する。逆比an -1は、遅
延レジスタ84と掛算器83で構成した累積乗算器85に与え
られ、掛算器83によって、前もって第3のメモリから読
み出して遅延レジスタ84にセットした1時点前の逆特性
係数An-1 -1と乗算して、n時点での逆特性係数An -1を計
算する。
【0084】他方、補間データxy cosφと-xy sinφ
は、掛算器86と87によって、それぞれ逆数平方根演算回
路81で得られた1/xyを乗算して、cosφと-sinφとす
る。求めた入出力位相差φの余弦値と正弦値は、三角関
数累算器818によって累算される。三角関数累算器818
は、式(18)の示す三角関数の加法公式を、回路により構
成したもので、累積位相差φnの余弦値と正弦値が遅延
レジスタ814と815とに格納される。 cosφn = cos(φn-1 + φ) 、 = cosφn-1cosφ - sinφn-1sinφ 、 -sinφn = -sin(φn-1 + φ) 、 = -sinφn-1cosφ - cosφn-1-sinφ 、 ‥‥‥式(18)
は、掛算器86と87によって、それぞれ逆数平方根演算回
路81で得られた1/xyを乗算して、cosφと-sinφとす
る。求めた入出力位相差φの余弦値と正弦値は、三角関
数累算器818によって累算される。三角関数累算器818
は、式(18)の示す三角関数の加法公式を、回路により構
成したもので、累積位相差φnの余弦値と正弦値が遅延
レジスタ814と815とに格納される。 cosφn = cos(φn-1 + φ) 、 = cosφn-1cosφ - sinφn-1sinφ 、 -sinφn = -sin(φn-1 + φ) 、 = -sinφn-1cosφ - cosφn-1-sinφ 、 ‥‥‥式(18)
【0085】三角関数累算器818の余弦値出力と正弦値
出力に、累積乗算器85で求めたAn -1を、掛算器816と817
とによってそれぞれ乗算し、逆特性データaI、aQとを求
め、出力端子230と231とからそれぞれ出力する。求めた
信号レベル代表値piでの逆特性データaI、aQは、第2の
メモリのi番目の記憶領域に格納する。また、累積乗算
器85の遅延レジスタ84のデータ及び三角関数累算器818
の遅延レジスタ814と815のデータを、第3のメモリのi
番目の記憶領域に格納する。こうして、全ての信号レベ
ル代表値pi(i = 0〜N)に対する逆特性データを計算す
る。
出力に、累積乗算器85で求めたAn -1を、掛算器816と817
とによってそれぞれ乗算し、逆特性データaI、aQとを求
め、出力端子230と231とからそれぞれ出力する。求めた
信号レベル代表値piでの逆特性データaI、aQは、第2の
メモリのi番目の記憶領域に格納する。また、累積乗算
器85の遅延レジスタ84のデータ及び三角関数累算器818
の遅延レジスタ814と815のデータを、第3のメモリのi
番目の記憶領域に格納する。こうして、全ての信号レベ
ル代表値pi(i = 0〜N)に対する逆特性データを計算す
る。
【0086】図8の実施例では、図4で説明したよう
に、帰還制御ループの中に累積乗算器があるので、n時
点での入出力信号レベル逆比an -1が多少の誤差を含んで
いても、An -1は求める値A-1(電力増幅器の入出力レベ
ル比の逆数)に収束する。図9によって、図8の実施例
で用いる逆数平方根演算回路の具体的な実施例を説明す
る。
に、帰還制御ループの中に累積乗算器があるので、n時
点での入出力信号レベル逆比an -1が多少の誤差を含んで
いても、An -1は求める値A-1(電力増幅器の入出力レベ
ル比の逆数)に収束する。図9によって、図8の実施例
で用いる逆数平方根演算回路の具体的な実施例を説明す
る。
【0087】図9は、本発明の逆数平方根演算回路の一
実施例の構成を示すブロック図である。232は入力端
子、90,91,93,94は掛算器、92は加算器、95は遅延レ
ジスタ、243は値3の入力端子、244は値1/2の入力端子、
233は出力端子である。逆数平方根を計算するために
は、図7の逆数演算と同様にニュートン法による漸化式
(式(12))を用いる。入力値aの逆数平方根x = a-1/2を
変形して、x-2 =aとし、f(x)= x-2 - a = 0を解く。f
(x)を微分すると、f′(x) = -2x-3となる。f(x)
およびf'(x)を式(12)に代入して、次式(19)がが求め
る漸化式である(ただし、繰り返し回数nは1ずらして
ある)。 xn = xn-1 + xn-1 3(xn-1 -2 - a )/2 = xn-1(3 - axn-1 2)/2 ‥‥‥式(19)
実施例の構成を示すブロック図である。232は入力端
子、90,91,93,94は掛算器、92は加算器、95は遅延レ
ジスタ、243は値3の入力端子、244は値1/2の入力端子、
233は出力端子である。逆数平方根を計算するために
は、図7の逆数演算と同様にニュートン法による漸化式
(式(12))を用いる。入力値aの逆数平方根x = a-1/2を
変形して、x-2 =aとし、f(x)= x-2 - a = 0を解く。f
(x)を微分すると、f′(x) = -2x-3となる。f(x)
およびf'(x)を式(12)に代入して、次式(19)がが求め
る漸化式である(ただし、繰り返し回数nは1ずらして
ある)。 xn = xn-1 + xn-1 3(xn-1 -2 - a )/2 = xn-1(3 - axn-1 2)/2 ‥‥‥式(19)
【0088】式(19)ように、まず、入力端子232を介し
て与えられた入力値xを掛算器90と91とよって遅延レジ
スタ95に格納したxn-1(前回の計算結果)と2度の乗算
を行う。次に乗算結果を加算器92により、入力端子243
から与えられる数値3から引く。そして、その減算結果
に、再び掛算器93によって遅延レジスタ95に格納したx
n-1を乗算する。更に掛算器94で、入力端子244から与え
られる数値1/2を掛けて、新しい計算結果xnを得出力端
子233を介して出力する。この新しい計算結果xnは、次
回の計算のために遅延レジスタ95に格納する。
て与えられた入力値xを掛算器90と91とよって遅延レジ
スタ95に格納したxn-1(前回の計算結果)と2度の乗算
を行う。次に乗算結果を加算器92により、入力端子243
から与えられる数値3から引く。そして、その減算結果
に、再び掛算器93によって遅延レジスタ95に格納したx
n-1を乗算する。更に掛算器94で、入力端子244から与え
られる数値1/2を掛けて、新しい計算結果xnを得出力端
子233を介して出力する。この新しい計算結果xnは、次
回の計算のために遅延レジスタ95に格納する。
【0089】以上の演算を繰り返し実行すると、xnはx
-1/2に(入力値xの逆数平方根)に収束する。遅延レジ
スタに設定する初期値は結果に近い値にすると、収束が
速い。図9に示した実施例では、掛算演算だけで割算を
用いないので、通常のDSPで精度良く逆数平方根を求め
ることができる。
-1/2に(入力値xの逆数平方根)に収束する。遅延レジ
スタに設定する初期値は結果に近い値にすると、収束が
速い。図9に示した実施例では、掛算演算だけで割算を
用いないので、通常のDSPで精度良く逆数平方根を求め
ることができる。
【0090】図10によって、本発明による第2の補間
部4の具体的実施例を説明する。図10は、本発明の第
2の補間部の一実施例の構成を示すブロック図である。
234は入力信号pの同相成分pI入力端子、235は入力信号p
の直交成分pQ入力端子、236はf0入力端子、237はf1入力
端子、238はf2入力端子、240は数値1/2入力端子、241
は1/Δp入力端子、100,102,104,107,108,1010は
加算器、101,102,105,106,109は掛算器である。
部4の具体的実施例を説明する。図10は、本発明の第
2の補間部の一実施例の構成を示すブロック図である。
234は入力信号pの同相成分pI入力端子、235は入力信号p
の直交成分pQ入力端子、236はf0入力端子、237はf1入力
端子、238はf2入力端子、240は数値1/2入力端子、241
は1/Δp入力端子、100,102,104,107,108,1010は
加算器、101,102,105,106,109は掛算器である。
【0091】補間には、図6の実施例と同様に、式(4)
に示したラグランジの補間多項式を用いる。しかしなが
ら、第2の補間部4では、補間係数の分母は信号レベル
代表値pi = iΔp(Δp = pmax/N )で与えられるの
で、補間係数をΔpで正規化すると、式(4)の補間多項式
は、p - pK(KはM/2(M = 偶数)、または(M - 1)/
2( M = 奇数))べき乗級数の形に変形できて、式(20)
で表される。 f(p) = Σj=0 Mfj[Πk≠j(p - pk)/(pj - pk)] = Σl=0 Mdl(fj、pj)(p - pK)l ‥‥‥式(20)
に示したラグランジの補間多項式を用いる。しかしなが
ら、第2の補間部4では、補間係数の分母は信号レベル
代表値pi = iΔp(Δp = pmax/N )で与えられるの
で、補間係数をΔpで正規化すると、式(4)の補間多項式
は、p - pK(KはM/2(M = 偶数)、または(M - 1)/
2( M = 奇数))べき乗級数の形に変形できて、式(20)
で表される。 f(p) = Σj=0 Mfj[Πk≠j(p - pk)/(pj - pk)] = Σl=0 Mdl(fj、pj)(p - pK)l ‥‥‥式(20)
【0092】式(20)における補間係数dl(fj、pj)を2
次補間多項式の場合について具体的に計算すると、式(1
4)において、p1 - p2 = -Δp 、 p2 - P0 = 2Δp
、 p0 - p1 = -Δpとなり、またp - p1 = Δ・Δpと
おくと、p - p0 = (Δ + 1)・Δp 、 p - p2 = (Δ
- 1)・Δpとなるので、補間多項式は式(21)で与えられ
る(ただし、Δ = (p - p1)/Δp)。 f(p)= f0/2(Δ2 - Δ)- f1(Δ2 - 1)+ f2/2(Δ2 + Δ) =(f0/2 - f1 + f2/2)Δ2 + (f2 - f0)/2Δ + f1 = d2Δ2 + d1Δ + d0 ‥‥‥式(21) 即ち式(21)では、補間係数の計算に割算演算が不要とな
る。
次補間多項式の場合について具体的に計算すると、式(1
4)において、p1 - p2 = -Δp 、 p2 - P0 = 2Δp
、 p0 - p1 = -Δpとなり、またp - p1 = Δ・Δpと
おくと、p - p0 = (Δ + 1)・Δp 、 p - p2 = (Δ
- 1)・Δpとなるので、補間多項式は式(21)で与えられ
る(ただし、Δ = (p - p1)/Δp)。 f(p)= f0/2(Δ2 - Δ)- f1(Δ2 - 1)+ f2/2(Δ2 + Δ) =(f0/2 - f1 + f2/2)Δ2 + (f2 - f0)/2Δ + f1 = d2Δ2 + d1Δ + d0 ‥‥‥式(21) 即ち式(21)では、補間係数の計算に割算演算が不要とな
る。
【0093】式(21)の演算を回路で実行する実施例を図
10によって説明する。入力信号レベルpに最も近い信
号レベル代表値pi = iΔpを求め、第2のメモリのi-1、
i、i+1番目の記憶領域から、逆特性データaIを読み出
し、p1 = pi、f0 = aI i-1 、 f1 = aI i 、 f2 = aI
i+1とする。加算器100によって、入力端子234と235とか
ら与えられる入力信号p0とp1とから同相成分p-p1を計算
し、更に、掛算器101で、入力端子241から与えられる1
/Δpを乗算し、Δ = (p - p1)/Δpを求める。求め
たΔは加算器105と掛算器109とにそれぞれ与えられる。
10によって説明する。入力信号レベルpに最も近い信
号レベル代表値pi = iΔpを求め、第2のメモリのi-1、
i、i+1番目の記憶領域から、逆特性データaIを読み出
し、p1 = pi、f0 = aI i-1 、 f1 = aI i 、 f2 = aI
i+1とする。加算器100によって、入力端子234と235とか
ら与えられる入力信号p0とp1とから同相成分p-p1を計算
し、更に、掛算器101で、入力端子241から与えられる1
/Δpを乗算し、Δ = (p - p1)/Δpを求める。求め
たΔは加算器105と掛算器109とにそれぞれ与えられる。
【0094】また、掛算器102と106により、入力端子23
6から与えられるf0と入力端子238から与えられるf2と
に、入力端子240から与えられる数値1/2を掛けて、そ
れぞれ1/2とする。そして、この1/2f0のデータと1/2
f2のデータとはそれぞれ加算器103と加算器107とに与え
られる。また、入力端子237を介して、f1が加算器104と
加算器1010とに系数d0として与えられる。
6から与えられるf0と入力端子238から与えられるf2と
に、入力端子240から与えられる数値1/2を掛けて、そ
れぞれ1/2とする。そして、この1/2f0のデータと1/2
f2のデータとはそれぞれ加算器103と加算器107とに与え
られる。また、入力端子237を介して、f1が加算器104と
加算器1010とに系数d0として与えられる。
【0095】加算器103では加算を行い、1/2f0 + 1/2
f2を出力し加算器104に与え、加算器107ではd0を減算し
て、1/2f2 - 1/2f0( = d1)を出力し加算器108に与
える。更に、加算器104では、f1を減算して、係数d2 =
f0/2 - f1 + f2/2を得、系数d2を加算器105に与え
る。
f2を出力し加算器104に与え、加算器107ではd0を減算し
て、1/2f2 - 1/2f0( = d1)を出力し加算器108に与
える。更に、加算器104では、f1を減算して、係数d2 =
f0/2 - f1 + f2/2を得、系数d2を加算器105に与え
る。
【0096】掛算器105では、入力したΔと系数d2とを
乗算して(d2Δ)加算器108に与え、加算器108はその値
を系数d1と加算して(d1+d2Δ)掛算器109に与える。掛
算器109では、加算器108の値(d1+d2Δ)にΔを乗算し
て(d1Δ+d2Δ2)、加算器1010に与える。加算器1010で
は、(d1Δ+d2Δ2)にd0( = f0)を加え、補間出力f
(Δ)、即ち、入力信号レベルpにおける逆特性データa
Iが得られる。同様の演算をaQについても行う。図10
の実施例では、aIとaQについて同じ構成で計算ができ、
しかも割算演算が不要であり、高速かつ高精度な補間演
算が可能となる。
乗算して(d2Δ)加算器108に与え、加算器108はその値
を系数d1と加算して(d1+d2Δ)掛算器109に与える。掛
算器109では、加算器108の値(d1+d2Δ)にΔを乗算し
て(d1Δ+d2Δ2)、加算器1010に与える。加算器1010で
は、(d1Δ+d2Δ2)にd0( = f0)を加え、補間出力f
(Δ)、即ち、入力信号レベルpにおける逆特性データa
Iが得られる。同様の演算をaQについても行う。図10
の実施例では、aIとaQについて同じ構成で計算ができ、
しかも割算演算が不要であり、高速かつ高精度な補間演
算が可能となる。
【0097】以上説明した本発明の実施例では、フィー
ドバックループを用いているので、電力増幅器の特性変
動、雑音の混入、演算の誤差等があっても自動的に収束
するように動作する。また、無線機の動作終了時に、第
2および第3のメモリに記憶してあるデータを、不揮発
性のメモリに退避し、動作再開時に退避したデータを第
2、第3のメモリに初期設定してから動作させることに
より、フィードバックループの収束を非常に短時間に行
える利点がある。
ドバックループを用いているので、電力増幅器の特性変
動、雑音の混入、演算の誤差等があっても自動的に収束
するように動作する。また、無線機の動作終了時に、第
2および第3のメモリに記憶してあるデータを、不揮発
性のメモリに退避し、動作再開時に退避したデータを第
2、第3のメモリに初期設定してから動作させることに
より、フィードバックループの収束を非常に短時間に行
える利点がある。
【0098】今回提案するアダプティブプリディストー
ションの概略構成を図11に示す。図11は本発明のリ
ニアライザの一実施例の構成を示すブロック図である。
HPA18′の非線形特性を取得するため、局部直交復調器1
1′を備えた構成としている。出力信号yと入力信号xの
比(y/x = an・ejΔφ)を計算し、これを逆特性計算
部で累積乗算(An = an・An-1)および累積加算(φn =
Δφ+ φn-1)して非線形特性(Anexp(jφn))を求める。
この逆数を入力信号に複素乗算して非線形補償信号x′
を得る。逆特性計算の構成を簡単化するため、HPA18′
入力信号変化範囲をN分割し、分割点でのみ逆特性を求
め、他の点では補間、および逆補間によって計算を行
う。以上の構成により本方式のアダプティブプリディス
トーションでは、累積乗算器247、および累積加算器248
のレジスタに初期値1、および0を与え漸近的に非線形特
性を求めていくので、HPA18′の非線形特性を測定する
ための特別なテスト信号は必要としない。
ションの概略構成を図11に示す。図11は本発明のリ
ニアライザの一実施例の構成を示すブロック図である。
HPA18′の非線形特性を取得するため、局部直交復調器1
1′を備えた構成としている。出力信号yと入力信号xの
比(y/x = an・ejΔφ)を計算し、これを逆特性計算
部で累積乗算(An = an・An-1)および累積加算(φn =
Δφ+ φn-1)して非線形特性(Anexp(jφn))を求める。
この逆数を入力信号に複素乗算して非線形補償信号x′
を得る。逆特性計算の構成を簡単化するため、HPA18′
入力信号変化範囲をN分割し、分割点でのみ逆特性を求
め、他の点では補間、および逆補間によって計算を行
う。以上の構成により本方式のアダプティブプリディス
トーションでは、累積乗算器247、および累積加算器248
のレジスタに初期値1、および0を与え漸近的に非線形特
性を求めていくので、HPA18′の非線形特性を測定する
ための特別なテスト信号は必要としない。
【0099】HPAの特性モデルとして、図12に示すよ
うな振幅位相非線形特性を与えて計算機シミュレーショ
ンを行った。変調方式、及びアダプティブプリディスト
ーション部の諸元を図13に示す。補償前後のHPA出力
信号スペクトルを図14に示す。
うな振幅位相非線形特性を与えて計算機シミュレーショ
ンを行った。変調方式、及びアダプティブプリディスト
ーション部の諸元を図13に示す。補償前後のHPA出力
信号スペクトルを図14に示す。
【0100】隣接チャネル(fC±12.5 KHz、fCは搬送周
波数)において約20 dBの改善が得られた。また、初期
値( A0 = 1、φ0 = 0 )からの収束時間は約400シンボ
ル(50 msec)であった。漏洩電力改善量はレベル分割
数Nに依存し、分割数Nを上げれば、更に改善が期待でき
る。また、収束時間は逆特性計算部の初期値を前回のデ
ータとすることで、大幅に短縮することが可能である。
波数)において約20 dBの改善が得られた。また、初期
値( A0 = 1、φ0 = 0 )からの収束時間は約400シンボ
ル(50 msec)であった。漏洩電力改善量はレベル分割
数Nに依存し、分割数Nを上げれば、更に改善が期待でき
る。また、収束時間は逆特性計算部の初期値を前回のデ
ータとすることで、大幅に短縮することが可能である。
【0101】
【発明の効果】以上のように、本発明によれば、電力増
幅器の非線形振幅特性をアダプティブプリディストーシ
ョン方式を用いて補償する非線形補償回路において、参
照テーブルや、近似式を用いずに直接非線形特性の逆特
性を計算することができ、高精度な非線形補償が可能で
ある。また、特別なテスト信号を必要としない。更に、
逆特性の計算は少数の信号レベル代表点についてのみ行
うので、非線形特性データの取得はレベル変動範囲全部
に亘って行う必要がない。そのため、非線形補償回路の
動作を高速にすることができ、また回路規模も小さくす
ることができる。また更には、参照テーブルを用いない
ので、製造時に初期特性を測定する必要が無く、特別な
テスト信号による検査も不必要であり、無線機の通常動
作時に自動的に非線形補償動作が実行することができ
る。また、電力増幅器の特性が温度、電源電圧、経時変
動などで変動しても、局部復調器を備えて常に特性監視
しているので、非線形補償特性が変わることは無い。ま
た、本発明の実施については、割算演算、関数計算等の
処理を必要としないので、通常用いられているディジタ
ル回路による方法でも、DSPによるソフトウェア処理で
も実現が可能であり、安価、小型、低消費電力の非線形
補償回路が得られる。
幅器の非線形振幅特性をアダプティブプリディストーシ
ョン方式を用いて補償する非線形補償回路において、参
照テーブルや、近似式を用いずに直接非線形特性の逆特
性を計算することができ、高精度な非線形補償が可能で
ある。また、特別なテスト信号を必要としない。更に、
逆特性の計算は少数の信号レベル代表点についてのみ行
うので、非線形特性データの取得はレベル変動範囲全部
に亘って行う必要がない。そのため、非線形補償回路の
動作を高速にすることができ、また回路規模も小さくす
ることができる。また更には、参照テーブルを用いない
ので、製造時に初期特性を測定する必要が無く、特別な
テスト信号による検査も不必要であり、無線機の通常動
作時に自動的に非線形補償動作が実行することができ
る。また、電力増幅器の特性が温度、電源電圧、経時変
動などで変動しても、局部復調器を備えて常に特性監視
しているので、非線形補償特性が変わることは無い。ま
た、本発明の実施については、割算演算、関数計算等の
処理を必要としないので、通常用いられているディジタ
ル回路による方法でも、DSPによるソフトウェア処理で
も実現が可能であり、安価、小型、低消費電力の非線形
補償回路が得られる。
【図1】 本発明のリニアライザの一実施例の構成を示
すブロック図。
すブロック図。
【図2】 プリディストーション方式の非線形補償の原
理を説明する図。
理を説明する図。
【図3】 従来のアダプティブプリディストーション方
式リニアライザの構成を示すブロック図。
式リニアライザの構成を示すブロック図。
【図4】 本発明の非線形補償ループの動作原理を説明
するための構成を示すブロック図。
するための構成を示すブロック図。
【図5】 本発明のデータ入力部とデータメモリ部の一
実施例の構成を示すブロック図。
実施例の構成を示すブロック図。
【図6】 本発明の第1の補間部の一実施例の構成を示
すブロック図。
すブロック図。
【図7】 本発明の逆数演算回路の一実施例の構成を示
すブロック図。
すブロック図。
【図8】 本発明の逆特性計算部の一実施例の構成を示
すブロック図。
すブロック図。
【図9】 本発明の逆数平方根演算回路の一実施例の構
成を示すブロック図。
成を示すブロック図。
【図10】 本発明の第2の補間部の一実施例の構成を
示すブロック図。
示すブロック図。
【図11】 本発明のリニアライザの一実施例の構成を
示すブロック図。
示すブロック図。
【図12】 本発明の計算機シミュレーションに使用し
たHPA特性モデルを示す図。
たHPA特性モデルを示す図。
【図13】 本発明の計算機シミュレーションに使用し
た変調方式及びアダプティブプリディストーション部の
諸元を表した図。
た変調方式及びアダプティブプリディストーション部の
諸元を表した図。
【図14】 本発明の計算機シミュレーションの結果を
表す図。
表す図。
1:非線形補償部、 2:レベル検出器、 3:複素掛算
器、 4,4′:第2の補間部、 5:第2のメモリ、
6,6′:逆特性計算部、 7,7′:第1の補間部、
8:第1のメモリ、 9,9′:データ入力部、 10,1
0′:直交変調器、11:直交復調器、 11′:局部直交
復調器、 12:D/A変換器、 13:A/D変換器、 14:周
波数混合器、 15,16:局部発振器、 17:周波数混合
器、 18:HPA、 19:方向性結合器、 20:複素掛算
器、 21:逆特性計算部、 22:直交変調器、 23:直
交復調器、 24:HPA、 30:参照テーブル、 31:パ
ワー計算部、 32:複素掛算器、 33:複素掛算器、
34:係数計算部、 35:データ更新部、 40:逆特性補
償器、 41:遅延レジスタ、 42:掛算器、 43:累積
乗算器、 44:加算器、 45:遅延レジスタ、 46:累
積加算器、 47:複素数比演算回路、 48:HPA、 5
0,51:遅延回路、 52:直交復調器、 53,54:掛算
器、 55:加算器、 56,57:掛算器、 58:加算器、
59:掛算器、60〜65:加算器、 66〜68:割算器、
69:掛算器、 70:掛算器、 71:遅延レジスタ、 7
2:加算器、 73:掛算器、 74:逆数演算回路、 7
5:掛算器、80:掛算器、 81:逆数平方根演算回路、
82,83:掛算器、 84:遅延レジスタ、 85:累積乗
算器、 86〜89:掛算器、 90,91:掛算器、 92:加
算器、 93,94:掛算器、 95:遅延レジスタ、 10
0:加算器、 101,102:掛算器、 103,104:加算
器、 105,106:掛算器、 107,108:加算器、 10
9:掛算器、 201,202:入力端子、 203:出力端子、
217,218:入力端子、 219:出力端子、 245,24
6:遅延レジスタ、 247:累積乗算器、 248:累積加
算器、 510:掛算器、 511:加算器、 512,513:掛
算器、 514:加算器、515:メモリ制御回路、 516:
メモリ、 610,611:掛算器、 612〜614:データレジ
スタ、 615〜618:掛算器、 619,620:加算器、 62
1〜623:平方根演算回路、 810:加算器、 811,81
2:掛算器、 813:加算器、 814:遅延レジスタ、 8
15:遅延レジスタ、 816,817:掛算器、 818:三角
関数累算器、 1010:加算器。
器、 4,4′:第2の補間部、 5:第2のメモリ、
6,6′:逆特性計算部、 7,7′:第1の補間部、
8:第1のメモリ、 9,9′:データ入力部、 10,1
0′:直交変調器、11:直交復調器、 11′:局部直交
復調器、 12:D/A変換器、 13:A/D変換器、 14:周
波数混合器、 15,16:局部発振器、 17:周波数混合
器、 18:HPA、 19:方向性結合器、 20:複素掛算
器、 21:逆特性計算部、 22:直交変調器、 23:直
交復調器、 24:HPA、 30:参照テーブル、 31:パ
ワー計算部、 32:複素掛算器、 33:複素掛算器、
34:係数計算部、 35:データ更新部、 40:逆特性補
償器、 41:遅延レジスタ、 42:掛算器、 43:累積
乗算器、 44:加算器、 45:遅延レジスタ、 46:累
積加算器、 47:複素数比演算回路、 48:HPA、 5
0,51:遅延回路、 52:直交復調器、 53,54:掛算
器、 55:加算器、 56,57:掛算器、 58:加算器、
59:掛算器、60〜65:加算器、 66〜68:割算器、
69:掛算器、 70:掛算器、 71:遅延レジスタ、 7
2:加算器、 73:掛算器、 74:逆数演算回路、 7
5:掛算器、80:掛算器、 81:逆数平方根演算回路、
82,83:掛算器、 84:遅延レジスタ、 85:累積乗
算器、 86〜89:掛算器、 90,91:掛算器、 92:加
算器、 93,94:掛算器、 95:遅延レジスタ、 10
0:加算器、 101,102:掛算器、 103,104:加算
器、 105,106:掛算器、 107,108:加算器、 10
9:掛算器、 201,202:入力端子、 203:出力端子、
217,218:入力端子、 219:出力端子、 245,24
6:遅延レジスタ、 247:累積乗算器、 248:累積加
算器、 510:掛算器、 511:加算器、 512,513:掛
算器、 514:加算器、515:メモリ制御回路、 516:
メモリ、 610,611:掛算器、 612〜614:データレジ
スタ、 615〜618:掛算器、 619,620:加算器、 62
1〜623:平方根演算回路、 810:加算器、 811,81
2:掛算器、 813:加算器、 814:遅延レジスタ、 8
15:遅延レジスタ、 816,817:掛算器、 818:三角
関数累算器、 1010:加算器。
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成12年4月19日(2000.4.1
9)
9)
【手続補正1】
【補正対象書類名】図面
【補正対象項目名】図1
【補正方法】変更
【補正内容】
【図1】
【手続補正2】
【補正対象書類名】図面
【補正対象項目名】図11
【補正方法】変更
【補正内容】
【図11】
───────────────────────────────────────────────────── フロントページの続き Fターム(参考) 5J090 AA01 AA41 CA21 CA65 CA92 FA08 FA17 GN03 KA00 KA15 KA26 KA32 KA33 KA34 KA46 KA53 KA55 KA68 MA11 NN11 SA14 TA01 TA02 TA03 5J091 AA01 AA41 CA21 CA65 CA92 FA08 FA17 KA00 KA15 KA26 KA32 KA33 KA34 KA46 KA53 KA55 KA68 MA11 SA14 TA01 TA02 TA03 5K004 AA05 FA09 FE10
Claims (18)
- 【請求項1】 ディジタル無線機に用いる電力増幅器の
非線形特性を補償するために、該電力増幅器で発生する
該非線形特性の逆特性を求め、求めた該逆特性を該電力
増幅器に入力する入力信号に与えるプリディストーショ
ン方式のリニアライザにおいて、 前記電力増幅器の出力信号と前記入力信号の比を計算す
るアダプティブプリディストーション部を設け、 計算した該比によって前記電力増幅器の前記非線形性の
逆特性を求めることを特徴とするリニアライザ。 - 【請求項2】 請求項1記載のリニアライザにおいて、
前記電力増幅器の前記入力信号の変化範囲をN分割し(N
は自然数)、該分割した点で前記非線形性の逆特性を求
め、他の点では補間及び逆補間を行うことによって前記
非線形性の逆特性を求めることを特徴とするリニアライ
ザ。 - 【請求項3】 局部直交復調器を用いて電力増幅器の出
力をモニタし、入力信号と出力信号の比を計算すること
によって、前記電力増幅器の非線形特性を取得し、該非
線形特性を補償する逆特性の非線形性を前記入力信号に
与えた非線形補償信号を、直交変調器を介して前記電力
増幅器に入力することにより、前記電力増幅器の非線形
特性を補償する非線形補償回路において、 前記入力信号xの同相成分(xI)及び直交成分(xQ)の
2乗和の平方根を計算して入力信号レベルp(p2 = xI 2
+ xQ 2)を求めるレベル検出器と、前記入力信号xの同相
成分(xI)と直交成分(xQ)データ、及び前記出力信号
yの同相成分(yI)と直交成分(yQ)データとを取得す
るデータ入力部と、取得した入力信号データ(xIとxQ)
及び出力信号データ(yIとyQ)を x2 = xI 2 + xQ 2、 y2 = yI 2 + yQ 2、 xy cosφ = xIyI + xQyyQ、 -xy sinφ = xQyI - xIyQ (ただし、φは前記入力信号xと前記出力信号yの位相
差)の4つの値からなる中間処理データに変換し、該中
間処理データを記憶しておく第1のメモリと、該第1の
メモリに記憶した中間処理データから補間により、予め
定めた複数の入力信号レベルに対する中間処理データ値
を計算する第1の補間部と、該第1の補間部によって求
めた中間処理データから、前記予め定めた複数の入力信
号レベルにおける前記電力増幅器の非線形特性の逆特性
を計算する逆特性計算部と、該逆特性計算部により求め
られた逆特性データを記憶しておく第2のメモリと、該
第2のメモリに記憶した逆特性データから補間により、
前記レベル検出器で検出された入力信号レベルに対応す
る逆特性データを計算する第2の補間部と、該第2の補
間部で求めた逆特性データを前記入力信号に乗積し非線
形補償動作を行う複素掛算器とから構成することを特徴
とするリニアライザ。 - 【請求項4】 請求項3記載のリニアライザにおいて、
前記予め定めた複数の入力信号レベルを、0から最大入
力レベル(pmax)までのレベル範囲を予め定めた間隔で
N分割(Nは自然数)した入力レベル代表値pi(i= 0〜
N、p0 = 0、pN= pmax )とし、前記逆特性計算部の計算
と、前記第2のメモリのデータ記憶を入力レベル代表値
でのみ行うことを特徴とするリニアライザ。 - 【請求項5】 請求項3または請求項4記載のリニアラ
イザのレベル検出器において、前記入力信号振幅の2乗
和の平方根(p =(xI 2 + xQ 2)1/2)の計算を、該振幅
2乗和の逆数平方根演算(aに対してa-1/2を求める)
と、更にその逆数演算(bに対してb-1を求める)により
行い、該逆数平方根演算を、ニュートン法により逆数平
方根を求める漸化式 zn+1 = zn( 3 - azn 2)/2(た
だし、aは入力値、nは自然数、z1は初期値、znはn回
目の計算結果)を用いて計算し、 前記逆数演算を、ニュートン法により逆数を求める漸化
式 zn+1 = z(ただし、bは入力値)を用いることで漸
近的に前記信号レベルの計算を行うことを特徴とする非
線形補償回路。 - 【請求項6】 請求項3〜請求項5のいずれかに記載の
リニアライザにおいて、前記データ入力部で取得した入
力信号データ、及びレベル検出器で求めた入力信号レベ
ルのデータを遅延して、データ入力部で得られる出力信
号データと時間合わせを行い、求めた入力信号レベルp
( =(xI 2 + xQ 2)1/2 )の値が前記入力レベル代表値
pi(pi = iΔp)に対して、 (pi-1 + pi)/2 <p<(pi + pi+1)/2の範囲にある
入力信号レベルp及び入出力信号データ(xI、xQ、yI、y
Q)から計算される前記中間処理データ値を前記第1の
メモリのi番目のデータとして書き込むと同時に、デー
タが書き込まれたことを示すフラグFiをセットし、全て
のフラグがセットされた時点で前記第1の補間部にデー
タを転送することを特徴とするリニアライザ。 - 【請求項7】 請求項6記載のリニアライザにおいて、
前記データ入力部で得られたデータが、既に前記第1の
メモリに書き込まれている場合、新たに得られたデータ
および既に書き込まれているデータの入力信号レベルを
比較し、入力レベル代表値 pi(pi = iΔp)に近いほ
うのデータを記憶保存することを特徴とするリニアライ
ザ。 - 【請求項8】 請求項3〜請求項7のいずれかに記載の
リニアライザにおいて、前記第1のメモリに書き込みが
終了した時点で、i番目の入力レベル代表値p iに対する
中間処理データ値( x2、y2、xy cosφ、-xy sinφ )
を、iを挟むM個(Mは正整数)のデータ(j = i−M/2
〜 i + M/2 )を補間して求め、前記N個の入力レベル
代表値pi(i = 0 〜N )に対する中間処理データ値が計
算終了した時点で前記第1のメモリのフラグをリセット
するように、前記第1の補間部を構成することを特徴と
するリニアライザ。 - 【請求項9】 請求項8記載のリニアライザの第1の補
間部において、補間にラグランジのM次補間多項式 (ただし、x = pi、fはxI、xQ、yI、yQ、また、Πk≠j
はk = j以外のすべての項の積を表わす) を用い、前記第1のメモリからpj及びfj(j = i - M/2
〜 i + M/2 )を読み出して補間係数ccj(x)を求
め、前記補間多項式によって入力レベル代表値pi( i =
0 〜 N )に対する中間処理データ値( x2、yy2、xy
cosφ、-xy sinφ)を計算することを特徴とするリニア
ライザ。 - 【請求項10】 請求項9記載のリニアライザの第1の
補間部において、前記補間係数 cj(x) = Πk≠j(x - pk)/( pj - pk) の計算における割算演算を、前記ニュートン法により逆
数を求める漸化式で求めた逆数を乗算することにより行
うことを特徴とするリニアライザ。 - 【請求項11】 請求項3〜請求項10記載のいずれか
のリニアライザにおいて、前記第1の補間部で求めた入
力レベル代表値piでの中間処理データ値( x 2、y2、xy
cosφ、-xy sinφ )を前記逆特性計算部に入力し、該
逆特性計算部において、入力レベル代表値pi( i = 0
〜 N )での逆特性データaI、aQを、 aI =(xIyI + xQyQ)/(yI 2 + yQ 2)= xy cosφ/y2 aQ =(xQyI - xIyQ)/(yI 2 + yQ 2)= -xy sinφ/y2 (ただし、xy =((xI 2 + xQ 2)(yI 2 + yQ 2))1/2 )
によって計算し、前記第2のメモリに格納することを特
徴とするリニアライザ。 - 【請求項12】 請求項11記載のリニアライザの逆特
性計算部において、前記逆特性データaIとaQを求める計
算における出力信号振幅2乗和(yI 2 + yQ 2)による割
算演算を、前記ニュートン法により逆数を求める漸化式
で求めた逆数を乗算することにより計算することを特徴
とするリニアライザ。 - 【請求項13】 請求項11記載のリニアライザの逆特
性計算部において、前記振幅補償係数の計算を前記入力
信号振幅の2乗値(x2 = xI 2 + xQ 2 )と、前記出力信
号振幅の2乗値(y2 = yI 2 + yQ 2)の積の逆数平方根
( 1/xy)を求め、該入出力信号振幅の積の逆数に前記
入力信号振幅2乗値xx2を掛けて、入出力信号振幅比
(an -1 = x/y =((xI 2 + xQ 2)/(yI 2 +
yQ 2))1/2)を求め、該入出力信号振幅比を、乗算器と
遅延レジスタから構成した累積乗算器により累積乗算し
て、振幅補償係数(An -1 = Πn(an -1)、Πnはn項ま
での累乗積を表わす)を求め、前記入力レベル代表値pi
でのデータ値のxy cosφ、及び-xy sinφに前記入出力
信号振幅の積の逆数( 1/xy)を掛けて、入出力信号位
相差の余弦値cosφ、及び正弦値-sinφを求め、該入出
力信号位相差の余弦値および正弦値を、乗算器と加算器
と遅延レジスタから構成した累積加算器により、 (ただしφnはn回目の累積位相差)を計算することに
より累積加算して、位相補償値の余弦値、および正弦値
を求め、該位相補償値の余弦値、および正弦値に前記振
幅補償係数An -1を掛けて、 aI = An -1cosφn aQ = -A-1sinφn を計算することによって前記逆特性データaIとaQを求め
ることを特徴とするリニアライザ。 - 【請求項14】 請求項13記載のリニアライザの逆特
性計算部において、前記入出力信号振幅の2乗値の積の
逆数平方根を求める演算を、前記ニュートン法により逆
数平方根を求める漸化式を用いて計算し、該漸化式の途
中結果および、前記累積乗算器と前記累積加算器の遅延
レジスタのデータを第3のメモリに格納し、前記漸化式
および累積乗算器、累積加算器の計算を漸近的に行うこ
とを特徴とするリニアライザ。 - 【請求項15】 請求項13または請求項14記載のリ
ニアライザにおいて、前記リニアライザの動作終了時に
前記第2および第3のメモリに記憶されているデータ
を、不揮発性のメモリに退避し、動作再開時に該退避し
たデータを前記第2および第3のメモリに初期設定して
から動作させることを特徴とするリニアライザ。 - 【請求項16】 請求項3〜請求項15のいずれかに記
載のリニアライザにおいて、前記レベル検出器により検
出された入力信号レベルpを挟むM個の入力レベル代表値
pi(i = 1 〜 M、p1<p<pi)に対する逆特性データを
前記第2のメモリから読み出し、補間によって入力信号
レベルpに対する逆特性データ(aI +jaQ)を計算するよ
うに前記第2の補間部を構成し、前記複素掛算器によっ
て前記逆特性データを入力信号(xI + jxQ)と複素乗積
することにより非線形補償動作を行うことを特徴とする
リニアライザ。 - 【請求項17】 請求項16記載のリニアライザの第2
の補間部において、補間にラグランジのM次補間多項式 f(p) = Σj=0 Mfj[Πk≠j(p - pk)/(pj - pk)] = Σj=0 Mfjcj(p)) (ただし、fはaI、aQ )を用い、前記第2のメモリから
pj及びfj(j = i - M/2 〜 i + M/2 )を読み出して
補間係数cj(p)を求め、前記補間多項式によって入力
信号レベルpに対する逆特性データ値(aI、aQ)を計算
することを特徴とするリニアライザ。 - 【請求項18】 請求項17記載のリニアライザの第2
の補間部において、前記ラグランジのM次補間多項式
を、前記入力信号レベルpと入力レベル代表値pKの差 (p - pK)(但し、KはM = 偶数のときM/2、M = 奇数
のとき(M - 1)/2) のべき乗級数の形式に変形した f(p) = Σj=0 Mfj[Πk≠j(p - pk)/(pj - pk)] = Σl=0 Mdl(fj、pj)(p - pK)l (ただし、fはaI、aQ ) を補間に用い、前記第2のメモリからpj及びfj(j = i
- M/2 〜 i +M/2 )を読み出して補間係数dl(fj、
pj))を求め、前記補間多項式によって入力信号レベルp
に対する逆特性データ値(aI、aaQ )を計算すること
を特徴とするリニアライザ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000077353A JP2001268150A (ja) | 2000-03-21 | 2000-03-21 | リニアライザ |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000077353A JP2001268150A (ja) | 2000-03-21 | 2000-03-21 | リニアライザ |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2001268150A true JP2001268150A (ja) | 2001-09-28 |
Family
ID=18594934
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000077353A Pending JP2001268150A (ja) | 2000-03-21 | 2000-03-21 | リニアライザ |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2001268150A (ja) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002344418A (ja) * | 2001-05-22 | 2002-11-29 | Sumitomo Electric Ind Ltd | 前置型非線形歪補償器 |
JP2003174332A (ja) * | 2001-12-05 | 2003-06-20 | Hitachi Ltd | プリディストーション型増幅装置 |
JP2004015310A (ja) * | 2002-06-05 | 2004-01-15 | Mitsubishi Electric Corp | 歪補償回路および送信機 |
JP2004032609A (ja) * | 2002-06-28 | 2004-01-29 | Nec Corp | 非線形歪み補償回路 |
JP2004088303A (ja) * | 2002-08-26 | 2004-03-18 | Ntt Docomo Inc | 送信機 |
JP2004112151A (ja) * | 2002-09-17 | 2004-04-08 | Hitachi Kokusai Electric Inc | 歪み補償方法及び送信機 |
JP2004363733A (ja) * | 2003-06-02 | 2004-12-24 | Sumitomo Electric Ind Ltd | 前置型非線形歪補償器 |
JP2005073032A (ja) * | 2003-08-26 | 2005-03-17 | Hitachi Kokusai Electric Inc | 歪補償増幅装置及び歪補償方法 |
JP2005072747A (ja) * | 2003-08-21 | 2005-03-17 | Nippon Telegr & Teleph Corp <Ntt> | 非線形歪補償装置 |
WO2006030919A1 (ja) * | 2004-09-15 | 2006-03-23 | Advantest Corporation | 逆特性測定装置および歪み補償装置、方法、プログラム、記録媒体 |
US7170342B2 (en) | 2002-12-10 | 2007-01-30 | Ntt Docomo, Inc. | Linear power amplification method and linear power amplifier |
JP2007518378A (ja) * | 2004-01-16 | 2007-07-05 | クゥアルコム・インコーポレイテッド | 増幅器予歪および自動較正の方法および装置 |
JP2008503169A (ja) * | 2004-06-30 | 2008-01-31 | インテル・コーポレーション | 周波数領域内で予歪を用いる電力増幅器の線形化方法および装置 |
JP2010041470A (ja) * | 2008-08-06 | 2010-02-18 | Mitsubishi Electric Corp | 送信機及び送受信装置 |
EP2202879A1 (en) | 2008-12-26 | 2010-06-30 | Fujitsu Limited | A predistortion apparatus and predistortion method |
WO2010073803A1 (ja) * | 2008-12-27 | 2010-07-01 | 住友電気工業株式会社 | 歪補償回路及び無線基地局 |
JP2010157937A (ja) * | 2008-12-27 | 2010-07-15 | Sumitomo Electric Ind Ltd | 歪補償回路及び無線基地局 |
US7864881B2 (en) | 2006-04-11 | 2011-01-04 | Hitachi, Ltd. | Digital predistortion transmitter |
US8179123B2 (en) | 2007-05-28 | 2012-05-15 | Advantest Corporation | Measurement apparatus and recording medium |
US8270530B2 (en) | 2008-09-10 | 2012-09-18 | Fujitsu Limited | Apparatus and method of compensating distortion |
JP2018182895A (ja) * | 2017-04-12 | 2018-11-15 | キヤノン株式会社 | モータ制御装置、方法及びプログラム |
JP2018196007A (ja) * | 2017-05-18 | 2018-12-06 | 日本電気株式会社 | ディジタル変調される信号の復調回路および変調回路 |
CN109959825A (zh) * | 2017-12-26 | 2019-07-02 | 东南大学 | 一种基于拉格朗日插值法的非线性压控衰减器拟合校准方法 |
US11251755B2 (en) | 2019-07-10 | 2022-02-15 | Samsung Electronics Co., Ltd. | Device and method for compensating for nonlinearity of power amplifier |
-
2000
- 2000-03-21 JP JP2000077353A patent/JP2001268150A/ja active Pending
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002344418A (ja) * | 2001-05-22 | 2002-11-29 | Sumitomo Electric Ind Ltd | 前置型非線形歪補償器 |
JP4644974B2 (ja) * | 2001-05-22 | 2011-03-09 | 住友電気工業株式会社 | 前置型非線形歪補償器 |
JP2003174332A (ja) * | 2001-12-05 | 2003-06-20 | Hitachi Ltd | プリディストーション型増幅装置 |
US7161990B2 (en) | 2001-12-05 | 2007-01-09 | Hitachi, Ltd. | Predistortion amplifier |
JP2004015310A (ja) * | 2002-06-05 | 2004-01-15 | Mitsubishi Electric Corp | 歪補償回路および送信機 |
JP2004032609A (ja) * | 2002-06-28 | 2004-01-29 | Nec Corp | 非線形歪み補償回路 |
US7333562B2 (en) | 2002-06-28 | 2008-02-19 | Nec Corporation | Nonlinear distortion compensating circuit |
JP2004088303A (ja) * | 2002-08-26 | 2004-03-18 | Ntt Docomo Inc | 送信機 |
JP2004112151A (ja) * | 2002-09-17 | 2004-04-08 | Hitachi Kokusai Electric Inc | 歪み補償方法及び送信機 |
US7170342B2 (en) | 2002-12-10 | 2007-01-30 | Ntt Docomo, Inc. | Linear power amplification method and linear power amplifier |
JP2004363733A (ja) * | 2003-06-02 | 2004-12-24 | Sumitomo Electric Ind Ltd | 前置型非線形歪補償器 |
JP2005072747A (ja) * | 2003-08-21 | 2005-03-17 | Nippon Telegr & Teleph Corp <Ntt> | 非線形歪補償装置 |
JP2005073032A (ja) * | 2003-08-26 | 2005-03-17 | Hitachi Kokusai Electric Inc | 歪補償増幅装置及び歪補償方法 |
JP2007518378A (ja) * | 2004-01-16 | 2007-07-05 | クゥアルコム・インコーポレイテッド | 増幅器予歪および自動較正の方法および装置 |
US7915954B2 (en) | 2004-01-16 | 2011-03-29 | Qualcomm, Incorporated | Amplifier predistortion and autocalibration method and apparatus |
JP2008503169A (ja) * | 2004-06-30 | 2008-01-31 | インテル・コーポレーション | 周波数領域内で予歪を用いる電力増幅器の線形化方法および装置 |
JP4846715B2 (ja) * | 2004-06-30 | 2011-12-28 | インテル・コーポレーション | 周波数領域内で予歪を用いる電力増幅器の線形化方法および装置 |
US7683631B2 (en) | 2004-09-15 | 2010-03-23 | Advantest Corporation | Inverse characteristic measuring apparatus, distortion compensation apparatus, method, program, and recording medium |
WO2006030919A1 (ja) * | 2004-09-15 | 2006-03-23 | Advantest Corporation | 逆特性測定装置および歪み補償装置、方法、プログラム、記録媒体 |
US7995674B2 (en) | 2006-04-11 | 2011-08-09 | Hitachi, Ltd. | Digital predistortion transmitter |
US7864881B2 (en) | 2006-04-11 | 2011-01-04 | Hitachi, Ltd. | Digital predistortion transmitter |
US8179123B2 (en) | 2007-05-28 | 2012-05-15 | Advantest Corporation | Measurement apparatus and recording medium |
JP2010041470A (ja) * | 2008-08-06 | 2010-02-18 | Mitsubishi Electric Corp | 送信機及び送受信装置 |
US8270530B2 (en) | 2008-09-10 | 2012-09-18 | Fujitsu Limited | Apparatus and method of compensating distortion |
US7920644B2 (en) | 2008-12-26 | 2011-04-05 | Fujitsu Limited | Predistortion apparatus and predistortion method |
EP2202879A1 (en) | 2008-12-26 | 2010-06-30 | Fujitsu Limited | A predistortion apparatus and predistortion method |
WO2010073803A1 (ja) * | 2008-12-27 | 2010-07-01 | 住友電気工業株式会社 | 歪補償回路及び無線基地局 |
JP2010157937A (ja) * | 2008-12-27 | 2010-07-15 | Sumitomo Electric Ind Ltd | 歪補償回路及び無線基地局 |
US8390377B2 (en) | 2008-12-27 | 2013-03-05 | Sumitomo Electric Industries, Ltd. | Distortion compensating circuit and radio base station |
JP2018182895A (ja) * | 2017-04-12 | 2018-11-15 | キヤノン株式会社 | モータ制御装置、方法及びプログラム |
JP2018196007A (ja) * | 2017-05-18 | 2018-12-06 | 日本電気株式会社 | ディジタル変調される信号の復調回路および変調回路 |
CN109959825A (zh) * | 2017-12-26 | 2019-07-02 | 东南大学 | 一种基于拉格朗日插值法的非线性压控衰减器拟合校准方法 |
US11251755B2 (en) | 2019-07-10 | 2022-02-15 | Samsung Electronics Co., Ltd. | Device and method for compensating for nonlinearity of power amplifier |
US11671059B2 (en) | 2019-07-10 | 2023-06-06 | Samsung Electronics Co., Ltd. | Device and method for compensating for nonlinearity of power amplifier |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2001268150A (ja) | リニアライザ | |
Cavers | The effect of quadrature modulator and demodulator errors on adaptive digital predistorters for amplifier linearization | |
US11303251B2 (en) | Digital predistortion adjustment based on determination of load condition characteristics | |
EP0658975B1 (en) | Baseband predistortion system for the adaptive linearization of power amplifiers | |
EP1723725B1 (en) | Digital predistortion system and method for linearizing an rf power amplifier with nonlinear gain characteristics and memory effects | |
US6489846B2 (en) | Distortion compensating device and distortion compensating method | |
JP3699602B2 (ja) | プリディストーション装置及びその方法 | |
FI105366B (fi) | Linearisointimenetelmä ja vahvistinjärjestely | |
EP1047185A2 (en) | On-line calibration of linear amplifiers | |
JP5659557B2 (ja) | 先行歪ませ器、先行歪ませ方法及び先行歪ませシステム | |
KR100429981B1 (ko) | Aqm의 에러보상장치 및 방법 | |
JP3994308B2 (ja) | プリディストーション型歪補償回路 | |
US7848455B2 (en) | Transmission circuit comprising multistage amplifier, and communication device | |
JPH10322137A (ja) | プリディストーション型歪補償回路付送信装置 | |
EP1450482B1 (en) | Circuit and method for compensating for nonlinear distortion of power amplifier | |
JP2006261952A (ja) | 歪補償装置および歪補償係数更新方法 | |
JP2001284980A (ja) | プリディストーション型非線形歪み補償回路およびこれを用いたディジタル送信機 | |
US20220393648A1 (en) | Amplitude and phase control device, amplitude and phase control method, amplification device, and radio transmitter | |
CN101606315B (zh) | 功率放大器时间延迟不变的预失真方法和装置 | |
KR101196584B1 (ko) | 무선 장치, 왜곡 보상 장치 및 왜곡 보상 방법 | |
JP3145015B2 (ja) | 送信装置 | |
US20060083330A1 (en) | Distortion compensation table creation method and distortion compensation method | |
JP2016127577A (ja) | 歪補償装置及び歪補償方法 | |
JP2001268151A (ja) | プリディストーション歪補償装置 | |
JP2003218970A (ja) | 非線形歪み補償装置及びその方法並びにプログラム |