发明名称 一种迭代学控制参考轨迹的优化匹配组合方法
摘要 迭代学控制参考轨迹的优化匹配组合方法,基于Kabsch算法设计NURBS曲线优化匹配方法,首先利用NURBS来描述任意曲线,曲线的分裂和组合就具有统一的数学形式;其次,利用Kabsch算法描述两段NURBS曲线的相似性;在满足相似度的情况下,求取使期望参考曲线分割次数最少时的匹配方法;最后,将所有匹配的基元依次组合起来,此组合曲线与原期望曲线相似。
申请公布号 CN105093926A 申请公布日期 2015.11.25
申请号 CN201510446408.6 申请日期 2015.07.27
申请人 浙江工业大学 发明人 徐建明;臧永灿;朱自立;孙明轩;俞立
分类号 G05B13/02(2006.01)I 主分类号 G05B13/02(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 一种迭代学习控制参考轨迹的优化匹配组合方法,包括对于已知的参考轨迹或期望轨迹l<sub>d</sub>(v)|<sub>v∈[0,1]</sub>(v为NURBS轨迹的节点矢量参数,为了便于优化匹配算法研究,整条轨迹的定义域为规范参数域v∈[0,1])和轨迹库CurveBank,轨迹库CurveBank其中含有n条NURBS轨迹l<sub>i</sub>(v)|<sub>v∈[0,1]</sub>(1≤i≤n,i∈N),分割参数<img file="FDA0000767918350000011.GIF" wi="47" he="73" />和<img file="FDA0000767918350000012.GIF" wi="353" he="76" />分割其中的第i条轨迹形成临时基元<img file="FDA0000767918350000013.GIF" wi="237" he="76" />迭代学习控制简称ILC;步骤如下:Step1:寻找与整条参考轨迹l<sub>d</sub>(v)|<sub>v∈[0,1]</sub>最相似的一条轨迹,令l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>=l<sub>d</sub>(v)|<sub>v∈[0,1]</sub>(j=1,l<sub>d1</sub>(v)|<sub>v∈[0,1]</sub>表示第1段参考轨迹),分裂参数<img file="FDA0000767918350000017.GIF" wi="35" he="67" />和搜索步长b初始化;Step2:假设L<sub>d</sub>和L<sub>l</sub>分别是在轨迹l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>(其中l<sub>dj</sub>(v)|<sub>v=0</sub>和l<sub>dj</sub>(v)|<sub>v=1</sub>分别为轨迹l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>的起止点)和l<sub>l</sub>(v)|<sub>v∈[0,1]</sub>(临时基元)上取出的在世界坐标系{B}下的坐标序列,k+1为此序列的坐标点个数,如式(1)、(2)所示;O<sub>d</sub>和O<sub>l</sub>分别是轨迹l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>和l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>的质心,w<sub>i</sub>表示L<sub>d</sub>和L<sub>l</sub>的第i点对应的权值<img file="FDA0000767918350000014.GIF" wi="448" he="105" />如式(3)、(4)所示;以轨迹l<sub>d</sub>(v)|<sub>v∈[0,1]</sub>和l<sub>l</sub>(v)|<sub>v∈[0,1]</sub>各自质心为原点的笛卡尔直角坐标系,暂称作质心坐标系,<sup>d</sup>L<sub>d</sub>和<sup>l</sup>L<sub>l</sub>分别是L<sub>d</sub>和L<sub>l</sub>相对各自质心坐标系(与世界坐标系的姿态一致)下的坐标序列,R为使L<sub>d</sub>和L<sub>l</sub>的均方根偏差最小时的最优旋转矩阵,如式(5)所示;<maths num="0001" id="cmaths0001"><img file="FDA0000767918350000015.GIF" wi="1234" he="148" /></maths><maths num="0002" id="cmaths0002"><img file="FDA0000767918350000016.GIF" wi="1202" he="150" /></maths><maths num="0003" id="cmaths0003"><img file="FDA0000767918350000021.GIF" wi="1014" he="139" /></maths><maths num="0004" id="cmaths0004"><img file="FDA0000767918350000022.GIF" wi="993" he="139" /></maths>D<sub>f</sub>=R<sup>l</sup>L<sub>l</sub>‑<sup>d</sup>L<sub>d</sub>   (5)<img file="FDA0000767918350000023.GIF" wi="1038" he="265" /><maths num="0005" id="cmaths0005"><img file="FDA0000767918350000024.GIF" wi="968" he="109" /></maths>其中,D<sub>f</sub>为偏差矩阵,D'<sub>f</sub>为D<sub>f</sub>的转置,tr(D<sub>t</sub>)表示矩阵D<sub>t</sub>的迹,通过公式(1)~(7)得到两条轨迹的最小均方根偏差lrm,给定一个相似度ε&gt;0,若lrm&lt;ε,则两条轨迹相似,跳入Step3;反之不相似,跳入Step4;Step3:若<img file="FDA0000767918350000025.GIF" wi="124" he="72" />跳入Step6;若<img file="FDA0000767918350000026.GIF" wi="125" he="73" />则判断搜索步长是否达到最小;若没达到,则令搜索步长b=b/2,分裂参数<img file="FDA0000767918350000027.GIF" wi="189" he="70" />跳入Step2;若达到则跳入Step5;Step4:令搜索步长b=b/2,分裂参数<img file="FDA0000767918350000028.GIF" wi="187" he="74" />跳入Step2;Step5:此时可得被分割的第j段期望轨迹<img file="FDA0000767918350000029.GIF" wi="198" he="86" />和与其对应的匹配基元l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>,这样,第j次匹配结束,期望轨迹l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>被分割成两段期望轨迹l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>与l<sub>d(j+1)</sub>(v)|<sub>v∈[0,1]</sub>;将l<sub>d(j+1)</sub>(v)|<sub>v∈[0,1]</sub>作为整条期望轨迹,令j=j+1,分裂参数<img file="FDA00007679183500000210.GIF" wi="40" he="71" />和搜索步长b初始化,跳入Step2;Step6:此时可得所有期望轨迹段l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>和与其对应的匹配基元l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>(j=1,……,m),优化匹配结束;Step7:对匹配基元l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>(j=1,……,m)(<sup>j</sup>l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>为匹配基元l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>在其质心坐标系{j}下的表示)进行旋转和平移变换:l"<sub>j</sub>(v)|<sub>v∈[0,1]</sub>=R<sub>j</sub>·<sup>j</sup>l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>+<sup>dj</sup>P<sub>BORG</sub>其中:R<sub>j</sub>和<sup>dj</sup>P<sub>BORG</sub>分别是基元<sup>j</sup>l<sub>j</sub>(v)|<sub>v∈[0,1]</sub>与对应的期望曲线段l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>之间最优的旋转矩阵和平移向量;Step8:组合匹配基元l<sub>dj</sub>(v)|<sub>v∈[0,1]</sub>(j=1,……,m),获得ILC参考轨迹<sub>ld(v)|v∈[0,1]</sub>的相似参考轨迹:<img file="FDA0000767918350000031.GIF" wi="553" he="102" />其中<img file="FDA0000767918350000032.GIF" wi="304" he="87" />表示将m条曲线l"<sub>j</sub>(v)|<sub>v∈[0,1]</sub>依次按各自质心在世界坐标系下的组合,<img file="FDA0000767918350000033.GIF" wi="42" he="79" />表示此组合轨迹与参考轨迹l<sub>d</sub>(v)|<sub>v∈[0,1]</sub>叠合。
地址 310014 浙江省杭州市下城区潮王路18号