文档介绍：该【matlab-实验报告 】是由【幸福人生】上传分享，文档一共【9】页，该文档可以免费在线阅读，需要了解更多关于【matlab-实验报告 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。matlab-实验报告数学试验报告一、实验问题:有一形状较为复杂,,将其放置于某一空间坐标系下,测得曲面上若干个点的坐标如下:Yx -5-4-3-2-1012345------4-----3-----2-----------:方法1:如图所示,曲面面积微元S近似等于以OP1和OP2为邻边的平行四边形的面积S=∥OP1×OP2∥,只要n取的值较大,d足够小,误差会较小。曲面微元:方法2:由第一型面积分公式:S=(σ)1+Fx2+Fy2dxdy,曲面面积微元S等于1+Fx2+,、y的偏导数?通过查找资料,发现Matlab软件中有求梯度的命令gradient,因此,曲面的面积微元S可以求得,而且与方法一类似,当d取足够小时,误差可以忽略不计。三、问题求解:(1)、使用Matlab软件编写如下程序:x=-5:1:5;y=-5:1:5;[xb,yb]=meshgrid(x,y);zb=[,-,-,-,,0,-,,,,0;-,-,-,,,0,,,,0,-;-,-,,,,0,,,0,-,-;-,,,,-,0,,0,-,-,-;,,,-,-,0,0,-,-,-,;0,0,0,0,0,0,0,0,0,0,0;-,,,,0,0,,,-,-,-;,,,0,-,0,,-,-,-,;,,0,-,-,0,-,-,-,,;,0,-,-,-,0,-,-,,,;0,,-,-,,0,-,,,,-];surf(xb,yb,zb)图1、由已知坐标点画出的工件表面图形(2)、使用Matlab软件编写如下程序:x=-5:1:5;y=-5:1:5;[xb,yb]=meshgrid(x,y);zb=[,-,-,-,,0,-,,,,0;-,-,-,,,0,,,,0,-;-,-,,,,0,,,0,-,-;-,,,,-,0,,0,-,-,-;,,,-,-,0,0,-,-,-,;0,0,0,0,0,0,0,0,0,0,0;-,,,,0,0,,,-,-,-;,,,0,-,0,,-,-,-,;,,0,-,-,0,-,-,-,,;,0,-,-,-,0,-,-,,,;0,,-,-,,0,-,,,,-];xc=-5::5;yc=-5::5;[xcb,ycb]=meshgrid(xc,yc);zcb=interp2(xb,yb,zb,xcb,ycb,'spline')surf(xcb,ycb,zcb)图2、插值计算出的所有点的纵坐标值(截取了一部分)图3、插值曲面(3)、求工件表面面积:①利用方法一求面积微元:x=-5:1:5;y=-5:1:5;[xb,yb]=meshgrid(x,y);zb=[,-,-,-,,0,-,,,,0;-,-,-,,,0,,,,0,-;-,-,,,,0,,,0,-,-;-,,,,-,0,,0,-,-,-;,,,-,-,0,0,-,-,-,;0,0,0,0,0,0,0,0,0,0,0;-,,,,0,0,,,-,-,-;,,,0,-,0,,-,-,-,;,,0,-,-,0,-,-,-,,;,0,-,-,-,0,-,-,,,;0,,-,-,,0,-,,,,-];d=;%改变d的值即可改变曲面四边形的个数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;whilep<10/d+1whileq<10/d+1s=sqrt(d*d*(zcb(p,q+1)-zcb(p,q))^2+d*d*(zcb(p+1,q)-zcb(p,q))^2+d^4);%s为每一个曲面四边形的面积ss=s+ss;%对s循环相加q=q+1;endq=1;p=p+1;endss(=)②使用方法二求面积微元:x=-5:1:5;y=-5:1:5;[xb,yb]=meshgrid(x,y);zb=[,-,-,-,,0,-,,,,0;-,-,-,,,0,,,,0,-;-,-,,,,0,,,0,-,-;-,,,,-,0,,0,-,-,-;,,,-,-,0,0,-,-,-,;0,0,0,0,0,0,0,0,0,0,0;-,,,,0,0,,,-,-,-;,,,0,-,0,,-,-,-,;,,0,-,-,0,-,-,-,,;,0,-,-,-,0,-,-,,,;0,,-,-,,0,-,,,,-];d=;%改变d的值即可改变曲面四边形的个数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;whileq<10/d+1whilep<10/d+1s=sqrt(1+Fx(p,q)^2+Fy(p,q)^2)*d*d;%s为每一个曲面四边形的面积ss=s+ss;%对s循环相加p=p+1;endp=1;q=q+1;endss(=)③两种方法的比较:,两种方法算得的工件表面面积最终都逼近到676左右,因此,可以认为工件表面面积约为676。还可以看出,当d取的值不是很大的情况下,方法一的结果更接近准确值,方法一误差相对较小。