matlab-實(shí)驗(yàn)報(bào)告

1、數(shù)學(xué)試驗(yàn)報(bào)告一、實(shí)驗(yàn)問題：有一形狀較為復(fù)雜，但表面很光滑的曲面工件.通過科學(xué)手段，將其放置于某一空間坐標(biāo)系下，測得曲面上若干個(gè)點(diǎn)的坐標(biāo)如下： Yx-5-4-3-2-1012345-513.6-8.2-14.8-6.61.40-3.81.413.616.80-4-8.2-15.8-7.92.23.800.67.310.10-16.8-3-14.8-7.92.55.82.302.75.10-10.1-13.7-2-6.62.25.93-0.301.90-5.1-7.3-1.4-11.43.82.3-0.3-0.900-1.7-2.7-0.63.80000000000001-3.80.62.71.7

2、000.90.3-2.3-3.8-1.421.47.35.10-1.700.3-3.1-5.8-2.26.6313.610.10-5.1-2.70-2.3-5.8-2.57.914.8416.80-10.1-7.3-0.60-3.8-2.27.915.88.25016.3-13.6-1.43.80-1.46.614.88.2-13.6要求：（1）、畫出該曲面工件的圖形.（2）、在已知相鄰的橫、縱坐標(biāo)之間插入三個(gè)分點(diǎn)，用interp2命令計(jì)算出所有點(diǎn)處的豎坐標(biāo)，畫出相應(yīng)的插值曲面.（3）、分別用不同的方法求出該曲面工件表面積的近似值.二、問題分析：（1）、由x=-5:1:5;y=-5:1:5產(chǎn)生

3、“格點(diǎn)”矩陣，再用已知的各“格點(diǎn)”的縱坐標(biāo)z，用surf（x,y,z）命令即可畫出未進(jìn)行插值的工件表面圖形。（2）、要在已知相鄰的橫、縱坐標(biāo)之間插入三個(gè)分點(diǎn)，則要用xc=-5:0.25:5;yc=-5:0.25:5產(chǎn)生另一個(gè)“格點(diǎn)”矩陣，然后使用interp2命令計(jì)算所有“格點(diǎn)”處的縱坐標(biāo)z，然后再用surf（x,y,z）命令即可畫出該插值曲面。（3）、該工件是一個(gè)曲面，要求其表面積的近似值，可以想到的用“分”，“勻”，“和”，“精”的思想，將工件曲面分成n個(gè)在xoy平面上的投影是邊長為d的矩形的曲面四邊形（從第（2）題畫出的曲面可以一目了然的看出），分別求出各個(gè)曲面四邊形的面積再求和，n取的

4、越大，d越小，越接近準(zhǔn)確值。對(duì)于如何求每一個(gè)曲面四邊形的面積:方法1：如圖所示，曲面面積微元S近似等于以O(shè)P1和OP2為鄰邊的平行四邊形的面積S=OP1×OP2，只要n取的值較大，d足夠小，誤差會(huì)較小。曲面微元：方法2：由第一型面積分公式：S=（）1+Fx2+Fy2dxdy ，曲面面積微元S等于1+Fx2+Fy2dxdy. 由上述分割工件曲面的思想，dxdy即等于d2. 問題是如何求得曲面方程的對(duì)x、y的偏導(dǎo)數(shù)？通過查找資料，發(fā)現(xiàn)Matlab軟件中有求梯度的命令gradient，因此，曲面的面積微元S可以求得，而且與方法一類似，當(dāng)d取足夠小時(shí)，誤差可以忽略不計(jì)。三、問題求解：（1）、

5、使用Matlab軟件編寫如下程序：x=-5:1:5;y=-5:1:5;xb,yb=meshgrid(x,y);zb=13.6,-8.2,-14.8,-6.6,1.4,0,-3.8,1.4,13.6,16.8,0; -8.2,-15.8,-7.9,2.2,3.8,0,0.6,7.3,10.1,0,-16.8; -14.8,-7.9,2.5,5.8,2.3,0,2.7,5.1,0,-10.1,-13.7; -6.6,2.2,5.9,3.0,-0.3,0,1.9,0,-5.1,-7.3,-1.4; 1.4,3.8,2.3,-0.3,-0.9,0,0,-1.7,-2.7,-0.6,3.8; 0,0,0

6、,0,0,0,0,0,0,0,0; -3.8,0.6,2.7,1.7,0,0,0.9,0.3,-2.3,-3.8,-1.4; 1.4,7.3,5.1,0,-1.7,0,0.3,-3.1,-5.8,-2.2,6.6; 13.6,10.1,0,-5.1,-2.7,0,-2.3,-5.8,-2.5,7.9,14.8; 16.8,0,-10.1,-7.3,-0.6,0,-3.8,-2.2,7.9,15.8,8.2; 0,16.3,-13.6,-1.4,3.8,0,-1.4,6.6,14.8,8.2,-13.6;surf(xb,yb,zb)圖1、由已知坐標(biāo)點(diǎn)畫出的工件表面圖形（2）、使用Matlab軟件

7、編寫如下程序：x=-5:1:5;y=-5:1:5;xb,yb=meshgrid(x,y);zb=13.6,-8.2,-14.8,-6.6,1.4,0,-3.8,1.4,13.6,16.8,0; -8.2,-15.8,-7.9,2.2,3.8,0,0.6,7.3,10.1,0,-16.8; -14.8,-7.9,2.5,5.8,2.3,0,2.7,5.1,0,-10.1,-13.7; -6.6,2.2,5.9,3.0,-0.3,0,1.9,0,-5.1,-7.3,-1.4; 1.4,3.8,2.3,-0.3,-0.9,0,0,-1.7,-2.7,-0.6,3.8; 0,0,0,0,0,0,0,0

8、,0,0,0; -3.8,0.6,2.7,1.7,0,0,0.9,0.3,-2.3,-3.8,-1.4; 1.4,7.3,5.1,0,-1.7,0,0.3,-3.1,-5.8,-2.2,6.6; 13.6,10.1,0,-5.1,-2.7,0,-2.3,-5.8,-2.5,7.9,14.8; 16.8,0,-10.1,-7.3,-0.6,0,-3.8,-2.2,7.9,15.8,8.2; 0,16.3,-13.6,-1.4,3.8,0,-1.4,6.6,14.8,8.2,-13.6;xc=-5:0.25:5;yc=-5:0.25:5;xcb,ycb=meshgrid(xc,yc);zcb=in

9、terp2(xb,yb,zb,xcb,ycb,'spline')surf(xcb,ycb,zcb)圖2、插值計(jì)算出的所有點(diǎn)的縱坐標(biāo)值（截取了一部分）圖3、插值曲面（3）、求工件表面面積：利用方法一求面積微元：x=-5:1:5;y=-5:1:5;xb,yb=meshgrid(x,y);zb=13.6,-8.2,-14.8,-6.6,1.4,0,-3.8,1.4,13.6,16.8,0; -8.2,-15.8,-7.9,2.2,3.8,0,0.6,7.3,10.1,0,-16.8; -14.8,-7.9,2.5,5.8,2.3,0,2.7,5.1,0,-10.1,-13.7; -6

10、.6,2.2,5.9,3.0,-0.3,0,1.9,0,-5.1,-7.3,-1.4; 1.4,3.8,2.3,-0.3,-0.9,0,0,-1.7,-2.7,-0.6,3.8; 0,0,0,0,0,0,0,0,0,0,0; -3.8,0.6,2.7,1.7,0,0,0.9,0.3,-2.3,-3.8,-1.4; 1.4,7.3,5.1,0,-1.7,0,0.3,-3.1,-5.8,-2.2,6.6; 13.6,10.1,0,-5.1,-2.7,0,-2.3,-5.8,-2.5,7.9,14.8; 16.8,0,-10.1,-7.3,-0.6,0,-3.8,-2.2,7.9,15.8,8.2;

11、 0,16.3,-13.6,-1.4,3.8,0,-1.4,6.6,14.8,8.2,-13.6;d=0.25； %改變d的值即可改變曲面四邊形的個(gè)數(shù)nxc=-5:d:5;yc=-5:d:5;xcb,ycb=meshgrid(xc,yc);zcb=interp2(xb,yb,zb,xcb,ycb,'spline');p=1;q=1;ss=0;while p<10/d+1 while q<10/d+1 s=sqrt(d*d*(zcb(p,q+1)-zcb(p,q)2+d*d*(zcb(p+1,q)-zcb(p,q)2+d4); %s為每一個(gè)曲面四邊形的面積 ss=s+

12、ss; %對(duì)s循環(huán)相加 q=q+1; end q=1; p=p+1;endss（步長d取0.25時(shí)得工件表面積約為ss =670.0184）使用方法二求面積微元：x=-5:1:5;y=-5:1:5;xb,yb=meshgrid(x,y);zb=13.6,-8.2,-14.8,-6.6,1.4,0,-3.8,1.4,13.6,16.8,0; -8.2,-15.8,-7.9,2.2,3.8,0,0.6,7.3,10.1,0,-16.8; -14.8,-7.9,2.5,5.8,2.3,0,2.7,5.1,0,-10.1,-13.7; -6.6,2.2,5.9,3.0,-0.3,0,1.9,0,-5.

13、1,-7.3,-1.4; 1.4,3.8,2.3,-0.3,-0.9,0,0,-1.7,-2.7,-0.6,3.8; 0,0,0,0,0,0,0,0,0,0,0; -3.8,0.6,2.7,1.7,0,0,0.9,0.3,-2.3,-3.8,-1.4; 1.4,7.3,5.1,0,-1.7,0,0.3,-3.1,-5.8,-2.2,6.6; 13.6,10.1,0,-5.1,-2.7,0,-2.3,-5.8,-2.5,7.9,14.8; 16.8,0,-10.1,-7.3,-0.6,0,-3.8,-2.2,7.9,15.8,8.2; 0,16.3,-13.6,-1.4,3.8,0,-1.4,6

14、.6,14.8,8.2,-13.6;d=0.25; %改變d的值即可改變曲面四邊形的個(gè)數(shù)nxc=-5:d:5;yc=-5:d:5;xcb,ycb=meshgrid(xc,yc);zcb=interp2(xb,yb,zb,xcb,ycb,'spline');Fx,Fy=gradient(zcb,d,d);p=1;q=1;ss=0;while q<10/d+1 while p<10/d+1 s=sqrt(1+Fx(p,q)2+Fy(p,q)2)*d*d; %s為每一個(gè)曲面四邊形的面積 ss=s+ss; %對(duì)s循環(huán)相加 p=p+1; end p=1; q=q+1;endss（步長d取0.25時(shí)得工件表面積約為ss = 661.1878）兩種方法的比較：d的取值方法一的ss方法二的ss0.250 670.0184661.1878