![智能计算:原理与实践](https://wfqqreader-1252317822.image.myqcloud.com/cover/961/45852961/b_45852961.jpg)
1.5.1 基于U-支持向量机的正交小波盲均衡算法
基于U-支持向量机的正交小波盲均衡算法(Wavelet Blind Equalization Algorithm Based onU-Support Vector Machine, USVMWTCMA)的原理如图1.5.1所示。
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_01.jpg?sign=1738834480-1aR9tiffYIuqTRaclC9CqZxL1YNJSLSW-0-cbc95cb48c96f1475b60314932f8ec12)
图1.5.1 基于U-支持向量机正交小波盲均衡算法原理
图中,a(n)为发射信号,h(n)为信道的冲激响应,v(n)是信道输出端的加性高斯白噪声。y(n)为均衡器的接收信号,z(n)为均衡器的输出信号,为判决器输出信号。由图1.5.1可得
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_03.jpg?sign=1738834480-dS22btBCfL15cAVctDMyqRYzUMgoaQR6-0-137f708cd752bd47db9c235eaec184cd)
式中,V是正交小波变换(WT)矩阵。
在U-支持向量机中,集合U起着至关重要的作用。参考文献[9]和[10]对集合U进行了研究。选取合适的集合U,得到的结果与标准支持向量机相比,预测精度会有不同程度的提高;如果随机产生集合U,结果对预测精度几乎没有改变。
本节采用Vapnik提出的通用构造集合U的方法,即集合U中的输入用正类点的输入和负类点的输入均值生成x∗(k)=(x+(i)+x-(j))/2,x+(i)为正类点的输入,x-(j)为负类点的输入,x∗(k)为集合U中的一个元素。此时,支持向量机的训练集为T∪U。T={(x(1),y(1)),(x(2),y(2)),…,(x(N),y(N))}∈(RN×R)N,U={x∗(1),…,x∗(Nu)}。
U-支持向量机的目标是求出超平面(w·x)+b=0,使对训练集T满足间隔最大,且使U中的输入尽可能地接近该超平面。
设y(n)为均衡器接收信号的前N组输入数据,利用支持向量机对均衡器权向量w(n)进行初始化,即可转化为求解代价函数的最小值。
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_04.jpg?sign=1738834480-rVvFJcmoFHkqDz3jKCETa4EMEh6NvCvt-0-25fb9ad2ae3e0401c66f9e402d65714c)
约束条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_05.jpg?sign=1738834480-34VRTn3qmGhJz5MbKkB390OfIwXkVnqc-0-79778f4917898428a26ecadebdcab194)
式中,φ∗=(φ(1),φ∗(1),…,φ(Nu),φ∗(Nu))T;C1,C2为惩罚因子;ξ为松弛变量。
式(1.5.1)的最优化问题可以转换为求解拉格朗日函数的鞍点问题,引入拉格朗日函数
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_06.jpg?sign=1738834480-wVfbcwtU63AXQOD8TWjZqikneb5o5ABb-0-053bdffd657b0e6f2a5b3e0e72563fd5)
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_01.jpg?sign=1738834480-HET71iifJL4boHRvtuaWqRRhcjSFtnUS-0-29c236d5e9db4d7d8d507f7a8047cca2)
式中,α(n)≥0,η(n)≥0,ν(s)≥0,μ(s)≥0,γ(s)≥0,γ∗(s)≥0;η=(η(1),…,η(N))T,α=(α(1),…,α(N))T;ν=(ν(1),…,ν(Nu))T;γ=(γ(1),…,γ(Nu))T;μ=(μ(1),…,μ(Nu))T为乘子向量。
将式(1.5.6)分别对w、b求导
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_02.jpg?sign=1738834480-Rf4YwMrYuS6GMCV7XP6KrYhRZBXqFOTm-0-3dd84a9386d8a7572960b70848fbadd5)
然后,将其结果代入式(1.5.7),则其最优化问题为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_03.jpg?sign=1738834480-nFY6ysV42g8vloXvRGnfo0UqMGf7yoHl-0-307bd9d6844f4f615d10a6ef97999da3)
约束条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_04.jpg?sign=1738834480-E5FpyOQ6HU00sYxjBt0uGntfSvQIeJi0-0-878b2636a098ccffe483714ff2e3ef67)
在约束条件下求解式(1.5.7),计算出均衡器权系数w(n),进行循环迭代,直到结果满足切换条件时,切换至小波盲均衡算法中。切换条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_05.jpg?sign=1738834480-ROqiEXUTdPTgWcWEr605jV8xazGr9z53-0-2b3e1134d1515e140fc11a9c0b13634c)
式中,AME(n)为平均调制误差,η为切换阈值。
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_06.jpg?sign=1738834480-fWiaaDzTyxsPVBrnELamMXdvuP02eLvG-0-f48ab93212d50d8fb0c71030ffb60989)
式中,R1为输入点到各个收敛点的距离。