数学建模社区-数学中国

标题: 数值分析模型 [打印本页]

作者: 总有以后 时间: 2014-11-8 22:16
标题: 数值分析模型
数值分析分析是研究用计算机求解各种数学计算问题的数值计算方法及其理论与软件实现的学科。也可以认为数值分析就是介绍如何用计算机来解决数学问题，以各种各样的程序语言来设计出数值计算程序，然后依靠计算机的强大计算能力来求解这些数学问题，数值分析对数学理论与程序设计并重。
运用数值分析解决问题的过程可分为如下几步：实际问题→数学模型→数值计算方法→程序设计→上机计算求出结果。
数值分析这门学科有如下特点：
1、       面向计算机。
2、       有可靠的理论分析
3、       要有好的计算复杂性。
4、       要有数值实验。
5、       要对算法进行误差分析。
数值分析包括的主要内容有插值法、函数逼近、曲线拟合、数值积分、数值微分、解线性方程组的直接方法，解线性方程组的迭代法，非线性方程球根，微分方程的数值解法。
1插值法
插值法是函数逼近的一种重要方法，是数值计算的基本内容。
（1）拉格朗日插值
拉格朗日插值是n次多项式插值，求n次多项式插值函数问题就是构造插值基函数。其基本思想是将求的n次多项式插值函数改写成另一种表示方式，再利用插值条件确定其中的待定函数，从而求出插值多项式。
两点和做一条直线。  为确定a，b，把两点带入方程，
  令
把叫做点的一次插值基函数。
（2）牛顿插值
基本思路：
给定插值点序列（。构造牛顿插值多项式。输入要计算的函数点并计算的值，利用牛顿插值公式，当增加一个节点时，只需在后面多计算一项，而前面的计算仍有用；另一方面的各项系数恰好又是各阶差商，而各阶差商可用差商公式来计算。
公式：

2非线性方程求根
   执行步骤
   1．计算f (x)在有解区间[a, b]端点处的值，f (a)，f (b)。
   2．计算f (x)在区间中点处的值f (x1)。
   3．判断若f (x1) = 0，则即是根，否则检验：
   （1）若f (x1)与f (a)异号，则知解位于区间[a, x1]，以x-1代替b；
   （2）若f (x1)与f (a)同号，则知解位于区间[x1, b]，x1代替a。
   反复执行步骤2、3，便可得到一系列有根区间：(a, b), (a1, b1), …, (ak, bk), …
其中每个区间都是前一个区间的一半，因此区间长度为这些区间最终必收敛于一点x*，该点就是所求的根。
3迭代法
迭代法是求方程近似根的一个重要方法，也是计算方法中的一种基本方法
（1）简单迭代法
简单迭代法的基本思想是：将方程f (x) = 0化为一个等价的方程  （1）
从而构成序列  （2）
给定一个初值x0，可算得x1 = j (x0)，再将x1又可得x2 =j (x1)，…。我们称{xk}为迭代序列，而称（1）中的j (x)为迭代函数，（2）为迭代格式。如果连续，迭代序列{ }收敛于x*，则x*就是方程（2）的解。
所以，如果迭代序列收敛，总能收敛于原方程的解。实际计算中，无穷过程不可能实现，只迭代到一定程度，取xk+1作为原方程的近似根。这种求根方法称为简单迭代法。
序列{ }的收敛速度，取决于曲线在根附近的斜率，由拉格朗日定理可知    其中x k在xk和xk-1之间。所以在根x*附近，|j’(x) |恒小于1，则此迭代序列收敛，若|j’(x) | 1，则此序列发散。
定理：不定点存在定理
如果满足下列条件（1）当xÎ[a, b]时，j(x)Î[a, b] ，（2）存在正常数L<1,当任意x,yÎ[a, b]，都有  。则 j (x)在[a, b]上有唯一的根x*。
（2）牛顿迭代法
设：已知方程f (x) = 0的一个近似根x0，把f (x)在x0处作泰勒展开，

若取前两项来近似代替f (x)（称为f (x)的线性化），则得近似的线性方程
  设，解之得
我们取x作为原方程f (x) = 0的近似根x1，即，一般地f (x) 1 0。
再重复用上述方法得    一般地，有迭代公式
  称为求解f (x) = 0的牛顿迭代公式。
（3）牛顿下山法
由牛顿法的收敛性定理知，牛顿法对初始值x0的选取要求是很高的，为保证收敛，常要利用条件f (x0) × f2 (x0) >0来选取初值x0，但对有些问题，往往很难检验满足条件的初x0，这时我们可利用所谓下山法来扩大初值的选取范围，将牛顿的迭代公式修改为
   其中l是一个参数，l的选取应使成立，当或时（其中e 1，e 2为事先给定的精度，称e 1为残量精确度，e 2为根的误差限），就停止迭代，且取x*» xk+1，否则再减小l，继续迭代。
按上述迭代过程计算，实际上得到一个以零为下界的严格单调下降的函数值序列{|f (xk)|}。这个方法就称为牛顿下山法。l称为下山因子，要求满足，el为下山因子下界，为了方便，一般开始时可简单地取l = 1，然后逐步分半减少，即可选取，且使
   牛顿下山法计算步骤可归纳如下：
   （1）选取初始近似值x0；
   （2）取下山因子l = 1；
   （3）计算
   （4）计算f (xk+1)，并与的大小，分以下二种情况：
         1）若，则当时，取x*» xk+1，计算过程结束；当时，则把xk+1作为新的xk值，并重复回到（3）。
         2）若，则当l≤el且，取x*» xk，计算过程结束；否则若l≤el，而时，则把xk+1加上一个适当选定的小正数，即取xk+1+d作为新的xk值，并转向（3）重复计算；当l＞el；且，则将下山因子缩小一半，取l/2代入，并转向（3）重复计算。
   牛顿下山法不但放宽了初值x0的选取，且有时对某一初值，虽然用牛顿法不收敛，但用牛顿下山法却可能收敛。
4弦截法和抛物线法
牛顿迭代法虽然有较高的收敛速度，但要计算导数值f’ (xk)，这对复杂的函数f (x)是不方便的。为避免导数的计算，用平均变化率  来替代迭代公式中的导数 f’ (xk)，于是得  按这个公式进行迭代计算就称弦截法。如果牛顿迭代公式中的导数f’(xk)改用来代替，就可以得到迭代公式：，由这个公式确定的迭代法称双点弦截法。

作者: lingxin179 时间: 2015-9-4 15:08
看看，谢谢分享

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)