本发明专利技术涉及一种基于变尺度混沌模拟退火算法的圆度误差评定方法，先将模拟退火算法作为寻优的粗搜索，然后利用混沌优化算法的高效性作为寻优的细搜索，确保寻优过程简洁高效，该模拟退火算法的优化结果不依赖于初始值，因此先采用模拟退火算法能保证寻找到实际全局最优区域，而混沌优化算法具有很好的执行效率，适合在全局最优区域内进行细搜索，从而搜索到全局最优解，因此在优化过程中具有更小的迭代次数，同时本发明专利技术提出的混沌映射xtin=cos(1/xtin2)能更好地实现混沌寻优，在执行过程中[a1,i+1，b1,i+1]、[a2,i+1，b2,i+1]根据解空间变化而不断减小，使得寻优过程不断逼近全局最优解。
【技术实现步骤摘要】
本专利技术涉及一种。
技术介绍
圆度误差是衡量加工零件是否是理想圆的重要指标。评定圆度误差的过程就是将理想圆和被测工件截面的实际形状进行比较的过程。测量工件的参数有圆度、直径、跳动、 同轴度等。测量时的形位参数、偏心误差、工件回转误差以及圆度误差将会使测量结果与计量室检测结果差异较大。其中尤其以圆度误差影响最大，因此圆度误差的测量和评定决定了加工零件的回转精度、加工精度以及装配精度。在工程实践中计算圆度常用的方法有罚函数法、单纯形法、遗传法、曲率法和区域搜索法等。缪立栋等所著的《基于多种群实数编码遗传算法的圆度误差评定》在实数编码的基础上引进多种群遗传算法，提高了全局搜索速度和精度。李秀明等所著的《圆度误差评定的曲率法研究》提出了基于曲率的圆度误差评定方法，根据曲率半径不断筛选候选点， 直至剩余满足判别准则的3点，筛选条件是不断去掉曲率最大或曲率最小的点。黄富贵等所著的《基于区域搜索的圆度误差评定方法》选取被测圆周上均匀分布的3个点作为初始圆，求取初始圆的圆心，然后以此圆心作为正方形搜索的中心确定搜索的步长和大小，进而得到圆度误差，该方法避免了最小二乘法带来的非线性运算，提高了圆度误差的计算精度。评定圆度误差常用的方法有最小二乘法、最小区域圆法、最小外接圆法和最大内切圆法。其中最小二乘法最简单，但是圆度误差较大，且不符合国家标准规定的圆度误差测量标准。最小外接圆法和最大内切圆法均要对数据采集进行线性化处理，同时为保证误差精度需要对计算附加额外条件，因此评定过程较为复杂，不易被工程技术人员掌握。混沌理论是20世纪物理学三大重大专利技术之一。作为一门新兴学科，其研究的理论成果在计算机科学、信息科学、最优化理论、物理学、经济学和气象学等多个领域得到广泛应用。混沌是非线性系统中普遍存在的现象，其行为表现为不确定性、不可预测、不可重复等特点，集中体现为貌似随机的不确定运动。混沌优化的核心是利用混沌的不确定性和遍历性产生相应的混沌映射，把该映射作为随机数发生器产生相应的扰动因子来搜索全局最优点。优化的目标函数为y=f(x)， X e R+，混沌优化的步骤为( I)算法初始化，生成自变量X相同维度的随机数；(2)按照混沌迭代发生器产生混沌变量；(3)根据混沌变量进行判断，若同时满足循环次数不大于最大混沌迭代次数和混沌现态值小于预设初始值，则当前最优值为现态值，与此同时初始值变为现态值，若不满足，则保持预设初始值不变，同时转至步骤(2);(4)判断循环次数是否大于最大混沌迭代次数，若大于则停止搜索，输出最优值。模拟退火算法是Kirkpartrick等人利用固体退火的物理过程建立相应的数学模型而得到的自启发优化方法。它是一种随机搜索方法，与普通的局部搜索方法相比，其在搜索过程中按一定概率接受比目标解更差的搜索点，因此能跳出局部极值点，最大限度的搜索到全局最优解。模拟退火过程中初始温度的大小影响优化问题的最终解，初始温度越高，退火过程经历时间越长，则寻优结果越接近于最优解；初始温度的选择决定了算法寻优的质量和效率。模拟退火算法的步骤为( I)算法初始化，设置初始退火温度T0、随机产生一组初始向量作为初始解和初始能量值；(2)对初始解在一定邻域范围内进行随机扰动，产生搜索解及其对应的现态能量值；(3)若现态能量值小于初始能量值，则判断搜索有效，继续执行步骤(2)，否则对邻域范围以概率为exp (能量差值/kT)接受,而以Ι-exp (能量差值/kT)的概率拒绝接受现态数值；(4)根据随机扰动状态转移进行判断，若同时满足循环次数不大于最大迭代次数和现态能量转移函数，则重复执行步骤⑵和步骤(3);(5)执行退温函数，降低退火温度；(6)重复执行步骤(2)和(5)，判断循环次数是否大于最大迭代次数，若大于则停止搜索，输出最优值。混沌优化过程中单纯增大迭代次数并不能提高搜索的效率，甚至会造成搜索局限于局部最优点附近。模拟退火利用马科夫链描述退火过程中由一个状态跳变到另一个状态转换关系，计算对应温度的马科夫链稳定状态的转移概率，然后再降低退火温度，直到固体内部达到稳定状态。当温度保持不变和温度趋近于O时，马科夫链的极限分布存在，并且在温度为O时，马科夫链的极限分布和的概率为I。
技术实现思路
本专利技术的目的在于提供一种， 克服了传统算法易收敛于局部极小值点和运算效率较低的缺点，具有较强的鲁棒性。一种，具体包括如下步骤(I)优化的目标为 F = min,其中权利要求1.一种，其特征在于包括如下步骤全文摘要本专利技术涉及一种，先将模拟退火算法作为寻优的粗搜索，然后利用混沌优化算法的高效性作为寻优的细搜索，确保寻优过程简洁高效，该模拟退火算法的优化结果不依赖于初始值，因此先采用模拟退火算法能保证寻找到实际全局最优区域，而混沌优化算法具有很好的执行效率，适合在全局最优区域内进行细搜索，从而搜索到全局最优解，因此在优化过程中具有更小的迭代次数，同时本专利技术提出的混沌映射xtin=cos(1/xtin2)能更好地实现混沌寻优，在执行过程中、根据解空间变化而不断减小，使得寻优过程不断逼近全局最优解。文档编号G06F19/00GK102982240SQ201210470650公开日2013年3月20日 申请日期2012年11月19日 优先权日2012年11月19日专利技术者傅文渊, 凌朝东, 李国刚, 莫冰, 黄瑞敏 申请人:华侨大学本文档来自技高网...
【技术保护点】
一种基于变尺度混沌模拟退火算法的圆度误差评定方法，其特征在于包括如下步骤：（1）优化的函数为F=min[Rer(x0,y0)]，其中Rer(x0,y0)=mini∈A1∪A2[maxA2((xi+x0)2+(yi+y0)2)-minA1((xi+x0)2+(yi+y0)2)]为优化函数，A1和A2分别表示包容被测工件轮廓的最小同心圆和最大同心圆，(x0,y0)为所求理想圆圆心，(xi,yi)为实际工件的采样坐标，i∈[1,n]，其中n为工件的坐标点数，x0,y0为优化自变量，[a1,i,b1,i]和[a2,i,b2,i]分别为为优化变量x0和y0的取值区间；设置初始退火温度T、随机产生一组初始向量作为初始解xintial和yintial、初始能量值E0、最大迭代次数M1、混沌迭代次数M2、M3，其中M3设置为0.2*M2~0.8*M2、相对误差因子e；（2）对初始解xintial和yintial在邻域范围内进行随机扰动，初始解扰动按照如下方式产生：xt=xintial+(rand?0.5)*exp(?2*M1/kT);yt=yintial+(rand?0.5)*exp(?2*M1/kT);其中rand是(0,1)内的随机数，Et为对应的现态能量值，k为波尔兹曼常数；（3）若现态能量值Et小于初始能量值E0，则搜索有效，继续执行步骤(2)，否则对邻域范围以概率为exp(ΔE/kT)接受，而以1?exp(ΔE/kT)的概率拒绝接受现态能量值，其中ΔE=Et?E0，k为波尔兹曼常数；（4）根据随机扰动状态转移进行判断，若同时满足循环次数不大于最大迭代次数M1和现态能量转移函数，则重复执行步骤(2)和步 骤(3)；（5）执行退温函数Tk+1=Tk/(M1+1)，降低退火温度；（6）重复执行步骤(2)至(5)，判断循环次数是否大于最大迭代次数M1，若大于则停止搜索，输出当前最优解xtin、ytin和相应的圆度误差值Rer(xtin,ytin)；（7）将xtin和ytin作为混沌迭代的初始值生成相应的混沌迭代变量，其混沌迭代的具体关系式为：xtin=cos(1/xtin2)，ytin=cos(1/ytin2)，设定混沌迭代次数为p；（8）将生成的混沌迭代变量执行M2次混沌搜索，现态混沌变量xtink=a1,i+(b1,i?a1,i)*xtin，ytink=a2,i+(b2,i?a2,i)*ytin，若现态混沌变量xtink，ytink满足Rer(xtink，ytink)bj,i，则bj,i+1=bj,i，其中j＝（1,2)；（10）对xtmin，ytmin进行混沌变换，同时迭代M3次重复执行步骤(7)；（11）重复执行步骤(7)至步骤(9)，直到迭代M3次后xtmin，ytmin对应的函数值偏差不超过相对误差因子e，该相对误差因子e的数值 根据实际问题变化而不同，一般为10?3~10?6；（12）重复执行步骤(11)，直到迭代次数为p时，算法寻优结束，输出最优解xtmin，ytmin。...
【技术特征摘要】
【专利技术属性】
技术研发人员：傅文渊，凌朝东，李国刚，莫冰，黄瑞敏，
申请(专利权)人：华侨大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。