一种基于改进的紧凑式教学优化算法的机械参数软测量方法，涉及制造领域的机械安装参数软测量。为了解决系统零件安装完成后用传统的传感器和测量方法常常难于测量其安装机械参数的问题。所述方法包括：建立待测系统模型，确定需测量的机械参数；测量系统模型中可以直接测得的参数，代入系统模型中，建立模型方程；建立优化算法目标函数；测量待测系统的n组输入和对应的n组输出，代入目标函数；采用改进的紧凑式教学优化算法迭代过程；将获得的结果对应到需测量的机械参数。本发明专利技术用于机械系统中测量机械参数。
【技术实现步骤摘要】
本专利技术涉及制造领域的机械安装参数软测量。
技术介绍
信息时代的到来，零部件越来越精细，加工的特征尺寸越来越精密，使得超精密测 量技术、超精密定位技术和更加先进的制造技术显得越来越重要。在机械系统的制造过程 中，机械系统的动作精度不仅取决于零件的制造精度，同时取决于零件在系统中的安装精 度。但在目前，受限于安装过程中复杂多变的环境因素，精密零件的安装精度通常远远小于 零件的制造精度，无法满足高精度系统的要求，因此需要对安装完成后的机械系统参数进 行测量，但在零件安装完成后用传统的传感器和测量方法常常难于测量其安装机械参数， 因而需要有新的测量方法来完成机械系统的安装参数及误差测量。
技术实现思路
本专利技术的目的是为了解决系统零件安装完成后用传统的传感器和测量方法常常 难于测量其安装机械参数的问题，本专利技术提供一种基于改进的紧凑式教学优化算法的机械 参数软测量方法。 本专利技术的，所述方法 包括如下步骤： 步骤一：建立待测系统模型，确定需测量的机械参数； 步骤二：测量系统模型中可以直接测得的参数，代入系统模型中，建立模型方程； 步骤三：建立优化算法目标函数； 步骤四：测量待测系统的η组输入和对应的η组输出，代入目标函数； 步骤五：初始化计数器t = 0,均值初始值yt = 0,方差初始值σ t = λ ;i =0,…η ;构成PV矩阵的初始值，PV的每一行包含高斯分布的一组均 值和方差；其中t为算法迭代次数计数器的计数，由PV矩阵生成随机向量Trt的初始值； 步骤六：由PV矩阵生成随机向量Stt，Stt中的每一个元素对应PV矩阵中一组均 值和方差决定的高斯函数的随机值； 步骤七：计算均值DMeant= randiX (Trt_round(l+rand2(0,1)) X μ t)，round 函数 接受一个参数返回与参数最近的整数，randJP rand 2均为随机函数； 步骤八：更新； 步骤九：将潑厂和Trt分别代入目标函数，将获得的函数值进行比较，若 CN 105160069 A 说明书 2/6 页 表示目标函数得到的较差解向量，winner表示目标函数得到的较优 解向量； 步骤十：轮流对winner和loser的每一项进行一次相互替换，然后比较替换后 的winner和替换前的winner两者在目标函数下的表现，如果替换后的表现较好，则替换 winner和loser对应的项； 步骤十一：更新均值和标准差：，更新PV矩阵； Np表示紧凑式教学优化算法虚拟人口数； 步骤十二：由步骤十一获得的PV矩阵随机生成; _8] 步骤十三：将汾广：和分别代入目标函数，将获得的函数值进行比较，若，生 步骤十四：将新的Stf1和Trt分别代入目标函数，将获得的函数值进行比较，若 步骤十五：轮流对winner和loser的每一项进行一次相互替换，然后比较替换后 的winner和替换前的winner两者在目标函数下的表现，如果替换后的表现较好，则替换 winner和loser对应的项； 步骤十六：更新均值和标准差：，更新PV矩阵； 步骤十七：更新 Trt+1= winner ; 步骤十八：t = t+1，判断t是否等于设定的值iteration maximum，若是，则转入 步骤十六，若否，则转入步骤五； 步骤十九：取Trt(t e )的最大值Trmax作为最优解向量 St〇pt; 步骤二十：将步骤十九所得的最优解Stcipt*的参数按需测量的机械参数的顺序 对应到所求的参数，获得系统模型中不能直接测得的机械参数。 步骤三：建立优化算法目标函数： 对于待测系统模型表达为Y = A ·Χ的线性系统，待测系统模型描述中X表示输入， Y表示输出，A为系统参数，建立其目标函数为： 待测系统中有η组输出，η为大于1的常数，i = 0，"·η ; 对于待测系统模型表达为Y = A Χ+Β的线性系统，建立其目标函数为： A和B均为系统参数； 对于待测系统模型表达为Y = A · X · C+B的线性系统，建立其目标函数为： A、B和C均为系统参数； 对于输入和输出为非线性关系的待测系统，可通过数学方法对其进行线性化，近 似为上述的线性系统，建立对应的目标函数； 对于输入输出不能表达为显式函数关系的待测系统，待测系统模型为G(R，X，Υ) =0, R为参数向量，建立其目标函数为： Ci. Np 等于 2〇。 λ 等于 10。 本专利技术的有益效果在于，改进的紧凑式教学优化算法的采用，使误差参数的测量 不受测量的参数和需求参数的维数所限制，不受传统方法解算方程时参数矩阵不能为奇异 矩阵的限制，误差参数的个数和需要测量的参数的个数不受限制；目标函数建立方法具有 很高的适应性，能应用于线性系统，可线性化的非线性系统，模块化的复杂系统等系统；算 法的性能良好，迭代速度快，根据需要可以满足各类精度需要。【附图说明】 图1为【具体实施方式】中的步骤五至步骤十一的原理示意图。 图2为【具体实施方式】中的步骤十二至步骤十六的原理示意图。【具体实施方式】 结合图1和图2说明本实施方式，本实施方式所述的一种基于改进的紧凑式教学 优化算法的机械参数软测量方法，所述方法包括如下步骤： 步骤一：建立待测系统模型，确定需测量的机械参数； 步骤二：测量系统模型中可以直接测得的参数，代入系统模型中，建立模型方程； 步骤三：建立优化算法目标函数： 对于待测系统模型表达为Y = A ·Χ的线性系统，待测系统模型描述中X表示输入， Y表示输出，A为系统参数，建立其目标函数为：，' 待测系统中有η组输出，η为大于1的常数，i = 0，"·η ; 对于待测系统模型表达为Y = A Χ+Β的线性系统，建立其目标函数为：_ A和B均为系统参数； 对于待测系统模型表达为Y = A · X · C+B的线性系统，建立其目标函数为：_ A、B和C均为系统参数； 对于输入和输出为非线性关系的待测系统，可通过数学方法对其进行线性化，近 似为上述的线性系统，建立对应的目标函数； 对于输入输出不能表达为显式函数关系的待测系统，待测系统模型为G(R，X，Υ) =0, R为参数向量，建立其目标函数为： 〇 对于复杂机械系统，通常可以将其模块化，各个子模块分别等价于一个简单模 型； 对于模型中含有误差参数的复杂函数的情况，直接套用上面的方法建立目标函数 进行计算可能会出现计算速度过慢的情况，也不符合紧凑式教学优化算法用于嵌入式和微 控制系统的初衷，可以考虑通过泰勒展开去除高次项或与之相近的方法在不影响最终结果 的情况下建立简化的近似模型加快迭代速度。 步骤四：测量待测系统的η组输入和对应的η组输出，代入目标函数； 步骤五：初始化计数器t = 0,均值初始值μ t = 0,方差初始值σ t = λ， λ等于10 ;i = 〇,…η ;当前第1页1 2 本文档来自技高网...
【技术保护点】
一种基于改进的紧凑式教学优化算法的机械参数软测量方法，其特征在于，所述方法包括如下步骤：步骤一：建立待测系统模型，确定需测量的机械参数；步骤二：测量系统模型中可以直接测得的参数，代入系统模型中，建立模型方程；步骤三：建立优化算法目标函数；步骤四：测量待测系统的n组输入和对应的n组输出，代入目标函数；步骤五：初始化计数器t＝0，均值初始值μt[i]＝0，方差初始值σt[i]＝λ；i＝0,…n；构成PV矩阵的初始值PV=μt[0]σt[0]......μt[n]σt[n],]]>PV的每一行包含高斯分布的一组均值和方差；其中t为算法迭代次数计数器的计数，由PV矩阵生成随机向量Trt的初始值；步骤六：由PV矩阵生成随机向量Stt，Stt中的每一个元素对应PV矩阵中一组均值和方差决定的高斯函数的随机值；步骤七：计算均值DMeant＝rand1×(Trt‑round(1+rand2(0,1))×μt)，round函数接受一个参数返回与参数最近的整数，rand1和rand2均为随机函数；步骤八：更新Sttnew=Stt+DMeant;]]>步骤九：将和Trt分别代入目标函数，将获得的函数值进行比较，若f(Trt)≤f(Sttnew),]]>则loser=Sttnew,]]>winner＝Trt，若f(Sttnew)<f(Trt),]]>则loser＝Trt，loser表示目标函数得到的较差解向量，winner表示目标函数得到的较优解向量；步骤十：轮流对winner和loser的每一项进行一次相互替换，然后比较替换后的winner和替换前的winner两者在目标函数下的表现，如果替换后的表现较好，则替换winner和loser对应的项；步骤十一：更新均值和标准差：μt+1[i]=μt[i]+1Np(winner[i]-loser[i])σt+1[i]=(σt[i])2+(μt+1[i])2-(μt[i])2+1Np(winner2[i]-loser2[i]),]]>更新PV矩阵；Np表示紧凑式教学优化算法虚拟人口数；步骤十二：由步骤十一获得的PV矩阵随机生成步骤十三：将和分别代入目标函数，将获得的函数值进行比较，若f(Sttnew2)<f(Sttnew),]]>生成新的Sttnew=Sttnew+rand3(Sttnew-Sttnew2),]]>若f(Sttnew)<f(Sttnew2),]]>生成新的Sttnew=Sttnew+rand3(Sttnew2-Sttnew);]]>步骤十四：将新的和Trt分别代入目标函数，将获得的函数值进行比较，若f(Trt)≤f(Sttnew),]]>则loser=Sttnew,]]>winner＝Trt，若f(Sttnew)<f(Trt),]]>则loser＝Trt，winner=Sttnew;]]>步骤十五：轮流对winner和loser的每一项进行一次相互替换，然后比较替换后的winner和替换前的winner两者在目标函数下的表现，如果替换后的表现较好，则替换winner和loser对应的项；步骤十六：更新均值和标准差：μt+1[i]=μt[i]+1Np(winner[i]-loser[i])σt+1[i]=(σt[i])2+(μt+1[i])2-(μt[i])2+1Np(winner2[i]-loser2[i]),]]>更新PV矩阵；步骤十七：更新Trt+1＝winner；步骤十八：t＝t+1，判断t是否等于设定的值iteration maximum，若是，则转入步骤十六，若否，则转入步骤五；步骤十九：取Trt(t∈[0,iteration maximum])的最大值Trmax作为最优解向量Stopt；步骤二十：将步骤十九所得的最优解Stopt中的参数按需测量的机械参数的顺序对应到所求的参数，获得系统模型中不能直接测得的机械参数。...
【技术特征摘要】
【专利技术属性】
技术研发人员：陈兴林，赵为志，董岳，刘杨，宋法质，王一光，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。