基于最小二乘均衡原理

基于最小二乘（LS ）估计的信道均衡原理

1、最小二乘估计原理

考虑短波数字通信系统的均衡问题，传输信道与均衡器的简化框图如图 1所示。

图 1 基于LS 估计的信道均衡

设均衡器采用图 2所示的N+M+1抽头的FIR 横向滤波器。均衡器工作在训练模式，d(n)为训练序列，也即均衡器的期望输出响应。均衡器的输入信号可以表示为卷积过程

u (n )=ℎ(n )∗d (n )+v (n )

式中，h(n)是信道的等效冲击响应，v(n)是加性白噪声过程。均衡器的作用即是降低短波信道h(n)造成的码间干扰等不利影响。

图 2 N+M+1抽头系数的横向均衡器（滤波器）

现有输入信号u(n)的L 个输入数据u(0), u(1), …, u(L-1)，期望响应（训练序列）d(n)的L 个样本d(0), d(1), …, d(L-1)。一般情况下要求L>2N+2M+1。

定义n 时刻的输入信号向量为

u ⃗ (n )=[u (n +M ) u (u +M −1) ⋯ u (n ) ⋯ u (n −N )]T

均衡器（滤波器）的权向量为（注意不是h(n)的权向量，无需计算信道h(n)的权向量）

w ⃗⃗ =

[w −M w −M+1 ⋯ w 0 ⋯ w N ]T

于是均衡器的估计输出为

d ̂(n )=w ⃗⃗ H u ⃗

(n )=u ⃗ T (n )w ⃗⃗ ∗=∑w k ∗u (n −k )N

k=−M

(1)

上式中，d

̂(n )为存在码间干扰的均衡器估计输出。第n 时刻的估计值需要输入第n 时刻之前的M 个样本值和第n 时刻之后的N 个样本值。注意M 可以取值为0。估计误差为

e (n )=d (n )−d ̂(n )=d (n )−u ⃗ T (n )w ⃗⃗ ∗

实际上进行滤波时，仅考虑当输入信号u(n)完全进入均衡器各节点的情况，即进入均衡器的样本数大于或等于N+M+1，此时各时刻的滤波估计误差为：

当n=N 时(根据图 2，第一个输出时刻是n=N,不是n=M)，有

e (N )=d (N )−d ̂(N )=d (N )−u ⃗ T (N )w ⃗⃗ ∗=d (N )−∑w k ∗u (N −k )N

k=−M

当n=N+1时，有

e (N +1)=d (N +1)−d ̂(N +1)=d (N +1)−u ⃗ T (N +1)w ⃗⃗ ∗

⋯⋯

当n=L-M-1时

e (L −M −1)=d (L −M −1)−∑w k ∗u (L −M −1−k )N

k=−M

将以上各式写成方程组的形式，并取共轭，有 ()()()()()()()

()()

()()()

()

()()*

*

*

110111111121M M N e N d N u N M u N M u w e N d N u N M u N M u w e L M d L M u L u L u L M N w --+++-⎡⎤⎡⎤⎡⎤

⎡⎤⎢⎥⎢⎥⎢

⎥⎢⎥+++++⎢⎥⎢⎥⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥

⎢⎥---------⎢⎥⎢⎥⎢⎥⎣⎦

⎣⎦⎣⎦⎣

⎦

（2）

定义误差向量e 和期望响应（训练序列）向量d

分别为 e =[e (N ) e (N +1) ⋯ e (L −M −1)]H d =[d (N ) d (N +1) ⋯d (L −M −1)]H 定义数据矩阵

()()()()()()

()()()

()

()()*

111121110H

u N M u N M u u N M u N M u u L u L u A u N u N u L M L M N ++-⎡⎤

⎢

⎥+++⎢⎥⎢⎥⎢

⎥

-----⎢⎥

⎣⎡⎤=+⎦--⎣⎦

= 由式（1）定义估计输出向量ˆd

为 ()()()H

ˆˆˆˆ11d d N d N d L M Aw ⎡⎤=+--=⎣

⎦

所以式（2）可以表示为

ˆe d d

d Aw =-=-（3） 均衡器的设计原则是，寻找权向量w ⃗⃗ 使得误差信号e(n)在某种意义下取得极小值。注意

到式（2）或者上式是由L-M-N 个方程构成的线性方程组，未知量w k ,k =−M,…,0,…,N 的个数是滤波器的抽头个数N+M+1。一般来说，输入数据的个数总是比滤波器权值的个数大得多，即L-M-N>N+M+1。此时对于式（2）或者式（3）来说，方程组的个数大于未知量的个数。此时，若令e =0，由式（3）可得方程组=Aw d ，但是该方程组在一般意义下无解，即找不到一个向量w ⃗⃗ 满足全部L-M-N 个方程。

尽管不能找到一个解w ⃗⃗ ，使得方程=Aw d 成立，但是可以找到一个w ⃗⃗ 的估计值ˆw

，使得

误差向量

ˆe d Aw

=- 在某种意义下取得极小值。

在最小二乘（LS ）意义下，使估计误差的模的平方和

()()ˆˆH

H

J e e d Aw

d Aw

==-- 取得极小值，所得到的解称为最小二乘解。要得到J 的极小值，首先求J 关于w ⃗⃗ 的梯度

*

2

22H H J

J A d A Aw w

∂∇==-+∂ 令0J ∇=，得

ˆH H A Aw

A d = （4） 上式是使J 取得极小值时,w ⃗⃗ 必须满足的条件，称为确定性正则方程。当L-M-N>N+M+1时，

如果方阵A H A 是非奇异的，那么在最小二乘意义下，确定性正则方程的解为

()1

ˆH H LS w

A A A d -= （5） 经常将ˆˆd

Aw =称为对期望响应向量d 的最小二乘估计，简称LS 估计。 2、最小二乘问题的求解

在工程实现时，如要获得均衡器的权向量估计值ˆw

，就要求解式4。如果直接求解确定性正则方程（式4）ˆH

H A Aw

A d =，需要进行矩阵求逆运算，因此必须考虑A H A 的奇异性。如果A H A 是非奇异的，矩阵求逆运算不仅计算量大，还有可能发散，工程实现上也不易实现。

式4的求解可以使用奇异值分解来求解，也可以使用递归最小二乘（RLS ）算法（求解算法很多，暂时不考虑其他算法）。本节将讨论RLS 算法。

RLS 算法使用迭代的方法求解最小二乘的确定性正则方程（式4），其基本思路是，已知

n -1时刻的均衡器权向量的最小二乘估计()ˆ1w

n -，利用当前n 时刻新得到的观测数据，用迭代的方法计算出n 时刻的滤波器权向量的最小二乘估计()ˆw

n 。 这里忽略推导过程，直接给出RLS 算法： 步骤1 初始化：()()()

111

0M N M N P I C

δ++⨯++-=∈，δ是小的正数（0.005左右），P 是

输入数据的时间相关矩阵的逆；()ˆ0=0w

，遗忘因子λ一般取值接近于1，如λ=0.98，λ=0.998等。遗忘因子使得离当前时刻近的观测值，对相关矩阵和互相关向量的影响较大，而较久远的值则影响较小。

步骤2 当n = 0, 1, 2, …, L-1时，完成如下的迭代运算：