matlab經(jīng)典圖像圖形程序

1、（平面擬合 m= a=m(:,1:3);a1=m(:,4:6);a2=m(:,7:9);a3=m(:,10:12);hold on;c=1,8;9,16;16,24;24,32;32,40;40,48;48,56;56,64;64,72;72,80;80,88;88,96;96,104;b=zeros(3,13); for k=1:13x=a(c(k,1):c(k,2),1);y=a(c(k,1):c(k,2),2);z=a(c(k,1):c(k,2),3);Xcolv=x(:);Ycolv=y(:);Zcolv=z(:);Const=ones(size(Xcolv);Coefficients

2、1(:,k) = Xcolv Ycolv ConstZcolv;XCoeff=Coefficients1(1);YCoeff=Coefficients1(2);CCoeff=Coefficients1(3);L=plot3(x,y,z,ro); %繪制三維圖形set(L,Markersize,2*get(L,Markersize) set(L,Markerfacecolor,r)hold on; grid on;xx, yy=meshgrid(0:0.1:1,0:0.1:1);zz = XCoeff * xx + YCoeff * yy + CCoeff; surf(xx,yy,zz)titl

3、e(sprintf(Plotting plane z=(%f)*x+(%f)*y+(%f),XCoeff, YCoeff, CCoeff)endCoefficients1%由程序計算所得的系數(shù)for k=1:13 x=a1(c(k,1):c(k,2),1);y=a1(c(k,1):c(k,2),2);z=a1(c(k,1):c(k,2),3);Xcolv=x(:);Ycolv=y(:);Zcolv=z(:);Const=ones(size(Xcolv);Coefficients2(:,k) = Xcolv Ycolv ConstZcolv;XCoeff=Coefficients2(1);YCo

4、eff=Coefficients2(2);CCoeff=Coefficients2(3);L=plot3(x,y,z,ro); %繪制三維圖形set(L,Markersize,2*get(L,Markersize) set(L,Markerfacecolor,r)hold on; grid on;xx, yy=meshgrid(0:0.1:1,0:0.1:1);zz = XCoeff * xx + YCoeff * yy + CCoeff; surf(xx,yy,zz)title(sprintf(Plotting plane z=(%f)*x+(%f)*y+(%f),XCoeff, YCoef

5、f, CCoeff)endCoefficients2%由程序計算所得的系數(shù)for k=1:13 x=a2(c(k,1):c(k,2),1);y=a2(c(k,1):c(k,2),2);z=a2(c(k,1):c(k,2),3);Xcolv = x(:);Ycolv = y(:);Zcolv = z(:);Const = ones(size(Xcolv);Coefficients3(:,k) = Xcolv Ycolv ConstZcolv;XCoeff=Coefficients3(1);YCoeff=Coefficients3(2);CCoeff=Coefficients3(3);L=plot

6、3(x,y,z,ro); %繪制三維圖形set(L,Markersize,2*get(L,Markersize) set(L,Markerfacecolor,r)hold on; grid on;xx, yy=meshgrid(0:0.1:1,0:0.1:1);zz = XCoeff * xx + YCoeff * yy + CCoeff; surf(xx,yy,zz)title(sprintf(Plotting plane z=(%f)*x+(%f)*y+(%f),XCoeff, YCoeff, CCoeff)endCoefficients3%由程序計算所得的系數(shù)m=a=m(:,1:3);%

7、1989觀測數(shù)據(jù)觀測數(shù)據(jù)for k=1:13觀測數(shù)據(jù)x=a3(c(k,1):c(k,2),1);y=a3(c(k,1):c(k,2),2);z=a3(c(k,1):c(k,2),3);Xcolv = x(:);Ycolv = y(:);Zcolv = z(:);Const = ones(size(Xcolv);Coefficients4(:,k) = Xcolv Ycolv XCoeff = Coefficients4(1);YCoeff = Coefficients4(2); CCoeff =Coefficients4(3);L=plot3(x,y,z,ro); %繪制三維圖形set(L,M

8、arkersize,2*get(L,Markersize) set(L,Markerfacecolor,r)hold on; grid on;xx, yy=meshgrid(0:0.1:1,0:0.1:1);zz = XCoeff * xx + YCoeff * yy + CCoeff; surf(xx,yy,zz)title(sprintf(Plotting plane z=(%f)*x+(%f)*y+(%f),XCoeff, YCoeff, CCoeff)endCoefficients4%由程序計算所得的系數(shù)程序二：（中心點坐標的確定）首先應用matlab 建立ff2.m 文件functi

9、onf=ff2(x,A,B) x1=A(:,1);y=A(:,2);z=A(:,3);f=sum(x1-x(1).2+(y-x(2).2+(z-(B(1)*x(1)+B(2)*x(2)+B(3).2);然后，再運行如下程序：m=觀測數(shù)據(jù)觀測數(shù)據(jù)觀測數(shù)據(jù)F4=m(:,10:12);%2011P1=-0.0008309780.003417388 0.471956284-0.0008177780.003628518 5.887189422-0.103531579-0.170598336160.1073406-0.084434924-0.139125499137.2555526-0.093172538

10、-0.156337182155.8157108-0.10905983 -0.155763638169.0503192-0.097249562-0.132905095154.0980362-0.100944809-0.138739994162.7653381-0.107075519-0.14903722 175.1283597-0.107132219-0.156237569182.2592083-0.135031497-0.217484619234.2576489-0.145292163-0.233311605252.6149174-0.147842995-0.253950265268.9965

11、036;P2=-0.0014879980.003714862 0.684399642-0.0004928270.003781979 5.617202869-0.103877083-0.170596425160.2976036-0.082898692-0.139164862136.4014982-0.093547416-0.156341154156.0256571-0.108695055-0.155620482168.7628486-0.096893482-0.132785851153.8288996-0.101347724-0.138656209162.9445585-0.106635181-

12、0.148858581174.7784003-0.107683182-0.156283925182.5910487-0.135533438-0.21745348 234.5201978-0.14590835 -0.233379092252.9945534-0.148170132-0.254010047269.2086747;P3=-0.002373299-0.0037697875.079985167-0.003971777-0.0001935489.661221920.171779992-0.104986237-30.246139390.140823091-0.082759129-19.892

13、823950.157627698-0.09089535 -20.558345520.168465935-0.118761803-7.6328559860.134087929-0.0952972043.2699959360.139156581-0.1007542596.7667143040.147528599-0.1048175387.6364023680.15272301-0.1034554557.3010547170.215152629-0.135879385-6.9951433480.232413579-0.147310618-6.5611672740.2413723-0.15772857

14、9-2.08825351;P4=-0.002333174-0.0037402875.040576042-0.002554701-0.00253047510.061014040.172474839-0.106260378-29.982248370.14091257-0.085030797-18.772167250.158721438-0.09191605 -20.653404770.168403679-0.119171259-7.3901692050.133926192-0.0949300623.1619691530.139361927-0.1007942716.6674070160.14823

15、9426-0.1072088578.4637523070.155180017-0.1071291347.8085213890.216632193-0.138845389-6.2951075630.231243928-0.142257642-8.5539260470.241592539-0.158407287-1.864415599;c=1,8;9,16;17,24;25,32;33,40;41,48;49,56;57,64;65,72;73,80;81,88;89,96;97,104;b1=zeros(3,13); for k=1:13A = a=P1(k,:);b1(1:2,k) = fmi

16、nsearch(x)ff2(x,A,a),1,2);b1(3,k)=b1(1,k)*P1(k,1)+b1(2,k)*P1(k,2)+P1(k,3);end b2=zeros(3,13); for k=1:13A = a=P2(k,:);b2(1:2,k) = fminsearch(x)ff2(x,A,a),1,2);b2(3,k)=b2(1,k)*P2(k,1)+b1(2,k)*P2(k,2)+P2(k,3);end b3=zeros(3,13); for k=1:13A = a=P3(k,:);b3(1:2,k) = fminsearch(x)ff2(x,A,a),1,2);b3(3,k)=

17、b3(1,k)*P3(k,1)+b3(2,k)*P3(k,2)+P3(k,3);end b4=zeros(3,13); for k=1:13A = a=P4(k,:);b4(1:2,k) = fminsearch(x)ff2(x,A,a),1,2);b4(3,k)=b4(1,k)*P4(k,1)+b4(2,k)*P4(k,2)+P4(k,3);endb1 ,b2,b3,b4（傾斜角的計算 clearclca=566.664741567.2640516522.7105282522.25409941.78737503755.10854517566.6649707567.2543522.710180

18、5522.23661.78300077755.11965447566.7268041567.336522.7014735522.21481.7644997955.091566.7269181567.3375522.7013548522.21351.76325044955.087;%各次底層中心點及塔尖中心點的坐標值b=zeros(4,1);for i=1:4b(i,1)=atan(sqrt(a(3*i-2,1)-a(3*i-2,2)2+(a(3*i-1,1)-a(3*i-1,2)2)/(a(3*i,2)-a(3*i,1); endb（空間曲線方程的擬合 clearclca1=566.66474

19、1 566.7196032 566.7251066 566.7841861 566.824563 566.8660307566.9166678 566.9538279 566.9897481 567.0267344 567.0564911 567.1007311567.148014 567.2640516522.7105282 522.668344 522.5474854 522.5418216 522.4961758 522.4639359522.4669521 522.4505947 522.4317649 522.4184393 522.3456014 522.3016699522.26

20、14503 522.25409941.787375037 7.320250093 12.28766395 16.70028124 21.31764794 25.8471228729.52711866 33.04949985 36.55586028 39.89098906 44.08502791 48.3605846452.5192089 55.10854517; x=a1(1,:);y=a1(2,:);z=a1(3,:);z 變量，za2=566.6649707566.7205451 566.7265463566.7869552566.827267566.8696053566.92055035

21、66.9579355 566.9946391567.0317229567.0620474567.1064981567.1540086567.2543522.7101805 522.6674562 522.5459168 522.5394603 522.4932815 522.4605572522.4630366 522.4463037 522.4268193 522.4132488 522.3398815 522.2955148522.2552224 522.23661.783000777 7.314628384 12.2829702 16.69614048 21.31277587 25.84

22、08361329.52177461 33.04364328 36.54837315 39.88566431 44.07846016 48.3546905252.51373486 55.11965447; x=a2(1,:);y=a2(2,:);z=a2(3,:);z 變量，za3=566.7268041566.763964 566.8798373566.8829498566.9237686567.0100817567.0214421567.0722482 567.1256312567.1797278567.2555827567.3031922567.3514709567.336522.7014

23、735 522.6692942 522.5896155 522.5816831 522.5500444 522.4897794 522.482311522.449441 522.4153258 522.3647064 522.3068783 522.2648641 522.2190398522.21481.76449979 7.309000247 12.2677573 16.68898035 21.30717388 25.8371993429.50962342 33.03954325 36.54536393 39.88097158 44.08064918 48.3526381152.48580

24、886 55.091; x=a3(1,:);y=a3(2,:);z=a3(3,:);z 變量，za4=566.7269181566.7641741 566.8809476566.8830262566.9252064567.0107413567.0222307567.0732454 567.1263284567.1816335567.257516567.3043501567.3528621567.3375522.7013548 522.6689947 522.5891489 522.5805209 522.5487839 522.4889408522.4816388 522.44856522.4

25、139989 522.3625187 522.3044631 522.2651696522.2174357 522.21351.763250449 7.290500213 12.2599313 16.67333868 21.29915907 25.8308607629.50188293 33.03600551 36.52682599 39.86353266 44.07156577 48.3355487252.48075565 55.087; x=a4(1,:);y=a4(2,:);z=a4(3,:);z 變量，z（曲率K 值的計算 clearclcsyms x1 y1 x2 y2 x3 y3

26、x4 y4 z;x1=-0.00007*z2+0.0064*z+566.65624; y1=-0.00001*z2-0.00863*z+522.71088; x2=-0.00006*z2+0.0069*z+566.65375; y2=0.00001*z2+0.00856*z+522.70036; x3=0.00002*z2+0.01078*z+566.70374; y3=-0.00004*z2-0.00682*z+522.70905; x4=0.00002*z2-0.01083*z+566.70353; y4=-0.00002*z2+0.01083*z+566.70353;dx1=diff(x

27、1,z,1) ddx1=diff(x1,z,2) dy1=diff(y1,z,1) ddy1=diff(y1,z,2) dx2=diff(x2,z,1) ddx2=diff(x2,z,2) dy2=diff(y2,z,1) ddy2=diff(y2,z,2) dx3=diff(x3,z,1) ddx3=diff(x3,z,2) dy3=diff(y3,z,1) ddy3=diff(y3,z,2) dx4=diff(x4,z,1) ddx4=diff(x4,z,2) dy4=diff(y4,z,1) ddy4=diff(y4,z,2)n=1.787375037 7.320250093 12.28

28、766395 16.70028124 21.31764794 25.8471228729.52711866 33.04949985 36.55586028 39.89098906 44.08502791 48.3605846452.5192089 55.108545171.783000777 7.314628384 12.2829702 16.69614048 21.31277587 25.8408361329.52177461 33.04364328 36.54837315 39.88566431 44.07846016 48.3546905252.51373486 55.119654471

29、.76449979 7.309000247 12.2677573 16.68898035 21.30717388 25.8371993429.50962342 33.03954325 36.54536393 39.88097158 44.08064918 48.3526381152.48580886 55.0911.763250449 7.290500213 12.2599313 16.67333868 21.29915907 25.8308607629.50188293 33.03600551 36.52682599 39.86353266 44.07156577 48.3355487252

30、.48075565 55.087;% 中心點的z 軸數(shù)據(jù)for k=1:length(n)a(k,1)=sqrt-(7*n(1,00)*(50000)-(-8100000-n(1,000)*(50000)2+(0000)2+(50000)2)/sqrt5-(7*n(1,k000)2+(-8100000-n(1,000)2+1)3);b(k,1)=sqrt(69/10000-(3*n*50000)-(112500-n(2,000)*(25000)2+(500 0)2+50000)2)/sqrt(69/10000-(3*n2+(112500-n(2,000)2+1)3);c(k,1)=sqrt(n

31、53950000)*(12500)-(-0000-n(3,k)/1250000)2+12500)2+(25000)2)/sqrt(n(0+53950000)2+(-350000-n(3,k)/12500)2+1)3);d(k,1)=sqrt(n-108/100000)*(5000)-(10100000-n(4,000)*(25000)2+(5 000)2+(25000)2)/sqrt(n(0-108/100000)2+(10100000-n(4,000)2+1)3);enda,b,c,d各% 次測量的的曲率值（扭轉(zhuǎn)量的計算 P=566.664741522.7105282 566.6649707

32、 522.7101805566.7268041522.7014735566.7269181522.7013548566.7196032522.668344566.7205451522.6674562566.763964522.6692942566.7641741522.6689947566.7251066522.5474854566.7265463522.5459168566.8798373522.5896155566.8809476522.5891489566.7841861522.5418216566.7869552522.5394603566.8829498522.5816831566.

33、8830262522.5805209566.824563522.4961758566.827267522.4932815566.9237686522.5500444566.9252064522.5487839566.8660307522.4639359566.8696053522.4605572567.0100817522.4897794567.0107413522.4889408566.9166678522.4669521566.9205503522.4630366567.0214421522.482311567.0222307522.4816388566.9538279522.450594

34、7566.9579355522.4463037567.0722482522.449441567.0732454522.44856566.9897481522.4317649566.9946391522.4268193567.1256312522.4153258567.1263284522.4139989567.0267344522.4184393567.0317229522.4132488567.1797278522.3647064567.1816335522.3625187567.0564911522.3456014567.0620474522.3398815567.2555827522.3

35、068783567.257516522.3044631567.1007311522.3016699567.1064981522.2955148567.3031922522.2648641567.3043501522.2651696567.148014522.2614503567.1540086522.2552224567.3514709522.2190398567.3528621522.2174357;A1=F(:,1:4);A2=F(:,3:6);A3=F(:,5:8);B1(:,1)=P(:,3)-P(:,1);B1(:,2)=P(:,4)-P(:,2);B2(:,1)=P(:,5)-P(

36、:,3);B2(:,2)=P(:,6)-P(:,4);B3(:,1)=P(:,7)-P(:,5);B3(:,2)=P(:,8)-P(:,6);for m=1:8x(m,1)=fminbnd(x)niuqv1(x,A1(m,:),B1(1,:),-0.005,0.005); endfor m=9:16x(m,2)=fminbnd(x)niuqv1(x,A1(m,:),B1(2,:),-0.005,0.005); endfor m=17:24x(m,3)=fminbnd(x)niuqv1(x,A1(m,:),B1(3,:),-0.005,0.005); endfor m=25:32x(m,4)=f

37、minbnd(x)niuqv1(x,A1(m,:),B1(4,:),-0.005,0.005); endfor m=33:40 x(m,5)=fminbnd(x)niuqv1(x,A1(m,:),B1(5,:),-0.005,0.005); endfor m=41:48x(m,6)=fminbnd(x)niuqv1(x,A1(m,:),B1(6,:),-0.005,0.005); endfor m=49:56x(m,7)=fminbnd(x)niuqv1(x,A1(m,:),B1(7,:),-0.005,0.005); endfor m=57:64x(m,8)=fminbnd(x)niuqv1

38、(x,A1(m,:),B1(8,:),-0.005,0.005); endfor m=65:72x(m,9)=fminbnd(x)niuqv1(x,A1(m,:),B1(9,:),-0.005,0.005); endfor m=73:80 x(m,10)=fminbnd(x)niuqv1(x,A1(m,:),B1(10,:),-0.005,0.005);endfor m=81:88x(m,11)=fminbnd(x)niuqv1(x,A1(m,:),B1(11,:),-0.005,0.005);endfor m=89:96x(m,12)=fminbnd(x)niuqv1(x,A1(m,:),B

39、1(12,:),-0.005,0.005);endfor m=97:104x(m,13)=fminbnd(x)niuqv1(x,A1(m,:),B1(13,:),-0.005,0.005);end a1(:,1)=asin(sum(x)/8); for m=1:8x(m,1)=fminbnd(x)niuqv1(x,A2(m,:),B2(1,:),-0.005,0.005); endfor m=9:16x(m,2)=fminbnd(x)niuqv1(x,A2(m,:),B2(2,:),-0.005,0.005); endfor m=17:24x(m,3)=fminbnd(x)niuqv1(x,A

40、2(m,:),B2(3,:),-0.005,0.005); endfor m=25:32x(m,4)=fminbnd(x)niuqv1(x,A2(m,:),B2(4,:),-0.005,0.005); endfor m=33:40 x(m,5)=fminbnd(x)niuqv1(x,A2(m,:),B2(5,:),-0.005,0.005); endfor m=41:48x(m,6)=fminbnd(x)niuqv1(x,A2(m,:),B2(6,:),-0.005,0.005); endfor m=49:56x(m,7)=fminbnd(x)niuqv1(x,A2(m,:),B2(7,:),-0.005,0.005); endfor m=57:64x(m,8)=fminbnd(x)niuqv1(x,A2(m,:),B2(8,:),-0.005,0.005); endfor m=65:72x(m,9)=fminbnd(x)niuqv1(x,A2(m,:),B2(9,:),-0.005,0.005); endfor m=73:80 x(m,10)=fminbnd(x)niuqv1(x,A2(m,:),B2(10,:),-0.005,0.005);endfor m=81:88x(m,11)=fminbnd(x)niuqv1(x,A2(m,:),B2