MATLAB--计算国土面积(2)

MATLAB--计算国土面积 数值分析 课程设计 历时两周 完美通过!!

x6=[542 555];

为了作出一个合理的图像，首先考虑用线性插值。记为yni,其中n=1,2,3……

以x1为例

x1i=17:1:72;

y1i=interp1(x1,y1,x1i);

y1i就是以区间[17 72]运用线性插值构造的函数值，x1,y1是插值节点，x1i是被插值点。所以我们对被分段的上下边疆的数据进行线性插值。

2．2求面积

对每段构造的函数值分别用复化梯形公式（trapz）求积分,记为sn,其中n=1,2,3……

以s1为例 s1=trapz(y1i)*1

设区间[17 72]之间的面积为s1。s1的面积是蓝色的部分，同理可求上、下边疆的面积。sx=s1+s2+s3+s4+s5-s6; ss=s7+s8+s9+s10+s11+s12-s13.