发明名称 基于伪杂交混合遗传算法的机构运动链同构识别方法
摘要 本发明公开了一种基于伪杂交混合遗传算法的机构运动链同构识别方法,先根据机构运动链的结构形成其对应的机构拓扑图;设定邻接矩阵的遗传编码;取待判定的两机构拓扑图邻接矩阵确定目标函数,计算每个个体的目标函数值和适应度;再采用伪杂交算子,在随机产生的群体中随机选取两个个体,随机产生两个杂交点进行空中扩展;从扩展生成的中间种群中选择两个个体进行初步局部寻优;最后结合局部搜索算子形成伪杂交混合遗传算法,进行机构运动链同构识别的运算和仿真,能对运动链进行高效、准确的同构判别。
申请公布号 CN101655928B 申请公布日期 2012.06.20
申请号 CN200910184138.0 申请日期 2009.08.25
申请人 江苏大学 发明人 杨平;曾科翰
分类号 G06N3/12(2006.01)I;G06F17/50(2006.01)I 主分类号 G06N3/12(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 1.一种基于伪杂交混合遗传算法的机构运动链同构识别方法,其特征是依次采用如下步骤:1)根据机构运动链的结构形成其对应的机构拓扑图;2)根据机构拓扑图设定邻接矩阵的遗传编码;3)取待判定的两机构拓扑图邻接矩阵确定目标函数,计算每个个体的目标函数值和适应度;目标函数为:<img file="FSB00000609317300011.GIF" wi="598" he="158" />C=TAT<sup>T</sup>,式中,d为个体长度,T=P<sub>1</sub>P<sub>2</sub>……P<sub>n</sub>,其中P<sub>i</sub>为矩阵,A、B为两个机构拓扑图的邻接矩阵,C为邻接矩阵A经过变换矩阵T进行行变换和列变换后得到的矩阵;4)采用伪杂交算子,在随机产生的父代群体中随机选取两个个体x、y,随机产生两个杂交点k<sub>1</sub>、k<sub>2</sub>进行空中扩展;从扩展生成的中间种群中选择两个个体进行初步局部寻优;所述两个个体x、y表示为:<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><msub><mi>k</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>x</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><msub><mi>k</mi><mn>2</mn></msub></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><msub><mi>k</mi><mn>1</mn></msub></msub><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>y</mi><msub><mi>k</mi><mn>2</mn></msub></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>y</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>(0≤k<sub>1</sub>≤k<sub>2</sub>≤d),d为个体长度;令tx=x,ty=y,对x进行如下操作:步骤1:令i=k<sub>1</sub>+1;步骤2:令j=0;步骤3:如果tx<sub>j</sub>≠ty<sub>i</sub>,则j=j+1,重复步骤3;步骤4;如果(i≠j)∩(k<sub>1</sub>+1≤j≤k<sub>2</sub>-1),则转步骤5;否则转步骤6;步骤5:交换x<sub>i</sub>和x<sub>j</sub>;步骤6:i=i+1,如果i=k<sub>2</sub>,结束,否则转步骤2;对y进行如下操作:步骤1:令i=k<sub>1</sub>+1;步骤2:令j=0;步骤3:如果ty<sub>j</sub>≠tx<sub>i</sub>,则j=j+1,重复步骤3;步骤4:如果(i≠j)∩(k<sub>1</sub>+1≤j≤k<sub>2</sub>-1),则转步骤5;否则转步骤6;步骤5:交换y<sub>i</sub>和y<sub>j</sub>;步骤6:i=i+1;如果i=k<sub>2</sub>,结束;否则转步骤2;所述初步局部寻优的步骤为:步骤1:令i=1;步骤2:如果y<sub>i</sub>≠i,则交换x<sub>i</sub>和<img file="FSB00000609317300021.GIF" wi="76" he="47" />得到<maths num="0003"><![CDATA[<math><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mo>=</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><msub><mi>y</mi><mi>i</mi></msub></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>否则转步骤4;步骤3:若F(x′)>F(x),则x=x′;步骤4:i=i+1,如果i>d,则结束;否则转步骤2;5)结合局部搜索算子形成伪杂交混合遗传算法,进行机构运动链同构识别的运算和仿真;局部搜索算子步骤如下:步骤1:令i=1;步骤2:令j=1;步骤3:如果j≠i,则交换x<sub>i</sub>和x<sub>j</sub>的位置,得到x′;否则转步骤5;步骤4:比较x和x′的适应度,如果F(x′)>F(x),则x=x′;步骤5:j=j+1;如果j≤d,则转步骤3;步骤6:i=i+1;如果i≤d,转步骤2;否则结束。
地址 212013 江苏省镇江市学府路301号