发明名称 识别机械加工表面轮廓分形特征长度尺度参数的小波方法
摘要 识别机械加工表面轮廓分形特征长度尺度参数的小波方法,属于机械加工表面质量技术领域,本发明从定性到定量确定了自相似过程的db2小波分解系数与分形特征长度尺度参数之间的关系,其主要步骤为:(一)、确定粗糙表面轮廓高度;(二)、对表面轮廓进行小波分解,得到各分解层的小波系数;(三)确定轮廓的分形特征长度尺度参数。该方法操作简单,思路简洁,只需将轮廓仪测得的轮廓高度数据导入到确定机中,再由MATLAB软件实现对所测轮廓高度数据序列进行小波分解,即可获得轮廓的分形特征长度尺度参数。本发明获得的表面轮廓的特征长度尺度参数,为进一步的机械结合部位的刚度、阻尼、热传导、电传导的分形建模,提供高精度的分形特征长度尺度参数。
申请公布号 CN105631120A 申请公布日期 2016.06.01
申请号 CN201510999583.8 申请日期 2015.12.28
申请人 太原科技大学 发明人 张学良;王余松;温淑花;兰国生;陈永会
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 太原高欣科创专利代理事务所(普通合伙) 14109 代理人 崔雪花
主权项 识别机械加工表面轮廓分形特征长度尺度参数的小波方法,其特征是:识别机械加工表面轮廓分形特征长度尺度参数的步骤如下:假设分形维数D已知,(一)、确定表面轮廓高度(1)、确定表面轮廓采样的点数A和采样间隔l,l的单位为μm;(2)、获取表面轮廓高度,将获得的粗糙表面的轮廓高度值序列记为z(i),(i=1,2,…,A),A为采样点的个数;(二)、对表面轮廓进行小波分解(1)、确定由计算机软件MATLAB对表面轮廓进行db2小波分解MATLAB计算机软件中将Daubechies小波记为dbN小波,小波的消失矩N值取为2,即为db2小波,db2小波分解时的最大分解层数M为:<img file="FDA0000893200770000011.GIF" wi="796" he="87" />其中,<img file="FDA0000893200770000012.GIF" wi="131" he="70" />向下取整函数;A—轮廓采样点数;lw—滤波器长度,lw=2×N=2×2;(2)、对确定的表面轮廓高度值序列进行M层小波分解,得到各分解层的小波系数<img file="FDA0000893200770000017.GIF" wi="79" he="80" />其中,m为分解层数,k为小波系数序号,其值为1,2,…,2<sup>m</sup>‑1;(3)、在半对数坐标系<img file="FDA0000893200770000013.GIF" wi="350" he="111" />下,即横坐标为分解层数m,纵坐标为小波系数平方取均值的对数<img file="FDA0000893200770000014.GIF" wi="302" he="117" />在分解层从1到M的M个数据点,找出其中相邻且接近于一条直线的点,将所述点所在分解层确定为有效分解尺度;(三)、确定表面轮廓的分形特征长度尺度参数G<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mi>G</mi><mo>=</mo><msup><mn>2</mn><mrow><msub><mi>Alog</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>E</mi><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msubsup><mi>d</mi><mi>k</mi><mi>m</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mo>)</mo></mrow><mo>+</mo><mi>B</mi></mrow></msup></mtd></mtr><mtr><mtd><mrow><mi>A</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mi>B</mi><mo>=</mo><mi>a</mi><mo>(</mo><mfrac><mn>1</mn><mrow><mi>D</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mn>1</mn><mo>)</mo><mo>-</mo><mfrac><mn>1.5</mn><mrow><mi>D</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mi>m</mi><mo>+</mo><mn>1.04</mn><mi>m</mi></mtd></mtr><mtr><mtd><mi>a</mi><mo>=</mo><mn>470.5</mn><msup><mi>l</mi><mrow><mo>-</mo><mn>0.003042</mn></mrow></msup><mo>-</mo><mn>451.4</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000893200770000015.GIF" wi="676" he="447" /></maths>其中,m为有效分解尺度中的任一值,<img file="FDA0000893200770000016.GIF" wi="65" he="62" />为各分解层的小波系数,l为采样间隔。
地址 030024 山西省太原市万柏林区窊流路66号