Question

大学matlab定积分求解

求地球质量

Answer 1

解题思路：将地球沿径向分割成球壳，将球壳的质量加总，即为地球的质量。

r = [0 1100 1500 2450 3400 3630 4500 5380 6060 6280 6380]*1000*100;

Rho = [13 12.4 12 11.2 9.7 5.7 5.2 4.7 3.6 3.4 3];

S = 4*pi*r.*r;  

RhoS = Rho.*S;

Z = trapz(r,RhoS)/1000

计算结果：

6.1087e+024 kg

Answer 2

在matlab中，对于表格形式定义的函数关系求定积分，可以用trapz(x,y)函数。这种思想是基于折边梯形法求面积。