Matlab曲面面积估计.docVIP

2013-2014第一学期 《数学软件与数学实验》 课程设计 2013年月日月日 实验所用软件及版本： Matlab（R2010b） 主要内容（要点）：（模型的求解原理、公式、推导、基本求解步骤、算法的流程图等） (1) 主要内容（要点）：（模型的求解原理、公式、推导、基本求解步骤、算法的流程图等）（接上页）： （1） 1、有题目数据表格给定，定义x，y范围以及步长，并将x、y所对应的值赋值给z x=2:0.1:3;y=1.0:0.1:2; 2、用linspace产生多个元素个数，并用产生的数据绘制网格 xi=linspace(2,3,80);yi=linspace(1,2,80); [XI,YI]=meshgrid(xi,yi); 3、利用三次样条（cubic）插值，并绘出曲面图像 ZI=interp2(X,Y,Z,XI,YI,cubic);surfc(XI,YI,ZI) 4、绘制原始网格曲面与数据点以及等高线 plot3(X,Y,Z,go,markeredgecolor,c) plot3(X,Y,Z,g) hold off rotate3d on figure contour(XI,YI,ZI,40) colorbar 5、利用sftool拟合工具箱产生拟合函数 sftool(X,Y,Z) （2） （转下页） 主要内容（要点）：（模型的求解原理、公式、推导、基本求解步骤、算法的流程图等）（接上页）： Y N Y N 1、将曲面分成若干个小三角形，由每个点算出三角形的边长，由海伦公式计算三角形的面积。a1=(x1-x2)^2+(y1-y2)^2+(z1-z2)^2;a1=sqrt(a1); b1=(x4-x2)^2+(y4-y2)^2+(z4-z2)^2;b1=sqrt(b1); c=(x1-x4)^2+(y1-y4)^2+(z1-z4)^2;c=sqrt(c); a2=(x1-x3)^2+(y1-y3)^2+(z1-z3)^2;a2=sqrt(a2); b2=(x4-x3)^2+(y4-y3)^2+(z4-z3)^2;b2=sqrt(b2); p1=(a1+b1+c)/2; p2=(a2+b2+c)/2; 主要内容（要点）：（模型的求解原理、公式、推导、基本求解步骤、算法的流程图等）（接上页）： 2、将每个小面积利用循环相加，得总面积。 sum=sum+sqrt(abs(p1*(p1-a1)*(p1-b1)*(p1-c))); sum=sum+sqrt(abs(p2*(p2-a2)*(p2-b2)*(p2-c))); （3） Y （转下页） 主要内容（要点）：（模型的求解原理、公式、推导、基本求解步骤、算法的流程图等）（接上页）： 1、定义function length=length(x1,y1,x2,y2) 2、步长取0.001，分别赋值 x11=x1+i*lengthx;x22=x11+lengthx*0.001; y11=y1+i*lengthy;y22=y11+lengthy*0.001; 3、带入方程，算出相应的z。 4、利用公式算出最终长度 leng=(lengthx*0.001)^2+(lengthy*0.001)^2+(z1-z2)^2; length=sqrt(leng)+length; 实验过程记录（含：主要程序清单及异常情况记录等）： （1） x=2:0.1:3; y=1.0:0.1:2; [X,Y]=meshgrid(x,y); Z=[ 5.11 5.13 5.14 5.13 5.09 5.04 4.98 4.93 4.89 4.85 4.85 5.39 5.49 5.51 5.46 5.32 5.14 4.94 4.74 4.59 4.49 4.48 5.61 5.77