首页 > 专利 > 南京信息工程大学 > 一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法专利详情

一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2015-07-20
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2015-12-16
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2018-01-05
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2035-07-20
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201510427082.2 申请日 2015-07-20
公开/公告号 CN105072064B 公开/公告日 2018-01-05
授权日 2018-01-05 预估到期日 2035-07-20
申请年 2015年 公开/公告年 2018年
缴费截止日
分类号 H04L25/03G06N3/12 主分类号 H04L25/03
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 3
权利要求数量 4 非专利引证数量 1
引用专利数量 4 被引证专利数量 0
非专利引证 1、黄伟.基于人工鱼群优化的小波盲均衡算法《.万方数据知识服务平台学术论文数据库》.2014,;
引用专利 CN102299875A、CN103346987A、CN103338170A、WO2005101655A1 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 南京信息工程大学 当前专利权人 南京信息工程大学
发明人 郭业才、吴华鹏、王惠 第一发明人 郭业才
地址 江苏省南京市宁六路219号 邮编 210044
申请人数量 1 发明人数量 3
申请人所在省 江苏省 申请人所在市 江苏省南京市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
南京经纬专利商标代理有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
许方
摘要
本发明公开了一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法(DNA‑GBA‑FS‑MMA),本发明将DNA遗传蝙蝠方法与分数间隔多模盲均衡方法(FS‑MMA)相结合,充分利用FS‑MMA和DNA遗传蝙蝠方法的优点,本发明将DNA遗传方法引入到蝙蝠方法中,得到一种DNA遗传蝙蝠方法(DNA‑GBA),利用这个新的优化方法来寻找蝙蝠群的全局最优位置向量,并作为多模盲均衡方法(MMA)初始化最优权向量的实部与虚部来提高收敛速度、减小剩余均方误差。与基于蝙蝠方法的分数间隔多模盲均衡方法(BA‑FS‑MMA)和基于蝙蝠方法的多模盲均衡方法(BA‑MMA)相比,本发明具有收敛速度最快、均方误差最小和全局最优的性能,在通信技术领域有很强的实用价值。
  • 摘要附图
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
  • 说明书附图:图1
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
  • 说明书附图:图2
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
  • 说明书附图:图3
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
  • 说明书附图:图4
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
  • 说明书附图:图5
    一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2018-01-05 授权
2 2015-12-16 实质审查的生效 IPC(主分类): H04L 25/03 专利申请号: 201510427082.2 申请日: 2015.07.20
3 2015-11-18 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,其特征在于,包括以下步骤:
步骤1、将发射信号a(k)经过第m个子信道cm(k)后加入第m个子信道噪声nm(k),得到第m个子信道多模盲均衡器输入信号ym(k),ym(k)=a(k)cm(k)+nm(k),其中,k为正整数且表示时间序列,m为整数;
步骤2、将步骤1所述的第m个子信道多模盲均衡器输入信号ym(k)的实部ymR(k)和虚部ymI(k)分别经过第m个实部均衡器得到输出信号zmR(k)、第m个虚部均衡器得到输出信号zmI(k);其中,zmR(k)=wmR(k)ymR(k),zmI(k)=wmI(k)ymI(k),wmR(k)为第m个盲均衡器权向量wm(k)的实部,wmI(k)为第m个盲均衡器权向量wm(k)的虚部,且wmR(k)和wmI(k)的更新公式为wmR(k+1)=wmR(k)-4μemR(k)ymR(k)
wmI(k+1)=wmI(k)-4μemI(k)ymI(k)
其中,emR(k)为误差信号的实部,emI(k)为误差信号的虚部,emR(k)和emI(k)是由误差生成函数生成,μ为迭代步长且是实数,0≤μ<1;
所述步骤2中盲均衡器权向量的实部wmR(k)和虚部wmI(k)的初始优化权向量的wmR(0)和wmI(0)是由DNA遗传蝙蝠方法获取,具体步骤如下:
步骤2-1、设置蝙蝠初始种群及参数:
设蝙蝠初始种群X=[X1,X2,…,XN],其中,Xi表示第i只蝙蝠的初始位置向量,1≤i≤N,N为蝙蝠个体数量,频率范围为[fmin,fmax],最大响度为A(0),最大频度为r(0),最大频度r(0)∈[0,1],最大响度A(0)∈[1,2],响度衰减系数为α,频度增加系数为γ,γ是大于0的常数,置换交叉概率pz,移位交叉概率py,变异概率pb,维数为d,搜索精度为tol,最大迭代次数为iter,fmin为最低频率,fmax为最高频率;
步骤2-2、确定适应度函数:
将多模盲均衡方法MMA的代价函数的倒数定义为DNA遗传蝙蝠方法的适应度函数,即其中,Xi(k)为DNA遗传蝙蝠方法中第i只蝙蝠k时刻的位置向量,用十进制数表示;JMMA(k)为MMA的代价函数,E{*}表示数学期望,zR(k)和zI(k)为盲均衡器输出信号的实部和虚部, 和 为统计模值的实部和虚部;
步骤2-3、计算每只蝙蝠位置向量的适应度函数值并将适应度函数值从大到小排列,其中,前一半对应的蝙蝠个体组成优质种群,后一半对应的蝙蝠个体组成劣质种群;适应度值最大的位置向量为当前全局最优位置向量X*;
步骤2-4、调整第i只蝙蝠超声波的频率fi,对第i只蝙蝠k-1时刻的速度Vi(k-1)和位置向量Xi(k-1)进行更新,得到第i只蝙蝠k时刻的位置向量Xi(k);
fi=fmin+(fmax-fmin)β
Vi(k)=Vi(k-1)+[Xi(k)-X*]fi
Xi(k)=Xi(k-1)+Vi(k)
其中,β为[0,1]上的随机数;
步骤2-5、产生一个服从均匀分布的随机频度rand1与第i只蝙蝠的频度ri进行比较,若rand1>ri,则对当前全局最优位置向量X*进行随机扰动,产生一个新的位置向量Xnew,用其代替第i只蝙蝠k时刻的位置向量Xi(k);
Xnew=X*+εA(k)
Xi(k)=Xnew
其中,ε为与Xi(k)维数相同的随机向量,随机向量中分量的取值为[-1,1]上的随机数,A(k)为k时刻蝙蝠群的平均响度;
步骤2-6、产生一个服从均匀分布的随机响度rand2与第i只蝙蝠的响度Ai进行比较,若*
rand2<Ai且JDNA-GBA(Xi(k))>JDNA-GBA(X),则用第i只蝙蝠k时刻的位置向量Xi(k)替代当前全局最优位置向量X*,并对第i只蝙蝠的响度Ai和频度ri分别进行更新;
Ai(k+1)=αAi(k)
ri(k+1)=r(0)[1-exp(-γk)]
其中,初始化时,所有蝙蝠的响度和频度都不相同;
步骤2-7、编码:
将各蝙蝠个体的十进制位置向量转换成四进制,得到DNA序列位置向量;
步骤2-8、交叉操作和变异操作:
产生一个随机数rand3∈(0,1),与置换交叉概率pz比较,若rand3步骤2-9、将经交叉、变异后得到的所有蝙蝠个体的DNA序列位置向量解码,用解码得到的位置向量计算适应度函数值,且从小到大排列并划分优质种群和劣质种群;
步骤2-10、选取适应度函数值最大的位置向量为当前全局最佳位置向量;
步骤2-11、达到最大迭代次数或搜索精度,则输出全局最优位置向量 否则转至步骤2-3,继续搜索;
步骤3、将步骤2所述的第m个实部均衡器输出信号zmR(k)和第m个虚部均衡器输出信号zmI(k)相加得到第m个子盲均衡器输出信号zm(k):zm(k)=zmR(k)+j·zmI(k),其中,j为虚数单位;
步骤4、根据步骤3依次求得所有子盲均衡器输出信号,求和得到盲均衡器的总输出信号z(k): 其中,M是子信道的个数。

2.根据权利要求1所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,其特征在于:所述步骤2-7中的编码,具体步骤如下:
步骤2-7-1、由第i只蝙蝠的位置向量Xi=[xi1,xi2,…,xiP]计算得到十进制位置过渡向量Bi=[bi1,bi2,…,biP],其中,xig表示第i只蝙蝠的位置向量Xi中第g个位置值,big表示十进制位置过渡向量中第g个位置值,g为正整数且1≤g≤P,P为十进制位置向量的维数,l为编码长度,dmaxg和dming分别为第i只蝙蝠的位置向量Xi中第g个位
置的最大值、最小值;
步骤2-7-2、将十进制位置过渡向量中第g个位置值big转换成一串四进制数sig,则第i只蝙蝠的DNA序列位置向量 由P串四进制数sig组成,其中,sig
表示第i只蝙蝠的DNA序列位置向量Si中第g个位置的整数串,长度为l, 表示第i只蝙蝠的DNA序列位置向量Si中第g个子整数串中第n位的数字,n为正整数且1≤n≤l。

3.根据权利要求1所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,其特征在于,所述步骤2-8中的交叉操作包括置换交叉操作和转位交叉操作;所述置换交叉操作、转位交叉操作、变异操作的具体如下:
步骤2-8-1、置换交叉操作:从DNA序列位置向量中任意选取两个作为父体,从两个父体中分别随机选取一段碱基数目相等的转座子序列位置向量进行交换,产生两个新的序列位置向量代替父体;选取父体时,一个从整个种群选取,一个从优质种群选取;
步骤2-8-2、转位交叉操作:从DNA序列位置向量中任意选取一个作为父体,从父体中截取一段任意碱基数目的转座子序列位置向量插入截取后父体的任意位置,产生一个新的序列位置向量代替父体;选取父体时,80%的概率从优质种群选取,20%概率从劣质种群选取;
步骤2-8-3、变异操作:从DNA序列位置向量中任意选取一个序列位置向量,该序列位置向量中任一碱基都有一定概率被其它碱基代替,从而产生一个新的序列位置向量。

4.根据权利要求2所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,其特征在于,所述步骤2-9中DNA序列位置向量解码,具体如下:
步骤2-9-1、将第i只蝙蝠的DNA序列位置向量 解码为
十进制位置过渡向量Bi=[bi1,bi2,…,biP],
步骤2-9-2、通过以下公式按比例将big转换成第i只蝙蝠的位置向量Xi中第g个位置值xig;
说明书

技术领域

[0001] 本发明涉及盲均衡技术领域,特别是一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法。

背景技术

[0002] 多模盲均衡方法(MMA,Multi-Modulus Algorithm)可以在不使用独立载波恢复系统的情况下同时实现盲均衡和载波相位恢复。与常模盲均衡方法(CMA,Constant Modulus Algorithm)相比,MMA收敛速度更快、稳态误差更小,可以有效地均衡多模信号。分数间隔均衡器(FSE,Fractionally Spaced Equalizer)的主要思想是对信号进行过采样,有效减小了盲均衡器的权长、获取了更多信道信息,有利于补偿信道失真,并恢复输入信号。将FSE与MMA相结合,形成的分数间隔多模盲均衡方法(FS-MMA,Fractionally Spaced Multi-Modulus Algorithm)可以减小稳态误差、减少计算量。蝙蝠方法(BA,Bat Algorithm)是一种基于种群的随机全局寻优方法,通过改变蝙蝠发出的超声波的频率、频度、响度搜寻全局最优位置,并利用自身特有的回波定位特性使搜索过程避免陷入局部搜索,提高了搜索全局最优位置的成功率,将搜索所得全局最优位置向量同时作为MMA初始权向量的实部与虚部,与普通MMA的中心抽头初始权向量相比,该全局最优位置向量能使MMA尽快达到收敛状态,且稳态误差受调制阶数的影响大大减小。因此,用搜索所得最优位置向量优化初始权向量,可以极大地加快收敛速度、减小稳态误差。
[0003] 如果能将DNA遗传蝙蝠方法应用到盲均衡方法的均衡器权向量优化中,则有利于获得更快的收敛速度和更低的稳态误差,但现有技术中尚缺乏将DNA遗传蝙蝠方法与盲均衡方法相结合的应用模型。

发明内容

[0004] 本发明所要解决的技术问题是克服现有技术的不足而提供一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,本发明将DNA遗传方法引入到蝙蝠方法中,得到一种DNA遗传蝙蝠方法(DNA-GBA),利用这个新的优化方法来寻找蝙蝠群的全局最优位置向量,并作为多模盲均衡方法(MMA)初始化最优权向量的实部与虚部来提高收敛速度、减小剩余均方误差。
[0005] 本发明为解决上述技术问题采用以下技术方案:
[0006] 根据本发明提出的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,包括以下步骤:
[0007] 步骤1、将发射信号a(k)经过第m个子信道cm(k)后加入第m个子信道噪声nm(k),得到第m个子信道多模盲均衡器输入信号ym(k),ym(k)=a(k)cm(k)+nm(k),其中,k为正整数且表示时间序列,m为整数;
[0008] 步骤2、将步骤1所述的第m个子信道多模盲均衡器输入信号ym(k)的实部ymR(k)和虚部ymI(k)分别经过第m个实部均衡器得到输出信号zmR(k)、第m个虚部均衡器得到输出信号zmI(k);其中,zmR(k)=wmR(k)ymR(k),zmI(k)=wmI(k)ymI(k),wmR(k)为第m个盲均衡器权向量wm(k)的实部,wmI(k)为第m个盲均衡器权向量wm(k)的虚部,且wmR(k)和wmI(k)的更新公式为
[0009] wmR(k+1)=wmR(k)-4μemR(k)ymR(k)
[0010] wmI(k+1)=wmI(k)-4μemI(k)ymI(k)
[0011] 其中,emR(k)为误差信号的实部,emI(k)为误差信号的虚部,emR(k)和emI(k)是由误差生成函数生成,μ为迭代步长且是实数,0≤μ<1;
[0012] 步骤3、将步骤2所述的第m个实部均衡器输出信号zmR(k)和第m个虚部均衡器输出信号zmI(k)相加得到第m个子盲均衡器输出信号zm(k):zm(k)=zmR(k)+j·zmI(k),其中,j为虚数单位;
[0013] 步骤4、根据步骤3依次求得所有子盲均衡器输出信号,求和得到盲均衡器的总输出信号z(k): 其中,M是子信道的个数。
[0014] 作为本发明所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法进一步优化方案,所述步骤2中盲均衡器权向量的实部wmR(k)和虚部wmI(k)的初始优化权向量的wmR(0)和wmI(0)是由DNA遗传蝙蝠方法获取,具体步骤如下:
[0015] 步骤2-1、设置蝙蝠初始种群及参数:
[0016] 设蝙蝠初始种群X=[X1,X2,…,XN],其中,Xi表示第i只蝙蝠的初始位置向量,1≤i≤N,N为蝙蝠个体数量,频率范围为[fmin,fmax],最大响度为A(0),最大频度为r(0),最大频度r(0)∈[0,1],最大响度A(0)∈[1,2],响度衰减系数为α,频度增加系数为γ,γ是大于0的常数,置换交叉概率pz,移位交叉概率py,变异概率pb,维数为d,搜索精度为tol,最大迭代次数为iter,fmin为最低频率,fmax为最高频率;
[0017] 步骤2-2、确定适应度函数:
[0018] 将多模盲均衡方法MMA的代价函数的倒数定义为DNA遗传蝙蝠方法的适应度函数,即
[0019]
[0020] 其中,Xi(k)为DNA遗传蝙蝠方法中第i只蝙蝠k时刻的位置向量,用十进制数表示;JMMA(k)为MMA的代价函数,E{*}表示数学期望,zR(k)和zI(k)为盲均衡器输出信号的实部和虚部, 和 为统计模值的实部和虚部;
[0021] 步骤2-3、计算每只蝙蝠位置向量的适应度函数值并将适应度函数值从大到小排列,其中,前一半对应的蝙蝠个体组成优质种群,后一半对应的蝙蝠个体组成劣质种群;适应度值最大的位置向量为当前全局最优位置向量X*;
[0022] 步骤2-4、调整第i只蝙蝠超声波的频率fi,对第i只蝙蝠k-1时刻的速度Vi(k-1)和位置向量Xi(k-1)进行更新,得到第i只蝙蝠k时刻的位置向量Xi(k);
[0023] fi=fmin+(fmax-fmin)β
[0024] Vi(k)=Vi(k-1)+[Xi(k)-X*]fi
[0025] Xi(k)=Xi(k-1)+Vi(k)
[0026] 其中,β为[0,1]上的随机数;
[0027] 步骤2-5、产生一个服从均匀分布的随机频度rand1与第i只蝙蝠的频度ri进行比较,若rand1>ri,则对当前全局最优位置向量X*进行随机扰动,产生一个新的位置向量Xnew,用其代替第i只蝙蝠k时刻的位置向量Xi(k);
[0028] Xnew=X*+εA(k)
[0029] Xi(k)=Xnew
[0030] 其中,ε为与Xi(k)维数相同的随机向量,随机向量中分量的取值为[-1,1]上的随机数,A(k)为k时刻蝙蝠群的平均响度;
[0031] 步骤2-6、产生一个服从均匀分布的随机响度rand2与第i只蝙蝠的响度Ai进行比较,若rand2<Ai且JDNA-GBA(Xi(k))>JDNA-GBA(X*),则用第i只蝙蝠k时刻的位置向量Xi(k)替代当前全局最优位置向量X*,并对第i只蝙蝠的响度Ai和频度ri分别进行更新;
[0032] Ai(k+1)=αAi(k)
[0033] ri(k+1)=r(0)[1-exp(-γk)]
[0034] 其中,初始化时,所有蝙蝠的响度和频度都不相同;
[0035] 步骤2-7、编码:
[0036] 将各蝙蝠个体的十进制位置向量转换成四进制,得到DNA序列位置向量;
[0037] 步骤2-8、交叉操作和变异操作:
[0038] 产生一个随机数rand3∈(0,1),与置换交叉概率pz比较,若rand3
[0039] 步骤2-9、将经交叉、变异后得到的所有蝙蝠个体的DNA序列位置向量解码,用解码得到的位置向量计算适应度函数值,且从小到大排列并划分优质种群和劣质种群;
[0040] 步骤2-10、选取适应度函数值最大的位置向量为当前全局最佳位置向量;
[0041] 步骤2-11、达到最大迭代次数或搜索精度,则输出全局最优位置向量 否则转至步骤2-3,继续搜索。
[0042] 作为本发明所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法进一步优化方案,所述步骤2-7中的编码,具体步骤如下:
[0043] 步骤2-7-1、由第i只蝙蝠的位置向量Xi=[xi1,xi2,…,xiP]计算得到十进制位置过渡向量Bi=[bi1,bi2,…,biP],其中,xig表示第i只蝙蝠的位置向量Xi中第g个位置值,big表示十进制位置过渡向量中第g个位置值,g为正整数且1≤g≤P,P为十进制位置向量的维数,l为编码长度,dmaxg和dming分别为第i只蝙蝠的位置向量Xi中第g个位置的最大值、最小值;
[0044] 步骤2-7-2、将十进制位置过渡向量中第g个位置值big转换成一串四进制数sig,则第i只蝙蝠的DNA序列位置向量 由P串四进制数sig组成,其中,sig表示第i只蝙蝠的DNA序列位置向量Si中第g个位置的整数串,长度为l, 表示第i只蝙蝠的DNA序列位置向量Si中第g个子整数串中第n位的数字,n为正整数且1≤n≤l。
[0045] 作为本发明所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法进一步优化方案,所述步骤2-8中的交叉操作包括置换交叉操作和转位交叉操作;所述置换交叉操作、转位交叉操作、变异操作的具体如下:
[0046] 步骤2-8-1、置换交叉操作:从DNA序列位置向量中任意选取两个作为父体,从两个父体中分别随机选取一段碱基数目相等的转座子序列位置向量进行交换,产生两个新的序列位置向量代替父体;选取父体时,一个从整个种群选取,一个从优质种群选取;
[0047] 步骤2-8-2、转位交叉操作:从DNA序列位置向量中任意选取一个作为父体,从父体中截取一段任意碱基数目的转座子序列位置向量插入截取后父体的任意位置,产生一个新的序列位置向量代替父体;选取父体时,80%的概率从优质种群选取,20%概率从劣质种群选取;
[0048] 步骤2-8-3、变异操作:从DNA序列位置向量中任意选取一个序列位置向量,将该序列位置向量中任一碱基都有一定概率被其它碱基代替,从而产生一个新的序列位置向量。
[0049] 作为本发明所述的一种基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法进一步优化方案,所述步骤2-9中DNA序列位置向量解码,具体如下:
[0050] 步骤2-9-1、将第i只蝙蝠的DNA序列位置向量 解码为十进制位置过渡向量Bi=[bi1,bi2,…,biP],
[0051] 步骤2-9-2、通过以下公式按比例将big转换成第i只蝙蝠的位置向量Xi中第g个位置值xig;
[0052]
[0053] 本发明采用以上技术方案与现有技术相比,具有以下技术效果:
[0054] (1)本发明首先将均衡器的输入信号作为DNA遗传蝙蝠方法的输入,并且把MMA的代价函数经适当变换后作为DNA遗传蝙蝠方法的适应度函数,利用DNA遗传蝙蝠方法的寻优能力来寻找盲均衡器最优初始权向量;
[0055] (2)本发明将DNA遗传蝙蝠方法与分数间隔多模盲均衡方法(FS-MMA)相结合,充分利用FS-MMA和DNA遗传蝙蝠方法的优点,本发明将DNA遗传方法引入到蝙蝠方法中,得到一种DNA遗传蝙蝠方法(DNA-GBA),利用这个新的优化方法来寻找蝙蝠群的全局最优位置向量,并作为多模盲均衡方法(MMA)初始化最优权向量的实部与虚部来提高收敛速度、减小剩余均方误差;
[0056] (3)本发明方法与基于蝙蝠方法的分数间隔多模盲均衡方法(BA-FS-MMA)和基于蝙蝠方法的多模盲均衡方法(BA-MMA)相比,本发明具有收敛速度最快、均方误差最小和全局最优的性能,在通信技术领域有很强的实用价值。

实施方案

[0062] 下面结合附图对本发明的技术方案做进一步的详细说明:
[0063] 将DNA遗传蝙蝠方法应用到分数间隔多模盲均衡方法中,进一步详述如下:
[0064] 将均衡器的输入信号作为DNA遗传蝙蝠方法的输入,并且把MMA的代价函数经适当变换后作为DNA遗传蝙蝠方法的适应度函数,利用DNA遗传蝙蝠方法的寻优能力来寻找均衡器的最佳权向量。分数间隔多模盲均衡方法原理框图如图1所示。
[0065] 图中a(k)为发射信号序列,cm(k)为第m个子信道,nm(k)为第m条支路的加性高斯白噪声,ymR(k)和ymI(k)分别为多模模块m输入信号ym(k)的实部与虚部,wmR(k)和wmI(k)分别为多模模块m权向量wm(k)的实部与虚部,zmR(k)和zmI(k)分别为多模模块m输出信号zm(k)的实部与虚部,emR(k)和emI(k)分别为多模模块m误差函数em(k)的实部与虚部,z(k)为整个分数间隔多模盲均衡系统的输出信号。
[0066] 对于图1中的(a)中第m条支路(m=0,1,…,M-1),各物理量的关系为:
[0067] cm(k)=c[(k+1)M-m-1]
[0068] 式中,c为整个系统的信道。
[0069] 盲均衡器输入信号为
[0070] ym(k)=a(k)cm(k)+nm(k)
[0071] 将输入信号ym(k)分为实部与虚部分别进行处理,得到均衡器输出信号的实部与虚部分别为
[0072] zmR(k)=wmR(k)ymR(k)
[0073] zmI(k)=wmI(k)ymI(k)
[0074] 输出信号为
[0075] zm(k)=zmR(k)+j·zmI(k)
[0076] 误差信号的实部与虚部分别为
[0077]
[0078]
[0079] 式中,TR和TI分别为发射信号a(k)实部和虚部的统计模值,分别定义为
[0080]
[0081]
[0082] 图1中的(b)第m路MMA的代价函数定义为
[0083]
[0084] 按照最速下降法,有
[0085]
[0086] 因此,多模模块m权向量实部wmR(k)和虚部wmI(k)的迭代公式分别为
[0087] wmR(k+1)=wmR(k)-4μemR(k)ymR(k)
[0088] wmI(k+1)=wmI(k)-4μemI(k)ymI(k)
[0089] 式中,μ∈(0,1)为步长。
[0090] 分数间隔多模盲均衡方法(FS-MMA)的输出信号为
[0091]
[0092] 以上构成了分数间隔多模盲均衡方法。
[0093] 本发明原理如图2所示,首先发射信号a(k)经过第m个子信道cm(k)后加入第m个子信道噪声nm(k),将得到多模盲均衡器输入信号作为DNA遗传蝙蝠方法的输入,并且把MMA的代价函数经适当变换后作为DNA遗传蝙蝠方法的适应度函数,利用DNA遗传蝙蝠方法的寻优能力来寻找盲均衡器初始最优权向量,具体地说,本发明包括如下步骤:
[0094] 步骤1、将发射信号a(k)经过第m个子信道cm(k)后加入第m个子信道噪声nm(k),得到第m个子信道多模盲均衡器输入信号ym(k),ym(k)=a(k)cm(k)+nm(k),其中,k为正整数且表示时间序列,m为整数;
[0095] 步骤2、将步骤1所述的第m个子信道多模盲均衡器输入信号ym(k)的实部ymR(k)和虚部ymI(k)分别经过第m个实部均衡器得到输出信号zmR(k)、第m个虚部均衡器得到输出信号zmI(k);其中,zmR(k)=wmR(k)ymR(k),zmI(k)=wmI(k)ymI(k),wmR(k)为第m个盲均衡器权向量wm(k)的实部,wmI(k)为第m个盲均衡器权向量wm(k)的虚部,且wmR(k)和wmI(k)的更新公式为
[0096] wmR(k+1)=wmR(k)-4μemR(k)ymR(k)
[0097] wmI(k+1)=wmI(k)-4μemI(k)ymI(k)
[0098] 其中,emR(k)为误差信号的实部,emI(k)为误差信号的虚部,emR(k)和emI(k)是由误差生成函数生成,μ为迭代步长且是实数,0≤μ<1;
[0099] 步骤3、将步骤2所述的第m个实部均衡器输出信号zmR(k)和第m个虚部均衡器输出信号zmI(k)相加得到第m个子盲均衡器输出信号zm(k):zm(k)=zmR(k)+j·zmI(k),其中,j为虚数单位;
[0100] 步骤4、根据步骤3依次求得所有子盲均衡器输出信号,求和得到盲均衡器的总输出信号z(k): 其中,M是子信道的个数;
[0101] 所述步骤2中盲均衡器权向量的实部wmR(k)和虚部wmI(k)的初始优化权向量的wmR(0)和wmI(0)是由DNA遗传蝙蝠方法获取,具体步骤如下:
[0102] 步骤2-1、设置蝙蝠初始种群及参数:
[0103] 设蝙蝠初始种群X=[X1,X2,…,XN],其中,Xi表示第i只蝙蝠的初始位置向量,1≤i≤N,N为蝙蝠个体数量,频率范围为[fmin,fmax],最大响度为A(0),最大频度为r(0),最大频度r(0)∈[0,1],最大响度A(0)∈[1,2],响度衰减系数为α,频度增加系数为γ,γ是大于0的常数,置换交叉概率pz,移位交叉概率py,变异概率pb,维数为d,搜索精度为tol,最大迭代次数为iter,fmin为最低频率,fmax为最高频率;
[0104] 步骤2-2、确定适应度函数:
[0105] 将多模盲均衡方法MMA的代价函数的倒数定义为DNA遗传蝙蝠方法的适应度函数,即
[0106]
[0107] 其中,Xi(k)为DNA遗传蝙蝠方法中第i只蝙蝠k时刻的位置向量,用十进制数表示;JMMA(k)为MMA的代价函数,E{*}表示数学期望,zR(k)和zI(k)为盲均衡器输出信号的实部和虚部, 和 为统计模值的实部和虚部;
[0108] 步骤2-3、计算每只蝙蝠位置向量的适应度函数值并将适应度函数值从大到小排列,其中,前一半对应的蝙蝠个体组成优质种群,后一半对应的蝙蝠个体组成劣质种群;适应度值最大的位置向量为当前全局最优位置向量X*;
[0109] 步骤2-4、调整第i只蝙蝠超声波的频率fi,对第i只蝙蝠k-1时刻的速度Vi(k-1)和位置向量Xi(k-1)进行更新,得到第i只蝙蝠k时刻的位置向量Xi(k);
[0110] fi=fmin+(fmax-fmin)β
[0111] Vi(k)=Vi(k-1)+[Xi(k)-X*]fi
[0112] Xi(k)=Xi(k-1)+Vi(k)
[0113] 其中,β为[0,1]上的随机数;
[0114] 步骤2-5、产生一个服从均匀分布的随机频度rand1与第i只蝙蝠的频度ri进行比较,若rand1>ri,则对当前全局最优位置向量X*进行随机扰动,产生一个新的位置向量Xnew,用其代替第i只蝙蝠k时刻的位置向量Xi(k);
[0115] Xnew=X*+εA(k)
[0116] Xi(k)=Xnew
[0117] 其中,ε为与Xi(k)维数相同的随机向量,随机向量中分量的取值为[-1,1]上的随机数,A(k)为k时刻蝙蝠群的平均响度;
[0118] 步骤2-6、产生一个服从均匀分布的随机响度rand2与第i只蝙蝠的响度Ai进行比较,若rand2<Ai且JDNA-GBA(Xi(k))>JDNA-GBA(X*),则用第i只蝙蝠k时刻的位置向量Xi(k)替代当前全局最优位置向量X*,并对第i只蝙蝠的响度Ai和频度ri分别进行更新;
[0119] Ai(k+1)=αAi(k)
[0120] ri(k+1)=r(0)[1-exp(-γk)]
[0121] 其中,初始化时,所有蝙蝠的响度和频度都不相同;
[0122] 步骤2-7、编码:
[0123] 将各蝙蝠个体的十进制位置向量转换成四进制,得到DNA序列位置向量;
[0124] 步骤2-8、交叉操作和变异操作:
[0125] 产生一个随机数rand3∈(0,1),与置换交叉概率pz比较,若rand3
[0126] 步骤2-9、将经交叉、变异后得到的所有蝙蝠个体的DNA序列位置向量解码,用解码得到的位置向量计算适应度函数值,且从小到大排列并划分优质种群和劣质种群;
[0127] 步骤2-10、选取适应度函数值最大的位置向量为当前全局最佳位置向量;
[0128] 步骤2-11、达到最大迭代次数或搜索精度,则输出全局最优位置向量Xo*pt;否则转至步骤2-3,继续搜索。
[0129] 所述步骤2-7中的编码,具体步骤如下:
[0130] 步骤2-7-1、由第i只蝙蝠的位置向量Xi=[xi1,xi2,…,xiP]计算得到十进制位置过渡向量Bi=[bi1,bi2,…,biP],其中,xig表示第i只蝙蝠的位置向量Xi中第g个位置值,big表示十进制位置过渡向量中第g个位置值,g为正整数且1≤g≤P,P为十进制位置向量的维数,l为编码长度,dmaxg和dming分别为第i只蝙蝠的位置向量Xi中第g个位置的最大值、最小值;
[0131] 步骤2-7-2、将十进制位置过渡向量中第g个位置值big转换成一串四进制数sig,则第i只蝙蝠的DNA序列位置向量 由P串四进制数sig组成,其中,sig表示第i只蝙蝠的DNA序列位置向量Si中第g个位置的整数串,长度为l, 表示第i只蝙蝠的DNA序列位置向量Si中第g个子整数串中第n位的数字,n为正整数且1≤n≤l。
[0132] 所述步骤2-8中的交叉操作包括置换交叉操作和转位交叉操作;所述置换交叉操作、转位交叉操作、变异操作的具体如下:
[0133] 步骤2-8-1、置换交叉操作:从DNA序列位置向量中任意选取两个作为父体,从两个父体中分别随机选取一段碱基数目相等的转座子序列位置向量进行交换,产生两个新的序列位置向量代替父体;选取父体时,一个从整个种群选取,一个从优质种群选取;
[0134] 步骤2-8-2、转位交叉操作:从DNA序列位置向量中任意选取一个作为父体,从父体中截取一段任意碱基数目的转座子序列位置向量插入截取后父体的任意位置,产生一个新的序列位置向量代替父体;选取父体时,80%的概率从优质种群选取,20%概率从劣质种群选取;
[0135] 步骤2-8-3、变异操作:从DNA序列位置向量中任意选取一个序列位置向量,将该序列位置向量中任一碱基都有一定概率被其它碱基代替,从而产生一个新的序列位置向量。
[0136] 所述步骤2-9中DNA序列位置向量解码,具体如下:
[0137] 步骤2-9-1、将第i只蝙蝠的DNA序列位置向量 解码为十进制位置过渡向量Bi=[bi1,bi2,…,biP],
[0138] 步骤2-9-2、通过以下公式按比例将big转换成第i只蝙蝠的位置向量Xi中第g个位置值xig;
[0139]
[0140] 为了验证本发明DNA-GBA-FS-MMA的性能,将MMA、FS-MMA、BA-MMA、GBA-MMA、BA-FS-MMA与本发明DNA-GBA-FS-MMA进行对比,以T/4分数间隔为例进行仿真实验。
[0141] 参数初始化:每个蝙蝠种群数中蝙蝠个体的数量N=20,频率范围[0,100],最大响度A(0)=1.5,最大频度r(0)=0.25,搜索精度tol=10-5,维数d=11,响度衰减系数α=0.9,频度增加系数γ=0.9,置换交叉概率pz=0.8,移位交叉概率py=0.3,变异概率pb=0.2,最大迭代次数iter=2000,运行次数runs=2000,信道c=[0.9556-0.09060.05780.2368],信噪比SNR=25,均衡器抽头个数L=11。
[0142] 实验1:采用16QAM调制信号,步长μMMA=μFS-MMA=0.02、μBA-MMA=μGBA-MMA=0.0005、μBA-FS-MMA=μDNA-GBA-FS-MMA=0.003。仿真结果如图3所示,图3中的(a)为均方误差曲线,图3中的(b)~(g)分别为MMA、BA-MMA、DNA-GBA-MMA、FS-MMA、BA-FS-MMA和本发明DNA-GBA-FS-MMA的星座图。图3是本发明16QAM信号仿真结果图;其中,图3中的(a)中的MMA为多模盲均衡方法,BA-MMA为基于蝙蝠方法的多模盲均衡方法,DNA-GBA-MMA为基于DNA遗传蝙蝠方法的多模盲均衡方法,FS-MMA为分数间隔多模盲均衡方法,BA-FS-MMA为基于蝙蝠方法的分数间隔多模盲均衡方法,本发明DNA-GBA-FS-MMA为基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法;图3中的(a)中的Iteration表示迭代次数,MSE表示均方误差,图3中的(b)~(g)中In-Phase和Quadrature分别表示同相和正交。
[0143] 实验2:采用16PSK调制信号,步长μMMA=μFS-MMA=0.02、μBA-MMA=μGBA-MMA=0.0018、μBA-FS-MMA=μDNA-GBA-FS-MMA=0.0035。仿真结果如图4所示,图4中的(a)为均方误差曲线,图4中的(b)~(g)分别为MMA、BA-MMA、FS-MMA、DNA-GBA-MMA、BA-FS-MMA和本发明DNA-GBA-FS-MMA的星座图。图4是本发明16PSK信号仿真结果图;图4中的(a)中的MMA为多模盲均衡方法,BA-MMA为基于蝙蝠方法的多模盲均衡方法,DNA-GBA-MMA为基于DNA遗传蝙蝠方法的多模盲均衡方法,FS-MMA为分数间隔多模盲均衡方法,BA-FS-MMA为基于蝙蝠方法的分数间隔多模盲均衡方法,本发明DNA-GBA-FS-MMA为基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法,图4中的(a)中Iteration表示迭代次数,MSE表示均方误差,(b)~(g)中In-Phase和
Quadrature表示同相和正交。
[0144] 实验3:采用16APSK调制信号,步长μMMA=μFS-MMA=0.01、μBA-MMA=μGBA-MMA=0.0008、μBA-FS-MMA=μDNA-GBA-FS-MMA=0.001。仿真结果如图5所示,图5中的(a)为均方误差曲线,图5中的(b)~(g)分别为MMA、BA-MMA、DNA-GBA-MMA、FS-MMA、BA-FS-MMA和本发明DNA-GBA-FS-MMA的星座图。图5中的(a)中的MMA为多模盲均衡方法,BA-MMA为基于蝙蝠方法的多模盲均衡方法,DNA-GBA-MMA为基于DNA遗传蝙蝠方法的多模盲均衡方法,FS-MMA为分数间隔多模盲均衡方法,BA-FS-MMA为基于蝙蝠方法的分数间隔多模盲均衡方法,本发明DNA-GBA-FS-MMA为基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法。图5中(a)中的Iteration表示迭代次数,MSE表示均方误差,(b)~(g)中In-Phase和Quadrature表示同相和正交。
[0145] 由图3可知,本发明DNA-GBA-FS-MMA和BA-FS-MMA迭代200次左右收敛,收敛速度比BA-MMA和DNA-GBA-MMA快了约100多次,比FSE-MMA和MMA快了约300多次;本发明DNA-GBA-FS-MMA的稳态误差达到约-24.5dB,比BA-FS-MMA降低了2.5dB,比FS-MMA降低了3.5dB,比DNA-GBA-MMA降低了4.5dB,比BA-MMA降低了6.5dB,比MMA降低了8dB,且本发明DNA-GBA-FS-MMA星座图的星座点最清晰、最紧凑。
[0146] 由图4可知,本发明DNA-GBA-FS-MMA、BA-FS-MMA、DNA-GBA-MMA和BA-MMA均迭代100次左右收敛,收敛速度比FSE-MMA和MMA快了约500次;本发明DNA-GBA-FS-MMA的稳态误差达到约-24dB,比BA-FS-MMA降低了1dB,比DNA-GBA-MMA降低了2dB,比FS-MMA降低了3dB,比BA-MMA降低了3.5dB,比MMA降低了7.5dB,且本发明DNA-GBA-FS-MMA星座图的星座点最清晰、最紧凑。
[0147] 由图5可知,本发明DNA-GBA-FS-MMA和BA-FS-MMA迭代200次左右收敛,收敛速度比BA-MMA和DNA-GBA-MMA慢了约150次,比FSE-MMA和MMA快了约400次;本发明DNA-GBA-FS-MMA的稳态误差达到约-25dB,比BA-FS-MMA降低了1dB,比FS-MMA降低了3dB,比DNA-GBA-MMA降低了4dB,比BA-MMA降低了5dB,比MMA降低了8dB,且本发明DNA-GBA-FS-MMA星座图的星座点最清晰、最紧凑。
[0148] 可见,使用本发明方法DNA-GBA-FS-MMA可以显著提高盲均衡方法的收敛速度和减少均方稳态误差,在通信技术领域有很强的实用价值。
[0149] 本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。

附图说明

[0057] 图1是分数间隔多模盲均衡方法FS-MMA原理框图;其中,(a)是分数间隔多模盲均衡方法FS-MMA原理框图,(b)是(a)中的多模模块。
[0058] 图2是本发明提供的第m个子信道中基于DNA遗传蝙蝠方法的分数间隔多模盲均衡方法DNA-GBA-FS-MMA。
[0059] 图3是本发明16QAM信号仿真结果图;其中,(a)是均方误差曲线,(b)为MMA输出星座图,(c)是BA-MMA输出星座图,(d)为DNA-GBA-MMA输出星座图,(e)为FS-MMA输出星座图,(f)为BA-FS-MMA输出星座图,(g)为DNA-GBA-FS-MMA输出星座图。
[0060] 图4是本发明16PSK信号仿真结果图;其中,(a)是均方误差曲线,(b)为MMA输出星座图,(c)是BA-MMA输出星座图,(d)为DNA-GBA-MMA输出星座图,(e)为FS-MMA输出星座图,(f)为BA-FS-MMA输出星座图,(g)为DNA-GBA-FS-MMA输出星座图。
[0061] 图5是本发明16APSK信号仿真结果图;其中,(a)是均方误差曲线,(b)为MMA输出星座图,(c)BA-MMA输出星座图,(d)为DNA-GBA-MMA输出星座图,(e)为FS-MMA输出星座图,(f)为BA-FS-MMA输出星座图,(g)为DNA-GBA-FS-MMA输出星座图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号