资源描述：
《数值计算论文-高斯消去法在电路分析中的应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、高斯消去法在电路分析中的应用—利用计算机求解一些复杂电路的参数电子信息工程学院自动化一班100401102摘要求解线性代数方程组的数值方法有很多，但归纳起来，可分为两类：一种是直接法，另一种为迭代法。直接法在不计运算过程的舍入误差时，经过有限次运算，可得到方程组的精确解。而本文将介绍这种方法之一的高斯消去法在求解一些电路分析问题时的应用，以及利用计算机更为方便的解出其参数。关键词：高斯消去法；电路分析；计算机13Gaussianeliminationcircuitanalysis-Usingthecomputer
2、tosolvesomecomplicatedcircuitAbstractManynumericalmethodsforsolvinglinearalgebraicequations,butsummedup,canbedividedintotwocategories:oneisthedirectmethod,anotheriterativemethod.DirectlyroundingerrorinexcludingOperationop,afterafinitenumberoftimes,theexactsol
3、utionoftheequationscanbeobtained.ThisarticlewillintroducethismethodisoneoftheGaussianeliminationinsolvingtheproblemofcircuitanalysis,aswellasusingthecomputermoreconvenienttosolveitsparametersKeywords:Gaussianelimination,circuitanalysis,computerprogramming13目录
4、摘要1Abstract2引言41高斯消去法51.1高斯消去法背景及定义51.2高斯消去法简单应用72利用计算机编程求解高斯消去法82.1计算步骤82.2编程步骤83高斯列主元消去法子程序93电路分析中的运用10结论12参考文献1213引言在传统的电路分析方法中，一些回路网孔比较简单的电路，使用一般方法可以较方便的解得结果，但在一些复杂的电路中，要求解其参数，将耗费很大的计算量，在这时就可用高斯消去法来简化运算，必要时也可借助计算机来完成本文介绍系统分析法。系统分析法是以基尔霍夫电流定律(Kirchhoff'scu
5、rrentlaw简称KCL)、基尔霍夫电压定律(Kirchhoff'svoltagelaw简称KVL)以及支路元件电压、电流约束关系(voltage-currentrelation简称VCR)为理论基础,以所选的电路分析变量为方程变量,列写电路方程的一种电路分析方法;根据所选的电路变量的不同,系统分析法包括:支路分析法(分支路电流分析法与支路电压分析法),回路电流分析法(当所选的独立回路为网孔时,也称网孔分析法),节点电压分析法。131高斯消去法1.1高斯消去法背景及定义高斯消去法是一个古老的求解线性方程组的直接
6、法，由它改进变形又得到了一些其他的方法，如选主元消去法，三角分解法等。仍然是目前计算机上解低阶稠密矩阵方程组的常用有效方法。高斯消去法是建立在逐次消去未知元的基础上，它的基本做法是把方程组（1）用逐次消去未知元的方法化为其等价的三角形方程组（1）（2）这个过程面的回代，这个过程称为回代过程，这就是高斯消去法。由线性代数知识知道，无论是消元过程还是回代过程都不需要对未知元或方程做真正的运算，只要对他们的系数和右端项做运算就足够了，换句话说，只要把方程组的系数和右端项从方程周分离出来，那么，消去法完全可通过增广矩阵的
7、行的初等变换来实现，其消去过程的基本步骤示意如下：13然后自下而上进行回代，按上述的步骤做，即便与编程又节约计算工作量。具体的计算方法常分为顺序高斯消去法、列主元素高斯消去法和全主元素高斯消去法。131.2高斯消去法简单应用例题：求解方程组的根。解：132利用计算机编程求解高斯消去法2.1计算步骤对于n元线性方程组：第一步：对方程组确定il,使为第一主元，交换第一个和i个方程，利用第一个方程将后n-1个方程中的消去。第二步：在第二列中寻找，重复以上过程，消去。n-1步后原方程组变为上三角形方程组，利用回代可求得结
8、果。2.2编程步骤对于增广矩阵：（1）对k=1,2,3,4,,n-1[1]选主元，确定r,使[2]交换中的r,k两行[3]对i=k+1,k+2,,n,计算[4]132.3高斯列主元消去法子程序subroutinegauss(N,A,B)dimensionA(N,N),B(N)do60K=1,N-1P=0do30I=K,NIf(abs(A(I,K)).LE.abs(P))g
