计算国体面积PPT-MATLAB.pptVIP

具体数据 下边疆采样坐标 x17182031415866727269576071104y299298288273262254234220207191175166160150x130146160163168179196223258282307315330352y13712111710683646356505246383221x377377392428462501524533555542550561574590y211614344346607595114138139133133x599610635644649669671677678696720723y139157162174188200207205206216218225 上边疆采样坐标 x723722710687676659647630619623626633608y225220240256256241245237245254273309308x596581558537511484464456449434425411394y315315290281270270272278290293301303308x368351332329312284281263251249244240247y297303311337342353358365356347346332314x233222217209189180169165165150138138132y297290297298301303307314325328332337336x12712210286656454322817y341338332328322316314314307299 模型假设 提出以下假设建模思路： 1. 国土的边界是光滑的连续曲线，能用3次多项式表示出来。 2. 国土边界曲线有连续二阶导数。 3. 题目所给数据是有选择性的，真实地反映了中国国土的大概轮廓。 4. 假设测量的地图和数据准确，由最西边界与最东边界分为n条连续的边界曲线，边界内所有的土地均为该国国土。 问题分析 1.从最西边界点到最东边界点，变量x∈[a,b],划分[a,b]为n小段[xi-1.xi],并由此将国土分成n小块，设每一块均为X型区域，即做垂直于x轴的直线穿过该区域，直线与边界曲线最多只有两个交点。 2.数值积分方法计算国土面积：根据测量数据利用matlab软件对上下边界进行三次多项式插值。数值积分法的基本思想是将上边界点与下边界点分别利用差值函数求出n条曲线，则曲线所围面积即为国土面积，然后根据比例缩放关系求出国土面积近似解。在求国土面积时，利用求平面图形面积的数值积分法--将该面积分成若干个小长方形，分别求出长方形的面积后相加即为该面积的近似解。 分割替代近似代替，在每个小区间[xi，xi-1]上任取一点ei（i=1,2,…n）,作以[xi，xi-1]为底，以f（ei）为高的小矩形，用其面积f（ei）ei近似代替同底小曲边梯形的面积。 模型建立 三次样条插值值法是一种分段插值法，由于在插值节点处具有二阶导数连续，从而具有更好的光滑性。三次样条插值用分段低次（不高于三次）的多项式去逼近，比较简单方便。 首先，我们要将原图形画出来（图3），观察图形的起伏程度，方便为后面求面积分段，例如，看下边疆，我们要在x中相同的点处将x分段，根据分段的x和对应的y值经行三次样条插值，画出图形（图1），与原图形经行比较（图2），观察哪里的差入比较大，在差入大的地方将x分段（根据个人观察，并不是分的越小越好）；直到画出的图形与原图差入不大。上边疆一样。 插值后的图形：（图4） 原图形与差值后的比较：（图5） 举例面积有重复的图形（浅灰色代表重复的面积）：（图6） 在图形上标出分段地方： 程序分析  1）根据差值函数，画出下边疆的图像（上边疆如同） 程序如下： x11=[17 18 20 31 41 58 66 72]; y11=[299 298 288 273 262 254 234 220]; x11i=17:1:72; y11i=interp1(x11,y11,x11i,spline); plot(x11i,y11i) hold on x=[72 72]; y=[220 207]; plot(x,y,b) x12=[72 69 57 60]; y12=[207 191 175 166]; x12i=57:1:72; y12i=interp1(x12,y12,x12i,spline); plot(x12i,y12i) hold on x13=[60 71 104 130 146 160 163 168