【正文】 所以=1，=，=，=，=。將模擬濾波器幅度平方的如下關(guān)系：，代入式(11)，則： (212)其在s平面的左半平面只有一個極點，由此可得巴特沃思濾波器系統(tǒng)函數(shù)為： (213)即為設(shè)計高斯逼近濾波器的模擬濾波器原型。Hd4(z)Hd4(z1)Hd4(z)Hd4(z1)— 1表面輪廓中線 8級表面粗糙度圖11 表面粗糙度測量的高斯濾波流程當或n=Nc時，要求傳輸系數(shù)為50%，所以 （218）由式(18)可以求得a，它是對逼近濾波器原型巴特沃思濾波器的極點進行修正，當Nc=1600時，可得a=。從上式可以看出，每次運算都只有少量的加減和乘法操作。對于16級零相移濾波器[22]，公式(22)。用高斯函數(shù)的逼近法設(shè)計出的高斯遞歸濾波算法，計算量小，計算速度快，濾波精度高，有效地解決了傳統(tǒng)方法計算效率和濾波精度兩者之間不可調(diào)和的矛盾。可見，適當增加濾波器節(jié)數(shù)，在不增加很多計算量的情況下，可以達到很高的精度。命令行運行方式可以通過直接在命令窗口輸入命令行來實現(xiàn)計算或作圖功能。Matlab的應(yīng)用范圍非常廣，包括[25]信號和圖像處理、通訊、控制系統(tǒng)設(shè)計、測試和測量、財務(wù)建模和分析以及計算生物學(xué)等眾多應(yīng)用領(lǐng)域。包括Matlab桌面和命令窗口、歷史命令窗口、編輯器和調(diào)試器、路徑搜索和用于用戶瀏覽幫助、工作空間、文件的瀏覽器。（2）簡單易用的程序語言Matlab一個高級的矩陣/陣列語言，它包含控制語句、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入和輸出和面向?qū)ο缶幊烫攸c。而且這種語言可移植性好、可拓展性極強，這也是Matlab能夠深入到科學(xué)研究及工程計算各個領(lǐng)域的重要原因。在通常情況下，可以用它來代替底層編程語言，如C和C++。可用于科學(xué)計算和工程繪圖。（6）實用的程序接口和發(fā)布平臺新版本的Matlab可以利用Matlab編譯器和C/C++數(shù)學(xué)庫和圖形庫，將自己的MMatlab程序自動轉(zhuǎn)換為獨立于Matlab運行的C和C++代碼。工具箱是Matlab函數(shù)的子程序庫，主要包括信號處理、控制系統(tǒng)、神經(jīng)網(wǎng)絡(luò)、模糊邏輯、小波分析和系統(tǒng)仿真等方面的應(yīng)用。F:\Matlab\work\39。 a=。 for i=1:7 for j=2:m r(j)=r(j1)+a*[ff(j)r(j1)]。 plot(j,z,39。)。4 Mat的運用Mat是mathTools公司推出的一個能將M文件轉(zhuǎn)化成相同功能C++代碼的工具。在安裝過程中出現(xiàn)選擇編譯器對話框，選擇“是”（圖41）。圖41 安裝圖示安裝完成后，啟動界面MIDEVA(圖42)。（如圖43）。 Wizard(exe)創(chuàng)建工程在VC中用MFC Wizard(exe)創(chuàng)建一個基于對話框名為Example的工程。圖46 拷貝文件圖示：include include 在void CExamleDlg::OnButton1()中分別添加一個初始化類庫調(diào)用函數(shù)和一個結(jié)束類庫調(diào)用函數(shù)（如圖47）。坐標系：確定表面機構(gòu)參數(shù)的坐標體系，其軸線形成一右旋笛卡爾坐標系，X軸與中線方向一致，Y軸也處于實際表面上，Z軸則在材料到周圍介質(zhì)的外延方向上。注：原始輪廓是評定原始輪廓參數(shù)的基礎(chǔ)。中線：具有幾何輪廓形狀并劃分輪廓的基準線。輪廓谷：連接兩相鄰交點向內(nèi)（從周圍介質(zhì)到材料）的輪廓部分。Rv：在一個取樣長度內(nèi)，最大的輪廓谷深Rz：在一個取樣長度內(nèi)，最大的輪廓峰高和最大的輪廓谷深之和的高度Rc：在一個取樣長度內(nèi)，輪廓單元高度的平均值Rt：在評定長度內(nèi)最大輪廓峰高和最大輪廓谷深之和平均值Ra：在一個取樣長度內(nèi)，縱坐標Z（x）絕對值的算術(shù)平均數(shù)Rq：在一個取樣長度內(nèi)，縱坐標值Z（x）的均方根Rsk：在一個取樣長度內(nèi)，縱坐標值Z（x）三次方的平均值與Rq的比值Rku：在一個取樣長度內(nèi)，縱坐標Z（x）四次方的平均值與與Rq的比值間距參數(shù)RSm：在一個取樣長度內(nèi),輪廓單元寬度的平均值混合參數(shù)R△ｑ：在一個取樣長度內(nèi)，縱坐標斜率的均方根值曲線和相關(guān)參數(shù)Rmr（c）：在一個取樣長度內(nèi)，在給定水平位置上輪廓實體材料長度與評定長度的比率：輪廓支承率隨水平位置而變的關(guān)系曲線Rmｒ：給定支承比率的兩個水平截面之間的垂直距離在計算出粗糙度的濾波中線后需要以下計算過程可算得國際標準中所規(guī)定的參數(shù)，其計算過程如下：%高度參數(shù)l=x(1601:3200)。c=ll(b)。c=ll(b)。%1個取評定長度yy=y(1601:9600)。RRp=pc。RRz=cp。 end Ra=mysum/1600。 Rq=sqrt(g)。 d=Rq^3。 end c= mysum/1600。由于測試結(jié)果有上百萬個數(shù)據(jù)量，且數(shù)據(jù)存儲在數(shù)據(jù)庫中，為了便于研究處理數(shù)據(jù)的方法，獲得可靠的測試結(jié)果，本文探討利用Microsoft Visual C++ 。 ADO數(shù)據(jù)庫訪問技術(shù)ADO(ActiveX Data Object)是一種基于COM(組件對象模型)的自動化接口(IDispatch)技術(shù)，并以O(shè)LE DB(對象連接和鑲?cè)氲臄?shù)據(jù)庫)為基礎(chǔ)，經(jīng)過OLE DB包裝后的數(shù)據(jù)庫訪問技術(shù)，利用它可以快速的創(chuàng)建數(shù)據(jù)庫應(yīng)用程序。(3)利用建立好的連接，通過三個智能指針執(zhí)行查詢、處理結(jié)果記錄集和SQL命令，使用完畢后關(guān)閉連接釋放對象。（2）MATLABMatlab與Visual C++的接口技術(shù)[29]Matlab系統(tǒng)提供一個非常重要的組件—Matlab應(yīng)用程序接口，它可以向MSATLAB環(huán)境傳送數(shù)據(jù)也可從Matlab環(huán)境接收數(shù)據(jù)；Visual C++可以調(diào)用Matlab編寫的程序，在這里數(shù)據(jù)處理系統(tǒng)采用Matlab引擎庫，Matlab引擎庫包含了若干個控制Matlab引擎的函數(shù)：engOpen ：啟動MATLAB引擎 engClose ：關(guān)閉MATLAB引擎 engGetArray ：從MATLAB引擎中獲取一個MATLAB矩陣 engPutArray ：向MATLAB引擎發(fā)送一個MATLAB矩陣 engEvalString ：執(zhí)行于一個MATLAB命令或自編的MATLAB函數(shù) engOutputBuffer ：創(chuàng)建一個存儲MATLAB文本輸出的緩沖區(qū) （3）Matlab與Visual C++的數(shù)據(jù)傳輸[30] Matlab與Visual C++的數(shù)據(jù)傳輸需要三個步驟，第一完成Visual C++讀取ACCESS數(shù)據(jù)庫的大量數(shù)據(jù)；第二把Visual C++讀取到的數(shù)據(jù)傳輸?shù)組atlab中，并在利用Matlab中自編的函數(shù)進行處理；第三把處理結(jié)果傳輸?shù)絍isual C++中。(6)合并測試數(shù)據(jù)在數(shù)據(jù)處理系統(tǒng)中，處理的測試數(shù)據(jù)往往對幾次測試的數(shù)據(jù)需要一起處理，系統(tǒng)可以把幾次測試代號所對應(yīng)數(shù)據(jù)按一定的順序進行合并。系統(tǒng)只需把數(shù)據(jù)處理方法的程序放入系統(tǒng)的處理函數(shù)里，就可以實現(xiàn)對測試數(shù)據(jù)的統(tǒng)計處理，本系統(tǒng)在數(shù)據(jù)存儲好后只需按“作圖求值”按鈕，就可以獲得我們所需的測試數(shù)據(jù)結(jié)果和統(tǒng)計圖。wt39。)。 fprintf(fid, 39。% %\n39。)。fprintf(fid, 39。% %\n39。)。fprintf(fid, 39。% %\n39。)。fprintf(fid, 39。% %\n39。 Visual C++應(yīng)用部分以下指出的此段程序為VC編程語句，可直接與Mat編譯語句聯(lián)合，這就是聯(lián)合編譯的方便之處。 pListResetContent()。 下圖為Visual C++所設(shè)計的工作界面（圖62）圖62 濾波器對話框點擊“粗糙度數(shù)據(jù)”按鈕時顯示所計算數(shù)據(jù)如下圖（圖63）圖62 粗糙度數(shù)據(jù)在點擊“濾波顯示”按鈕時顯示如下濾波結(jié)果（圖64） 圖64 濾波圖形顯示7結(jié)束語在論文中分別介紹了數(shù)字濾波器的研究背景，相關(guān)理論，以及對章節(jié)的安排，在第二章中著重介紹了高斯濾波理論，這也是此軟件設(shè)計的理論支持。對于Mat在聯(lián)合編譯中的作用我在第四章中作了簡要的介紹，以一個簡單的示例來介紹了它如何實現(xiàn)編譯。實際上VC與Matlab的接口實現(xiàn)方法有很多種，Mat只是其中一種，本設(shè)計主要運用的接口方式就是Mat,其它接口方式就不作介紹。由于在此之前對于濾波器的相關(guān)知識不足，在張老師的指導(dǎo)下，補充了對于濾波器的認識。作為同學(xué)她為我提供了力所能及的幫助。 dMm(x)。 dMm(m)。 dMm(i_test1_v0)。 line 1 f:/ call_stack_begin。 line 4 f:/_ z = x(colon(,))。 line 8 f:/_ y() = 。i_test1_i0()。 for (j_test1_i1=0。 line 12 f:/ } line 13 f:/_ j_test1_v2 = colon(,)。j_test1_i2++) { forelem(j_,j_test1_v2,j_test1_i2)。 line 17 f:/ } line 18 f:/_ j_ = colon(,m+)。}void CTestDlg::OnShuju() { // TODO: Add your control notification handler code here initM(MATCOM_VERSION)。 dMm(a)。 dMm(i_)。 dMm(j_gaoshilvbo_v2)。 dMm(b)。 dMm(Rz)。 dMm(RRz)。 dMm(Ra)。 dMm( \ i_gaoshilvbo_v5)。 dMm(Rku)。 line 2 f:/_ x = load(TM(F:\\Matlab\\work\\),Mc())。 line 6 f:/_ m = 。 int i_gaoshilvbo_i0。 line 10 f:/_ j_gaoshilvbo_v1 = colon(,m)。j_gaoshilvbo_i1++) { forelem(j_,j_gaoshilvbo_v1,j_gaoshilvbo_i1)。 for (j_gaoshilvbo_i2=0。 line 15 f:/ } line 16 f:/_ ff = y。 line 24 f:/_ b = find(l==p)。 line 29 f:/_ b = find(l==p)。 line 33 f:/ //Rz line 35 f:/ //Rc line 37 f:/_ xx = x(colon(,))。 line 41 f:/_ c = yy(b)。 line 45 f:/_ c = yy(b)。 line 51 f:/_ i_gaoshilvbo_v3 = colon(,)。i_gaoshilvbo_i3++) { forelem(i_,i_gaoshilvbo_v3,i_gaoshilvbo_i3)。 line 57 f:/_ i_gaoshilvbo_v4 = colon(,)。i_gaoshilvbo_i4++) { forelem(i_,i_gaoshilvbo_v4,i_gaoshilvbo_i4)。 line 61 f:/ //Rq line 63 f:/_ mysum = 。i_gaoshilvbo_i5()。 line 68 f:/_ d = mpower(Rq,)。 int i_gaoshilvbo_i6。 line 73 f:/_ mysum = mysum+mpower(x(),)。 line 77 f:/ //Rku line 79 f:/_ fid = fopen(TM(),TM(wt))。 line 83 f:/_ fprintf(fid,TM(% %\\n),(CL(Rv)))。 line 87 f:/_ fprintf(fid,TM(% %\\n),(CL(Rt)))。 line 91 f:/_ fprintf(fid,TM(% %\\n),(CL(Rq)))。 line 95 f:/_ fprintf(fid,TM(% %\\n),(CL。 line 93 f:/_ fprintf(fid,TM(% %\\n),(CL(Rsk)))。 line 89 f:/_ fprintf(fid,TM(% %\\n),(CL(Ra)))。 line 85 f:/_ fprintf(fid,TM(% %\\n),(CL(Rz)))。 line 81 f:/_ fprintf(fid,TM(% %\\n),(CL(Rp)))。 li