(数值分析简述及求解应用.docVIP

数值分析简述及求解应用 摘要：数值分析是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科它以数字计算机求解数学问题的理论和方法为研究对象。有可靠的理论分析要有数值实验主要内容插值法，函数逼近，曲线拟和，数值积分，数值微分，解线性方程组的直接方法，解线性方程组的迭代法，非线性方程求根，常微分方程的数值解法。运用数值分析解决问题的过程：实际问题→→数值计算方法→程序设计→上机计算求出结果 在自然科学研究和工程技术中有许多问题可归结为求解方程组的问题,方程组求解是科学计算中最常遇到的问题如在应力分析、电路分析、分子结构、测量学中都会遇到解方程组问题在很多广泛应用的数学问题的数值方法中，如三次样条、最小二乘法、微分方程边值问题的差分法与有限元法也都涉及到求解方程组 1、插值法 许多实际问题都用y=f(x)来表示，有的函数虽然有解析式，但由于计算复杂实用不方便，为了找一个既能反映函数的特性又便于计算的函数，我们利用插值法可以得到这个简单函数，插值法包括拉格朗日插值，插值，Hermite插值 拉格朗日插值是n次多项式插值，其成功地用构造插值基函数的 方法解决了求n次多项式插值函数问题。插值也是n次多项式插值，它提出另一种构造插值多项式的方法，与Lagrange插值相比，具有承袭性和易于变动节点的特点。Hermite插值是利用未知函数f(x)在插值节点上的函数值及导数值来构造插值多项式的,起其提法为:给定n+1个互异的节点x0,x1,……,xn上的函数值和导数值通过一系列的加减消元运算，也就是代数中的加减消去法，A对角线以下的元素化为零，将化为上三角矩阵；然后，再逐一回代求解出x向量。现举例说明如下： 第一步：将/3使的系数化为1再将、式中的系数都化为零，即由2×(1)(1)得 由(3)-4×得将除以2/3，使系数化为1得再将式中系数化为零，由-(-14/3)*(2)(2)得 将除以18/3，使系数化为1，得 经消元后，得到三角代数方程组 第二步：回代过程 由得将代入得将代入得所以，本题解为第步： 即原方程组被等价转化成为上三角方程组，然后，逐步回代得原方程组的解即可。 3、解非线性方程组的方法 解非线性方程组的方法包括牛顿法，迭代收敛的加速解法，弦解法和抛物线法等 牛顿法实质是一种非线性方程逐步归结为线性方程来求解的，牛顿迭代法原理如下： 设已知方程f(x)=0的近似根X0则在X0附近f(x)可用一阶泰勒多项式近似代替.因此, 方程f(x)=0可近似地表示为P(x)=0.用X1表示P(x)=0的根,它与的根差异不大. 设,由于满足解得 重复这一过程,得到迭代格式： 这就是著名的牛顿迭代公式,它相应的不动点方程为： 三、数值分析的一些应用 为了更好地说明说明数值分析在物理领域的应用,以及如何使用数值分析进行求解,本文将简单介绍用直接法和逼近法求解一些实际问题. 1、在电网中的应用 确定下图电网中的回路电流。 分析：在回路1中,电流I1流过三个电阻,且电压降为；在回路2中的电流也流经回路1的一部分,即从D到A的分支,对应的电压降为 6I2伏特．然而,回路1中电流在DA段的方向与回路2中选定的方向相反,因此,回路1中所有电压降的代数和为．由于回路1中的电压为+60伏特,由基尔霍夫电压定律,可得回路1的方程为 , 同理,可得回路2的方程为 , 其中, -6I1是回路1中流经DA分支的电流(因为电流与回路2中的电流方向相反,所以电压为负)；12I2是回路2中所有的电阻乘上回路电流的和； -2I3是回路3中流经CB分支上2欧姆电阻的电流,方向与回路2中该段的电流方向相反。回路3的方程为: 注意,在CB分支上10伏特的电池被当作是回路2和回路3中的一部分,但是由于回路3中电流方向,电池在回路3中为-10伏特．出于同样的道理,40伏特的电池也应取负值.综合上述讨论,上述电网的回路电流满足下列线性方程组: 写成矩阵形式为 (*) 对增广矩阵进行行变换,得 从而解得I1=3安培, I2=1安培, I3=-8安培．I3取负值说明回路3中的实际电流与图中显示的电流方向相反. 在方程组(*)中,如果将其系数矩阵记为R,右端列向量记为u,i =(I1,I2,I3)T, 则可得到以矩阵形式表示的欧姆定律:. 2.在提拉单晶硅中的应用 在提拉单晶硅的过程中，假定晶体导电，熔化物和晶体的