利用直接发求解线性方程组及其应用.docVIP

利用直接法解线性方程组及其应用 摘要: 在科学研究和工程技术中有许多问题可归结为求解线性代数方程组的问题和常微分方程的定解问题。本文主要讨论了解线性方程组的各种直接解法，常微分方程的边值问题及其应用，最后介绍了病态方程和及其求解方法。 关键字: 线性方程组 直接法 应用举例 一、引言 随着科学技术的发展，提出了大量复杂的数值计算问题，在建立电子计算机成为数值计算的主要工具以后，它以数字计算机求解数学问题的理论和方法为研究对象。有可靠的理论分析要有数值实验主要内容插值法，函数逼近，曲线拟和，数值积分，数值微分，解线性方程组的直接方法，解线性方程组的迭代法，非线性方程求根，常微分方程的数值解法。运用数值分析解决问题的过程：实际问题→→数值计算方法→程序设计→上机计算求出结果 在自然科学研究和工程技术中有许多问题可归结为求解线性代数方程组的问题,线性方程组求解是科学计算中最常遇到的问题如在应力分析、电路分析、分子结构、测量学中都会遇到解线性方程组问题在很多广泛应用的数学问题的数值方法中，如三次样条、最小二乘法、微分方程边值问题的差分法与有限元法也都涉及到求解线性方程组Gauss消元法 通过一系列的加减消元运算，也就是代数中的加减消去法，A对角线以下的元素化为零，将化为上三角矩阵；然后，再逐一回代求解出x向量。现举例说明如下：消元过程 第一步：将/3使的系数化为1再将、式中的系数都化为零，即由2×(1)(1)得 由(3)-4×得 第二步：将除以2/3，使系数化为1，得再将式中系数化为零，由-(-14/3)*(2)(2) ，得 第三步：将除以18/3，使系数化为1，得 经消元后，得到三角代数方程组 1.2 回代过程 由得将代入得将代入得所以，本题解为第步：第步：第步： 即原方程组被等价转化成为上三角方程组，然后，逐步回代得原方程组的解即可。 1.4高斯消元的公式 综合以上讨论，不难看出，高斯消元法解方程组的公式为 第一步，消元 令 aij(1) = aij , （i,j=1,2,3,…,n） bi(1) =bi , （i=1,2,3,…,n） 对k=1到n-1，若akk(k)≠0，进行 lik = aik (k) / akk(k) , （i=k+1,k+2,…,n） aij(k+1) = aij(k) - lik * akj(k), （i,j= k+1,k+2,…,n） bi(k+1) = bi(k) - lik * bk(k), （i= k+1,k+2,…,n） 第二步，回代 若ann(n) ≠ 0 xn = bn(n) / ann(n) xi = (bi(i) – sgm(aij(i) * xj )/- aii(i) ,（i = n-1,n-2,…,1）,( j = i+1,i+2,…,n ) 1.5 解的判断 设方程组的增广矩阵记为，则经过初等行变换可化为如下的阶梯形矩阵(必要是可重新排列未知量的顺序)： 其中cii(0(i=1，2，…，r)．于是可知： 当dr+1=0，且r=n时，原方程组有唯一解． 当dr+1=0，且rn时，原方程组有无穷多解． 当dr+1(0，原方程组无解． 2、LU分解法求解线性代数方程组除了高斯消元法外，还常用LU分解法（三角形分解法）。LU分解法的优点是当方程组左端系数矩阵不变，仅仅是方程组右端列向量改变，即外加激励信号变化时，能够方便地求解方程组设n阶线性方程组Ax=b 假设能将方程组左端系数矩阵A，分解成两个三角阵的乘积，即A=LU 式中，L为主对角线以上的元素均为零的下三角矩阵 且主对角线元素均为1的上三角矩阵；U为主对角线以下的元素均为零 所以有，LUx=b 令 Ux=y则 Ly=b…,n ai1 = li1u11 , i=1,2,…,n 推出 u1j = a1j, j=1,2,…,n li1 = ai1/u11, i=1,2,…,n 这样就定出了U的第一行元素和L的第一列元素。 设已定出了U的前k-1行和L的前k-1列，现在确定U的第k行和L的第k列。由矩阵乘法： 当rk时，lkr=0, 且lkk=1，因为 所以， 同理可推出计算L的第k列的公式： 因此得到如下算法——杜利特（Doolittle）算法: 将矩阵分解为A=LU,对k=1,2,…,n；j=k,k+1,…n； i=k,k+1,…n； 公式1 解Ly=b 解Ux=y 对大规模稀疏问题，如果能够通过调整方程及未知量的顺序使得方程组的系数矩阵成带状结构，则对系数矩阵使用通常的LU分解，可以保障单位下三角矩阵L及上三角矩阵U仍为带状结构. 3. 直接法解线性方程组的应用 例1：确定下图电网中的回路电流。 解 在回路1中