建模中常用的数值计算演示文件修改版.pptVIP

建模中常用的数值计算方法及应用实例 吉林大学数学科学院 刘金英 2005年9月 引言 由于数学模型可以描述为：对于现实世界的一个特定对象，为了一个特定目的，根据特有的内在规律，做出一些必要的简化和假设，运用适当的数学工具，得到的一个数学结构。所以一般来说建立一个数学模型的过程就可以分为表述、解释、验证这几个阶段. 表述（Formulation）是指根据建模的目的和掌握的信息（如数据、现象）将实际问题翻译成数学问题，用数学语言确切地表述出来. 解释（Interpretation）是指把数学语言表述的解答翻译回现实对象，给出实际问题的解答. 验证(Vertification)是指用现实对象的信息检验得到的解答，以确认结果的正确性. 从这个过程来看建模过程主要是双向翻译过程，如何将一个用普通语言文字表达的问题通过适当的简化和假设表述为数学问题是至关重要的.这里谈到的简化和假设要服从于数学概念的抽象、数学定理的逻辑，适于用数学工具处理，而不是合情合理的简化和假设。记住一点模型是种近似描述，而不是拷贝.所以不要试图将所有的因素都考虑进去，数学处理的可行性与实际问题的真实性是相对的.越逼真的模型常常越复杂，即使数学上能处理，这样的模型应用时所需要的费用也会相当高，而高费用不一定与复杂模型取得的效益相匹配，所以建模时往往需要在模型的逼真性与可行性，费用与效益之间做出折中的选择。 主要内容 1 多项式插值 2 曲线拟合 3 线性代数方程组的数值解法 4 非线性方程求根 5 Matlab中微分方程数值解的常用命令 求解的基本思路：构造一个相对简单的函数y=f(x)，（通常取多项式）使f通过全部节点，即f(xj)=yj(j=0,1,…,n)，再用f(x)近似代替插值点的函数值，即y*=f(x*)． 最简单的插值方法是：先根据基准数据，调用MATLAB的绘图指令（如plot,mesh），获得数据图形表现.然后，再估计所需点处的值．当然这样会比较粗略，下面介绍应用较为广泛的三种插值方法及其指令． 1.1 Lagrange多项式插值 Lagrange插值多项式的构造思想：对给定的n+1节点，先构造一组Lagrange基函数： 　 是n次多项式，满足 令 1.2分段线性插值 分段线性插值函数的构造思想：是将两个相邻的节点用直线连起来．如此形成的一条折线就是分段线性插值函数，也记作，它满足且在每个小区间[xi,xi+1]上是线性函数．可以表示为 (1.1.6) 对于二元函数的插值问题同样具有命令函数: z=interp2(x0,y0,z0,x,y，‘method’) 这里注意输入x0,y0的维数如果分别为n,m,那么z0的维数就是m*n.，x0,y0,z0是已知的数据，x，y是待求点即所求插值点，输出为插值点的函数值.例如94年赛题逢山开路，如果能将地形图画处理将十分有助于寻找搜索路径.我们利用二维插值给出了连续型地形图.原始数据地形图. Method有四种 (1)、’linear’双线性型；(2)、’nearest’最近邻域插值； (3)、‘spline’三次样条插值；(4)、’cubic双三次插值 * * 因此我们的原则是能用初等的方法不用高等的方法；能用线性的不用非线性的，也就是工程师原则-----线性化. 建模需要的能力是多方面的，但主要的一点是会将复杂的问题简单化、层次化、结构化。具体地说会将一个复杂的问题肢解成若干个相对独立的小问题，每个小问题都能借用类比的方法给出相应的模型，然后再能将小问题组装成一个综合问题。而每个小的模型都可以通过已知的算法进行求解。 MATLAB有强大的数值计算功能，由于时间所限，只介绍在数学模型中常用的几个数值计算问题．对于有关的数学知识，只做简单的介绍或不作介绍，具体可参阅计算方法的有关内容． 1 多项式插值 例1 已知温度和电阻的一组实测数据如下： t(0C) 20.5 32.7 51.0 73.0 82.3 95.7 100 R()765 826 873 942 ？ 1032 ？ 如何获得温度t＝82.3时的电阻以及t＝100时的电