《数值分析》课程设计-国土面积计算.docVIP

国土面积计算 摘 要：数学建模方法是处理科学理论的一种经典方法，也是解决各类实际问题的常用方法。本文采用插值、复化梯形公式的方法，并利用数学软件MATLAB对国土面积进行计算，使之计算结果与实际记录基本吻合。 关键词 ：插值 复化梯形公式 MATLAB 一、问题提出 图3.8是某国的地图，为了计算它的国土面积，首先对地图作如下测量：以由西向东方向为轴，由南到北方向为轴，选择方便的原点，得到了表3.6、表3.7的地图测量数据，比例尺为30（数据单位）：100公里（实际单位）。试由测量数据采用插值的方法产生一张需要的地图，计算该国国土的近似面积，与它的精确值156.6500万平方公里进行比较。 二、模型假设 题中所给数据是真实可用的，假设国土面积在一定时期内没有变化，各种因素对此题也无影响，可以进行实地测量。 三、问题分析与模型建立 在工程建设和地籍管理中，会经常遇到面积的测算和计算工作，传统的方法是在纸上，利用求积仪进行计算，存在绘图等操做的误差，精度较低。在现在工作中，全站仪的广泛使用我们能够容易得到一系列离散数据坐标，降低了对使用者的基础和计算机语言的要求，使计算不在是测算人员的负担。MATLAB提供了非常方便的绘图功能，越来越受测量人员的青睐。现根据题目给出的数据，我们可以使用MATLAB软件中的plot函数先大概画出该图形轮廓，从该图可知本题主要是求不规则平面图形的面积，题中已经把图形分为上下边疆线，所以可以采用线性插值与复化梯形公式分别求出上下边疆与横坐标所围的面积，最后求出图形的总面积。 根据题目中给出的已知条件，应用机理分析进行建模具体分析如下： 根据对某地地图测量的数据计算出某地地图的近似面积。 将测算数据与其精确值156.6500万平方公里进行比较。 四、模型求解 1.插值作图 在将数据进行插值时，由于数据中出现重复点，故不能直接插值，所以考虑分段处理，分别将上下边疆线合理分段，每段记为xi(i=1,2,……),而y的值也与之对应分段记为yi(i=1,2,……)，本文将下边疆线分为六段，将上边疆线分为七段，现以下边疆线为例分段： 第一段：x1=[17 18 20 31 41 58 66 72]；y1=[299 298 288 273 262 254 234 220]；第二段：x2=[72 69 57]；y2=[207 191 175]；第三段：x3=[57 60 71 104 130 146 160 163 168 179 196 223 258 282 307 315 330 352 377]；y3=[175 166 160 150 137 121 117 106 83 64 63 56 50 52 46 38 32 21 21]；第四段：x4=[377 392 428 462 501 524 533 555];y4=[16 14 34 43 46 60 75 95]；第五段：x5=[555 542];y5=[95 114]；第六段：x6=[542 550 561 574 590 599 610 635 644 649 669 671 677 678 696 720 723]；y6=[114 138 139 133 133 139 157 162 174 188 200 207 205 206 216 218 225]； 下面先将每段进行一维线性插值如第一段: x1=[17 18 20 31 41 58 66 72]; y1=[299 298 288 273 262 254 234 220]; x1i=17:1:72; y1i=interp1(x1,y1,x1i); %线性插值 plot(x1i,y1i) hold on; 再考虑将每一段进行一维三次样条插值(1)进行比较，同样以第一段为例：y1i=interp1(x1,y1,x1i,’spline’) %三次样条插值 plot(x1i,y1i) hold on; 其中x1,y1,x1i,y1i的含义同上，而’spline’为插值的类型是三次样条插值，再次得出插值后的图如下：（蓝线为插值后的线) 根据以上两图的比较可以看出线性插值与三次样条插值存在着一定的差别，但是只根据图形无法判断它们的精确性，因此下面再分别用它们插值后的线来计算面积。 2.求面积 对以上每一段所构造的函数值用复化梯形公式（I=） （h=1）求积分，记为An，n=（1，2，……6），下以A1为例： s=0; for i=1:54 运用for循环语句s=s+y1i(i); end A1=0.5*(299+220+2