機械精度實驗

1、Calculating flatness with least square method1. Matlab code:clcclearA=xlsread(flatness,1; for i=1:size(A,1 for j=1:size(A,2 if isnan(A(i,j A(i,j=0;endendend x,y,z=find(A; t=x,y; b=regress(z,t; X,Y=meshgrid(1:size(A,1,1:size(A,2; mesh(X,Y,b(1*X+b(2*Y; hold on plot3(x,y,z,k.;maxval,maxind=max(z-(b(1*x

2、+b(2*y; e=abs(maxval+abs(minval; 2.Graph of regressed plant and related points: 3 The numerical value of flatness (The conclusion calculate by matlab, the flatness is 11.5120 m The conclusion with least square method matches the result in jpg (The conclusion calculate with software, which is also 11

3、.5120 m Calculating flatness with Maximum straightness method clear clc A=xlsread(flatness,1;a=;b=;for i=1:size(A,1flag1=1;flag2=1;for j=2:size(A,2-1ifany(isnan(A(i,j=0 & any(isnan(A(i,j-1a=a,j; flag1=0; end if any(isnan(A(i,j=0 & any(isnan(A(i,j+1b=b,j; flag2=0; end end if flag1a=a 1; endif flag2b=

4、b,size(A,2;end endk=;f=;j=1;for i=1:size(A,1 k=k,(A(i,b(j-A(i,a(j/(b(j-a(j; maxval,maxind=max(A(i,a(j:b(j-myfun(a(j:b(j,k(j,b(j,A(i,b(j; minval,minind=min(A(i,a(j:b(j-myfun(a(j:b(j,k(j,b(j,A(i,b(j; f=f,abs(maxval+abs(minval; j=j+1;endf_final=max(f; The conclusion with this method is: (The mallab giv

5、es a result which is 7.215 umCalculating straightness with least square method1. Codeclcclear A=xlsread( e=; for i=1:2:3 a=A(i,:; b=A(i+1,:; d=regress(b,a;figure plot(a,c,r,a,b,b,100,k. .;enda=A(10,1:16;b=A(11,1:16; d=regress(b,a;maxval,maxind=max(b-c; minval,minind=min(b-c; figure plot(a,c,r,a,b,b,

6、100,k. . ; e=e abs(maxval+abs(minval;2. The numerical value of straightnessAccording to Matlab:e =0.0337 0.0312 0.0375 for three lines respectively.3.Theplot of three lines with related zone and curve:(The plot of first line (The plot of second line (The plot of third line II AB Straightness=AB+0=0.

7、0312um II AB Straightness=AB+0=0.0375um IABI BCStraightness=AB+BC=0.0337umCalculate straightness with two end points method1. Codeclearclc A=xlsread(straightness.xls,1; f=;for i=1:2:3 a=A(i,:; b=A(i+1,:; k=(b(1-b(length(b/(a(1-a(length(b; c=k*(a-a(1*ones(1,length(b+b(1*ones(1,length(b; maxval,maxind

8、=max(b-c;minval,minind=min(b-c;figureplot(a,c,r,a,b,b,a(maxind*ones(1,100,linspace(c(maxind,b(maxind,100,k. ,a(minind*ones(1,100,linspace(b(minind,c(minind,100,c. ,a,(k*(a-a(maxind*ones(1,length(b+b(maxind*ones(1,length(b,m. ,a,(k*(a-a(minind*ones(1,length(b+b(minind*ones(1,length(b,y;f=f abs(maxval

9、+abs(minval; enda=A(10,1:16;b=A(11,1:16;k=(b(1-b(length(b/(a(1-a(length(b;c=k*(a-a(1*ones(1,length(b+b(1*ones(1,length(b;maxval,maxind=max(b-c;minval,minind=min(b-c; figure plot(a,c,r,a,b,b,a(maxind*ones(1,100,linspace(c(maxind,b(maxind,100,k. ,a(minind*ones(1,100,linspace(b(minind,c(minind,100,c.,a

10、,(k*(a-a(maxind*ones(1,length(b+b(maxind*ones(1,length(b,m.,a,(k*(a-a(minind*ones(1,length(b+b(minind*ones(1,length(b,y;f=f abs(maxval+abs(minval;2.The numerical value of straightnessAccording to matlab: f =0.0226 0.0287 0.0333 for three lines separately.3.Theplot of three lines with related zone an

11、d curve:Import data from excel, and save it as A Calculate the inclines with two and points. And define lines with then. Finally calculate the vertical intercepts. Plot graph containing all required information. Save the straightness in F Save the straightness in F These process is a repeating of pr

12、eceding steps. ( No need to payattention toIIIIII ABStraightness=AB+0=0.0226um ABABStraightness=AB+0=0.0287umStraightness=AB+0=0.0333umCalculating cylindricity with regressive circle method clc,clear;A=xlsread(cgroup4_1.CSV;B=xlsread(cgroup4_2.CSV;C=xlsread(cgroup4_3.CSV;D=xlsread(cgroup4_4.CSV; n=l

13、ength(X(:,1;y=ones(n,1;b=rand(1*1000 rand(1 rand(1;fun=inline(X(:,1.2+X(:,2.2+a(1*X(:,1+a(2*X(:,2+a(3,a,X;a,r,J = nlinfit(X,y,fun,b ;vpa(a,10 ;X1 =-a(1/2;Y1 =-a(2/2;R =sqrt(a(12+a(22-4*a(3/2;X1=ones(n,1*X1;Y1=ones(n,1*Y1;maxval=max(sqrt(X(:,1-X1.2+(X(:,2-Y1.2;minval=min(sqrt(X(:,1-X1.2+(X(:,2-Y1.2;R

14、1=maxval-minval;The graph of cylinder:Calculating roughness with least squareRead the x,y coordinatefrom excels for differentheight and save then onA,B,C,D respectivelyProject all data to oneplane and regress all datato a specific circle.Find the Radius, centralcoordinates of regressed circle.Find t

15、he difference ofdistance between largestand smallest points fromcentral pointmethod1.Programing code clc;clear; A=xlsread(roughness.xlsx;z=find(isnan(A(:,3; for i=length(z:-1:1 A(z(i,:=;end Z=A(:,3; t=A(:,1,A(:,2,ones(size(A,1,1;b=regress(Z,t; f=abs(b(1*A(:,1+b(2*A(:,2-Z+b(3*ones(size(A,1,1/sqrt(b(1

16、2+b(22+(-12;Ra=mean(f;plot3(A(:,1,A(:,2,Z,k.;hold on X,Y=meshgrid(0:700,0:700; mesh(X,Y,b(1*X+b(2*Y+b(3*ones(size(X,1,size(X,2The conclusion:The Ra found by this method is:Ra= 0.1634 um, which matches the conclusion calculated by software: 3.The graph of regressed plane and measured points: Import the data from