本发明的预失真器可以在任何种类的无线通信中使用,例如,蜂 窝电话,
数字视频广播,数字音频广播,或任何种类的无线通信,例 如,数字用户线(DSL),以增强具有最小非线性失真的高功率放大 器所发送的功率。本发明能够具有在手持无线通信设备和数字卫星通 信中的直接未来应用。
本发明为预失真器。该预失真器是
电子非线性
信号处理设备,它 被放置在依次连接到无线通信系统的传输天线的高功率放大器之前。 高功率放大器的目的是给由该高功率放大器发送至传输天线的 OFDM信号提供尽可能大的功率。然而,功率中大的增加会迫使高功 率放大器中的信号超出该高功率放大器的线性范围。为了实现高功率 放大器输出功率中的这种增加,同时使失真最小化,在放大器之前插 入了预失真器。该预失真器反转放大器的非线性,以便预失真器和高 功率放大器的组合在超出高功率放大器的正常线性范围之外展现线性 特性。这个过程被称为线性化。
所示例发明的特定特征在于,预失真器的设计是基于用于预失真 器输入-输出特性描述的精确解析表达,该输入输出特性是基于用于高 功率放大器的解析模型。这允许了由OFDM
信号传输系统所完成的以 上线性化任务的执行中的准确性和效率。
指导本
申请的基本原理是正交频分复用具有几个所希望的特征, 这使得它成为许多新兴无线通信标准,例如,IEEE802.11a和g WLAM和ETSI地面广播的主要候选。然而,OFDM信号引发的主要 问题之一是其高峰值对平均功率比,这由于高峰值对平均功率比所产 生的非线性失真而严重限制了高功率放大器的效率。
示例
实施例提供了新的混合计算-解析方法,用于在该高功率放 大器为具有随时间变化特性的行波管放大器(TWTA)或固态功率放 大器(SSPA)的情况下,对这种非线性失真的补偿。在数字卫星频道 的情况下,当要求高传输功率时在无线通信系统中使用行波管放大器, 而固态功率放大器用于基于地面的移动无线通信系统。与前面基于查 询表或自适应方案的预失真器技术相比较,该示例实施例依靠与非线 性参数估计算法组合的Saleh行波管放大器模型和Rapp的固态功率 放大器模型的解析反转。这导致预失真器稀疏但准确的表示,具有有 效
跟踪高功率放大器的任何快速随时间变化行为的能力。计算机仿真 结果举例说明并验证了所介绍的方法。
在示例实施例中,我们通过使用Saleh行波管放大器模型和Rapp 的固态功率放大器模型用于这些设备,而对预失真器描述了新的方法 用于高功率放大器,并诉诸依靠仅有几个参数表示的严格封闭式表达 用于其反转。这种方法避免了普通近似表达(类似于多项式近似法) 会要求精确表示的大量参数。
在示例方法中,我们利用了用于固态功率放大器和行波管放大器 的解析模型,以给分别标注为预失真器I和预失真器II的两个预失真 器推到出令人信服的算法。预失真器I算法适用于固态功率放大器, 而预失真器I算法适用于行波管放大器。
我们使用这两种类型的高功率放大器的理由是,这两种类型对于 今天的无线通信系统非常重要。行波管放大器通常用于卫星通信,而 固态功率放大器用于移动通信系统。因为行波管放大器严重的非线性, 已经在失真补偿方面对这种类型放大器做了相当可观的工作。然而, 希望OFDM以与码分多址(CDMA)组合的形式,即,多载波码分 多址(MC-CDMA)或多载波直接序列码分多址(MC-DS-CDMA), 成为用于下一代蜂窝系统的标准。码分多址是一种使用扩展
频谱技术 的数字蜂窝技术。与竞争系统不同,CDMA不为每个用户分配
指定的 频率。而是,每个信道都使用完整的可用频谱。单独会话使用伪随机 数字序列编码。CDMA一贯提供比其他的商业移动技术更好的性能用 于语音和数据通信,允许更多的用户在任何给定时间连接。多载波 (MC)CDMA是直接序列(DS)CDMA(码分多址)与OFDM技 术的组合技术。它在频域内使用了扩频序列。
因此,固态功率放大器的重要性就会比现在大得多。由于这个原 因,我们也使用固态功率放大器作为高功率放大器。虽然用于Saleh 模型反转的封闭式表达已知,在其中高功率放大器特性是随时间变化 的示例实施例中,这种反转没有在其预失真器的实现中使用。我们将 用于高功率放大器特性反转的封闭式表达与连续非线性参数估计算法 相结合,它允许预失真器的稀疏实施和对高功率放大器随时间变化行 为的准确跟踪或者适应。
与以上提及的其它
现有技术方法比较,通过以下描述的计算机仿 真示范并验证,我们的算法是快速的,准确的,并且具有低复杂度。
虽然已经或者将为具有功能性解释的合乎语法的变通性而描述 该装置与方法,应当清楚理解,除非是在35USC112下明确阐述,权 利要求书不被理解为必须以任何方式由“方法”或“步骤”限制的结构所 限制,而要符合在等同物的公正宗旨之下根据由
权利要求书所提供定 义的方法或等价物的整个范围,并且在权利要求书在35USC112下明 确阐述的情况下,要符合35USC112之下的完整法定等价物。本发明 可以通过现在求助于以下图示而更好地形象化,其中相同的单元用相 同的
附图标记表示。
图1为具有本发明的预失真器和高功率放大器的简化OFDM通 信发射器。
图2为显示作为归一化输入的函数的归一化输出的Saleh的行波 管放大器模型的非线性幅度和
相位转移函数的图示。
图3为显示作为归一化输入的函数的归一化输出的Rapp的固态 功率放大器模型的非线性幅度转移函数的图示。
图4为显示作为归一化输入的函数的归一化输出的具有预失真器 的Saleh的行波管放大器模型的幅度补偿效应图示。
图5为与随时间变化的高功率放大器结合的预失真器的简化框 图。
图6a为显示作为归一化输入的函数的归一化输出的使用预失真 器的Rapp的固态功率放大器模型的补偿效应图示。
图6b为显示作为归一化输入的函数的归一化输出的使用预失真 器的Rapp的固态功率放大器模型的补偿与削波效应图示。
图7a为显示I信道对比Q信道,有行波管放大器,无预失真器 的所接收OFDM信号构象的图示。
图7b为显示I信道对比Q信道,有行波管放大器,有预失真器 的所接收OFDM信号构象的图示。
图8为显示在具有随时间变化的行波管放大器的OFDM系统中, 有和没有预失真器的误比特率(BER)输出性能的图示,它显示BER 作为以db为单位的输入Eb/N0比率的函数,其中Eb是每比特的信号
能量,而N0是噪声
功率谱密度。那就是Eb/N0=SNR(
信噪比)。
图9a为显示作为归一化输入的函数的归一化输出的饱和条件下 的信号幅度图示,其中超过1的归一化信号被削波。
图9b为饱和条件下信号相位的图示。这个图显示了归一化输入 幅度对比输出相位失真,因为输出相位失真是归一化输入幅度的函数。
图10为显示具有在参数均匀分布,IBO(输入补偿)=6dB的随 时间变化的行波管放大器的OFDM系统中,有和无预失真器的BER 输出性能的图示,其中给和不给预失真器提供跟踪,图示显示BER作 为以db为单位的输入Eb/N0比率的函数,其中Eb是每比特的信号能 量,而N0是噪声功率谱密度。那就是Eb/N0=SNR(信噪比)。
图11为显示具有在参数均匀分布,IBO=7dB的随时间变化的行 波管放大器的OFDM系统中,有和无预失真器的BER输出性能的图 示,其中给和不给预失真器提供跟踪,图示显示BER作为以db为单 位的输入Eb/N0比率的函数,其中Eb是每比特的信号能量,而N0是 噪声功率谱密度。那就是Eb/N0=SNR(信噪比)。
图12a为显示I信道对比Q信道,有固态功率放大器,无预失真 器的所接收OFDM信号构象的图示。
图12b为显示I信道对比Q信道,有固态功率放大器,有预失真 器的所接收OFDM信号构象的图示。
图13为具有不随时间变化的固态功率放大器的OFDM系统中预 失真器的BER性能的图示,当A0=p=1时显示BER作为以db为单位 的输入Eb/N0的函数,其中Eb是每比特的信号能量,而N0是噪声功 率谱密度。那就是Eb/N0=SNR(信噪比)。
图14为预失真器的BER性能的图示,当参数在1≤A0≤1.5,1≤p≤1.5 的范围内均匀分布,IBO=6dB时,显示BER作为输入Et/N0比率(以 db为单位)的函数,其中Eb为比特差错数而N0为输入比特的总数。
图15为预失真器的BER性能的图示,当参数在1≤A0≤2,1≤p≤2 的范围内均匀分布,IBO=6dB时,显示BER作为输入Eb/N0的函数, 其中Eb是每比特的信号能量,而N0是噪声功率谱密度。那就是 Eb/N0=SNR(信噪比)。
图16为预失真器的BER性能的图示,当参数在1≤A0≤2,1≤p≤2 的范围内均匀分布,IBO=7dB时,显示BER作为输入Eb/N0的函数, 其中Eb是每比特的信号能量,而N0是噪声功率谱密度。那就是 Eb/N0=SNR(信噪比)。
图17显示了Saleh的TWTA模型中具有高斯和均匀分布的两个 变化参数,β,ε的收敛情况。
图18为显示具有在参数为高斯和均匀两种分布,IBO(输入补 偿)=6dB的随时间变化的行波管放大器的OFDM系统中,有和无预 失真器的BER输出性能的图示,其中给和不给预失真器提供跟踪,图 示显示BER作为以db为单位的输入Eb/N0比率的函数,其中Eb是每 比特的信号能量,而N0是噪声功率谱密度。那就是Eb/N0=SNR(信 噪比)。
图19为显示具有在参数为高斯和均匀两种分布,IBO(输入补 偿)=7dB的随时间变化的行波管放大器的OFDM系统中,有和无预 失真器的BER输出性能的图示,其中给和不给预失真器提供跟踪,图 示显示BER作为以db为单位的输入Eb/N0比率的函数,其中Eb是每 比特的信号能量,而N0是噪声功率谱密度。那就是Eb/N0=SNR(信 噪比)。
图20显示了Rapp的SSPA模型中具有高斯分布的A0,p的两个 变化参数的收敛情况。
图21为当参数为高斯分布时预失真器的BER性能图示,方差 =0.1,IBO=6dB,图示显示BER作为以db为单位的输入Eb/N0比率 的函数,其中Eb是每比特的信号能量,而N0是噪声功率谱密度。那 就是Eb/N0=SNR(信噪比)。
图22为当参数为高斯分布时预失真器的BER性能图示,方差 =0.1,IBO=7dB,图示显示BER作为以db为单位的输入Eb/N0比率 的函数,其中Eb是每比特的信号能量,而N0是噪声功率谱密度。那 就是Eb/N0=SNR(信噪比)。
本发明及其各种实施例现在可以通过求助于以下优选实施例的 详细描述而得到更好的理解,优选实施例作为在权利要求书中定义的 本发明示例实例提出。应当清楚的理解,如权利要求书所定义的本发 明可以比以下所描述的示例实施例更为广泛。
系统描述
图1为显示系统结构的本发明的简化
框图,总体表示为附图标记 10,用于OFDM系统高功率放大器非线性的补偿。OFDM基带模
块 12生成OFDM格式的信号到预失真器14,其数字输出被
数模转换器 16转换为模拟格式以产生
相移QAM输出到乘法器18和20,它们在 加法器22处组合并相加,并且输入到功率放大器24用于发送到无线 或有线通信系统。必须理解,图1中的
硬件可以按照很多种等价方式 实现。例如,预失真器14是数字
电路,它可以是使用硬件和/或软硬 件的组合的专用
数字信号处理器,或者可以是具有适合信号
接口的计 算机,它由
软件进行计算机排列或配置,以如本发明教导的那样处理 数字信息。在可以实现预失真器14的指定技术方面没有限制,而且明 确预期所有现在已知的或以后设计的方法在本发明的范围之内。
通常,OFDM信号x(t)可以如下解析表示:
x ( t ) = 1 N Σ k = 0 N - 1 X [ k ] e j 2 π f k t - - - ( 1 ) 其中X[k]代表正交幅度调制(QAM)符号,N为
子载波的数量, 而fk为第k个子
载波频率,它可以如下表示:
f k = k · 1 N T s - - - ( 2 ) 其中,Ts为x(t)的
采样周期。QAM是将两种幅度调制(AM)信 号组合到信号信道中的方法,从而使有效带宽加倍。QAM在数字系 统中,尤其在无线应用中用于脉冲幅度调制(PAM)。在QAM信号 中具有两种载波,每种具有相同频率但是相位相差90度(四分之一周 期,从其中出现求
积分项)。一种信号被称为I信号,而另外一种被 称为Q信号。从数学上说,信号中的一个可以用
正弦波表示,而另外 一个用余弦波表示。这两个调制载波在发送的源头处被组合起来。在 目的地处,该载波分离开,数据从彼此中抽出,并且随后该数据被组 合为原始调制信息。
通过在t=nTs处将x(t)离散化,我们有以下公式:
x ( n ) ≡ x ( n T s ) = 1 N Σ k = 0 N - 1 X [ k ] e j 2 πkn N - - - ( 3 ) 预失真器14是预先计算并消除跟随预失真器14之后的零记忆高 功率放大器24中出现的非线性失真的非线性零记忆设备。
行波管放大器模型
作为高功率放大器模型,我们展示了Saleh的建立得很好的行波 管放大器模型。在这个模型中,行波管放大器的AM/FM和AM/PM 转换可以如下表示:
u [ r ] = αr 1 + β r 2 - - - ( 4 ) Φ [ r ] = γ r 2 1 + ϵ r 2 - - - ( 5 ) 其中,u为幅度响应,Φ为相位响应,r为行波管放大器的输入 幅度,而α,β,γ和ε是四个可调节参数。公式(4)和(5)的行为 在图2的图示中举例说明,其中行波管放大器的归一化输出作为归一 化输入的函数显示。在图2中,我们使用α=1.9638;β=0.9945;γ=2.5293; 和ε=2.8168,如同在Saleh的原著中一样。无预失真器14的行波管放 大器24的输出z(t)可以如下表示:
z(t)=u[r]cos(ωct+φ(t)+Φ[r])
(6)
其中(t)为
输入信号的相位,而ωc为载波频率。
固态功率放大器模型
对于固态功率放大器模型24,我们使用归一化的Rapp的模型。 在这个模型中,我们假设AM/PM转换足够小,以便它可以被忽略。 接着,固态功率放大器的AM/AM和AM/PM转换可以如下表示:
u [ r ] = r ( 1 + ( r A 0 ) 2 p ) 1 2 p Φ[r]≈0 (7),(8)
其中,r为固态功率放大器24的输入幅度,A0为最大输出幅度, 而p为影响该转移的平滑度的参数。公式(7)的行为在图3的图示中 举例说明,其中归一化输出作为归一化输入的函数显示。无预失真器 14的固态功率放大器24的输出z(t)可以如下表示:
z(t)=u[r]cos(ωct+φ(t))
(9)
其中(t)为输入信号的相位。
预失真器
现在根据本发明对行波管放大器24和固态功率放大器24二者都 考虑预失真器14。使q和u分别代表预失真器14和高功率放大器24 的非线性零记忆输入和输出映射,而xl(n)代表预失真器14的输入, yl(n)代表预失真器14的输出,它也是到高功率放大器24的输入,而 zl(n)代表如图1中所示的高功率放大器24的输出。那么对任何给定的 高功率放大器24,根据本发明的理想预失真器14是输入输出映射满 足以下公式的一个:
u[q(xl(n))]=kxl(n) (10)
其中,k为期望的预先指定的线性放大常数。在这个例子中,我 们假设k=1。
用于行波管放大器的预失真器
不随时间变化的情况
在行波管放大器24中,用于预失真器14的输入xl(n)和输出yl(n) 的通用基带(等效的低通信号)表达式为:
xl(n)=r(n)ejφ(n)
yl(n)=q[r(n)]ej(φ(n)+0[r(n)])
(11),(12)
其中函数q和要通过要求满足公式(10)确定。根据公式(4) 和(5),行波管放大器24的输入和输出为:
y(t)=q[r(t)]cos(ωct+φ(t)+θ[r(t)])
z(t)=u[q[r(t))]cos(ωct+φ(t)+θ[r(t)]+φ{q(t)])
(13),(14)
其中
u [ q ( r ) ] = αq 1 + β q 2 Φ [ q ( r ) ] = γ q 2 1 + ϵ q 2 (15),(16)
为了满足(10),必须具备以下条件:
αq 1 + β q 2 = r γ q 2 1 + ϵ q 2 = - θ (17),(18)
从公式(17)
rβq2-αq+r=0
(19)
可以求解这个公式得到q以产生:
q ( r ) = α - α 2 - 4 r 2 β 2 rβ , r ≤ 1 - - - ( 20 ) 同时对零相位失真,我们必须要:
0(r)+Φ(q)=0
(21)
或
θ ( r ) = - Φ ( q ) = - γ ( q ( r ) ) 2 1 + ϵ ( q ( r ) ) 2 - - - ( 22 ) 如果r>1,公式(20)无解。根据图4的描述,这对应于信号的 削波,图4中归一化输出显示为对于具有预失真器14的行波管放大器 24的归一化输入的函数。公式(20),(22)的解析解以前已经由 Brajal和Chouly获得。
随时间变化的自适应情况
我们如下将这种解法扩展到随时间变化的情况。作为随时间变化 的模型,我们假设四个参数α,β,γ和ε随时间变化。我们表达:
J ( α , β ) = E ( αq 1 + β q 2 - u ) 2 - - - ( 23 ) 其中J为应当被最小化的成本函数,E为关于α,β的期望值。 对于α部分求微分并使结果等于0,我们得到:
∂ J ( α , β ) ∂ α = E [ 2 ( αq 1 + β q 2 - u ) q 1 + + β q 2 ] = 0 αE ( q 2 ( 1 + β q 2 ) 2 ) = E ( qu 1 + β q 2 ) (24),(25)
对于β与之类似继续进行,我们得到:
∂ J ( α , β ) ∂ β = E [ 2 ( αq 1 + β q 2 - u ) ( - αq ( 1 + β q 2 ) 2 ) q 2 ] = 0 - - - ( 26 ) 或
αE ( q 4 ( 1 + β q 2 ) 3 ) = E ( q 3 u ( 1 + β q 2 ) 2 ) - - - ( 27 ) 为简单起见,让我们定义以下公式:
A ( β ) = E ( q 2 ( 1 + β q 2 ) 2 ) B ( β ) = E ( qu 1 + β q 2 ) C ( β ) = E ( q 4 ( 1 + β q 2 ) 3 ) D ( β ) = E ( q 3 u ( 1 + β q 2 ) 2 ) (28),(29),(30),(31)
根据公式(25),(28)和(29)
α = B ( β ) A ( β ) - - - ( 32 ) 而根据公式(27),(30),(31),(32)
B ( β ) A ( β ) C ( β ) = D ( β ) - - - ( 33 ) 这样,我们的方法为:在图5中所示的估计器26中数值解算公 式(33)以求解它是β的估计值,并且随后在公式(32)中代替 以获得α的估计值公式(28),(29),(30),(31)中的期 望值可以使用以下公式估计:
A ^ ( β ) = 1 N Σ n = 1 N q n 2 ( 1 + β q n 2 ) 2 B ^ ( β ) = 1 N Σ n = 1 N q n u n 1 + β q n 2 C ^ ( β ) = 1 N Σ n = 1 N q n 4 ( 1 + β q n 2 ) 3 D ^ ( β ) = 1 N Σ n = 1 N q n 3 u n ( 1 + β q n 2 ) 2 (34),(35),(36),(37)
γ和ε也可以按照如上所述的相同方式准确估计。这种方法在图 5的框图中举例说明,图5显示了用于随时间变化的高功率放大器的 预失真器14,其中提供参数估计器26以从高功率放大器24取得参数, 并将它们提供给估计器26以产生参数估计用于预失真器14。
为了从(33)获取β的最优估计值,我们使用以下公式:
β ^ opt = min β | B ( β ) C ( β ) - A ( β ) D ( β ) | 2 - - - ( 38 ) 为了最小化由以下公式定义的均方差(MSE),定义了满足(38) 的最优系数
J ( β ) = E [ B ^ ( β ) C ^ ( β ) - A ^ ( β ) D ^ ( β ) ] 2 - - - ( 39 ) 其中J为要被最小化的成本函数,而E为关于β的期望值。
然后,J对于β求导
∂ J ( β ) ∂ β = 2 E ( B ^ ( β ) C ^ ( β ) - A ^ ( β ) D ^ ( β ) ) · ( ∂ B ^ ( β ) ∂ β C ^ ( β ) + B ^ ( β ) ∂ C ^ ( β ) ∂ β - ∂ A ^ ( β ) ∂ β D ^ ( β ) - A ^ ( β ) ∂ D ^ ( β ) ∂ β ) - - - ( 40 ) 其中
∂ A ^ ( β ) ∂ β = - 2 N Σ n = 1 N q n 4 ( 1 + β q n 2 ) 3 ∂ B ^ ( β ) ∂ β = - 1 N Σ n = 1 N q n 3 u n ( 1 + β q n 2 ) 2 ∂ C ^ ( β ) ∂ β = - 3 N Σ n = 1 N q n 6 ( 1 + β q n 2 ) 4 ∂ D ^ ( β ) ∂ β = - 2 N Σ n = 1 N q n 3 u n ( 1 + β q n 2 ) 3 (41)(42)(43)(44)
此后,LMS(最小均方)算法可以如下表示:
β ^ ( n + 1 ) = β ^ ( n ) - μ β ^ · ( B ^ ( β ^ ( n ) ) C ^ ( β ^ ( n ) ) - A ^ ( β ^ ( n ) ) D ^ ( β ^ ( n ) ) ) ) · ( ∂ B ^ ( β ^ ( n ) ) ∂ β ^ ( n ) C ^ ( β ^ ( n ) ) + B ^ ( β ^ ( n ) ) ∂ C ^ ( β ^ ( n ) ) ∂ β ^ ( n ) - ∂ A ^ ( β ^ ( n ) ) ∂ β ^ ( n ) D ^ ( β ^ ( n ) ) - A ^ ( β ^ ( n ) ) ∂ D ^ ( β ^ ( n ) ) ∂ β ^ ( n ) ) - - - ( 45 ) 其中为LMS算法的步长。
一旦我们得到β的估计值,我们就很容易从(32)得到α的估计 值。γ和ε可以用如上所述的相同方式获得。
用于固态功率放大器的预失真器
不随时间变化的情况
正如在行波管放大器24中那样,用于固态功率放大器24的预失 真器14的输入xl(n)和输出yl(n)的通用基带(等同于低通信号)表达 式为:
x l ( n ) = r ( n ) e jφ ( n ) y l ( n ) = q [ r ( n ) ] e jφ ( n ) (46),(47)
其中函数q和要通过要求满足公式(10)确定。由于我们假设 忽略相位失真,我们不必考虑相位前置补偿。根据公式(7)和(8), 固态功率放大器24的输入和输出为:
yc(t)=q[r(t)]cos(ωct+φ(t))
z(t)=u[q[r(t)]]cos(ωct+φ(t))
(48),(49)
其中
u [ q ( r ) ] = q ( r ) ( 1 + ( q ( r ) A 0 ) 2 p ) 1 2 p - - - ( 50 ) 根据公式(50),公式(10)隐含:
q ( r ) ( 1 + ( q ( r ) A 0 ) 2 p ) 1 2 p = r - - - ( 51 ) 因此,在一些代数处理之后,我们能够找到用于预失真器特性q(r) 的准确表达:
q ( r ) = r ( 1 - ( r A 0 ) 2 p ) 1 2 p , r < A 0 - - - ( 52 ) 补偿效应的示例在图6中展示。当r>A0时,公式(52)无解。 在这种情况下,我们如图6中那样削减输入信号。
随时间变化可调节的情况
由于高功率放大器24是随时间变化的系统,作为随时间变化的 模型,我们假设固态功率放大器模型中的参数A0和p随时间变化。为 了跟踪两个参数A0和p,我们使用了训练符号。通过使用训练符号, 我们获得预失真器14的输入q(n)和预失真器14的输出u(n)。在训练 阶段内,我们假设预失真器14是关闭的。那就是说,预失真器14的 输入和输出会相同(r(n)=q(n))。
为了估计参数A0和p,首先,我们如下改变公式(50):
A 0 = q · u ( q 2 p - u 2 p ) 1 2 p - - - ( 53 ) 为了归纳算法,如果我们知道p的话,我们可以轻易地从公式(53) 获得A0。然而,我们假设A0和p二者都随时间变化。首先,发送两 个训练符号,然后我们知道高功率放大器24的输入幅度q和输出幅度 u。随后从公式(53),对应两个不同的训练符号,我们可以获得A0 的两个不同估计值,如以下公式(54)和(55)所给出的名为A01和 A02。如果我们选择了正确的p,它在训练时间内对高功率放大器24 相同,即A01和A02的A0的两个不同值具有几乎相同的值,或者由于 步长,具有非常接近的值。我们可以找到用于这一点的p,它具有两 个估计值A0之间的最小距离,即Dmin=|A01-A02|2。那么,从公式(53) 和p的估计值,我们从最小距离Dmin=|A01-A02|2得到
A ^ 0 = A 01 ≈ A 02 . 这个 算法在计算方面非常简单。我们使用仅两个训练符号并且没有反复, 因此导致非常小的延迟。
该算法的简要描述
作为一种更实用的方式,如果我们知道了p,我们可以很容易从 公式(53)得到A0。然而,我们假设A0和p二者都随时间变化。在 这种情况下,我们建议遵循算法。首先,发送两个训练符号,然后我 们知道高功率放大器24的输入幅度q和高功率放大器的输出幅度u。 此后,从公式(53),对应于两个不同训练符号,我们得到A0的两个 不同估计值A01和A02。
A 01 = q 1 · u 1 ( q 1 2 p - u 1 2 p ) 1 2 p A 02 = q 2 · u 2 ( q 2 2 p - q 2 2 p ) 1 2 p (54),(55)
其中,q1,u1分别为用于第一训练符号的预失真器14和高功率 放大器24的输出幅度,而q2,u2分别为用于第二训练符号的预失真 器14和高功率放大器24的输出幅度。训练符号不受如我们前面所声 明的预失真器14的函数的影响。在训练期间内,我们将q1和q2替代 为r1和r2,这是训练符号的原始幅度。我们可以使用以下公式估计未 知的A0和p。
p ^ opt = min p | A 01 ( p ) - A 02 ( p ) | 2 A ^ 0 = A 01 ( p ^ opt ) ≈ A 02 ( p ^ opt ) (56),(57)
其中为A0的估计值,而为我们可以从公式(56)得到的的 最优值。
仿真结果与讨论
现在考虑示例预失真器技术的测试用于补偿如计算机仿真所证 明的高功率放大器非线性失真。假设加性高斯白噪声(AWGN)信道 明显观察到非线性和由所示例预失真器14所做的性能提升的效果。考 虑具有128个子载波和16个QAM的OFDM系统。如果输入幅度非 常高的话,高功率放大器24以高非线性状态运行。如果输入幅度非常 小的话,高功率放大器24以非常小的失真运行。在高功率放大器24 的操作中,需要相对的功率补偿
水平以减少失真。然而,并不希望这 种功率补偿,因为它减少了功效。在我们的算法中,在范围r<A0内总 是存在补偿解决方案,其中A0为最大输出幅度。因此,如果输入平均 功率与A0 2相同的话,我们得到最大的功效,但是得到高度非线性结果。 从而,我们需要一种标准以显示多么需要来自最优功效的功率补偿。 在仿真中,我们定义IBO(输入补偿)为:
IBO = 10 log 10 ( A 0 2 P in ) - - - ( 58 ) 其中Pin为输入平均功率(OFDM信号的平均功率)。与之类似, 我们也可以定义OBO(输出补偿)为:
OBO = 10 log 10 ( A 0 2 P out ) - - - ( 59 ) 其中Pout为输出平均功率(高功率放大器24的平均输出功率)。
用于行波管放大器的预失真器
不随时间变化的情况
现在假设参数α,β,γ和ε不随时间变化的条件下考虑OFDM 仿真结果。图7a和7b为将α作为I的函数描述,并且显示分别没有 和有预失真器14的信号构象之间区别的图示。在图7a和7b中,我们 使用IBO=6dB。在图8的图示中显示的误码比率或误码率(BER)性 能曲线显示BER作为Eb/N0的函数,其中Eb为每比特的信号能量, 而N0为噪声功率谱密度,并且显示了预失真器14能够明显地降低 OFDM系统10中的非线性失真。BER是用某个规定时间内的错误比 特数除以所传输,接收,或处理比特的总数。误码率的例子为(a)传 输BER,即,用错误比特数除以接收到的传输比特的总数;(b)信 息BER,即,用错误解码(纠正)比特的数除以解码(纠正)比特的 总数。BER通常作为系数和10的幂表达;例如,100,000个所发送比 特中的2.5错误比特会是105分之2.5或2.5×10-5。
均匀分布的随时间变化的自适应情况
如上所述,高功率放大器24是随时间变化的系统。假设四个参 数α,β,γ和ε现在是随时间变化的,那么我们应该跟踪α,β,γ和 ε的变化。我们假设这四个参数根据以下条件平均分布变化:
(1)这四个参数在以下范围内变化
1.01≤α≤2 (60)
0.01≤β≤1 (61)
1.5≤γ,ε·≤ 3 (62)
(2)输入和输出归一化条件,β=α-1。
(3)
饱和度条件,信号超过1就被削波,如图9a和9b的图示 中所示。
我们为什么在幅度和相位方面选择以上条件的原因是,要在输入 和输出两个方面维持归一化约束和以上范围(r>A0)内的约束条件, 即使是幅度被改变了。这些约束仅是为了表达方便,因此在真实系统 中,即使没有满足以上条件,我们的算法也运行得很好。下面的表1 显示了使用我们的算法跟踪α,β,γ和ε之后的错误。我们使用了以 下公式获得表1的结果。
Error ( α ) = | α - α ^ | | α max - α min | Error ( β ) = | β - β ^ | | β max - β min | Error ( γ ) = | γ - γ ^ | | γ max - γ min | Error ( ϵ ) = | ϵ - ϵ ^ | | ϵ max - ϵ min | (63),(64),(65),(66)
我们仅使用两个训练符号、计算1000次并将结果平均得到了表1 的结果。
表1参数的错误
步长 错误(α) 错误(β) 错误(γ) 错误(ε) 0.1 1.02×10-2 2.74×10-2 6.3×10-3 1.72×10-2 0.01 9.47×10-4 2.5×10-3 6.04×10-4 1.7×10-3 0.001 9.49×10-5 2.54×10-4 6.18×10-5 1.69×10-4
表1的结果显示仅有两个训练符号对我们的算法就足够了。这表 示我们的算法非常快,并且几乎没有延迟。在具有随时间变化的高功 率放大器24的OFDM 10中预失真器14的BER性能在图10和图11 的图示所示。在这些曲线中,我们假设步长=0.01。正如从图10和图 11中清楚地看到的那样,如果高功率放大器24的变化没有被跟踪的 话,与跟踪的情况比较性能就会差很多。这个仿真结果因此显示了跟 踪参数中变化的这种能力给系统性能增值。
按照高斯分布和LMS算法的随时间变化可调节的情况
我们再次仿真了我们的PD,但是参数分布不同。我们假设4个 参数α,β,γ,ε是高斯和均匀分布地随时间变化,并且使用LMS(最 小均方)算法跟踪参数的变化。首先我们在图17中显示我们算法的收 敛情况。我们仅显示两个参数β和ε的原因是,正如我们以前所显示 的,一旦我们得到了β和ε二者,其他的参数α和γ就可以很容易得 到。在这次仿真中,我们假设β是均匀分布的,而ε是如在Saleh的 原始模型中具有均值E(ε)=2.8168的高斯分布和方差0.01。为了快速 收敛我们使用步长μβ=6000000和με=600000000000。
现在我们显示BER性能在具有和没有跟踪之间的比较。在这些 仿真中,我们假设四个参数根据以下条件变化。
(1)这两个参数在以下范围内变化
1.01≤α≤2 (67)
0.01≤β≤1 (68)
(2)相位参数γ和ε每个都按照平均值E(γ)=2.5293,E(ε)=2.8168 和方差σ=0.1的高斯分布变化。
(3)输入和输出归一化条件,β=α-1。
(4)饱和度条件,信号超过1就被削波,如图9a和9b的图示 中所示。
正如在前面部分中我们所解释的,这些约束仅是为了表达方便。 具有随时间变化HPA的OFDM中PD的BER性能在图18(IBO=6dB) 和图19(IBO=6dB)中显示。在这些BER性能仿真中,我们假设步 长μβ=50000000和με=10000000000。我们使用两个训练符号并
迭代 1000次。甚至通常PD需要少得多的迭代,我们使用足够多的迭代以 确保所有的参数都收敛。正如从图18和图19中所显见的那样,如果 HPA的变化没有被跟踪的话,与跟踪的情况比较性能就差得多。仿真 结果因此显示跟踪参数中变化的能力给系统性能增值。
用于固态功率放大器的预失真器
不随时间变化的情况
假设固态功率放大器24为不随时间变化的系统而考虑OFDM仿 真结果。在这次仿真中,使用了16QAM作为调制方案并且使用了128 个子载波。因为高的峰值对平均功率比,OFDM比单载波系统需要多 得多的IBO。图12a和12b分别显示了有和无预失真器14的信号构 象输出。与行波管放大器的情况相比较,幅度失真不是那么严重并且 没有相位失真存在。然而,没有预失真器14,即使IBO=6dB,幅度 失真也非常高。在图13中,BER性能曲线显示我们的预失真器14能 够明显降低OFDM系统10中非线性失真的效应。在图13中,我们 使用A0=p=1。
均匀分布的随时间变化可调节的情况
如我们在前面所提及,高功率放大器14是随时间变化的系统。 假设两个参数A0和p是随时间变化的,这样我们应当跟踪A0和p的 变化。正如在行波管放大器24中的情况,两个参数A0和p均匀分布。 仿真使用了简单的搜寻算法。表2显示了在使用我们的算法跟踪A0 和p之后的错误。我们使用了以下公式得到表2的结果。
Error ( A 0 ) = | A 0 - A ^ 0 | | A max - A min | Error ( p ) = | p - p ^ | | p max - p min | (69),(70)
其中和是使用简单搜寻算法跟踪的参数,而|Amax-Amin|和 |pmax-pmin|为变化范围。我们计算公式(69)和(70)1000次并平均每 次的错误。根据表2,即使步长为0.1,错误也是非常小的。
表2和的错误
步长 1≤A0,p≤1.5 1≤A0,p≤2 1≤A0,p≤3 错误(A0) 错误(p) 错误(A0) 错误(p) 错误(A0) 错误(p) 0.1 3.86×10-2 5.1×10-2 2.29×10-2 2.56×10-2 1.40×10-2 1.23×10-2 0.01 3.7×10-3 5.0×10-3 2.22×10-3 2.5×10-3 1.5×10-3 1.3×10-3 0.001 3.66×10-4 4.8875×10-4 2.1718×10-4 2.4870×10-4 1.4934×10-4 1.2852×10-4
我们现在了解了用于随时间变化的固态功率放大器24的预失真 器14的BER性能。我们在下面BER性能仿真中使用步长0.01。在 图14中,我们假设两个参数在范围1≤A0,p≤1.5内均匀分布,每个具有 平均值=1.25,IBO=6dB。在没有跟踪的情况下,我们对两个参数使 用平均值1.25。在图15和图16中,我们显示了当两个参数在宽度范 围1≤A0,p≤2内平均分布,每个具有平均值=1.5,IBO=6dB和7dB时, 用于随时间变化的固态功率放大器24的预失真器14的BER性能。在 没有跟踪的情况下,我们对两个参数使用平均值1.5。
按照高斯分布随时间变化可调节的情况
现在,我们假设两个参数A0和p是按照高斯分布随时间变化, 并使用LMS算法跟踪该变化。首先,我们在图20中仿真我们算法的 收敛。在这次仿真中,我们假设两个参数A0和p按照高斯分布(均值 E(A0)=1.5,E(p)=1.5,方差
σ A 0 = 0.01 , σp=0.01)连续变化。为了快 速收敛,我们使用步长
μ p ^ ( n ) = 10000 . 作为MSE(均方差),我们每次 计算错误100次并将其平均。由于A0的MSE依赖于p的MSE,它们 的MSE显示出类似的特性。在图21(IBO=6dB)和图22(IBO=7dB) 中,我们比较了跟踪参数p和A0变化的情况和没有跟踪参数p和A0 变化的情况。在这些仿真中,我们假设两个参数p和A0按照方差0.1 的高斯分布。由于在真实系统中,HPA的特征不是变化得那么快,我 们假设两个参数p和A0变化每768符号变化,并且我们知道参数何时 可能会变化。如果参数变化更快的话,那么我们只要减少训练阶段的 时间以及时跟踪两个参数的变化。为了快速收敛,我们使用步长
μ p ^ ( n ) = 5000 . 在没有跟踪的情况下,我们使用两个参数p和A0的平均值, 每个为1.5。另外一件事,我们应当提到的是关于选择训练符号,我们 应当从HPA函数中的足够非线性
位置选择符号。如果输入非常小的 话,HPA运算非常接近线性情况。那就是说,这种情况下输入=输出。 那么从公式(53),A0趋于无穷大,而且我们不能够找到两个参数p 和A0。然而,HPA总是具有非线性区域。(如果它不具有非线性部 分的话,我们就不需要使用预失真器)。我们总是能够找到两个适当 参数p和A0。
现在可以理解基于OFDM的无线通信10中使用的随时间变化的 高功率放大器24中用于消除或减轻非线性失真的上述基于模型的预 失真器方法的优点了。该方法使用行波管放大器的Saleh模型和固态 功率放大器的Rapp的模型的封闭式反转,在反转的表达中要求非常 少的参数。这种稀疏但是准确的表示实现了对高功率放大器24的随时 间变化行为的快速跟踪。这些特性已经由简单的计算机仿真所证明。
在不背离本发明的思想和范围下可以由本领域的普通技术人员 做出很多变更与修正。因此,必须理解,所示例实施例仅出于举例的 目的而提出,并且不应当被当作如以下发明及其各种实施例所定义的 那样限制本发明。
因此,必须理解,所示例实施例仅出于举例的目的而提出,并且 不应当当作限制如以下权利要求书所定义的本发明。例如,尽管权利 要求书的元素以下以某种组合提出的事实,必须清楚地理解,本发明 包括更少,更多或不同元素的其他组合,即使是当它们没有一开始就 以这种组合
声明时也在以上中公开。还要将两个元素以所声明组合方 式组合的教导也理解成允许所声明的组合,其中两个元素相互之间不 组合,但是可以单独使用或在其他组合中组合使用。本发明任何公开 元素的删除明确期望在本发明的范围之内。
在要描述本发明及其各种实施例的
说明书中使用的词不仅要从 其通常定义的含义的意义上理解,还要包括超出通常定义的含义范围 的在该说明书中的结构,材料或行为中的特别定义。这样如果元素能 够在该说明书的上下文中理解为包括多于一种含义,那么其在权利要 求书中的用途必须要被理解为对由该说明书和由该词本身所支持的所 有可能含义。
因此,以下权利要求书的词或元素的定义就在该说明书中定义, 以不仅包括按照字面意思提出的元素的组合,还包括以基本相同的方 式执行基本相同功能以获得基本相同结果的所有等价结构,材料或行 为。在此意义上,因此期望对以下权利要求书中的任何一个元素可以 进行两个或更多元素的等价替代,或者权利要求书中的单个元素可以 替代为两个或更多元素。尽管元素可以在以上描述为以某种组合中运 行并且甚至最初是这样声明的,应当清楚地理解,来自所声明组合的 一个或更多元素在一些情况下可以从该组合中拆出,并且所声明的组 合可以针对子组合或子组合的变更。
明确期望现在已知的或以后设计的,本领域的普通技术人员所观 察到的,来自所声明主题的非实质性变化等价地在权利要求书的范围 之内。因此,现在或以后对本领域普通技术人员已知的各种替代被定 为在所定义元素的范围之内。
这样权利要求书被理解为包括以上特别例举和描述的内容,概念 上等价的内容,能够明显替代的内容和基本包含本发明的基本观点的 内容。
相关申请
本申请涉及2004年8月19日申请的美国临时
专利申请 No.60/602,905,它在此作为参考,并且其优先权声明依据35USC119。