【2017年整理】计算方法在数学建模中的应用.pptVIP

计算方法在数学建模中的应用; 计算方法又称“数值分析”，是为各种数学问题的数值解答研究提供最有效的算法。计算方法在数学建模中有很广泛的应用，特别是处理具有大规模数据的数学建模问题。其中应用最多的是误差分析，数值逼近以及数值积分等。在历年的全国大学生数学建模竞赛(CUMCM)赛题之中，有很多都用到了计算方法，其中以数据拟合方法用得最多，其次就是数值积分方法。下面分别对这些方法进行介绍，并通过几个赛题来熟悉这些方法的具体应用。;一、数据拟合; 应用最小二乘法，使 取最小。对上式各参变量求偏导数后，令其偏导数为0，就得到了以下线性方程组 从该线性方程组中求解出系数aj(j=0,1,…,k)之后，就可以写出拟合的多项式。; 有时还会采用三角函数，指数函数等来拟合。一般来说，多项式函数拟合的性态更易描述，而且误差分析较为简单。 ;例1：CUMCM1992A题施肥效果分析; ;问题的分析和求解; 然后应用最小二乘法，令 取最小，其中为xi施肥量，yi为产量，n为样本个数。对上式各参变量求偏导数后，令其偏导数为0，就得到了关于a,b,c的线性方程组 从中求出这三个系数，进而写出其近似关系式。可以用matlab软件的polyfit命令或者mathematics软件的fit命令来求出对应的系数。; ; ; ; 以上结果从数据的分布情况得出的，但根据实际情况，并非施肥量越多，作物的产量就越高。施肥量过高时会导致土壤中溶液浓度过高而导致“烧苗”现象。N肥的关系式体现出了这一点，但是P肥和K肥的关系式则没有体现出这一现象，因此对于P肥和K肥也可以利用二次函数进行数据拟合，其拟合结果为： 在实际工作中，三种肥料除了与产量有直接关系外，还有彼此之间的交互关系，因此要准确研究施肥量与产量的关系，还需进一步考虑到几种肥料之间的交互作用。;例2：CUMCM2004B题电力市场的输电阻塞管理; 电网中每条线路上的有功潮流(输电功率和方向)取决于电网结构和各发电机组的出力。电网每条线路上的有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限） 。如果各机组出力分配方案使某条线路上发生输电阻塞时，需要研究如何制订既安全又经济的调度计划。某电网有8台发电机组，6条主要线路，表1和表2中的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案1~32给出了围绕方案0的一些实验数据，试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。; ; ;问题的分析和求解; 根据表一和表二中给出的0-32号方案中各条线路上有功潮流与各个发电机组的出力，用最小二乘法，求拟合函数中的54个系数，归结为以下最小化问题： 对上式各参变量求偏导数后，令其偏导数为0，就得到了以下线性方程组：; 求解以上方程组，得到系数aij的值。进而就得到了各线路上的潮流与各机组出力的近似线性关系式为：; ;二、数值积分; ; ; ; ;例3：CUMCM2010A题储油罐的变位识别与罐容表标定;(1) 为了掌握罐体变位后对罐容表的影响，利用小椭圆型储油罐(图4)，分别对罐体无变位和倾斜角为α=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 (2) 对于实际储油罐(图1)，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数(纵向倾斜角度α和横向偏转角度β)之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据(附件2)，根据所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。; ;问题的分析和求解; 通过积分可以得到储油罐内油的体积随高度h的变化关系式为 可以通过积分法求出来 将油位高度h间隔为1cm时的体积值计算出来，即就是所求的罐容表。为了和测量值作比较，通过画图给出不同油位高度h下储油量V的理论值和测量值的比较情况。; ; 从上图的结果可以看出，理论储油量与实测储油量之间存在较大误差, 而且测量值普遍高于理论值，这是由于理论模型忽略了油量探针浸入油内的部分体积，从而导致同样高度下的理论储油量较多。为了简单起见，不妨将该偏差 用关于高度h的三次多项式函数来近似，即假设 然后，通过不同高度下的 以及其对应的高度h的数据，应用最小二乘数据拟合的方法，就可以求出对应的系数，结果为 ; 通过