一种盲均衡方法和一种盲均衡系统的制作方法

一种盲均衡方法和一种盲均衡系统的制作方法
【技术领域】
[0001] 本发明涉及无线和数字通信领域，尤其涉及一种盲均衡方法和一种盲均衡系统。
【背景技术】
[0002] 在遥感卫星通信中，发送信号在传输时由于受到信道、发送器件和接收器件等的 不理想因素的影响，会造成发送信号的畸变或相互干扰，即产生失真，其中最为显著的即为 线性失真和非线性失真。线性失真表现为幅频特性不平坦和群时延失真，会使得接收到的 信号产生码元串扰（ISI，IntersymbolInterference);非线性失真一般由发送端功率放大 器，如行波管放大器或固态功率放大器，工作于近饱和或饱和区域造成，会使得接收到的信 号产生畸变。这些由线性或非线性失真导致的ISI或畸变均会造成极大的误码，使得接收 端无法辨认出发送信号，因此，在接收端，为了抵消掉这些失真的影响，还原出原始发送信 号，就需要加入一个均衡器，对接收到的信号进行失真补偿；同时，实际卫星通信中，基于传 输效率的考虑，卫星是不会发送训练序列供均衡器进行先期训练的，因此，必须采用盲均衡 的方式，即在期望信号未知的情况下进行失真补偿操作。
[0003] 目前，失真补偿技术可以分为两大类，一类是发送端补偿，也就是预失真；一类是 接收端补偿，也就是均衡。预失真可分为固定和自适应两种，固定预失真无法跟踪器件的特 性变化，而自适应预失真则需要在卫星中加入模拟域到数字域的反馈电路，硬件实现复杂 度过高，不利于卫星的稳定可靠工作。
[0004] 现在常用的均衡器为线性均衡器，如线性横向均衡器（LTE，LinearTransversal Equalizer)等。但是线性均衡器只能补偿线性失真，对非线性失真的补偿效果则很差。因 此，在接收端补偿时，为了能够补偿更多的失真，更好地完成信号的校正，需要使用非线性 均衡器。
[0005] 非线性均衡器现有两类，一类是基于多项式滤波的，如volterra均衡器；一类是 基于神经网络的均衡器。基于多项式滤波的均衡器一般不支持联想记忆，结构较为固定， 需要使用高阶多项式才能保证好的均衡性能，实现复杂度较高，不利于硬件实现；基于神 经网络的均衡器结构灵活，其特有的并行处理特性能够支持高码率时的并行操作要求，支 持联想记忆、特征提取和知识处理，具备很强的分类和识别能力，其中的双层前向神经网络 (TLFNN，Two-layerFeedforwardNeuralNetwork)实现较简单，便于硬件实现。而基于 TLFNN的盲均衡器正是本发明研究的重点。
[0006] 在使用TLFNN盲均衡器的过程中，最为重要的一点就是需要避免均衡陷入局部收 敛状态。因为一旦局部收敛，很可能使得均衡的输出误码率（SER，SymbolErrorRate)比 均衡的输入SER还高，即在接收端加入盲均衡器之后反而使得信号质量降低了，均衡器起 到了"反效果"。但是，TLFNN的非线性结构会使得均衡性能曲面上的局部收敛点急剧增多， 从而使得均衡更易陷入局部收敛状态，同时，不合理的初始递归更新参数值也可能导致均 衡过程陷入局部收敛；而且，由于期望信号是未知的，无法有效探知盲均衡器是否处于局部 收敛状态，也就无法采取手段来跳出局部收敛。而盲均衡器是否会陷入局部收敛状态与其 所采用的盲均衡算法是息息相关的，因此，需要找出适用于TLFNN盲均衡器的能够有效避 免局部收敛的盲均衡算法。
[0007]基于最小均方误差（LMS，LeastMeanSquare)的面向判决（DD，Decision-Directed)算法是比较经典的盲均衡算法，以其实现简单而广泛应用。但是，DD 算法需要极小的初始SER和合适的初始递归更新参数才能确保较好收敛，而在实际通信应 用中，并不是总能保证初始SER极小的，而且由于信道特性往往是未知的，也不可能在初始 时就配置出合适的递归更新参数。针对这种情况，现阶段通常采用以下两种方法进行改进， 一种是基于额外信息辅助的方法，一种是基于双模式切换的方法。
[0008] 基于额外信息辅助的方法通过发掘更多的有别于DD算法中采用的输出判决信 息的辅助信息来提升盲均衡的收敛效果。M.Lazaro和I.Santamaria等人于2005年提出 的SQD(Stochasticquadraticdistance)算法，就使用了均衡输出信号的概率密度函数 (pdf,probabilitydensityfunction)信息，并将其与标准星座的pdf进行匹配，达到提升 收敛效果的目的。D.Erdogmus和J.C.Principe则于2002年以均衡输出误差的熵作为辅助 信息，提出一种最小误差熵算法（MEE，Error-EntropyMinimizationAlgorithm)。而著名 的Sato算法及Godard算法则可以看做使用了输出信号的高维统计量作为辅助信息。基于 双模式切换的方法则通过模式切换的方式，结合盲均衡捕获算法和跟踪算法的优点，达到 避免局部收敛同时获得较好性能的目的。在盲均衡初期的捕获阶段，采用捕获能力较强的 算法(如恒模算法（CMA，ConstantModulusAlgorithm))保证收敛效果，避免局部收敛；之 后切换到跟踪阶段，采用跟踪能力较强的算法(如DD)，获得较好的均衡性能。几个经典的例 子有BG(Benveniste-Goursat)算法、DM-CMA(DualModeCMA)算法及其扩展SAG-DM-CMA (Stop-and-GoDM-CMA)算法。
[0009] 上述两种改进方法的收敛效果均优于DD算法，但是，上述算法存在两点不足：一 是受递归更新参数的初始配置值的影响较大，不合理的配置值会使得这些算法仍然较易陷 入局部收敛；二是针对TLFNN盲均衡器所采用的非线性结构，这些算法较难避免局部收敛。 这是因为，上述算法中采用的辅助信息均来自于均衡输出信号，而输出信号由均衡器生成， 直接受均衡器所采用前馈结构及递归更新参数初始值的影响。因此，对于TLFNN盲均衡器 的非线性结构以及无法合理设置递归更新参数初始值的情况，这些算法仍然无法避免局部 收敛。

【发明内容】

[0010] (一）要解决的技术问题
[0011] 本发明提供一种盲均衡方法和一种盲均衡系统，以解决TLFNN盲均衡器中由于非 线性结构和无法合理设置递归更新参数初始值而使得盲均衡较易陷入局部收敛状态的技 术问题。
[0012] (二）技术方案
[0013] 为解决上述技术问题，本发明提供一种盲均衡方法，包括：
[0014] 将输入信号经输入变换后生成输入变换后矢量；
[0015] 利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信号；
[0016] 在期望信号未知的情况下，利用所述输入信号获得输入判决信息，利用所述输出 信号获得输出判决信息，利用所述输入判决信息和所述输出判决信息联合设计盲均衡算 法，获得反馈误差；
[0017] 利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参数。
[0018] 进一步地，
[0019] 所述将输入信号经输入变换后生成输入变换后矢量包括：
[0020] 将N维输入信号矢量表不为：xk，2，…，xk，N)T，令N维输入M维输出的输 入变换函数为F(X)，则输入信号Xk经输入变换后生成的M维输入变换后矢量表示为： ①k_(小k，1，小k，2，? ??，小k，M);
[0021] 所述利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信号包 括：
[0022] 将输出信号yk表示为：凡=g((丨< ，
[0023] 其中8(1)为激活函数，扎=〇115，1，1115,2，"、1115，11)1^为递归更新参数。
[0024] 进一步地，
[0025] 所述利用所述输入信号获得输入判决信息包括：检测输入信号Xk得到中间元素判 决值，即为输入判决信息ik，ik=dec(xk，J，其中L=(N+l)/2,dec( ?)为判决函数；
[0026] 所述利用所述输出信号获得输出判决信息包括：检测输出信号yk得到判决值，即 为输出判决信息dk，dk=dec(yk)。
[0027] 进一步地，
[0028] 所述在期望信号未知的情况下，利用所述输入判决信息和所述输出判决信息联合 设计盲均衡算法，获得反馈误差，包括：利用下述公式获得反馈误差ek:
[0029]ek= (dk_yk) +ak |dk_yk | (ik_yk)，
[0030] 式中a k为自适应系数，随收敛过程逐渐变化：
[0031]ak=yEk
[0032]Ek+1= 3Ek+ (1_ 3 ) (dk_yk)2
[0033] 其中线性权重Y和遗忘参数P为经验正值；
[0034] 或，利用下述公式获得反馈误差ek:
[0035]
[0036] 进一步地，
[0037] 所述利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参数，包 括：利用下述公式更新递归更新参数：
[0038]
[0039]其中U为更新步进。
[0040] 另一方面，本发明还提供一种盲均衡系统，包括：输入变换单兀、输出信号单兀、输 入判决单元、输出判决单元、盲均衡单元和反馈更新单元，其中 ：
[0041] 输入变换单元，用于将输入信号经输入变换后生成输入变换后矢量，输送至输出 信号单元；
[0042] 输出信号单元，用于利用所述输入变换后矢量和递归更新参数，通过前馈网络获 取输出信号，输送至输出判决单兀；
[0043] 输入判决单元，用于利用所述输入信号获得输入判决信息，并输送至盲均衡单 元；
[0044] 输出判决单元，用于利用所述输出信号获得输出判决信息，并输送至盲均衡单 元；
[0045] 盲均衡单元，用于在期望信号未知的情况下，利用所述输入判决信息和输出判决 信息联合设计盲均衡算法，获得反馈误差，输送至反馈更新单元；
[0046] 反馈更新单元，用于利用所述反馈误差和输入变换后矢量，通过反馈过程更新递 归更新参数，并回传至所述输出信号单元。
[0047] 进一步地，
[0048] 所述输入变换单元包括： 变换子单元，用于将N维输入信号矢量Xk表示为： xk=(Xk，i，xk，2，…，xk，N)T，经N维输入M维输出的输入变换函数F(X)进行输入变换之后生成M 维输入变换后矢量(1\表示为：(1\=〇1 5，1，<1)15,2，"％<1 )15，11)'
[0049] 所述输出信号单元包括：计算子单元，用于通过前馈网络计算获得输出信号yk，表 示为如下公式=g(巾丨//，），其中g(x)为激活函数，^=0^,11!^，…，hk，M)T为递归更新参 数。
[0050] 进一步地，
[0051] 所述输入判决单元包括：输入判决子单元，用于检测输入信号Xk得到中间元素判 决值，即为输入判决信息ik，公式为：ik=dec(Xk』)，其中L=(N+l)/2,dec( ?)为判决函数；
[0052] 所述输出判决单元包括：输出判决子单元，用于检测输出信号yk得到判决值，即为 输出判决信息dk，dk=dec(yk)。
[0053] 进一步地，
[0054] 所述盲均衡单元包括：第一盲均衡子单元，或第二盲均衡子单元；
[0055] 所述第一盲均衡子单元用于在期望信号未知的情况下，利用所述输入判决信息ik 和所述输出判决信息dk联合设计盲均衡算法，获得反馈误差ek :
[0056] ek= (dk_yk) +ak |dk_yk | (ik_yk)，
[0057] 式中ak为自适应系数，随收敛过程逐渐变化：
[0058] ak=yEk
[0059] Ek+1= 3Ek+ (1- 3 ) (dk_yk)2
[0060] 其中线性权重Y和遗忘参数0为经验正值；
[0061] 第二盲均衡子单元用于在期望信号未知的情况下，利用所述输入判决信息ik和所 述输出判决信息dk联合设计盲均衡算法，获得反馈误差ek:
[0062]
[0063] 进一步地，
[0064] 所述反馈更新单元包括：更新子单元，用于利用下述公式更新递归更新参数：
[0065]
[0066] 其中U为更新步进。
[0067](三）有益效果
[0068] 在本发明提出的盲均衡方法和系统中，加入了输入判决信息来辅助反馈过程，由 于输入判决信息是不受均衡器本身影响的，所以能够有效避免盲均衡中受前馈网络结构 及递归更新参数初始设置值影响而常出现的局部收敛的问题，保证了盲均衡结果的稳健收 敛；同时，联合输入判决信息和输出判决信息的盲均衡算法的采用，也保证了盲均衡的性 能。
【附图说明】
[0069] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根 据这些附图获得其他的附图。
[0070] 图1是盲均衡基础算法模型示意图；
[0071] 图2是本发明实施例盲均衡方法的基本流程示意图；
[0072] 图3是本发明一个优选实施例盲均衡方法的流程示意图；
[0073] 图4是本发明实施例基于输入输出判决辅助的盲均衡结构示意图；
[0074] 图5是本发明实施例ID算法输出SER与输入SER结果对比图；
[0075]图6是本发明实施例四种算法确定初始状态线性均衡收敛性仿真结果图；
[0076]图7是本发明实施例四种算法随机初始状态线性均衡收敛性仿真结果图；
[0077] 图8是本发明实施例四种算法确定初始状态非线性均衡收敛性仿真结果图；
[0078] 图9是本发明实施例四种算法随机初始状态非线性均衡收敛性仿真结果图；
[0079] 图10是本发明实施例四种算法线性均衡SER性能对比图；
[0080] 图11是本发明实施例四种算法非线性均衡SER性能对比图；
[0081] 图12是本发明实施例盲均衡系统的基本结构示意图；
[0082] 图13是本发明一个优选实施例盲均衡系统的结构示意图。
【具体实施方式】
[0083] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0084]图1为盲均衡的基础算法模型，图中Xk=(Xlu，xk，2，…，xk，N)T为N维输入信号矢 量，①，…，为输入信号经过输入变换后生成的M维输入变换后矢量， 扎二队，"hk2,…，hkM)T为递归更新参数，yk为输出信号，dk=dec(yk)为判决函数dec( ?)的 输出判决值。图1中F(X)为N维输入M维输出的输入变换函数，g(x)为激活函数。
[0085] 盲均衡的目的就是在期望信号ak未知的情况下通过Xk和递归更新参数Hk获取yk， 以使yk与ak误差最小甚至无误差。实现这个目标有两个过程，一是通过前馈网络获取yk， 前馈方程如下式：
[0086]
[0087] 二是通过反馈过程更新扎。反馈方程如下式：
[0088]
[0089] 式中U为更新步进。
[0090] 本模型同时适用于线性横向均衡器（LTE，LinearTransversalequalizer)和基 于双层前向神经网络（TLFNN，Two-LayerFeedforwardNeuralNetwork)的非线性均衡器。
[0091] 当F (X) =X且g (x) =x时，该模型为LTE模型，则：
[0092]
[0093]
[0094] 当F(X)或g(x)为非线性函数时，该模型转化为非线性盲均衡模型。对于TLFNN， 其F(X)可以是切比雪夫多项式、勒让德多项式或其它扩展函数；其g(x)通常采用著名的双 曲正切函数g (x) =tanh (x)。
[0095] 当使用4阶切比雪夫多项式扩展且g (x) =tanh (x)时，则：
[0096]
[0097]
[0098]其中：
[0099] 〇k=(l. 0,T〇(xk,D.T^x,,,l),T2(xk,l),T3(xk, 1) ???,!〇 (xk,N)^! (xk,N),T2 (xk,N),T3 (xk,N))tTq(x)=1. 0,I\(x)=x，Tn+1WzSxTnWUx)，C>k 中常量 1. 0 用来补偿可能存在的 直流偏置。
[0100] 对于DD算法，反馈方程中的ek=dk_yk。
[0101]对于DM-CMA算法：
[0102]
[0103] (
[0104] 上还DD和DM-CMA算法1 乂使用湔出判決信息进行自适应，输出判决信息来源于输 出信号，而输出信号由均衡器生成，直接受均衡器所采用前馈结构及递归更新参数初始值 的影响，仿真发现该类基于输出判决信息的盲均衡方法较易局部收敛，使得均衡效果很差， 为了避免局部收敛，需要考虑输入判决信息，由于输入判决信息来自于输入信号，不受均衡 器本身结构和参数的影响，因而能够有效地避免局部收敛的产生。
[0105] 因此，本发明实施例首先提出一种盲均衡方法，参见图2,包括：
[0106] 步骤201 :将输入信号经输入变换后生成输入变换后矢量。
[0107] 步骤202:利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信 号。
[0108] 步骤203:在期望信号未知的情况下，利用所述输入信号获得输入判决信息，利用 所述输出信号获得输出判决信息，利用所述输入判决信息和所述输出判决信息联合设计盲 均衡算法，获得反馈误差。
[0109] 步骤204:利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参 数。
[0110] 可见，在本发明实施例提出的盲均衡方法中，加入了输入判决信息来辅助反馈过 程，由于输入判决信息是不受均衡器本身影响的，所以能够有效避免盲均衡中受前馈网络 结构及递归更新参数初始设置值影响而常出现的局部收敛的问题，保证了盲均衡结果的稳 健收敛；同时，联合输入判决信息和输出判决信息的盲均衡算法的采用，也保证了盲均衡的 性能。
[0111] 在本发明的一个实施例中，优选地，将输入信号经输入变换后生成输入变换后 矢量可以包括：将N维输入信号矢量表不为：Xk=(xtl，xk，2，…，xk，N)T，令N维输入M维输 出的输入变换函数为F(X)，则输入信号Xk经输入变换后生成的M维输入变换后矢量表 示为(小u，，…，^M)T。利用输入变换后矢量和递归更新参数，通过前馈网络 获取输出信号可以包括：将输出信号yk表示为：.、其中g(x)为激活函数， ，…，hk，M)T为递归更新参数。利用输入信号获得输入判决信息可以包括：检测输 入信号Xk得到中间元素判决值，即为输入判决信息ik，ik=dec(Xk』)，其中L=(N+l)/2。
[0112] 在本发明的另一个实施例中，当单纯考虑输入判决信息，所采用的输 入判决算法（ID,Inpute-Directed)可以为：在更新递归更新参数时，利用公式
々=U?…进行计算，其中u为更新步进。
[0113] 可见ID算法将DD算法的输出判决信息dk替换为输入判决信息ik，定义输入SER 为：
[0114] pin=p(ak 关ik)
[0115] 输出SER为：
[0116] P〇ut=P (ak ^ dk)
[0117] 经过仿真发现如果输入SERpinS够小，ID算法得到的输出SERP()Ut总能小于pin，即 ID算法总能保证盲均衡器的稳健收敛。
[0118] ID算法虽然能够保证稳健收敛性，但是SER性能却较DD算法差，因此我们结合输 入判决信息和输出判决信息两方面的优点，在保证SER性能与其他算法相当的情况下完成 均衡过程的稳健收敛，这样得到的盲均衡既不受递归更新参数初始值的影响，也不受所采 用前馈结 构的影响，保证了可靠性和稳定性。
[0119] 在本发明的一个实施例中，结合输入判决信息和输出判决信息的盲均衡方法中， 利用输入信号获得输入判决信息可以包括：检测输入信号Xk得到中间元素判决值，即为输 入判决信息ik，ik=dec(xk，J，其中L= (N+1)/2,dec( ?)为判决函数；利用输出信号获得输出 判决信息可以包括：检测输出信号yk得到判决值，即为输出判决信息dk，dk=dec(yk)。
[0120] 其中，基于输入判决信息和输出判决信息的盲均衡方法可以分别利用 BG-IOD(Benveniste_GoursatInput-Output-Decision)方法和SAG-IOD(Stop-and-Go Input-Output-Decision)方法进行计算，在更新递归更新参数时，可以利用公式
，其中u为更新步进。
[0121] 在本发明的一个实施例中，对于BG-I0D算法，可以在期望信号未知的情况下，利 用输入判决信息和输出判决信息联合设计盲均衡算法，获得反馈误差：
[0122] ek= (dk-yk) +akIdk-ykI(ik-yk)
[0123] 式中的ak为自适应系数，在初始捕获状态有一个较大的值，并在收敛过程中逐渐 减小。这里使用输出误差对ak进行调节：
[0124]ak=yEk
[0125]Ek+1= @Ek+ (1_ @ ) (dk_yk)2
[0126] 其中线性权重Y和遗忘参数P为经验正值。
[0127] 当输出误差dk_yk比较大的时候，ek表达式中的后半部分起主导作用，此时均衡器 工作类似于ID;随着输出误差的逐渐减小，ek表达式中的前半部分的权重逐渐增大，均衡器 工作趋向于DD。这个过程即为从输入判决信息到输出判决信息的软切换过程。
[0128] 对于SAG-I0D算法，也可以在期望信号未知的情况下，利用输入判决信息和输出 判决信息联合设计盲均衡算法，获得反馈误差
参数矢量Hk仅在输 出判决信息dk等于输入判决信息ik时才进行更新。
[0129] 下面以利用BG-I0D算法进行盲均衡为例，来具体说明本发明实施例的一个实现 过程，参见图3。
[0130] 步骤301:将输入信号经输入变换后生成输入变换后矢量。
[0131] 本发明实施例的基于输入判决信息和输出判决信息的盲均衡结构参见图4。本 步骤中，将N维输入信号矢量表示为：Xk=(xk;1，xk，2，…，xk，N)T，令N维输入M维输出的输 入变换函数为F(X)，则输入信号Xk经输入变换后生成的M维输入变换后矢量表示为： ①k_ (小k，1，小k，2，? ? ?，小k，M)。
[0132] 步骤302:利用输入变换后矢量和递归更新参数，通过前馈网络获取输出信号。
[0133] 本步骤中，将输出信号yk表示为：凡=g(叫乂），
[0134] 其中8(1)为激活函数，扎=〇115，1，11 15,2，"、1115，11)1为递归更新参数。
[0135] 步骤303:利用输出信号获得输出判决信息。
[0136] 本步骤中，利用输出信号yk通过判决函数dec( ?)获得输出判决信息：dk=dec(yk)。
[0137] 步骤304:利用输入信号获得输入判决信息。
[0138] 本步骤中，检测输入信号Xk得到中间元素判决值，即为输入判决信息ik， ；[15=(16(3(叉15， 1)，其中1=(奸1)/2。
[0139] 步骤305 :在期望信号ak未知的情况下，利用输入判决信息和输出判决信息联合 设计盲均衡算法，获得反馈误差。
[0140] 对于136-100算法，反馈误差：615=((115-：715) + (115|(115-：715|(;[15-：7 15)
[0141] 上式中ak为自适应系数，随收敛过程逐渐减小：
[0142]ak=yEk
[0143][email protected]+([email protected])(dk-yk)2
[0144] 其中线性权重Y和遗忘参数P为经验正值。
[0145] 步骤306:利用输入判决信息、输出判决信息、输入变换后矢量和输出信号，通过 反馈过程更新递归更新参数。
[0146] 本步骤中，利用如下公式更新递归更新参数：
[0147]
[0148] u为更新步进。
[0149] 步骤307:通过更新后的递归更新参数更新输出信号。
[0150] 此步骤中，利用更新后的递归更新参数Hk+1更新输出信号yk，并将更新的输出信号 代入下次盲均衡算法的计算中。
[0151] 至此，则完成了本发明实施例利用BG-I0D算法进行盲均衡的全过程。
[0152] 另外，需要说明的是，上述基于图3的所有流程描述是本发明盲均衡方法的一种 优选的实现过程，在本发明盲均衡方法的实际实现中，可以根据需要在图3所示流程的基 础上进行任意变形，可以是选择图3中的任意步骤来实现，各步骤的先后顺序也可以根据 需要调整等。
[0153] 为验证本发明实施例方法的效果，采用如下方式对上述各项算法进行了仿真。仿 真传输信道包括线性和非线性两个部分，线性部分的脉冲响应传递函数为
[0154]
[0155] 式中A设为2. 9,这样获得特征值比（EVR，EigenValueRatio)为6. 08,从而获 得Z域归一化转换函数：
[0156] 0. 209+0. 995z_1+0. 209z-2
[0157] 为了检测均衡器的非线性特性，本发明实施例在线性部分之后加入下面的非线性 部分：
[0158]NL=0:Nonk=Link
[0159]NL=1:Nonk=tanh(Link)
[0160] 其中NL=0对应线性信道模型，NL=1对应非线性信道模型。
[0161]传输信号为4-PAM信号，其格式为{-0. 9,-0. 3,+0. 3,+0. 9}，信道输出端叠加一个 零均值的高斯白噪声。
[0162] 对于线性模型，信道NL=0,均衡器采用11阶LTE;对于非线性模型，信道NL=1，均 衡器采用3阶TLFNN，应用4阶切比雪夫多项式扩展及双曲正切激活函数。在本TLFNN中，输 入信号通过切比雪夫多项式扩展后变换为13维输入变换后矢量(包括一个常量偏置1. 0)。
[0163] 本发明实施例采用这种线性和非线性均衡模型，对各种算法进行仿真验证，包括 ID、BG-IOD、SAG-IOD、DM-CMA和SQD。
[0164] 通过大量的实验，得到这些算法的最好的参数配置值。最终LTE和TLFNN的更新 步进分别取 1E-3 和 1E-4 ;对BG-I0D取Y=1. 0, 0 =0. 999995,EF30. 0 ;对DM-CMA，经验值d 取 0? 1;对SQD，参数取值为：a=l. 0,b=0. 0,a=〇? 000005,Epl. 0。
[0165] 首先对ID算法的输出SER与输入SER进行比较以证明ID算法的收敛效果；然后 针对后四种算法检验其在线性信道和非线性信道中的收敛性及SER性能。对于每一个信噪 比（SNR)值进行100次独立运算，每次运算包含107次迭代，为了获得较为稳定的SER值，对 每次仿真取最后1〇6次迭代进行运算获得最终结果。另外为了对比SER性能，还对有监督 学习算法进行了仿真，其公式为：
[0166]
[0167] 其中期望信号ak已知。
[0168] 在ID输出SER收敛性验证方面，图5中给出了 100次仿真中最大输出SER和最 小输入SER的对比曲线，线性信道时SNR范围设为[12, 24]dB，非线性信道时SNR范围设为 [12, 30]dB，左图为线性均衡结果，右图为非线性均衡结果。可以看到对于线性或非线性均 衡器，如果输入SER小于2E-1，ID算法所获得的输出SER总能小于输入SER。
[0169] 在收敛性方面，图6-9给出了四种盲均衡算法：86-100、5八6-100、01^嫩、5〇0收敛 效果的仿真对比：（3)01-〇麻；（13)5〇0;((3)86-100 ;((1)5六6-100。图6和图7为线性均 衡下的收敛性结果，其中图6为中心点初始递归更新参数（即递归更新参数的中间点为1， 其他点均为〇);图7为随机初始递归更新参数，SNR均为20dB。
[0170] 从图6可以看到，在中心点初始递归更新参数下四种均衡器均能达到好的收敛 效果，但是在图7中，当初始递归更新参数随机选取时，DM-CMA和SQD出现局部收敛，而 BG-I0D和SAG-I0D仍然工作良好。
[0171] 图8和图9给出了非线性均衡下的收敛性结果，其中图8为中心点初始递归更新 参数下的仿真结果；图9为随机初始递归更新参数下的仿真结果，SNR均为27dB。从以上两 个图中可以看到对于DM-CMA和SQD，在非线性均衡情况下无论初始递归更新参数采用中心 点方式还是随机选取，都会进入局部收敛状态，而BG-I0D和SAG-I0D则在上述两种初始递 归更新参数设置情况下均工作良好，其中SAG-I0D偶尔会进入局部收敛，表明SAG-I0D的稳 健收敛性要略逊于BG-I0D。
[0172] 四种算法的SER性能如图10和图11所示，图10为线性均衡，SNR范围设为 [12,23]dB，图11为非线性均衡，SNR范围设为[ll，33]dB。这里我们采用中心点初始递归 更新参数进行仿真，从图10可以看到对于线性均衡情况，四种盲均衡算法的SER性能与有 监督学习算法的SER性能相当。从图11可以看到BG-I0D的SER性能要略好于SAG-I0D，而 DM-CMA和SQD在非线性情况下已完全陷入局部收敛状态。
[0173] 本发明一个实施例还提出了一种盲均衡系统，见图12,包括：
[0174] 输入变换单元1201，用于将输入信号经输入变换后生成输入变换后矢量，输送至 输出信号单元1202;
[0175] 输出信号单元1202,用于利用所述输入变换后矢量和递归更新参数，通过前馈网 络获取输出信号，输送至输出判决单元1204;
[0176] 输入判决单元1203,用于利用所 述输入信号获得输入判决信息，并输送至盲均衡 单元1205 ;
[0177] 输出判决单元1204,用于利用所述输出信号获得输出判决信息，并输送至盲均衡 单元1205 ;
[0178] 盲均衡单元1205,用于在期望信号未知的情况下，利用所述输入判决信息和输出 判决信息联合设计盲均衡算法，获得反馈误差，输送至反馈更新单元1206 ;
[0179] 反馈更新单元1206,用于利用所述反馈误差和输入变换后矢量，通过反馈过程更 新递归更新参数，并回传至输出信号单元1202。
[0180] 在本发明的一个实施例中，优选地，输入变换单元1201可以包括：变换子单元 1301，见图13,用于将N维输入信号矢量表不为：xk，2，…，xk，N)T，令N维输入M维输 出的输入变换函数为F(X)，将输入信号Xk经输入变换后生成的M维输入变换后矢量表示 为(小tl，叭2，…，WT。输出信号单元1202可以包括：计算子单元1302,用于将输 出信号yk表示为如下公式=以叫仏），其中g(x)为激活函数，，…，hk,M)T为 递归更新参数。
[0181] ID算法虽然能够保证稳健收敛性，但是SER性能却较DD算法差，因此我们结合输 入判决信息和输出判决信息两方面的优点，在保证SER性能与其他算法相当的情况下完成 均衡过程的稳健收敛，这样得到的盲均衡既不受递归更新参数初始值的影响，也不受所采 用前馈结构的影响，保证了可靠性和稳定性。在本发明的一个实施例中，优选地，输入判决 单元1203可以包括：输入判决子单元1303,用于检测输入信号Xk得到中间元素判决值，即 为输入判决信息ik，公式为：ik=dec(Xk』)，其中L=(N+l)/2,dec( ?)为判决函数；输出判决 单元1204可以包括：输出判决子单元1304,用于检测输出信号yk得到判决值，即为输出判 决信息dk，dk=dec(yk)。
[0182] 在本发明的另一个实施例中，优选地，盲均衡单元1205可以包括：第一盲均衡子 单兀1305,或第二盲均衡子单兀1306。其中，第一盲均衡子单兀1305用于在期望信号未知 的情况下，利用所述输入判决信息ik和所述输出判决信息dk联合设计盲均衡算法，获得反 馈误差ek :
[0183]ek= (dk_yk) +ak |dk_yk | (ik_yk)，
[0184] 式中ak为自适应系数，随收敛过程逐渐变化：
[0185]ak=yEk
[0186]Ek+1= @Ek+ (1_ @ ) (dk_yk)2
[0187] 其中线性权重Y和遗忘参数P为经验正值；
[0188] 第二盲均衡子单元1306用于在期望信号未知的情况下，利用所述输入判决信息ik 和所述输出判决信息dk联合设计盲均衡算法，获得反馈误差ek:
[0189]
[0190] 在本发明的一个实施例中，优选地，反馈更新单元1206可以包括：更新子单元 1307,用于利用下述公式更新递归更新参数：
[0191]
[0192] 其中u为更新步进。
[0193] 需要说明的是，上述图13所示的盲均衡系统的各个实施例的结构可以进行任意 组合使用。
[0194] 可见，本发明实施例具有如下有益效果：
[0195] 在本发明实施例提出的盲均衡方法和系统中，加入了输入判决信息来辅助反馈过 程，由于输入判决信息是不受均衡器本身影响的，所以能够有效避免盲均衡中受前馈网络 结构及递归更新参数初始设置值影响而常出现的局部收敛的问题，保证了盲均衡结果的稳 健收敛；同时，联合输入判决信息和输出判决信息的盲均衡算法的采用，也保证了盲均衡的 性能。
[0196] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽 管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然 可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替 换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精 神和范围。
【主权项】
1. 一种盲均衡方法，其特征在于，包括： 将输入信号经输入变换后生成输入变换后矢量； 利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信号； 在期望信号未知的情况下，利用所述输入信号获得输入判决信息，利用所述输出信号 获得输出判决信息，利用所述输入判决信息和所述输出判决信息联合设计盲均衡算法，获 得反馈误差； 利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参数。2. 根据权利要求1所述的盲均衡方法，其特征在于： 所述将输入信号经输入变换后生成输入变换后矢量包括： 将N维输入信号矢量表不为Ak=(X)u, xk，2,…，xk，Ν)τ，令N维输入M维输出的输入 变换函数为F(X)，则输入信号Xk经输入变换后生成的M维输入变换后矢量表示为： ① k_ (小 k，1，小 k，2，· · ·，小 k，M); 所述利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信号包括： 将输出信号yk表不为：=#(Φ丨"J， 其中g(x)为激活函数，Hk=(hk l, hk 2,…，hk M)T为递归更新参数。3. 根据权利要求2所述的盲均衡方法，其特征在于： 所述利用所述输入信号获得输入判决信息包括：检测输入信号Xk得到中间元素判决 值，即为输入判决信息ik，ik=dec(xk，J，其中L=(N+l)/2, dec( ·)为判决函数； 所述利用所述输出信号获得输出判决信息包括：检测输出信号yk得到判决值，即为输 出判决信息dk，dk=dec (yk)。4. 根据权利要求3所述的盲均衡方法，其特征在于： 所述在期望信号未知的情况下，利用所述输入判决信息和所述输出判决信息联合设计 盲均衡算法，获得反馈误差，包括：利用下述公式获得反馈误差ek : ek= (dk-yk) + a k I dk-yk | (ik-yk), 式中ak为自适应系数，随收敛过程逐渐变化： α k= Y Ek Ek+i= ^ Ek+ (1- β ) (dk-yk)2 其中线性权重Y和遗忘参数β为经验正值； 或，利用下述公式获得反馈误差ek :5. 根据权利要求4所述的盲均衡方法，其特征在于： 所述利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参数，包括：利 用下述公式更新递归更新参数： 凡+1=巧+"ΚΧΦ?",,)%，& = 1，2，·· 其中u为更新步进。6. -种盲均衡系统，其特征在于，包括：输入变换单兀、输出信号单兀、输入判决单兀、 输出判决单元、盲均衡单元和反馈更新单元，其中： 输入变换单兀，用于将输入信号经输入变换后生成输入变换后矢量，输送至输出信号 单元； 输出信号单元，用于利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输 出信号，输送至输出判决单元； 输入判决单元，用于利用所述输入信号获得输入判决信息，并输送至盲均衡单元； 输出判决单元，用于利用所述输出信号获得输出判决信息，并输送至盲均衡单元； 盲均衡单元，用于在期望信号未知的情况下，利用所述输入判决信息和输出判决信息 联合设计盲均衡算法，获得反馈误差，输送至反馈更新单元； 反馈更新单元，用于利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更 新参数，并回传至所述输出信号单元。7. 根据权利要求6所述的盲均衡系统，其特征在于： 所述输入变换单元包括：变换子单元，用于将N维输入信号矢量Xk表示为： Xk= (xtl，xt2，…，xtN)T，经N维输入M维输出的输入变换函数F(X)进行输入变换之后生成M 维输入变换后矢量O k表示为ιΦ?ΗΦ?^，<K，2,…，<K，M)T ; 所述输出信号单元包括：计算子单元，用于通过前馈网络计算获得输出信号yk，表示为 如下公式：Ji 丨巧），其中g(x)为激活函数，Hk=Qilu, hk，2,…，hk，M)T为递归更新参数。8. 根据权利要求7所述的盲均衡系统，其特征在于： 所述输入判决单元包括：输入判决子单元，用于检测输入信号Xk得到中间元素判决值， 即为输入判决信息ik，公式为Jk=Clec(Xtl),其中L=(N+l)/2, dec( ·)为判决函数； 所述输出判决单元包括：输出判决子单元，用于检测输出信号yk得到判决值，即为输出 判决信息 dk，dk=dec (yk)。9. 根据权利要求8所述的盲均衡系统，其特征在于： 所述盲均衡单元包括：第一盲均衡子单元，或第二盲均衡子单元； 所述第一盲均衡子单元用于在期望信号未知的情况下，利用所述输入判决信息ik和所 述输出判决信息dk联合设计盲均衡算法，获得反馈误差ek : ek= (dk-yk) + a k I dk-yk | (ik-yk), 式中ak为自适应系数，随收敛过程逐渐变化： α k= Y Ek Ek+i=3 Ek+(l-3 ) (dk_yk)2 其中线性权重Y和遗忘参数β为经验正值； 第二盲均衡子单元用于在期望信号未知的情况下，利用所述输入判决信息ik和所述输 出判决信息dk联合设计盲均衡算法，获得反馈误差ek :10. 根据权利要求9所述的盲均衡系统，其特征在于： 所述反馈更新单元包括：更新子单元，用于利用下述公式更新递归更新参数： 巾々=12.… 其中U为更新步进。
【专利摘要】本发明提供一种盲均衡方法和一种盲均衡系统，包括：将输入信号经输入变换后生成输入变换后矢量；利用所述输入变换后矢量和递归更新参数，通过前馈网络获取输出信号；在期望信号未知的情况下，利用所述输入信号获得输入判决信息，利用所述输出信号获得输出判决信息，利用所述输入判决信息和所述输出判决信息联合设计盲均衡算法，获得反馈误差；利用所述反馈误差和输入变换后矢量，通过反馈过程更新递归更新参数。本发明加入了输入判决信息来辅助反馈过程，能够有效避免盲均衡中受前馈网络结构及递归更新参数初始设置值影响而常出现的局部收敛的问题，保证了盲均衡结果的稳健收敛；同时也保证了盲均衡的性能。
【IPC分类】H04L25/03
【公开号】CN104883330
【申请号】CN201410069372
【发明人】徐璐, 陈金树, 陆建华
【申请人】清华大学
【公开日】2015年9月2日
【申请日】2014年2月27日