![智能计算:原理与实践](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=1739279666-8wnqVDK0ZQTJl9oyfPDAmug00hXoVa0P-0-4b3f3d985441ed91c4b4ecef876c83df)
图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=1739279666-xBtu8wET7q02FGDg0VUiibVQXEAH1HaE-0-a43d4037b96d59105e7f2334254d9e6c)
式中,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=1739279666-nNbxjp6l3scwVpNrEKeuwNUfJ9Z4XO7k-0-bb748fd2766467ef0ab06b0d3aa0c8c5)
约束条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_05.jpg?sign=1739279666-9HsuqgvwojI7tiQDYZFCWuLA3HeQJ9IR-0-5d3ec3fcd5c4ccba68a7aabcbd24efa9)
式中,φ∗=(φ(1),φ∗(1),…,φ(Nu),φ∗(Nu))T;C1,C2为惩罚因子;ξ为松弛变量。
式(1.5.1)的最优化问题可以转换为求解拉格朗日函数的鞍点问题,引入拉格朗日函数
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/35_06.jpg?sign=1739279666-ljSOmp7Bz6qd3UdJzd0ujNbXR63tgIbP-0-bb598f58999b233fb70c2eb300179a8f)
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_01.jpg?sign=1739279666-4gNdqF9XRczvpPmzEdXZe4qbeR2CCnPl-0-cf0aeacba238d0cb6b5ad59010e19777)
式中,α(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=1739279666-OXPcjhNdWrIWxAXJrFguj1ttZFAhCJqF-0-d383b2b538a5bb091e1a2ebcbe50771c)
然后,将其结果代入式(1.5.7),则其最优化问题为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_03.jpg?sign=1739279666-PK4DVhW7l6M0oo9qjoxTys04Qmaj0068-0-4b604828cdc3d5abca9c00e914bca4a6)
约束条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_04.jpg?sign=1739279666-JMgVzYgDFFoSMt7j87LTxXEofSlueUu6-0-400a8fedec0edf036987f7cf0e892532)
在约束条件下求解式(1.5.7),计算出均衡器权系数w(n),进行循环迭代,直到结果满足切换条件时,切换至小波盲均衡算法中。切换条件为
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_05.jpg?sign=1739279666-7ez5gm51YJuzJEeNoJMIYcJxtommV7pW-0-1dfc3fc5aae5baa554106cef5c332411)
式中,AME(n)为平均调制误差,η为切换阈值。
![](https://epubservercos.yuewen.com/2C2AFB/24975619809585106/epubprivate/OEBPS/Images/36_06.jpg?sign=1739279666-yMZ2nPtnR54s6nA6zbaGQ1FWtXHpuO0X-0-8988469743fb90598e1d7ea08ee6ea71)
式中,R1为输入点到各个收敛点的距离。