【正文】 )==rmin B=b。%調(diào)用數(shù)據(jù)輸入函數(shù)讀取數(shù)據(jù) [h,i,k,m,n,l]=sjsrpm()。%求Y軸最大值b=max(n)。subplot(39。,x0)。subplot(39。mesh(m,n,l4)。,[,])。rmin=min(r)。LB=l(B)。z=inv(A)*D。%繪制被測曲面示意圖subplot(39。,[,])。mesh(m,n,l2)。subplot(39。)。%調(diào)用子函數(shù)求系數(shù)x=djpmht1(m,n,l)。LA=l(A)。A=[m(1),n(1),l(1)。對角線法繪圖程序function H=djpmht(m,n,l)。,[,])。hold on。l3=(MA*x(1)+x(2)*NA+x(3)*LAx(1)*mx(2)*n)/x(3)。,[,])。position39。一致性——不與其它已有的界面風(fēng)格截然向左，而與已經(jīng)開發(fā)的界面風(fēng)格盡量一致。計(jì)算機(jī)屏幕由系統(tǒng)默認(rèn)，無需設(shè)計(jì)，其它結(jié)點(diǎn)按設(shè)計(jì)原則根據(jù)實(shí)際需要設(shè)計(jì)。（2）在軸屬性編輯框中，設(shè)置如下屬性值：Box off %軸不封閉；Units normalized %采用相對度量單位，縮放時(shí)保持比例；Tag axes1 %；（3）在圖形區(qū)上方的靜態(tài)文本的屬性編輯框中，作如下屬性設(shè)置；FontSize 12 %字體大??；FontUnits normalized %采用相對度量單位，縮放時(shí)保持字體比例；HorzontalAlignment Center %字體中心對齊；String 測量平面模擬圖 %圖形名稱； Units normalized %采用相對度量單位，縮放時(shí)保持該區(qū)比例。x4,x5,x6。（17）在圖形窗口右邊上平板里的右列上靜態(tài)文本的屬性編輯框中，作如下屬性設(shè)置；FontSize 12 %字體大??；FontUnits normalized %采用相對度量單位，縮放時(shí)保持字體比例；HorzontalAlignment Center %字體中心對齊；String /mm %相應(yīng)編輯框的名稱； Units normalized %采用相對度量單位，縮放時(shí)保持該區(qū)比例。（15）在圖形窗口右邊上平板里的左列中靜態(tài)文本的屬性編輯框中，作如下屬性設(shè)置；FontSize 12 %字體大小；FontUnits normalized %采用相對度量單位，縮放時(shí)保持字體比例；HorzontalAlignment Center %字體中心對齊；String y: %相應(yīng)編輯框的名稱； Units normalized %采用相對度量單位，縮放時(shí)保持該區(qū)比例。（9）在按鍵右下方平板里的靜態(tài)文本的屬性編輯框中，作如下屬性設(shè)置；FontSize 12 %字體大??；FontUnits normalized %采用相對度量單位，縮放時(shí)保持字體比例；HorzontalAlignment Center %字體中心對齊；String μm %表示單位； Units normalized %采用相對度量單位，縮放時(shí)保持該區(qū)比例。圖43拖拉出控件后的打表法界面2）圖形窗口和控件參數(shù)的設(shè)置雙擊工作區(qū)或控件可引出圖形窗和相應(yīng)控件的“屬性編輯框Property Inspector”，如圖44所示。MATLAB系統(tǒng)內(nèi)部使用對象語言描述各種圖形單元，并將這些圖形單元按照樹形結(jié)構(gòu)組織起來進(jìn)行管理和實(shí)施各種操作，計(jì)算機(jī)屏幕作為該結(jié)構(gòu)的根，它的一級樹結(jié)點(diǎn)是圖形窗口對象；二級結(jié)點(diǎn)即圖形窗口的子對象可以是坐標(biāo)系對象、用戶界面控制元（有按鈕pushbutton、彈出式菜單popupmenu、文字說明text、編輯框edit等）和用戶界面菜單等；第三級樹結(jié)點(diǎn)是坐標(biāo)系對象的子對象，它們是線段對象、區(qū)域片對象、曲面對象、圖形對象、文字對象、光源對象等。所以設(shè)計(jì)用戶界面應(yīng)該遵循以下三個(gè)原則：簡單性——保持界面整潔、直觀，可有可無的功能和控件盡量刪除。)。subplot(39。mesh(m,n,l2)。,[,])。%繪制被測曲面示意圖subplot(39。z=inv(A)*D。%建立目標(biāo)函數(shù)A=zeros(3,3)。 endend%找出距離最大點(diǎn)和最小點(diǎn)的坐標(biāo)值MA=m(A)。（4）對角線法繪圖程序的編譯對角線法繪圖主程序[h,i,k,m,n,l]=sjsrpm()。l1=interp2(m,n,l,x2,y2,39。mesh(m,n,l3)。position39。subplot(39。%求Y軸最大值b=max(n)。m(k),n(k),l(k)]。MB=m(B)。%計(jì)算各測點(diǎn)到基準(zhǔn)平面的距離for k=1:N r(k)=(x(1)*m(k)+x(2)*n(k)+x(3)*l(k)+1)/sqrt(x(1)^2+x(2)^2+x(3)^2)。subplot(39。position39。%繪制上平面示意圖[MA,NA,LA,MB,NB,LB]=zxechtzcx(x)。[x,fval,exitflag]=fminsearch(39。x1=0:50:a。%求X軸最大值a=max(m)。LB=l(B)。rmin=min(r)。%調(diào)用數(shù)據(jù)輸入函數(shù)讀取數(shù)據(jù) [h,i,k,m,n,l]=sjsrpm()。cubic39。%繪制下平面示意圖l4=(MB*x(1)+x(2)*NB+x(3)*LBx(1)*mx(2)*n)/x(3)。,[,])。hold on。,[,])。 end最小區(qū)域法繪圖程序function H=zxqyht(m,n,l)。NA=n(A)。zxqyht139。%求X軸坐標(biāo)for j=1:i for k=1:h n(j,k)=(j1)*L。x(1,1)=0。),39。scz39。,39。Tag39。n=reshape(n39。%求Z軸坐標(biāo)l=c*x。%求橫向采點(diǎn)數(shù)a=xm(1,:)。%輸入返測值fc=str2num(get(findobj(39。))。string39。 水平儀測量方法的程序編譯原理 采用水平儀測量時(shí)，每一點(diǎn)的坐標(biāo)都是該點(diǎn)相對它前一點(diǎn)的讀數(shù)，所以該點(diǎn)的坐標(biāo)讀數(shù)也都是相對坐標(biāo)，它的實(shí)際坐標(biāo)與測量的橋板的跨距、分度值以測量的點(diǎn)數(shù)有關(guān)。%求被測點(diǎn)總采點(diǎn)數(shù)k=h*i。string39。),39。xsr39。cubic39。subplot(2,3,1)。[x2,y2]=meshgrid(x1,y1)。a=max(a)。0。對角線法繪圖子程序；function x=djpmht1(m,n,l)；%調(diào)用數(shù)據(jù)輸入函數(shù)讀取數(shù)據(jù)；[h,i,k,m,n,l]=sjsrpm()。z=rmaxrmin；%標(biāo)記距離最大的點(diǎn)；for a=1:N if r(a)==rmax A=a。z=inv(A)*D。%建立目標(biāo)函數(shù)求解系數(shù)矩陣；A=zeros(3,3)。N=length(m)。%繪制被測曲面光順示意圖[x2,y2]=meshgrid(x1,y1)。%繪制上平面示意圖[MA,NA,LA,MB,NB,LB]=sydpmhtzcx(x)。y1=0:100:b。n=n+1。%建立目標(biāo)函數(shù)求解系數(shù)矩陣；A=[m(1),n(1),l(1)。NA=n(A)。N=length(m)。%建立目標(biāo)函數(shù)求解系數(shù)；A=zeros(3,3)。N=length(m)。)。mesh(m,n,l3)。zxecht139。mesh(m,n,l)。 end最小二乘法繪圖程序；function H=zxecht(m,n,l)。MB=m(B)。endrmax=max(r)。 end最小二乘法繪圖主程序；function [MA,NA,LA,MB,NB,LB]=zxechtzcx(x)。 end %找出最大距離； dmax=max(d)。最小二乘法是以最小二乘中心面作為為理想要素，然后求各測點(diǎn)相對此理想要素的最大距離和最小距離，兩者的代數(shù)差即為為所測零件的平面度誤差。l1=interp2(m,n,l,x2,y2,39。mesh(m,n,l3)。position39。[x2,y2]=meshgrid(x1,y1)。position39。 end最小區(qū)域法繪圖程序；function H=zxqyht(m,n,l)。NA=n(A)。,x0)。%調(diào)用數(shù)據(jù)輸入函數(shù)讀取測得數(shù)據(jù)；[h,i,k,m,n,l]=sjsrpm()。%采用最小優(yōu)化函數(shù)解系數(shù)；[x,fval,exitflag]=fminsearch(39。（6）與其他軟件和語言有良好的對接性MATLAB與Maple、Fortran、C和Basic語言之間都可以實(shí)現(xiàn)很方便的連接，用戶只需將已有的EXE文件轉(zhuǎn)換成MEX文件即可。對于初學(xué)者只要有C語言的基礎(chǔ)，就能很容易掌握MATLAB語言。K矩陣可以通過MATLAB中矩陣除法運(yùn)算求得。平面在空間直角坐標(biāo)系中它的方程一般為：，若令=A/C，=B/C，=D/C， 則平面一般方程可寫成為: (6) 式(6)便是平面方程線性回歸的數(shù)學(xué)模型稱之為回歸函數(shù)，采用最小二乘法來估計(jì)模型式(6)中的參數(shù)，令： （7）= (8)其中（，）為實(shí)測點(diǎn)的坐標(biāo)值，則系數(shù)的最小二乘估計(jì)應(yīng)滿足以下條件，即最小二乘法的目標(biāo)函數(shù)：F()= (9)再用MATLAB優(yōu)化工具箱中合適的優(yōu)化函數(shù)對目標(biāo)函數(shù)方程（9）進(jìn)行優(yōu)化求解便可得出基準(zhǔn)平面方程的系數(shù)，至此已經(jīng)獲得了一個(gè)平面方程： (10)該方程所代表的是一個(gè)理想平面它反映了所有測點(diǎn)的平面趨勢這些測點(diǎn)有的位于該平面上有的位于其上下兩側(cè)在每一側(cè)總有一個(gè)測點(diǎn)至該平面距離為最大那么平面兩側(cè)最大的兩距離之和即為平面度誤差. 三遠(yuǎn)點(diǎn)法求平面度誤差的數(shù)學(xué)模型三遠(yuǎn)點(diǎn)法求平面度誤差的原理是按被測實(shí)際表面上相距最遠(yuǎn)的三點(diǎn)建立理想平面，以各測點(diǎn)對此平面的偏差中最大值與最小值的代數(shù)差作為被測實(shí)際表面的平面度誤差值。因此需建立平面度誤差評定的數(shù)學(xué)模型。 圖21平面度檢測示意圖最小區(qū)域法是以符合最小條件的理想平面作為平面度誤差的評定基準(zhǔn)平面的。用于平面度誤差的測量方法很多，常用的有打表法、水平儀和三坐標(biāo)機(jī)等。對角線法雖然評定結(jié)果有唯一性，但當(dāng)評定結(jié)果大于規(guī)定的平面度公差時(shí)，不能作出平面度合格與否的判斷，存在較大的局限性。在生產(chǎn)實(shí)踐中，測定平面度主要采用三坐標(biāo)機(jī)等儀器。形位公差就是研究這些要素在形狀及其相互間方向和位置方面的精度問題。形位誤差對零件的使用性能有及其重要的影響，它會(huì)影響到零件的裝配性、工作精度、配合性質(zhì)，密封性等。這一系列的的形位誤差確立、評定的過程在產(chǎn)品的設(shè)計(jì)加工過程中十分重要，它是決定產(chǎn)品是否合格的關(guān)鍵。隨著MATLAB軟件功能的日益完善，其擁有600 多個(gè)工程中要用的數(shù)學(xué)運(yùn)算函數(shù)，函數(shù)中所使用的算法都是科研和工程計(jì)算中的最新研究成果，而且經(jīng)過了各種優(yōu)化及容錯(cuò)處理，因此使用起來可靠性非常高。由于用最小區(qū)域法評定的平面度誤差接近于理想誤差值，且符合ISO標(biāo)準(zhǔn)，為此多年來有許多學(xué)者致力于這方面的研究。手工計(jì)算可以獲得相對精確的結(jié)果，但效率太低。如圖23所示。采用最小包容區(qū)域法評定時(shí)產(chǎn)生的誤差最小，精確度最高。因此便可得出基準(zhǔn)平面的方程表達(dá)式。矩陣和數(shù)組是MATLAB的核心，其提供的M語言可以與 Fortran 和C語言混合編程。（3）豐富的內(nèi)部函數(shù)MATLAB的內(nèi)部函數(shù)庫提供了相當(dāng)豐富的函數(shù)，這些函數(shù)可以方便地解決許多基本問題。針對這以上三種不同測量方法的數(shù)據(jù)特點(diǎn)可以將其分為兩類，一類是打表法和三坐標(biāo)機(jī)測量，這兩種方法所測量的數(shù)據(jù)是絕對的三維直角坐標(biāo)；另一類就是采用水平儀等工具測量的相對的三維直角坐標(biāo)。%計(jì)算各點(diǎn)到基準(zhǔn)平面的距離； for j=1:k d(j)=(x(1)*m(j)+x(2)*n(j)+x(3)*l(j)+1)/sqrt(x(1)^2+x(2)^2+x(3)^2)。%調(diào)用數(shù)據(jù)輸入函數(shù)讀取測得數(shù)據(jù)； [h,i,k,m,n,l]=sjsrpm()。rmin=min(r)。NB=n(B)。%求x軸最大值；a=max(m)。x1=0:100:a。[x,fval]=fminsearch(39。%繪制上平面示意圖；[MA,NA,LA,MB,NB,LB]=zxqyhtzcx(x)。position39。subplot(39。 x0=[0,0,0]。最小二乘法計(jì)算的子程序；function z=zxecfpmd(x)。 x0=[0,0,0]。%標(biāo)記距離最大的點(diǎn)；for a=