1、4.4.1 仿真球曲面面积计算step1=0.01;step2=0.01;R=1;mj=0;for z=0:step1:1-step1t=0:step2:2*pi-step2;r=sqrt(R2-z2);x1=r*cos(t);y1=r*sin(t);x3=x1(2:end) x1(1);y3=y1(2:end) y1(1);z=z+step1;r=sqrt(R2-z2);x2=r*cos(t);y2=r*sin(t);x4=x2(2:end) x2(1);y4=y2(2:end) y2(1);a=sqrt(x1-x2).2+(y1-y2).2+step1.2); b=sqrt(x1-x4).
2、2+(y1-y4).2+step1.2);c=sqrt(x2-x4).2+(y2-y4).2);s=(a+b+c)/2;s=sqrt(s.*(s-a).*(s-b).*(s-c);mj=mj+s;a=sqrt(x4-x1).2+(y4-y1).2+step1.2);b=sqrt(x4-x3).2+(y4-y3).2+step1.2);c=sqrt(x1-x3).2+(y1-y3).2);s=(a+b+c)/2;s=sqrt(s.*(s-a).*(s-b).*(s-c);mj=mj+s;enddisp('数值法计算的面积=',num2str(sum(mj)disp('直接法计算的面积=',num2str(2*pi*R2)
