||
利用MATLAB计算如下:
[K,E] = ellipke(1/2)
K =1.8541
E =1.3506
[K,E] = ellipke(1/4)
K = 1.6858
E =1.4675
由此能看出,利用MATLAB求解第一类和第二类完全椭球积分需要将输入参数平方才能得到正确结果;
千万注意:MATLAB求解完全椭圆积分时,需要特别注意如下:m=k^2
在求解Diskload荷载引起的地表位移时,公式如下:
注意,上式中E(R/r)和K(R/r),利用MATLAB求解时,应该输入(R/r)^2
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2025-1-15 19:42
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社