• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 13
    • 下载费用:30 金币  

    重庆时时彩的网站地址: 一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法.pdf

    关 键 词:
    一种 基于 量子 遗传 算法 雷达 海杂波 反演 蒸发 波导 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610877172.6

    申请日:

    2016.10.08

    公开号:

    CN106354979A

    公开日:

    2017.01.25

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20161008|||公开
    IPC分类号: G06F17/50; G06N3/12 主分类号: G06F17/50
    申请人: 西安电子科技大学
    发明人: 胡荣旭; 周伟豪; 张金鹏; 吴振森
    地址: 710071 陕西省西安市太白南路2号
    优先权:
    专利代理机构: 北京挺立专利事务所(普通合伙) 11265 代理人: 叶树明
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610877172.6

    授权公告号:

    |||

    法律状态公告日:

    2017.03.01|||2017.01.25

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法,该方法将QGA算法应用于RFC技术的实施,并将这种结合称为RFC?based?QGA?RFC?based技术。本发明RFC?based?QGA在蒸发波导折射率剖面反演中具有可行性和优良性能。RFC?based?QGA具有精确度高,耗时少,稳定好的特点,适合推广应用。

    权利要求书

    1.一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法,其特征在于,包括以下
    步骤:
    开始t←0
    (1)用Q-bit初始化蒸发波导高度种群其中表示种群第t
    代中的第j个表示蒸发波导高度的个体,N是初始种群大??;
    (2)对初始种群Q(t)中每个个体实施一次测量,得到相应的蒸发波导高度二进制确
    定解获得种群的二进制串表达
    (3)根据适应度目标函数Φ,计算B(t)中的每个确定解的适应度值获得适应
    度值集
    (4)保存F(t)中最佳适应度值并保存相应最佳个体作为当前所有个体下
    一代演化的最佳目标;
    (5)while(为满足终止条件)do
    开始
    (a)t←t+1
    (b)同步骤(2),一次测量种群Q(t-1)中的每个个体,产生二进制确定解群X(t);
    (c)同步骤(3),对(b)产生的X(t)进行适应值评估,得到其适应值集F(t);
    (d)通过量子旋转门更新Q(t-1)获得新一代种群Q(t);
    (e)同步骤(4),存储F(t)中最佳适应度值并保存相应最佳个体作为当前
    所有个体下一代演化的最佳目标;
    结束
    结束。

    说明书

    一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法

    技术领域

    本发明属于优化算法QGA应用技术领域,涉及一种基于量子遗传算法的雷达海杂
    波反演蒸发波导的方法。

    背景技术

    海面上低空对流层折射率常随着气象条件发生着时空的变化。当对流层的垂直温
    度和湿度发生异常变化时,就会产生大气波导。大气波导可以使处于其中的电磁波产生超
    视距和雷达空洞现象,这些现象极大的影响着海上雷达和通信系统的性能。因此实时获得
    大气波导的垂直折射率剖面信息对雷达性能评估和预测有着重要的作用。RFC
    (Refractivity from clutter)是利用雷达海杂波信息反演大气波导垂直折射率剖面信息
    的一种新兴技术,它实质上是将实际测量海杂波与模拟海杂波进行对比拟合的过程,当两
    组数据吻合最好时的模拟海杂波所对应的大气折射率垂直剖面即为反演的最佳折射率剖
    面。RFC技术的实施要用到优化算法,合适的优化算法可以让RFC技术的性能得到极大的提
    升,为此业界很多工作者都在致力于寻找更好的优化算法。量子遗传算法(QGA,Quantum
    Genetic Algorithm)是一个正处于发展中的优化算法,正式提出于2000年,它结合量子计
    算和遗传算法,主要用量子比特编码基因和个体,并用量子旋转门更新个体,具有种群小、
    多样性丰富和收敛快的优点,已经在很多领域被成功尝试。目前,现有技术中还没有将QGA
    算法应用于RFC技术的实施技术。

    发明内容

    为了弥补现有技术中的空白,本发明提出了一种基于量子遗传算法的雷达海杂波
    反演蒸发波导的方法,该方法将QGA算法应用于RFC技术的实施,并将这种结合称为RFC-
    based-QGA技术。其技术方案如下:

    一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法,包括以下步骤:

    开始t←0

    (1)用Q-bit初始化蒸发波导高度种群其中表示种
    群第t代中的第j个表示蒸发波导高度的个体,N是初始种群大小。

    (2)对初始种群Q(t)中每个个体实施一次测量,得到相应的蒸发波导高度二进
    制确定解获得种群的二进制串表达

    (3)根据适应度目标函数Φ,计算B(t)中的每个确定解的适应度值获得
    适应度值集

    (4)保存F(t)中最佳适应度值并保存相应最佳个体作为当前所有个体
    下一代演化的最佳目标。

    (5)while(为满足终止条件)do

    开始

    (a)t←t+1

    (b)同步骤(2),一次测量种群Q(t-1)中的每个个体,产生二进制确定解群X(t)。

    (c)同步骤(3),对(b)产生的X(t)进行适应值评估,得到其适应值集F(t)。

    (d)通过量子旋转门更新Q(t-1)获得新一代种群Q(t)。

    (e)同步骤(4),存储F(t)中最佳适应度值并保存相应最佳个体作为
    当前所有个体下一代演化的最佳目标。

    结束

    结束。

    本发明的有益效果为:

    本发明RFC-based-QGA在蒸发波导折射率剖面反演中具有可行性和优良性能。
    RFC-based-QGA具有精确度高,耗时少,稳定好的特点,适合推广应用。

    附图说明

    图1是蒸发波导修正折射率垂直剖面;

    图2是蒸发波导环境下的标准化海杂波功率群(2.84GHz,30.78m)

    图3是蒸发波导环境下的标准化海杂波功率群(10GHz,10m)

    图4是实际观测的标准化海杂波功率(10GHz,13m);

    图5是基于海杂波群图2的蒸发波导高度反演结果频率直方图和散点图,其中:

    图5(a)是基于图2中海杂波曲线(a)的蒸发波导高度PSO反演结果频率直方图。

    图5(b)是基于图2中海杂波曲线(b)的蒸发波导高度PSO反演结果频率直方图。

    图5(c)是基于图2中海杂波曲线(c)的蒸发波导高度PSO反演结果频率直方图。

    图5(d)是基于图2中海杂波曲线(d)的蒸发波导高度PSO反演结果频率直方图。

    图5(e)是基于图2中海杂波曲线(e)的蒸发波导高度PSO反演结果频率直方图。

    图5(f)是基于图2中海杂波曲线(a)的蒸发波导高度PSO反演结果散点图。

    图5(g)是基于图2中海杂波曲线(b)的蒸发波导高度PSO反演结果散点图。

    图5(h)是基于图2中海杂波曲线(c)的蒸发波导高度PSO反演结果散点图。

    图5(i)是基于图2中海杂波曲线(d)的蒸发波导高度PSO反演结果散点图。

    图5(j)是基于图2中海杂波曲线(e)的蒸发波导高度PSO反演结果散点图。

    图5(k)是基于图2中海杂波曲线(a)的蒸发波导高度QGA反演结果频率直方图。

    图5(l)是基于图2中海杂波曲线(b)的蒸发波导高度QGA反演结果频率直方图。

    图5(m)是基于图2中海杂波曲线(c)的蒸发波导高度QGA反演结果频率直方图。

    图5(n)是基于图2中海杂波曲线(d)的蒸发波导高度QGA反演结果频率直方图。

    图5(o)是基于图2中海杂波曲线(e)的蒸发波导高度QGA反演结果频率直方图。

    图5(p)是基于图2中海杂波曲线(a)的蒸发波导高度QGA反演结果散点图。

    图5(q)是基于图2中海杂波曲线(b)的蒸发波导高度QGA反演结果散点图。

    图5(r)是基于图2中海杂波曲线(c)的蒸发波导高度QGA反演结果散点图。

    图5(s)是基于图2中海杂波曲线(d)的蒸发波导高度QGA反演结果散点图。

    图5(t)是基于图2中海杂波曲线(e)的蒸发波导高度QGA反演结果散点图。

    图6是基于海杂波群图3的蒸发波导高度反演结果频率直方图和散点图,其中:

    图6(a)是基于图3中海杂波曲线(a)的蒸发波导高度PSO反演结果频率直方图。

    图6(b)是基于图3中海杂波曲线(b)的蒸发波导高度PSO反演结果频率直方图。

    图6(c)是基于图3中海杂波曲线(c)的蒸发波导高度PSO反演结果频率直方图。

    图6(d)是基于图3中海杂波曲线(d)的蒸发波导高度PSO反演结果频率直方图。

    图6(e)是基于图3中海杂波曲线(e)的蒸发波导高度PSO反演结果频率直方图。

    图6(f)是基于图3中海杂波曲线(a)的蒸发波导高度PSO反演结果散点图。

    图6(g)是基于图3中海杂波曲线(b)的蒸发波导高度PSO反演结果散点图。

    图6(h)是基于图3中海杂波曲线(c)的蒸发波导高度PSO反演结果散点图。

    图6(i)是基于图3中海杂波曲线(d)的蒸发波导高度PSO反演结果散点图。

    图6(j)是基于图3中海杂波曲线(e)的蒸发波导高度PSO反演结果散点图。

    图6(k)是基于图3中海杂波曲线(a)的蒸发波导高度QGA反演结果频率直方图。

    图6(l)是基于图3中海杂波曲线(b)的蒸发波导高度QGA反演结果频率直方图。

    图6(m)是基于图3中海杂波曲线(c)的蒸发波导高度QGA反演结果频率直方图。

    图6(n)是基于图3中海杂波曲线(d)的蒸发波导高度QGA反演结果频率直方图。

    图6(o)是基于图3中海杂波曲线(e)的蒸发波导高度QGA反演结果频率直方图。

    图6(p)是基于图3中海杂波曲线(a)的蒸发波导高度QGA反演结果散点图。

    图6(q)是基于图3中海杂波曲线(b)的蒸发波导高度QGA反演结果散点图。

    图6(r)是基于图3中海杂波曲线(c)的蒸发波导高度QGA反演结果散点图。

    图6(s)是基于图3中海杂波曲线(d)的蒸发波导高度QGA反演结果散点图。

    图6(t)是基于图3中海杂波曲线(e)的蒸发波导高度QGA反演结果散点图。

    图7是基于实测杂波图4的蒸发波导高度反演结果频率直方图和散点图,其中:

    图7(a)是基于图4所示实测海杂波曲线的蒸发波导高度PSO反演结果频率直方图。

    图7(b)是基于图4所示实测海杂波曲线的蒸发波导高度PSO反演结果散点图。

    图7(c)是基于图4所示实测海杂波曲线的蒸发波导高度QGA反演结果频率直方图。

    图7(d)是基于图4所示实测海杂波曲线的蒸发波导高度QGA反演结果散点图。

    具体实施方式

    下面结合附图和具体实施方式对本发明的技术方案作进一步详细地说明。

    理论和模型

    蒸发波导模型

    在RFC技术中,需要用到表达海洋低空对流层波导垂直剖面的参数模型,不同类型
    的对流层波导用不同的参数向量表示。常用的蒸发波导参数模型如下:

    <mrow> <mi>M</mi> <mrow> <mo>(</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>M</mi> <mn>0</mn> </msub> <mo>+</mo> <mn>0.125</mn> <mrow> <mo>(</mo> <mi>z</mi> <mo>-</mo> <mi>h</mi> <mi> </mi> <mi>l</mi> <mi>n</mi> <mfrac> <mrow> <mi>z</mi> <mo>+</mo> <msub> <mi>z</mi> <mn>0</mn> </msub> </mrow> <msub> <mi>z</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>

    式中,z为海面以上垂直高度,单位是米;z0为空气动力学粗糙因子,通常取1.5×
    10-4m;M0是海面高度处的大气修正折射率;h为蒸发波导高度,单位是米。在蒸发波导RFC技
    术中,h就是要反演的参数,只要反演得到了h,根据式(1)就得到了蒸发波导折射率垂直剖
    面M(z)。蒸发波导折射率垂直剖面如图1所示。

    海杂波模型

    RFC技术利用接收到的雷达海杂波反演海洋低空对流层波导环境,是个已知测量
    海杂波与可能波导环境下雷达接收海杂波的对比寻优过程。寻优过程的实施要基于可能环
    境中雷达海杂波功率的计算,在不考虑接雷达收机噪声时,雷达接收到来自距离x(km)处的
    海杂波功率Pc(x,Menv)可表示为:

    Pc(x,Menv)=-2L(x,Menv)+10log10(x)+σ0(x,θ(x))+C (2)

    其中L(x,Menv)为海面有效散射高度处单程传播损耗,可以由抛物方程计算出;σ0
    (x,θ(x))是海面后向散射系数,θ(x)是掠射角,蒸发波导情况下,在远距离处,θ(x)可看做
    常数。C是与雷达功率、增益等雷达参数有关的常数项。所以,只要计算出L(x,Menv),σ0(x,θ
    (x))和C,就可以模拟相应环境中的海杂波,如图2和图3所示的海杂波功率曲线。

    RFC框架

    RFC技术的寻优匹配的过程,从数学上来看,就是最小化目标函数
    的过程,即

    <mrow> <mover> <mi>h</mi> <mo>^</mo> </mover> <mo>=</mo> <munder> <mrow> <mi>arg</mi> <mi>min</mi> </mrow> <mi>h</mi> </munder> <mi>&Phi;</mi> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mi>c</mi> <mrow> <mi>o</mi> <mi>b</mi> <mi>s</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>P</mi> <mi>c</mi> <mrow> <mi>c</mi> <mi>a</mi> <mi>n</mi> </mrow> </msubsup> <mo>(</mo> <mi>h</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>

    其中,为反演出的蒸发波导修正折射率剖面式(1)中的蒸发波导高度;
    分别是实际观测海杂波和可能波导环境中雷达接收到的海杂波,当目标函数
    最小时,说明和是最佳匹配,此时所在的环境最有可能是
    所在的环境,此时所在环境的蒸发波导高度h就是要反演的结果

    反演方法

    粒子群优化算法(PSO)

    PSO(Particle swarm optimization)提出于1995年,是一种很受欢迎的随机搜索
    优化算法,它简单且收敛快,在很多领域都有成功的应用。在RFC技术中也常被用到,也是目
    前国内RFC技术所使用的优化算法中性能最好的一种算法。PSO算法及其在RFC技术中的应
    用(说明中有时会把这种结合称为RFC-based-PSO)可参考现有技术,此不赘述。

    量子遗传算法(QGA)

    下面简单介绍QGA算法及其在RFC技术中应用设计。

    量子比特(Q-bit)及决策变量的染色体表达

    一个Q-bit可以表示为,


    其中α和β是复数,分别表示基本态|0>和|1>的概率幅,|α|2和|β|2表示量子位
    处于状态|0>和状态|1>的概率,且|α|2+|β|2=1。

    一个染色体可由多个Q-bit编码为:

    (j=1,2,…,N;i=1,2,…,k;k,N是整数) (5)

    其中,表示第t代种群中的第j个染色体个体,N是种群大??;k表示构成染色体的
    Q-bit个数;为染色体中第i个量子比特的概率幅对;这里(5)表达的有2k个叠加
    态。

    更新策略

    染色体中的Q-bit可由量子旋转门更新上代Q-bit得到,其数学表达如下:

    <mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msubsup> <mi>&alpha;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&beta;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>U</mi> <mrow> <mo>(</mo> <msubsup> <mi>&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msubsup> <mi>&alpha;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&beta;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>cos&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>sin&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>sin&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>cos&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msubsup> <mi>&alpha;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>&beta;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

    其中,为t-1代染色体中第i个量子比特的概率幅对;表示量子旋
    转门;为旋转角度,的值根据如下调整策略确定:

    <mrow> <msubsup> <mi>&theta;</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>=</mo> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <mo>(</mo> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <msubsup> <mi>B</mi> <mi>j</mi> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mi>f</mi> <mrow> <mo>(</mo> <msubsup> <mi>B</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>&times;</mo> <mo>(</mo> <mrow> <msubsup> <mi>b</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>-</mo> <msubsup> <mi>b</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> <mo>.</mo> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> </mrow> <mo>)</mo> <mo>&times;</mo> <msubsup> <mi>&alpha;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&times;</mo> <msubsup> <mi>&beta;</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>&times;</mo> <mi>&Delta;</mi> <mi>&theta;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

    其中,sign(x)为符号函数,当x≥0时,sign(x)=1;当x<0时,sign(x)=-1。在量
    子门更新过程中,Δθ的大小和符号起关键作用,Δθ的幅度影响收敛速度,如果太大,会导
    致早熟,一般推荐使用0.005π~0.1π。是当前代染色体一次测量的二进制串,可被表示
    为:

    <mrow> <msubsup> <mi>B</mi> <mi>j</mi> <mi>t</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>b</mi> <mrow> <mi>j</mi> <mn>1</mn> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>j</mi> <mn>2</mn> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>j</mi> <mi>k</mi> </mrow> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>

    是二进制值“0”或“1”。是当前代染色体的适应度值。是种群当前最佳
    二进制串:

    <mrow> <msubsup> <mi>B</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> </mrow> <mi>t</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>b</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> <mn>.1</mn> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> <mn>.2</mn> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> <mo>.</mo> <mi>i</mi> </mrow> <mi>t</mi> </msubsup> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>b</mi> <mrow> <mi>b</mi> <mi>e</mi> <mi>s</mi> <mi>t</mi> <mo>.</mo> <mi>k</mi> </mrow> <mi>t</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>

    和相似,是二进制“0”或“1”。是种群当前最佳的适应度值。

    QGA在RFC中的实施

    开始t←0

    (1)用Q-bit初始化蒸发波导高度种群其中表示种
    群第t代中的第j个表示蒸发波导高度的个体,N是初始种群大小。

    (2)对初始种群Q(t)中每个个体实施一次测量,得到相应的蒸发波导高度二进
    制确定解获得种群的二进制串表达

    (3)根据适应度目标函数Φ,计算B(t)中的每个确定解的适应度值获得
    适应度值集

    (4)保存F(t)中最佳适应度值并保存相应最佳个体作为当前所有个
    体下一代演化的最佳目标。

    (5)while(为满足终止条件)do

    开始

    (a)t←t+1

    (b)同步骤(2),一次测量种群Q(t-1)中的每个个体,产生二进制确定解群X(t)。

    (c)同步骤(3),对(b)产生的X(t)进行适应值评估,得到其适应值集F(t)。

    (d)通过量子旋转门更新Q(t-1)获得新一代种群Q(t)。

    (e)同步骤(4),存储F(t)中最佳适应度值并保存相应最佳个体作为
    当前所有个体下一代演化的最佳目标。

    结束

    结束

    RFC-based-QGA应用举例

    RFC-based-QGA应用例子安排

    因为RFC的性能受雷达频率和天线高度的影响,所以文章用三组不同的应用例子
    来表明蒸发波导RFC-based-QGA技术的优良性和稳定性。第一组应用中所用到的观测海杂
    波群图2是频率为2.84GHz,天线高度为30.78米的雷达在不同蒸发波导高度环境中模拟出
    来的。第二组应用中用到的观测海杂波图3是频率为10GHz,天线高度为10米的雷达在不同
    蒸发波导高度环境中模拟出来的;第三组应用中所用到的实测海杂波曲线图4为2007年夏
    季在我国东南沿海某岛屿实际实验中测量到的,实验所用的雷达频率为10GHz,天线高度为
    13米。由于前两组应用中的观测海杂波都是模拟出来的,所以前两个应用例子实际上是模
    拟应用,后一个才为实际应用。另外,两组模拟应用中模拟观测海杂波群图2和图3所用天线
    的极化方式都为垂直极化,波束宽度都为0.39度,天线仰角都为0度,海杂波曲线(a)、(b)、
    (c)、(d)和(e)所对应的真实蒸发波导高度分别为6、12、18、24和30米。第三组观测海杂波所
    用天线的极化方式是水平极化,波束宽度为0.7度,天线仰角为0度。

    为了充分说明蒸发波导环境下RFC-based-QGA的系统性能,PSO算法被应用到上述
    三组例子中,只是为了和QGA算法做对比。算法中参数选取说明:种群大小为20,参数搜索范
    围为0-40米。QGA中采用19个Q-bit编码一个基因,量子旋转角步进Δθ=0.03π;PSO中粒子
    最大和最小搜索速度分别为0和40。

    蒸发波导RFC-based-QGA反演应用举例

    应用例1

    此例对图2所示的每条海杂波曲线分别实施了200次RFC-based-QGA和200次RFC-
    based-PSO反演,反演结果如图5所示。图5第一行,(a)、(b)、(c)、(d)和(e),和第二行,(f)、
    (g)、(h)、(i)和(j),分别是PSO算法反演结果的频率分布图和散点图;第三行,(k)、(l)、
    (m)、(n)和(o),和第四行,(p)、(q)、(r)、(s)和(t),分别是QGA算法反演结果的频率分布图
    和散点图,其中第一行和第三行频率分布图中的垂直虚线表示蒸发波导的真实高度。图5从
    左向右每一列所反演的真实蒸发波导高度分别是6,12,18,24和30米;从图5可以看出,QGA
    算法反演结果的精确度明显高于PSO算法的结果。

    为了进一步量化比较,表格1给出了两种算法200次反演结果的均值、标准差和分
    布区间以及计算机CPU耗时。从表格1可以看出,两种方法的均值都能真实反演真实波导高
    度,但是QGA的标准差和反演结果的分布区间比PSO更小更窄,这点和图5是相符的。同时从
    计算量看,QGA的平均每代CPU耗时比PSO小,虽然少的不多,但是对于反演代数比较多的情
    况下,还是比较节时的。譬如,对高度为12米的蒸发波导的反演,应用QGA反演结果均值为
    12.00,标准差为0.015,反演结果分布区间为[11.945 12.047],所耗平均CPU时间为17.602
    秒/代;而应用PSO反演结果均值为11.999,标准差为0.028,反演结果分布区间为[11.909
    12.082],所耗平均CPU时间为19.537秒/代。所以,RFC-based-QGA比RFC-based-PSO有更高
    的精度,更小的偏差和更好的稳定性以及更加节时。由此可见,对于实际反演次数往往比较
    少的情况下,RFC-based-QGA每次反演结果的可信度会更高一些。

    表1不同蒸发波导高度反演结果统计比较(2.84GHz,30.78m)


    应用例2

    类似于应用例1,图6和表格2是对图3每条海杂波曲线200个反演结果的分布图与
    统计结果。图6第一行和第二行是粒子群(PSO)方法反演结果,第三行和第四行是量子遗传
    算法(QGA)方法反演的结果。根据图6和表格2,明显RFC-based-QGA表现出比RFC-based-PSO
    更优良的性能,有更高的精度,更小的偏差、更好的稳定性和更节时以及更高的可信度,和
    应用例1的结论相同。

    表2不同蒸发波导高度反演结果统计比较(10GHz,10m)



    应用例3

    为了进一步说明RFC-based-QGA技术的优良性能,此例对东南沿海一次蒸发波导
    现象中实际测量的海杂波数据,如图4所示,分别实施200次RFC-based-QGA和200次RFC-
    based-PSO反演。经过3δ检验,PSO算法200个反演结果中190个为有效结果;QGA算法200个反
    演结果中199个为有效结果,这些有效反演结果的分布图和统计特性如图7和表3所示,图7
    中第一列为PSO反演结果,第二列为QGA反演结果。从图7可看出应用QGA算法得出的反演结
    果的精确度明显高于应用PSO算法的结果。再仔细观察表格3的统计特性,显然可得到:RFC-
    based-QGA比RFC-based-PSO有更高的精度,更小的偏差和更好的稳定性、更节时和更可信。

    这里要说明的是,在应用1和应用2中,每条海杂波曲线的每种算法的200个反演结
    果都经过3δ检验,只是检验结果表明所有结果都为有效值,所以没有进行反演结果有效性
    特别讨论。

    表3实测海杂波反演蒸发波导高度结果统计比较

    优化算法
    均值
    标准差
    反演结果分布区间
    200次反演结果的有效值个数
    平均CPU耗时
    PSO
    27.336
    0.017
    [27.314 27.386]
    190
    18.378秒/代
    QGA
    27.334
    0.010
    [27.314 27.363]
    199
    17.843秒/代

    以上所述,仅为本发明较佳的具体实施方式,本发明的?;し段Р幌抻诖?,任何熟
    悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简
    单变化或等效替换均落入本发明的?;し段?。

    关于本文
    本文标题:一种基于量子遗传算法的雷达海杂波反演蒸发波导的方法.pdf
    链接地址://www.4mum.com.cn/p-6027116.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • pk10走势图判断方法 江苏时时组三的几率 福建时时11选5下载 11选五财富计划软件 重庆时时官网投注 彩票双面盘1.995即可享受1.5返水 微信猜大小单双群软件 重庆时时是不是真的 光大彩票网手机投注站 天津时时几点开奖 歪理邪说平特一肖 北京pk拾全天赛车计划 好运来计划软件下载 云南时时中奖规则 pk106码怎么跟稳定 七乐彩专家免费预测