数值分等析k模型.docVIP

数值分析模型 数值分析分析是研究用计算机求解各种数学计算问题的数值计算方法及其理论与软件实现的学科。也可以认为数值分析就是介绍如何用计算机来解决数学问题，以各种各样的程序语言来设计出数值计算程序，然后依靠计算机的强大计算能力来求解这些数学问题，数值分析对数学理论与程序设计并重。 运用数值分析解决问题的过程可分为如下几步：实际问题→数学模型→数值计算方法→程序设计→上机计算求出结果。 数值分析这门学科有如下特点： 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]，以x1代替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 = ( (x0)，再将x1又可得x2 =( (x1)，…。我们称{xk}为迭代序列，而称（1）中的( (x)为迭代函数，（2）为迭代格式。如果 连续，迭代序列{ }收敛于x*，则x*就是方程（2）的解。 所以，如果迭代序列收敛，总能收敛于原方程的解。实际计算中，无穷过程不可能实现，只迭代到一定程度，取xk+1作为原方程的近似根。这种求根方法称为简单迭代法。 序列{ }的收敛速度，取决于曲线 在根附近的斜率 ，由拉格朗日定理可知 其中( k在xk和xk-1之间。所以在根x*附近，((’(x) (恒小于1，则此迭代序列收敛，若((’(x) ( 1，则此序列发散。 定理：不定点存在定理 如果 满足下列条件 （1）当x([a, b]时，((x)([a, b] ，（2）存在正常数L1,当任意x,y([a, b]，都有 。则 ( (x)在[a, b]上有唯一的根x*。 （2）牛顿迭代法 设：已知方程f (x) = 0的一个近似根x0，把f (x)在x0处作泰勒展开， 若取前两项来近似代替f (x)（称为f (x)的线性化），则得近似的线性方程 设 ，解之得 我们取x作为原方程f (x) = 0的近似根x1，即 ，一般地f (x) ( 0。 再重复用上述方法得 一般地，有迭代公式 称为求解f (x) = 0的牛顿迭代公式。 （3）牛顿下山法 由牛顿法的收敛性定理知，牛顿法对初始值x0的选取要求是很高的，为保证收敛，常要利用条件f (x0) ( f( (x0) 0来选取初值x0，但对有些问题，往往很难检验满足条件的初x0，这时我们可利用所谓下山法来扩大初值的选取范围，将牛顿的迭代公式修改为 其中(是一个参数，(的选取应使 成立 ，当 或 时（其中( 1，( 2为事先给定的精度，称( 1为残量精确度，( 2为根的误差限），就停止迭代，且取x*( xk+1，否则再减小(，继续迭代。 按上述迭代过程计算，实际上得到一个以零为下界的严格单调下降的函数值序列{(f (xk)(}。这个方