表面粗糙度测量中的高斯递归滤波算法摘要基于中心极限定理，提出了一个高斯滤波器逼近模型，并利用冲激响应不变法设计了一个数字高斯滤波器。通过级联方法，该滤波器的幅度传输特性偏差可控制在(-0.28%+1.14%)以内。该滤波器采用递归实现算法，算法简便，易于实现。在递归算法的基础上，给出了一种预测边缘初始值的方法，有效抑制了递归算法的边缘效应。实验结果表明，该滤波器能够稳定可靠地提取工程表面中线，在普通计算机上提取11200点数据表面中线仅耗时170ms，效率很高。关键词计量学表面粗糙度高斯滤波器边缘效应TP2文献标识码A国家标准学科分类与代码520.604GaussianRecursiveFilteringalgorithmtoSurfaceRoughnessMeasurementAbstractBasedonthecentrallimittheorem,anapproximationmodelofGaussianfilterisputforward.Then,adigitalGaussianfilterisdesignedbycombinationwithimpulseinvarianttransformationandthismodel.Throughthecascadeuse,theamplitudetransmissiondeviationofthisfiltercanbecontrolledwithin(-0.28%+1.14%).Itsrecursivealgorithmissosimplethatitisimplementedeasily.Amethodofforecastingtheinitialvalueoftheedgeisproposedtorestraintheendeffectcausedbythisrecursivealgorithm.Experimentalresultsshowthatthisfiltercanextractthemeanlineofengineeringsurfacestablyandreliably,onceextractionof11,200pointdataonageneralcomputeronlyconsumes170ms,efficiencyisdemonstrated.KeywordsMetrologySurfaceroughnessGaussianfilterEndeffect1引言在表面粗糙度测量中，表面中线的确定是一个关键问题，它对评定表面特征参数有着重要的作用[1]。中线的提取是依靠轮廓滤波器完成的。国际标准曾经规定2RC滤波中线作为表面粗糙度评定基准线，但由于2RC中线本身的非线性相移会扭曲滤波结果，最终被淘汰。1996年，ISO11562推荐高斯滤波器中线作为新的国际标准，时至今日，仍然被广泛采用[2]。高斯滤波的幅度特性是一个高斯函数，既要精确又要快速地进行高斯滤波并不是一件容易实现的事情。围绕着如何准确快速地实现国际标准中规定的传输特性，国际上众多学者进行了广泛的研究，继而出现了许多高斯滤波器的逼近算法和快速算法[3-5]。本文叙述了一种用高斯函数的逼近法和冲激响应不变法设计出的数字高斯滤波器，并给出了其零相移的递归滤波算法，算法简单，易于实现。但是由于采用递归的计算方法，其运算结果难免会出现边缘数据失真，这种失真也被称为边缘效应。给出了一种预测边缘初始值的方法，在最小二乘的原则下抑制边缘失真，有效保留了边缘特征。2高斯滤波器的逼近方法根据ISO11562[2]，高斯滤波器的权函数和传输特性分别为：(1)(2)---本文来源于网络，仅供参考，勿照抄，如有侵权请联系删除---其中是常数，为波长，是截止波长。，为采样间距，为采样点数，，为截止波长内采样点数。高斯滤波器具有线性相位，而且当，高斯滤波器的幅频特性为50%，所以。对于高斯函数，可以由泰勒级数展开为。当时，(3)定义函数(4)由中心极限定理可知，的级联可以逼近高斯函数，表示如下：(5)其中，为级联标识符。例如，时，式(5)为二阶级联，时，式(5)为十六阶级联。因为时的幅度传输特性为50%，所以(6)对于单个逼近单元式(4)，由于，，所以(7)将代入上式，得到(8)它的极点为，只有在左平面，所以(9)即为设计高斯逼近滤波器的模拟滤波器原型。3IIR型高斯逼近滤波器对式(9)应用冲激响应不变法，设计IIR型数字滤波器[6,7]。得到如下：(10)其中，同样对应的z变换为(11)由式(10)(11)的两个数字滤波器就可以构成零相移滤波器，有：(12)由，可得(13)以式(12)(13)这样的零相移滤波器，经2n-1阶级联，就可构成IIR型高斯逼近滤波器。(14)图1给出了一阶，四阶，八阶，十六阶高斯逼近滤波器的幅度传输特性图。可以看出，逼近滤波器的幅度传输特性随着级联阶数的增加逐渐逼于高斯滤波器。图2给出了逼近滤波器相对高斯滤波器的传输特性偏差。由图2可知，阶数越高，偏差越---本文来源于网络，仅供参考，勿照抄，如有侵权请联系删除---...