数值分析法建模.doc

第五章 数值分析法建模 在化学工程中经常会遇到计算高温状态下蒸汽压力和温度的问题，但是考虑到测量设备等的限制，我们希望利用低温状态下的压力等有关数据进行外推。表5.1给出了氨蒸汽的一组温度和压力数据。那么能否从所列的数据中计算75氨蒸汽的压力(kN/m2)？ 表5.1氨蒸汽的压力和温度 温度 20 25 30 35 40 45 50 55 60 压力 805 985 1170 1365 1570 1790 2030 2300 2610 类似的问题在科技研究或解决实际问题过程中会经常碰到。许多情况下，我们除了要进行一定的理论分析外，实验观测或相应的数据处理也是必不可少的步骤。由于实验测定实际系统的数据具有一定的代表性，所以在处理时必须充分利用这些信息。由于测定过程中不可避免会产生误差，所以在分析经验公式时又必须考虑这些误差的影响，这两者之间是相互制约的。据此合理建立实际系统数学模型的方法称为数值逼近法。本章介绍插值法与曲线拟合法建模。 5.1曲线拟合法建模 工程实践和科学实验中，常常需要从一组实验观测数据（）中，求自变量x与因变量y的一个近似的函数关系式y=f（x）。比如，程序控制的铣床加工精密工件时，加工时走刀方向是直线的，要加工得到光滑的曲线外廓，必须计算出外形曲线足够密的点的坐标，用它来控制走刀方向，以便加工出达到精度要求的工件（见图5.1）,曲线在各点的切线方向为走刀方向。一般来说的形状不是非常确定的，我们只能根据 （）,=1，2，… 的数据，给出 的足够精确的近似函数—经验公式。 从几何角度解释，就是找一条“最佳” 的曲线，使之与靠得最近，这就是所谓的曲线拟合。 通过找经验公式或曲线拟合将离散的数据条件变成连续的函数条件，从而可以使用连续函数的分析方法进行建模讨论。 曲线拟合法的一般步骤是先根据实验数据，结合相关定律，将要寻求的最恰当的拟合曲线方程形式预测出来，再用其他的数学方法将经验公式中的参数确定下来。 对于事先给定的一组数据，确定经验公式一般可分三步进行： （ = 1 \* roman i）确定经验公式的形式 根据系统和测定数据的特点，并参照已知图形的特点确定经验公式的形式。 （ = 2 \* roman ii）确定经验公式中的待定系数 计算待定系数的方法有许多，常用的方法有图示法、均值法、差分法、最小二乘法、插值法等。 （ = 3 \* roman iii）检验 求出经验公式后，还要将测定的数据与用经验公式求出的理论数据作比较，验证经验公式的正确性,必要时还要修正经验公式。 关于确定经验公式的形式，我们可以从下面几个方面着手： ( = 1 \* roman i) 利用已知的结论确定经验公式形式， 比如由已知的胡克定律可以确定在一定的条件下，弹性体的应变与应力呈线性关系等。 ( = 2 \* roman ii) 从分析实验数据的特点入手，将之与已知形式的函数图形相对照，确定经验公式的形式。 ( = 3 \* roman iii)用描点作图法，将已知的点用光滑曲线连接起来，寻找曲线的形式。 ( = 4 \* roman iv)多项式近似、线性插值或样条插值等。 多项式近似是工程中十分常见的方法，它首先需要我们确定多项式的次数，一般可以用差分法、差商法来估计。 差分与差商概念 一阶前向差分 二阶前向差分 m阶前向差分 一阶差商 二阶差商 m阶差商 差商与导数的联系 （微分中值定理）若y=f(x)在[a，b]上m次可导，且 ，则 ， 据此，若结点为等距分割点时，有，h为结点距，且 ，。 因此对n阶多项式有 为常数。据此,我们可以根据数据的差分或差商来确定多项式的次数。 考虑本章开始提出的问题，自变量取等距节点，可考虑差分，根据表5.1的数据可以得到差分表5.2。 表5.2 差分表 p 20 805 25 985 180 30 1170 185 5 35 1365 195 10 5 40 1570 205 10 0 -5 45 1790 220 15 5 5 50 2030 240 20 5 0 55 2300 270 30