[发明专利]基于改进遗传算法的三维成像声纳接收平面阵阵元稀疏优化方法有效
申请号: | 201610033074.4 | 申请日: | 2016-01-19 |
公开(公告)号: | CN105572658B | 公开(公告)日: | 2018-06-08 |
发明(设计)人: | 陈飞虎;黄勇;王朋;黄海宁;赵兴奋;陈苏广 | 申请(专利权)人: | 苏州桑泰海洋仪器研发有限责任公司 |
主分类号: | G01S7/527 | 分类号: | G01S7/527;G01S15/89;G06F17/50 |
代理公司: | 苏州市新苏专利事务所有限公司 32221 | 代理人: | 徐鸣 |
地址: | 215125 江苏省苏州*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 稀疏 三维成像 声纳 阵元 改进遗传算法 接收平面 优化 波束方向图 适应度函数 初始种群 问题转化 声纳阵 染色体 旁瓣 | ||
1.一种基于改进遗传算法的三维成像声纳接收平面阵阵元稀疏优化方法,其特征在于:所述方法包括下述步骤:
1)将三维成像声纳接收平面阵阵元稀疏优化的问题转化为染色体串,即将三维成像声纳的接收平面阵阵元位置的开启情况进行编码,若开启阵元,则对该阵元对应的位置进行置1,否则置0,此时通过round(rand(M,N)+0.5)进行生成阵元个体开启情况;
2)初始种群确定:
对阵元数为M*N维的接收平面阵进行阵元稀疏布阵,获得阵元位置个数为M*N,其中每个个体表示为1*M*N的行向量,该向量中的每个元素值非0即1,表示该接收平面阵中的阵元位置处是否有无阵元开启,此时随机生成初始父代,并确定初始父代种群的数目,假设种群数目为popsize,初始父代的种群通过MATLAB的rand()和round()函数实现,round(rand(popsize,M,N)+0.5);
3)三维成像声纳阵元稀疏的适应度函数确定:
所述三维成像声纳阵元稀疏的适应度函数表达式为
其中,BPLL(ux,vy)表示最大旁瓣峰值的期望值,B(ux,vy)为波束输出能量值,ux和vy分别为波束方向矢量值,对一个固定大小的平面阵来说,BPLL(ux,vy)为一个常数,Ω为满足以下两个条件的(ux,vy)集合:
其中,λ表示为发射信号的波长,L表示接收平面阵的阵列边长,且主瓣波束方向图的能量不在(ux,vy)集合范围内,R为阵元稀疏率,k1为加权因子,表示稀疏率与波束方向图的最大旁瓣峰值相比在优化适应度函数中的重要程度;
4)选择确定:
假设群体的目标平均适应度值为个体x的适应度值为fitness(x),若则需保留当前个体;否则,需用rand函数产生一个0到1之间的随机数,记为a,若则将当前的个体保留下来,否则剔除该个体;比较该个体x对应的适应度值与当前找到的所有的最好的个体对应的适应度值,两者中适应度值较大的对应的个体作为新的最好的个体x;
5)交叉确定:
当上述步骤3)中所述的适应度函数中的稀疏率为一随机值时,通过采用均匀交叉法,即随机生成一个1*M*N的行向量,且该向量由0和1组成,1所在的位置表示两个父代样本需要基因交换的位置,而0的位置是保持两个父代样本不变的位置;
当上述步骤3)中所述的适应度函数中的稀疏率为一固定值时,首先在父代中选择两个个体进行交叉,交叉运算如下:
其中,F1和F2表示两个交叉的个体,γ为两个交叉个体的交叉参数,通过对上述计算得到的新的个体S1中含有1的个数为N1,则在新的个体S2中随机选取[γ(N0-N1)]个数值等于γ的向上取整的数值,同样的选择[(1-γ)(N0-N1)]个数值等于(1-γ)的向上取整的数值,余下的则为γ或者(1-γ)个向下取整的整数,这样就保证了新得到的个体S1中1的个数为N0,同理,对S2进行处理,其中,N0表示稀疏平面阵中阵元的个数;接着对个体S1、S2、F1、F2分别按照适应度函数进行求解,选取其中适应度值较大的两个个体,将其作为交叉后的新个体进行输出,即为本次交叉获得的子代个体;
6)变异确定:
变异过程即为个体中出现基因突变的元素,所述的突变过程即由0变1或由1变0的过程,也即是开启还是关闭阵元;
为了避免遗传算法出现早熟现象,加快收敛速度,采用自适应小波变换对上述步骤5)中得到的个体进行变异操作得到新的个体,采用自适应小波变换的变异运算如下:
其中,为新的子代个体,为旧的子代个体,Sk为第k次迭代之后适应度值,Sfit为适应度值中最坏值,κ的表达式遵循Morlet变换表达式:
其中,较大的|κ|产生的变异较大,较小的|κ|产生较小的变异,另外,当κ为正时,此时会产生较好的适应度值,负的κ会使适应度值向坏的方向发展,这样就会使得收敛速度加快:
ψ=exp(In(kmax)*(1-(1-k/kmax)ε)) (6)
θ=5ψ(2rand-1)/2 (7)
其中,rand为0到1之间的随机数,k为当前迭代次数,kmax为预设的最大迭代次数,εmin和εmax是预设的最小值和最大值,ε的选取影响收敛的结果,ψ是一个随迭代次数不断增大的变量,刚开始,ψ的值很小使得ε的值足够大以产生较大的搜索空间,当ψ的值很大时,这样会得到很小的ε,从而使得搜索空间变小。
2.根据权利要求1所述的基于改进遗传算法的三维成像声纳接收平面阵阵元稀疏优化方法,其特征在于:所述步骤3)所述的B(ux,vy)的表达式为:
其中:
式(10)中,π为圆周率,α为方位角,β为俯仰角,α0和β0分别对应初始方位角和俯仰角。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于苏州桑泰海洋仪器研发有限责任公司,未经苏州桑泰海洋仪器研发有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201610033074.4/1.html,转载请声明来源钻瓜专利网。
- 上一篇:包裹安全检测系统和方法
- 下一篇:一种雷达坐标精度鉴定试验实时误差指示方法