[0078] 下面结合附图对本发明的技术方案做进一步的详细说明:
[0079] 将复数神经多项式网络和模糊神经网络控制器相结合,应用到卫星信道盲均衡方法中,本发明实现了一种卫星信道复数神经多项式网络盲均衡方法。
[0080] 本发明在充分利用CNP和FNN优点的基础上,利用复数神经多项式网络中单层神经网络实现对卫星信道线性特性进行补偿、用非线性处理器对卫星信道非线性特性进行处理;由模糊神经网络控制器自动控制复数神经多项式网络盲均衡方法权向量更新时的迭代步长,从而提高盲均衡方法的性能。
[0081] 本发明提供一种卫星信道复数神经多项式网络盲均衡系统,如图1所示。该系统由信道模块、复数神经多项式网络盲均衡模块和模糊神经网络控制器模块构成,其中,复数神经多项式网络盲均衡模块由单层神经网络、活化约束模块和非线性无记忆处理器构成。单层神经网络用于对卫星信道线性特性进行补偿、非线性无记忆处理器对于卫星信道的非线性特性进行补偿、活化约束利用非线性约束函数将非线性无记忆处理器输入限定在一定区域内,以保证非线性无记忆处理器的收敛性和稳定性。本发明中,信道模块采用Volterra卫星信道模型,一定区域为非线性约束函数的最小值与最大值区域。
[0082] 该系统中,卫星信号通过信道模块后作为复数神经多项式网络盲均衡模块中单层神经网络的输入,由单层神经网络得到的输出信号经活化约束模块约束在一定区域内,将活化约束模块约束后的信号作为非线性无记忆处理器输入信号,该信号经非线性无记忆处理器处理后得到复数神经多项式网络盲均衡模块输出信号。卫星信号在系统传输过程中,模糊神经网络控制器自动控制复数神经多项式网络盲均衡模块中单层神经网络权向量和非线性无记忆处理器权向量更新时的迭代步长,从而提高盲均衡方法的性能。
[0083] 图1中,a(n)为发射信号,c(n)为卫星信道、w(n)为信道噪声;Z-1表示延迟一个单位,x(n)表示单层神经网络的输入信号,fm(n)表示单层神经网络的第m个权系数,m=1,2,…M,M表示单层神经网络权向量长度;Y(n)为单层神经网络的输出信号,ψ(·)表示活化约束函数,YR(n)和YI(n)分别是Y(n)的实部与虚部,y(n)是活化约束函数的映射或称为非线性无记忆处理器的输入信号;hl(n)表示非线性无记忆处理器的第l个权系数,0≤l≤L-1,L为正整数,是非线性无记忆处理器权向量长度;e(n)是误差函数,S(n)是非线性无记忆处理器的输出信号或称为复数神经多项式网络盲均衡模块的输出信号。
[0084] 基于上述的一种卫星信道复数神经多项式网络盲均衡系统,本发明还提供一种卫星信道复数神经多项式网络盲均衡方法,具体如下:
[0085] 由图1得,单层神经网络的输入信号为
[0086] x(n)=c(n)a(n)+w(n) (1)
[0087] 单层神经网络的权向量为
[0088] F(n)=[f1(n),f2(n),…fM(n)] (2)
[0089] 其中,fm(n)为单层神经网络的第m个权系数。
[0090] 单层神经网络的输入向量为
[0091] X(n)=[x(n),x(n-1),…,x(n-M)] (3)
[0092] 单层神经网络的输出信号为
[0093] Y(n)=f1(n)x(n)+f2(n)x(n-1)+…fM(n)x(n-M)=F(n)XT(n) (4)
[0094] 其中,T表示转置,F(n)=[f1(n),f2(n),…fM(n)]。
[0095] 单层神经网络的输出信号Y(n)经活化约束函数ψ(·)映射后,得到非线性无记忆处理器的输入信号为
[0096] y(n)=ψ(Y(n))=ψ(YR(n))+jψ(YI(n)) (5)
[0097] 其中, 为虚数单位,活化约束函数ψ(·)为
[0098]
[0099] 其中,γ代表YR(n)或YI(n)
[0100] 非线性无记忆处理器的输入信号y(n)经非线性无记忆处理器g处理后,得其输出信号,也就是复数神经多项式网络盲均衡模块的输出信号为
[0101] S(n)=g(y(n))=h1(n)+h2(n)y(n)+…+hL(n)yL-1(n)=H(n)ZT(n) (7)[0102] 其中,H(n)是非线性无记忆处理器权向量,Z(n)是非线性无记忆处理器的输入向量,且
[0103] H(n)=[h1(n),h2(n),…,hL(n)] (8)
[0104] Z(n)=[1,y(n),y2(n),…,yL-1(n)] (9)
[0105] 为了获得式(2)或式(8)的更新方法,误差函数为
[0106] e(n)=|S(n)|2-R (10)
[0107] 式中,R为模值,定义为
[0108]
[0109] 式中,E表示数学期望。
[0110] 由误差函数定义常模代价函数为
[0111]
[0112] 根据复值原则,将权向量F(n)表示为
[0113] F(n)=FR(n)+jFI(n) (13)
[0114] 式中,FR(n)为F(n)的实部,FI(n)为F(n)的虚部。
[0115] J(n)分别对FR(n)和FI(n)求偏导,得
[0116]
[0117] 同理,得
[0118]
[0119]
[0120] 式中,Z′(n)为Z(n)的导数,Z′(n)=[0,1,2y(n),…,(L-1)yL-2(n)];*表示取共轭;FR(n)和FI(n)分别是F(n)的实部与虚部;ψ′(·)为ψ(·)的导数。
[0121] 采用随机梯度法,得FR(n)和FI(n)的更新方法分别为:
[0122]
[0123] 因此,得到单层神经网络的权向量更新方法为
[0124]
[0125] 式中,Z′(n)为Z(n)的导数,FR(0)、FI(0)与F(0)均采用中心抽头初始化。
[0126] 同理,对H(n)求偏导得,
[0127]
[0128] 因此,得到非线性处理器的权系数迭代公式为
[0129] H(n+1)=H(n)+μ2(n)e(n)Z*(n)S(n) (19)
[0130] 式中,H(0)采用中心抽头初始化,μ1(n)和μ2(n)分别为F(n)和H(n)的迭代步长,均由模糊神经网络控制器进行控制。
[0131] 在模糊神经网络控制器的设计过程中,先将模糊神经网络控制器的模糊控制规则定义如下:
[0132] 规则1:如果ΔMSE(n)为正,且MSE(n)大,则Δα(n)大;
[0133] 规则2:如果ΔMSE(n)为正,且MSE(n)中,则Δα(n)中大;
[0134] 规则3:如果ΔMSE(n)为正,且MSE(n)小,则Δα(n)中;
[0135] 规则4:如果ΔMSE(n)为零,且MSE(n)大,则Δα(n)中;
[0136] 规则5:如果ΔMSE(n)为零,且MSE(n)中,则Δα(n)中;
[0137] 规则6:如果ΔMSE(n)为零,且MSE(n)小,则Δα(n)中;
[0138] 规则7:如果ΔMSE(n)为负,且MSE(n)大,则Δα(n)中;
[0139] 规则8:如果ΔMSE(n)为负,且MSE(n)中,则Δα(n)中小;
[0140] 规则9:如果ΔMSE(n)为负,且MSE(n)小,则Δα(n)小;
[0141] 其中,MSE(n)是均方误差,ΔMSE(n)是均方误差增量,α(n)表示μ1(n)或μ2(n),Δα(n)是α(n)的变化量。
[0142] 根据上述模糊控制规则,设计了如图2所示的模糊神经网络控制器,该模糊神经网络控制器共有六层:第一层为输入层,有2个节点;第二层为模糊化层,有6个节点,分成两组,每组3个节点;第三层为规则层,有9个节点;第四层为模糊运算层,有5节点;第五层为误差引入层,有5节点;第六层为解模糊输出层,有1节点。
[0143] 模糊神经网络控制器各层的处理过程如下:
[0144] 第一层即输入层,将ΔMSE(n)和MSE(n)作为模糊神经网络控制器的输入变量,α1(n,1)=ΔMSE(n),α1(n,2)=MSE(n),O1(n,i1,j1)=α1(n,i1),i1=1,2,j1=1,2,3,其中,α1(n,1)表示n时刻第一层第一个节点的输入,α1(n,2)表示n时刻第一层第二个节点的输入,O1(n,i1,j1)为n时刻第一层第i1个节点第j1个输出,ΔMSE(n)=MSE(n)-MSE(n-1)。
[0145] 第二层即模糊化层,将ΔMSE(n)和MSE(n)划分为三个模糊等级,分别为:
[0146] α2(n,i2,j2)=O1(n,i1,j1) (20)
[0147]
[0148] 式中,m(i2,j2)表示n时刻第2层第i2个节点输入和第j2个模糊域的高斯函数中心,和q(i2,j2)表示n时刻第2层第i2个节点输入和第j2个模糊域的宽带;O2(n,i2,j2)表示n时刻第2层第i2个节点输入和第j2个输出;,i2=1,2,3,j2=1,2,…,9,α2(n,i2,j2)为第2层第i2个节点输入和第j2个输入。
[0149] 第三层即规则层,此层节点输出为第2层节点输出之积,其中:
[0150] O3(n,1)=O2(n,1,1)*O2(n,2,1);
[0151] O3(n,2)=O2(n,1,2)*O2(n,2,2);
[0152] O3(n,3)=O2(n,1,3)*O2(n,2,3);
[0153] O3(n,4)=O2(n,1,2)*O2(n,2,1);
[0154] O3(n,5)=O2(n,1,2)*O2(n,2,2);
[0155] O3(n,6)=O2(n,1,2)*O2(n,2,3);
[0156] O3(n,7)=O2(n,1,3)*O2(n,2,1);
[0157] O3(n,8)=O2(n,1,3)*O2(n,2,2);
[0158] O3(n,9)=O2(n,1,3)*O2(n,2,3);
[0159] 式中,O3(n,j2)表示n时刻第3层第j2个节点的输出。
[0160] 第四层即模糊运算层,该层节点输出由第三层节点输出按模糊控制规则运算获得,即
[0161] O4(n,1)=O3(n,1);
[0162] O4(n,2)=O3(n,2);
[0163]
[0164] O4(n,4)=O3(n,8);
[0165] O4(n,5)=O3(n,9);
[0166] 式中,O4(n,p)表示第n时刻第4层第p个的输出,p=1,2,3,4,5。
[0167] 第五层即误差引入层,引入误差MSE(n)后,得
[0168] O5(n,p)=β(p)*O4(n,p)*MSE(n) (22)
[0169] 式中,O5(n,p)表示第n时刻第5层第p个的输出;第5层的所有输出构成向量O5(n)={O5(n,1),O5(n,2),…,O5(n,5)};β(p)为控制量,主要用来调整该层的输出,完成规则的后部分。
[0170] 第六层即解模糊输出层:
[0171] Δα(n)=O6(n)=w6(n)O5(n) (23)
[0172] 式中,O6(n)为第n时刻第6层的输出;w6(n)为第n时刻第六层的向量,w6(n)=[w61(n),w62(n),…,w65(n)],其中,w6p(n)表示从第5层第p个输出到第6层节点支路上的权系数,w6(n)的更新方法为固定步长常模方法,具体为w6(n+1)=w6(n)-η(|Δα(n)|2-R1)Δα(n)O5(n),η为w6(n)的迭代步长,初始权向量w6(0)为中心抽头初始化; 为均方误差MSE(n)的模值;
[0173] 模糊神经网络控制器的第六层输出Δα(n),用来对α(n)进行更新,更新方法为[0174] α(n+1)=α(n)-|Δα(n)|
[0175] 式中,α(n)的初值α(0)=1,Δα(n)的初值Δα(0)=0。
[0176] 这样构成由模糊神经网络控制器控制单层神经网络和非线性无记忆处理器迭代步长的方法,从而发明了一种卫星信道复数神经多项式网络盲均衡盲均衡方法FNN-CNP-CMA(fuzzy neural network and complex-valued neural polynomial based constant modulus algorithm)。
[0177] 本发明方法采用复数神经多项式网络结构的均衡器,计算复度得到有效降低。以处理一帧数据所需要的加法和乘法次数来表示的计算复杂度。M和L分别表示单层神经网络权向量长度和非线性无记忆处理器权向量长度。由于进行一次乘法比进行一次加法所需运算量大得多,所以以每产生一次输出信号所需进行的乘法次数作为计算量的衡量标准。单层神经网络进行一次权系数更新所需的乘法次数为4L2+8M+8L;非线性无记忆处理器进行2
一次权系数更新所需的乘法次数为L +L+1。所以CNP-CMA进行一次权系数更新所需要的总的乘法次数为5L2+9L+8M+1。同样可得本发明方法进行一次权系数更新所需要的总的乘法次数为5L2+10L+9M+1。与CNP-CMA相比,本发明方法FNN-CNP-CMA每次迭代过程仅增加了M+L次乘法运算,计算复杂度没有变化。
[0178] 为了验证本发明方法FNN-CNP-CMA的有效性,采用最典型的非线性Volterra卫星信道进行仿真实验,并与常规的常模方法CMA、复数神经多项式网络盲均衡CNP-CMA(complex-valued neural polynomial based constant modulus algorithm)进行比较。
[0179] 实验中测试信道为
[0180]
[0181] 式(24)中的信道参数,如表1所示。发射信号为32-APSK,信噪比为25,50次蒙特卡罗仿真结果,如图3和图4所示。
[0182] 表1信道Volterra核
[0183]
[0184] 图3为均方误差收敛曲线,图4为发射信号为4+12+16+32-APSK时的均衡效果。图3和图4表明,与CMA和CNP-CMA相比,本发明方法FNN-CNP-CMA收敛速度分别加快了5000步和800步;而本发明方法FNN-CNP-CMA的稳态误差分别降低了11dB和7dB;本发明方法FNN-CNP-CMA的输出星座图比CMA和FNN-CMA的输出星座图更加清晰、紧凑,眼图张开更加清晰,均衡器输出信号更接近发射信号,均衡效果更好。
[0185] 以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。