发明名称 一种基于遗传算法的室内水下目标定位方法
摘要 本发明公开了一种基于遗传算法的室内水下目标定位方法,将水下目标的定位问题转化为优化问题,并应用混合遗传算法求解优化问题,得出定位结果。遗传算法中涉及的目标函数、适应度函数以及约束条件根据具体的室内水下目标定位系统确定。与现有技术相比,本发明极大的提高了目标定位精度,且其具有很高的收敛可靠性和较高的收敛速度,也可使测量误差对定位精度的影响较小。本发明亦可扩展应用于其它基于测距信息对目标进行定位的系统中。
申请公布号 CN105388460A 申请公布日期 2016.03.09
申请号 CN201510677240.X 申请日期 2015.10.19
申请人 东南大学 发明人 陈熙源;臧云歌
分类号 G01S5/30(2006.01)I 主分类号 G01S5/30(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 徐激波
主权项 一种基于遗传算法的室内水下目标定位方法,其特征在于:该方法具体步骤如下:步骤一:明确室内水下目标定位系统结构;室内水下定位系统应用短基线原理,其测量部分主要包括一个水声换能器A,三个应答器B、C、D以及一个深度传感器E,水声换能器A以及深度传感器E安装在水下目标上,三个应答器B、C、D以水域中心为中心点,以三角阵的方式分布在定位水域内,两两距离不小于三分之一水域宽度且不大于二分之一水域长度,水听器与水域边界保持一定距离,不小于四分之一水域宽度,且水听器安置深度大致为水域深度的一半,不接触水域底部;步骤二:设定已知量和自变量;在水域内以水域中心位置为原点建立三维坐标系,水域纵向为x轴,横向为y轴,深度方向为z轴,目标即换能器A与三个应答器的距离通过水声测距原理得到,深度信息由深度传感器E测得,故已知量为水域纵向长度a、水域横向宽度b、水域深度c,三个应答器B、C、D的坐标,分别为(x<sub>0</sub>,y<sub>0</sub>,z<sub>0</sub>)、(x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>)、(x<sub>2</sub>,y<sub>2</sub>,z<sub>2</sub>),以及目标与应答器B、C、D之间的距离s<sub>0</sub>,s<sub>1</sub>,s<sub>2</sub>和目标的深度h;设定的自变量为目标的xy平面坐标(x,y);步骤三:确立约束条件和目标函数;目标位于水域内,故约束条件为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mo>-</mo><mfrac><mi>a</mi><mn>2</mn></mfrac><mo>&le;</mo><mi>x</mi><mo>&le;</mo><mfrac><mi>a</mi><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mi>b</mi><mn>2</mn></mfrac><mo>&le;</mo><mi>y</mi><mo>&le;</mo><mfrac><mi>b</mi><mn>2</mn></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000824555190000011.GIF" wi="278" he="271" /></maths>根据几何定律,自变量与已知量之间存在以下定位关系:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><msqrt><mrow><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><msub><mi>s</mi><mn>0</mn></msub><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mrow><msqrt><mrow><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><msub><mi>s</mi><mn>1</mn></msub><mo>=</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msqrt><mrow><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><msub><mi>s</mi><mn>2</mn></msub><mo>=</mo><msub><mi>n</mi><mn>3</mn></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000824555190000012.GIF" wi="902" he="301" /></maths>n<sub>i</sub>为测量时引入的噪声,将定位问题转化为极值优化问题:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>f</mi><mi>i</mi><mi>n</mi><mi>d</mi><mo>:</mo></mrow></mtd><mtd><mrow><mi>X</mi><mo>=</mo><mo>&lsqb;</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>&rsqb;</mo><mo>,</mo><mi>X</mi><mo>&Element;</mo><mi>&Phi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>min</mi><mo>:</mo></mrow></mtd><mtd><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>s</mi><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>s</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>h</mi><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>s</mi><mn>2</mn><mn>2</mn></msubsup></mrow></msqrt></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000824555190000021.GIF" wi="1744" he="174" /></maths>式中Φ为定位方程组的解空间,f(x,y)为优化目标函数,当f(x,y)最小时,所对应的X即为目标定位结果;步骤四:应用混合遗传算法对f(x,y)进行优化,得到优化结果及对应的X,其具体步骤为:1)编码,确定变量上下限,生成初始群体;2)根据目标函数f(x)建立适应度函数<img file="FDA0000824555190000022.GIF" wi="447" he="135" />使用适应度函数区分群体中的好坏;3)将交叉算子作用于群体;4)将变异算子作用与群体;5)依据自适应混合算子概率p<sub>n</sub>,对群体进行经典算法局部搜索,自适应概率为p<sub>n</sub>(t)=p<sub>0</sub>e<sup>‑a(1‑t/T)</sup>,T为遗传算法中最大代数,t为当前进化的代数,常数p<sub>0</sub>∈(0,1],a取1;6)将选择算子作用与群体;7)终止条件判断;终止条件为是否超过最大进化代数;得到定位方程组解(x,y),目标定位结果即为(x,y,h)。
地址 211189 江苏省南京市江宁区东南大学路2号