1701008080
1701008081
1701008082
当且仅当H＝R，S＝2πR2＝2π×302＝5654.9。
1701008083
1701008084
实际用量将会比清真寺顶部面积多1.5％。
1701008085
1701008086
1701008087
1701008088
1701008089
具体的MATLAB程序如下：
1701008090
1701008091
clc,clear,close all %清屏和清除变量 warning off %消除警告 syms x y z R r H z=sqrt(R^2-x^2-y^2); dz_dx=diff(z); dz_dy=diff(z,‘y’); z1=sqrt((dz_dx)^2+(dz_dy)^2+1); z2=R/sqrt(R^2-r^2); z3=z2*r; Intxy=int(int(z3,‘r’,0,‘R’),‘y’,0,2*pi); >> pretty(Intxy) 2 1/2 2 (R ) R pi >> R=30; >> Intxy=2*(R^2)^(1/2)*R*pi Intxy = 5.6549e+003 >> Intxy*1.015 ans = 5.7397e+003
1701008092
1701008093
（2）对于半椭球而言，其方程如下：
1701008094
1701008095
1701008096
1701008097
1701008098
其偏导数为：
1701008099
1701008100
1701008101
1701008102
1701008103
则椭圆球半盖的表面积积分编程如下：
1701008104
1701008105
clc,clear,close all %清屏和清除变量 warning off %消除警告 >> syms x y z a b c >> z=sqrt(c^2*(1-y^2/b^2-x^2/a^2)); >> dz_dx=diff(z,‘x’) dz_dx = -1/(c^2*(1-y^2/b^2-x^2/a^2))^(1/2)*c^2*x/a^2 >> dz_dy=diff(z,‘y’) dz_dy = -1/(c^2*(1-y^2/b^2-x^2/a^2))^(1/2)*c^2*y/b^2 >> z1=sqrt((dz_dx)^2+(dz_dy)^2+1); >> pretty(z1) / 2 2 2 2 \1/2 | c x c y | |––––––— + ––––––— + 1| |/ 2 2 \ / 2 2 \ | || y x | 4 | y x | 4 | ||1 - –- - –-| a |1 - –- - –-| b | || 2 2 | | 2 2 | | \ b a / \ b a / / >> %Intxy=int(int(z1,‘x’,-sqrt(a^2*(1-y^2/b^2)),sqrt(a^2*(1-y^2/b^2))),‘y’,-b,b)
1701008106
1701008107
整理结果如下：
1701008108
1701008109
1701008110
1701008111
1701008112
从上式中可知，该方程为二次椭圆积分问题方程无相应的原函数与之对应求解，故运用MATLAB较难无法求解，采用半椭圆体表面积计算公式（GB）有：
1701008113
1701008114
1701008115
1701008116
1701008117
1701008118
其中，为离心率，a为长半轴，b为短半轴。
1701008119
1701008120
则椭圆表面积编程计算如下：
1701008121
1701008122
clc,clear,close all %清屏和清除变量 warning off %消除警告 a=30.6; %长半轴 b=29.6; %短半轴 %离心率 e=sqrt((a^2-b^2))/a %表面积 A1=pi*b*b+pi*a*b/e*asin(e) A2=pi*b*(a+b+(a^2-b^2)/6/a) A3=A1*1.015 A4=5800-A3
1701008123
1701008124
运行程序输出结果如下：
1701008125
1701008126
e = 0.2536 A1 = 5.6295e+03 A2 = 5.6286e+03 A3 = 5.7139e+03 A4 = 86.0831
1701008127
1701008128
整理结果得：
1701008129