表面纹理方向性研究 - 图文

[18] P Podsiadlo，G W stachowiak，Fractal-Wavelet based classification of tribological

surfaces．Wear，2003，254：1189－1198．

[19] 李丽伟，董申，陈凯．表面微观形貌定量表征中几种新方法的应用．中国机

械工程，2004，13（9）：1701－1705．

[20] 陈庆虎，李柱．表面粗糙度评定的小波基准线．计量学报，1998，19 （4）：

254－257．

[21] 陈庆虎，李柱．表面粗糙度提取的小波频谱法．机械工程学报，1999，35（3）：

41－43．

[22] 宋康，蒋庄德，李兵．基于小波函数迭代的外圆磨削表面分形分析．计量学

报，2004，7：211－214．

[23] 王安良，杨春信．评价机械加工表面形貌的小波变换方法．机械工程学报，

2001，37（8）：65－74．

[24] S H Lee，H Zahouani，R Caterini，T G.mathia．Morphological characterization of

engineered surfaces by wavelet transform．Tools Manufact, 1998，38 (5)：581－589．

[25] 章毓晋．图象工程[M]．清华大学出版社，2004

[26] 刘小君．粗糙表面形貌的表征及应用．合肥工业大学学报（自然科学版），1997，

20（3）：116－122．

[27] 余成波．数字图像处理及MATLAB实现[M]．重庆大学出版社，2003 [28] 董长红，赖志国，余啸海．Matlab图象处理与应用[M]．国防工业出版社，

2004

[29] Nadir Patir，H S Cheng．An Average Flow Model for Determining Effects of

Three-Dimensional Roughness on Partial Hydrodynamic Lubrication．ASME， 1978，100：12－14．

[30] Nadir Patir，H S Cheng．Application of Average Flow Model to Lubrication Between Rough Sliding Surfces．ASME，1979，101：220－230．

[31] 陈虎，王守尊，周朝辉．基于数学形态学的图像边缘检测方法研究．工程图 学学报，2004 ，2：112－115．

- 23 -

致 谢

值此论文完成之际，我的心情很激动，在这里首先向我的指导老师方新燕（讲师）致以衷心的感谢和崇高的敬意！感谢她在这段时间里所给予的帮助和指导，本人能够按时保质完成这篇论文，完全得益于方老师。从论文的选题、初稿、定稿的过程中，方老师倾注了大量的时间和精力，感谢方老师这三个多月以来的悉心关怀和细心指导！方老师学识渊博，每次与她交流都有很大的收获。她严谨细致、一丝不苟的作风一直是我工作、学习中的榜样。在此再次表示感谢。

其次要向我们毕业设计的专业负责人张少伍（副教授）表示衷心的感谢，感谢他在百忙中对我们进行指导和为我们解决设计期间遇到的一些问题。此外，在做毕业设计的过程中，有多少可敬的师长、同学、朋友给了我无尽的帮助，在这里请接受我诚挚的谢意！

最后对在论文中所引用文献的作者表示感谢！

- 24 -

附录Ⅰ MATLAB的简介

MATLAB语言的开发者是Cleve Moler教授，他在数值分析，特别是在数值线性代数的领域中很著名。由于当时用其他高级程序语言编程极为不便，便构思并开发了MALAB,这一软件了利用了当时数值线性代数领域最高水平的EISPACK和LINPCK两大软件包的Fortran语言编写了集指令翻译、科学计算与一体的交互式软件系统。

目前，MATLAB语言的设计思想代表了当前计算机语言的发展方向，它已经成为了一种具有广泛应用前景的全新的计算机高级编程语言，有人称它为“第四代”计算机语言，MATLAB已经成为国际上最流行的科学与工程计算的软件工具。 MATLAB即是一种直观、高效的计算机语言，同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数字和工程函数，工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。

MATLAB语言有一下特点：编程效率高、便于用户使用、扩充能力强、语句语法简单、高效方便的矩阵和数组运算、方便的绘图功能。

- 25 -

附录Ⅱ 部分源程序

1、 随机表面产生的MATLAB程序为：

设表面的粗糙度高度矩阵为B，B（i，j）是相应的x=i?l,y?j?l 处粗糙度幅度值，?l 为采样长度。

Rx(?)?E(B(i,j)*B(i??*,j)自相关函数的定义 ：

Ry(?)?E(B(i,j)*B(i,j??*)

*???*?l,?0.5 定义为轮廓自相关长度减小到 ??0时值的50%的?的值。

自相关函数的数字化方法：

NM?m1Rx(?)???B(i,j)*B(i?m,j)N(M?m)j?1i?1 N?nM1Ry(?)?B(i,j)*B(i,j?n)??M(N?n)j?1i?1

I=imread('p9c.tif'); I=double(I);

I1=-4.155+6.551/256*I; Rz=0;a=1;b=1; for i=1:100 for j=1:100

Rz=Rz+I1(i,j)*I1(i,j)/(100*100); end end for m=1:50 Rx=0; for i=1:100 for j=1:100

Rx=Rx+I1(i,j)*I1(i+m,j)/(100*(100-m)); end end

x=abs(Rx-Rz/2.0); if x

- 26 -