线性方程组的求解方法与应用

湖北民族学院理学院2016届本科毕业论文(设计)线性方程组的求解 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 及应用学生姓名：付世辉学号：021240712专业： 数学 数学高考答题卡模板高考数学答题卡模板三年级数学混合运算测试卷数学作业设计案例新人教版八年级上数学教学计划 与应用数学指导老师：刘先平答辩时间：2016年5月22日装订时间：2016年5月28日AGraduationThesis(Project)SubmittedtoSchoolofScience,HubeiUniversityforNationalitiesInPartialFulfillmentoftheRequiringforBSDegreeIntheYearof2016ThecalculationmethodandapplicationofthesystemoflinearequationsStudentName:FuShihuiStudentNo.:021240712Specialty:MathematicsAndAppliedMathematicsSupervisor:LiuXianpingDateofThesisDefense:2016.05.22DateofBookbinding:2016.05.28摘要线性方程组在数学领域中的应用非常广泛，是线性代数的主要内容之一.矩阵及其基本理论是学习线性代数的一种基本工具，矩阵的初等变换则是线性方程组求解的工具.线性方程组常用的求解方法有一般消元法、克拉默法则、LU分解法等一系列方法，根据问题的不同，我们在求解的过程中选择的方法也就多种多样.这些方法可以很好地解决线性方程组的求解问题，在求解过程中，向量和矩阵起着一个不可或缺的作用.在线性方程组的应用方面，除了跟数学理论知识有着密不可分的联系，还和我们的实际生活联系的极其紧密.关键词：线性方程组，矩阵，初等变换，克拉默法则，LU分解法IAbstractLinearequationsarewidelyusedinthefieldofmathematicsandtheyarethemaincontentsoflinearalgebra.TheMatrixanditsbasictheoryarebasictoolforlearninglinearalgebra,theelementarytransformationofthematrixisthetoolofthesolutionofthelinearequations,thecommonlyusedmethodsofsolvinglinearequationshavethegeneraleliminationmethod,Gramer,theLUdecompositionmethodandsoon,isaccordingtotheproblem,wechooseonefromavarietyofmethodintheprocessofsolving.Thesemethodscansolvetheproblemsolvinglinearequations,vectorsandmatricesplayintegralrolesintheprocessofsolving.Intheapplicationoflinearequations,ithasnotonlyacloselinktotheknowledgeofmathematicaltheory,butalsoveryclosetoourreallife.Keywords:linearequations,matrix,elementarytransformation,Gramer,theLUdecompositionmethodII 目录 工贸企业有限空间作业目录特种设备作业人员作业种类与目录特种设备作业人员目录1类医疗器械目录高值医用耗材参考目录 2.3.3一般线性方程组的解的结构......................................................................59III33455 总结 初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf 与展望....................................................................................................................1601绪言本课题阐述与线性方程组有关的求解方法及其广泛应用，线性方程组是贯穿大学线性代数的一个重要工具,它是贯穿向量、矩阵的桥梁.国内外许多著名的数学学家对线性方程组也做了不少的研究，并且取得了显著的科研成果.1.1课题背景线性代数是大学数学代数学科的一个重要分支，早在中世纪就开始了对线性代数的研究.而方程组理论则是代数学发展的一个重要方向，也是代数学的核心内容之一.关于线性方程组的求解，在中国历史上很早以前就进行了研究.对线性方程组的研究，最早记录在公元初《九章算术》中，远远早于欧洲.大学所学高等代数中求解线性方程组是用一些比较基本的方法，在解决一些比较复杂的问题上有一定的局限性.本文主要运用了一般消元法、克拉默法则、分解法等解法.针对不同的问题，我LU们解决这些问题所选择的方法也不尽相同，这些相关的问题都需要我们去解决.在现代科学计算中的许多问题，例如生活中的营养搭配问题、电路问题均与线性方程组的求解有关.11.2课题研究的目的和意义课题研究的意义：(1)给出线性方程组的一些求解方法，使读者对线性方程组有更深层次的了解；(2)线性方程组的应用与我们的生活息息相关，特别是与我们的饮食健康、经济平衡联系的比较紧密，我们可用它解决生活中的一些基本问题.1.3国内外概况对于线性方程组求解方法的研究，国内外许多著名的数学学家对此作出了不少的贡献.随着科学技术的进步，数学已经渗入到各学科之中，甚至渗入到我们的日常生活中.而我们所学线性方程组的理论知识，则是源自许多著名的国内外学家的著作.在实际生活中，我们需要确定所需目标，并对此目标作出一系列的决策，这些决策中的关键要素就是我们重点研究的对象.22预备知识2.1线性方程组的定义形如.qypypyp,qypypyp,qypypypsnsn22s11s2nn22221211nn1212111（2.1）的方程组的叫做线性方程组.其中，代 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 个未知量，是方程的个数，nyyy21,ns称为方程组的系数，称为常数项；系数njsipij2,1,2,1sjqj2,1ijp的第一个指标表示它在第个方程，第二个指标表示它是的系数.iijjx常记为.QPy其中，snssnnpppppppppP212222111211sqqqQ21当线性方程组的右端全为零时，该线性方程组就称为齐次线性方程3组；当线性方程组的右端不全为零时，该线性方程组就称为非齐次线性方程组.2.2线性方程组有解判别定理定理2.2.1线性方程组.qypypyp,qypypyp,qypypypsnsn22s11s2nn22221211nn1212111有解的充分必要条件是它的系数矩阵snssnnpppppppppP212222111211与它的增广矩阵ssnssnnqqqpppppppppP21212222111211有相同的秩.证明：不妨先引入向量ssnnnnssqqqQpppPpppPpppP2121222122121111,,,,4（2.2）于是定理中的线性方程组就可以写成QyPyPyPnn2211（2.3）很明显，该线性方程组有解得充要条件是向量可以由向量组Q线性表出，定理证明过程如下：nPPP,,,21必要性假设该线性方程组有解，那么就是说向量是向量组Q的线性组合，从而可以得到与向量组等价，nPPP,,,21nPPP,,,21QPPPn,,,,21又已知等价的向量组有相同的秩，故这两个向量组有相同的秩.并且这两个向量组分别是系数矩阵与增广矩阵的列向量组.所以，系数矩阵PP和增广矩阵有相同的秩.PP充分性假设系数矩阵与增广矩阵有相同的秩，那么它们的列向PP量组与有相同的秩，不妨让它们的秩都等于.nPPP,,,21QPPPn,,,,21r中的极大线性无关组是由个向量组成的，可以假设nPPP,,,21rrPPP,,,21是它的一个极大线性无关组，故向量可以由线性表出，再加QPr,,,21PP上等个线性无关的向量，可以知道向量可以经nrrPPP,,,21rnQ线性表出.所以，该线性方程组有解.nPPP,,,21证毕2.3线性方程组解的结构在上一节，我们解决了线性方程组有解的判别条件之后，我们这一节还需要探讨一下线性方程组解的结构.52.3.1齐次线性方程组的性质对于齐次线性方程组（2.4）.0ypypyp,0ypypyp,0ypypypnsn22s11snn2222121nn1212111它的解所构成的集合有以下两个重要性质：性质1:两个解的和还是方程组的解.假设与是该线性方程组的两个解.即将这两n21a,,a,an21b,,b,b个解代入到方程组中，每个方程都变成了恒等式.,01njjijap),,2,1(si,01njjijbp),,2,1(si将这两个解的和nnbababa,,,2211代入到该线性方程组中，可以得到.（2.5）n1jjjijbapn1jjijap000bpn1jjij)s,,2,1i(这就说明了两个解的和还是方程组的解.证毕6性质2：一个解的倍数还是方程组的解.设是该线性方程组的一个解，则有n21a,,a,a.（2.6）0apn1jjij)s,,2,1i(将这个解的倍代入到方程组中，就可以得到c.（2.7）00capccapn1jjijn1jjij)s,,2,1i(这就说明了一个解的倍数还是方程组的解.证毕2.3.2基础解系及其存在性（1）如果满足a.方程组的任意一个解可以表示成的线性组合；nPPP,,,21b.线性无关.nPPP,,,21那么这组解就称为齐次线性方程组的一个基础解系.n21P,,P,P（2）在齐次线性方程组有非零解的情况下，它有基础解系，并且基础解系所含有的解的个数等于，这里表示的是系数矩阵的秩.rnr2.3.3一般线性方程组的解的结构（1）线性方程组的两个解的差是它的导出组的解；7（2）线性方程组的一个解与它的导出组的一个解之和还是这个线性方程组的一个解；（3）在方程组有解的情况下，解是唯一的充要条件是它的导出组只有零解.83线性方程组的求解方法3.1一般消元法例3.2.1求消元法求解线性方程组.10y4y2y4,8y10y4y8,2y6y2y4321321321解：下面对这三个方程进行加减运算从而达到消元的目的.第二个方程减去第一个方程的2倍，得,4y2y832第三个方程减去第一个方程得,8y2y432将第一个方程、第四个方程、第五个方程综合起来就得到一个新的方程组.8y2y4,4y2y8,2y6y2y43232321再分别对这个方程组中的第二个和第三个方程进行加减运算，即第二个方程减去第三个方程的2倍就可以得到.122,824,2624332321yyyyyy可以解出9.6y,1y,9y321从而原方程组的解为（9，-1，-6）.一般消元法用来计算一些比较简单的线性方程组，是最简单最直接最有效的方法，它的基本思想就是将方程进行加减和代入运算，要转换成矩阵的行初等变换来求解.3.2克拉默法则3.2.1克拉默法则求解具备的条件利用克拉默法则求解线性方程组时需要具备两个条件：（1）线性方程组的方程个数必须与未知量的个数相等；（2）线性方程组的系数行列式不等于零.3.2.2克拉默法则设含有n个未知数的线性方程组的系数行列式（3.1）Qnnnnnnppppppppp2122221112110则该线性方程组有解，且只有唯一解，其解可以表示为QQy,,QQy,QQynn221110（3.2）其中是把系数行列式中第列的元素用常数项),,2,1(njQjQj代替后所得到的阶行列式，即nbbb,,21n（3.3）nnjnnjnnnjjnjjjppbppppbppppbppQ1,1,121,221,22111,111,111nj,2,1定理中包含着三个结论：方程组有解；)1(解是唯一的；)2(解由公式(3.2)给出.)3(下面来证明一下克拉默法则：证明：1.该线性方程组简写成.nibypnjijij2,1,1（3.4）首先需要证明（3.2）是（3.3）的解.将（3.2）代入到第i个方程，那么左端就为.（3.5）njjijnjjijQpQQQp11111由于.nssjsnjnjjjAbPbPbPbQ12211故有n1jn1ssjsijn1jjijPbpQ1QpQ1=n1ssn1jsjijb)Pp(Q1=.iibQbQ1这就与第i个方程的右边是一致的，故（3.2）是该线性方程组的解.2.假设是该线性方程组的一个解，那么就有那个恒等式)cc,c(n21.n,,2,1i,bcpn1jijij（3.6）下面证明，不妨取矩阵中第k列元素的代数余子式QQckk，再用他们乘上面这n个恒等式，就可以得到nkk2k1P,,P,P.n1jikijijikn,,2,1i,PbcpP将它们加起来，就可得.n1in1jn1iikijijikPbcpP（3.7）再来看（3.7）的左边，12n1in1jn1in1jjikijjijikcPpcpP=.n1jjn1iikijc)Pp(已知n1jikij.kj,0,kj,QPp当当故有.QccPpkjn1jn1iikij就可得.n,,2,1k,QQckk即方程组的一个解为，它必为)c,,c,c(n21QQ,,QQ,QQn21故方程组最多有一组解.证毕例3.2用克拉默法则解线性方程组：.4yy3yy3,3y2yyy3,5y2y3y3y3,6y2y3yy24321432143214321解：该线性方程组的系数行列式为13，0701313211323332312Q，0701314211323352316Q1，0701343213323532362Q2，0701413231325332612Q3，0704313311353336312Q4从而该线性方程组的唯一解为.（3.8）1,1,1,144332211QQyQQyQQyQQy克拉默法则只适用于系数行列式不等于0的线性方程组，那么对于系数行列式为0的情况就不能用克拉默法则解题了，克拉默法则主要在于理论上的应用.定理3.2.1如果齐次线性方程组14（3.9）.0ypypyp,0ypypyp,0ypypypnnnn2nn1nnn2222221nn1212111的系数矩阵的行列式，那么它只有零解.即，如果该齐次线性方程组有非零解，那么必然有.0Q下面来证明一下定理3.2.1：证：对该齐次线性方程组应用克拉默法则，因为行列式中有一列全jQ为零，则有（3.10）.0QQQn21故它的唯一解为.0,0,0从而定理3.2.1得证.例3.2.2求在什么条件下，方程组.0xx1,0x1x2121有非零解.解:根据定理3.2.2，若该齐次线性方程组有非零解，那么系数矩阵的行列式,011112故.1,02115可以验证当时，该齐次线性方程组确有非零解.1021或克拉默法则在计算一些基本的线性方程组是可行的，但是对于个未n知量n个方程的线性方程组就必须得计算个级行列式，相对而言计算量非1nn常大.3.3分解法LU分解法，又称三角形分解法，是求解线性方程组的重要方法之一.LU当方程组左边的系数矩阵不变，仅仅是方程组右边列向量发生改变，能够很好地求解方程组.设n阶线性方程组为.QPy若可以将方程组左边的系数矩阵分解成两个三角阵的乘积，即可P以写成，其中，为主对角线以上的元素均为零并且主对角线元素LUPL均为1的下三角矩阵，为主对角线以下的元素均为零的上三角矩阵.UnnnnnnnnnnnnUUUUUULLLPPPPPPPPP222112112121212222111211111故有，QLUy16令，xUy则有.QLx由，由矩阵的乘法公式：LUP,n,,2,1j,UPj1j1,n,,2,1i,ULP111i1i推出,n,2,1j,PUj1j1,n,,2,1i,UPL111i1i这样就可以确定出的第一行元素和的第一列元素.以此类推，设UL已经确定了的第行元素和的第列元素，由矩阵乘法：U1kL1k.n1rrjkrkjULP当时，有，因为kr1,0rrkrLL.nrrjkrkjkjULUP1所以,nrrjkrkjkjULPU1n,,1k,kj同理可以求出的第列的公式：Lk17,kknrrjkrikikUULPL1nkki,,1,所以可以得到如下算法——杜利特（Doolittle）算法：（1）将矩阵分解为n,2,1k,LUP对nrrjkrkjkjULPU1,n,,1k,kjkknrrjkrikikUULPL1,n,,1k,ki.1Lkk（2）解QLx.n,,2,1k,xLQx1k1rrkrkk（3）解xUy.1,,1n,nk,UxUxykkn1krrkrkkLU分解法的证明过程如下：证明:对，当时上述结论显然成立.nnP1n假设当上述结论成立，现证明时结论成立.1knkn对进行分块得，令，易得Ann1npCQPP1CPEL11n1n1其中.已知对成立,其中.QCPp0QPLP11nnn1n1QQCPP11nnn1kP1kP22UL18为单位下三角矩阵为非奇异上三角矩阵，于是2L2U.qLU100LLq0QULLP12221221令,，显然易见为单位下三角矩阵，为非L10021LLUqLU122LU奇异上三角矩阵.例3.3.1求解方程组.4722239yyyy40115618962569262424321解：由上面的分解法第一步的公式可得出，1233121121L1633216242U于是可以化为两个方程组，472223912331211214321xxxx.432143211633216242xxxxyyyy利用第二步和第三步的公式可得出该线性方程组的解为19.TTyx1,3,2,5,1,3,5,9204线性方程组的应用线性方程组的应用包括理论和实际方面的应用，下面我们来介绍一下它的具体应用.线性方程组在数学理论的应用还是很广泛的，特别在解析几何、代数上有重要的应用.我们可以用齐次线性方程组来解决初等数学的一些基本问题，应用起来可以适当避免一些繁琐的过程.4.1线性方程组在几何学中的应用定理4.1.1任意三点不在同一直线上的不同四点iQiiba,4,3,2,1i共圆的充分必要条件是.01111442424332323222222112121babababababababa证明：先证必要性：我们假设在圆iiibaQ,4,3,2,1i022FEbDabaA0A上，那么我们就可以将这四个点带入到圆方程中去,就可以得到.0FEbDabaA,0FEbDabaA,0FEbDabaA,0FEbDabaA442424332323222222112121我们已知，故此方程组是关于的齐次线性方程组，它必0AFEDA,,,21有非零解，从而可以知道它的系数行列式等于零，即.01111442424332323222222112121babababababababa从而必要性得证.再证充分性：如果，01111442424332323222222112121babababababababa那么关于的齐次线性方程组FEDA,,,.0FEbDabaA,0FEbDabaA,0FEbDabaA,0FEbDabaA442424332323222222112121就会存在非零解.我们可以假设是该齐次线性方程组的1111,,,FEDA一组非零解.又已知四点中任意三点都不共线，那么可以假设，否则01A的话不同时为零,从而导致三点共线.即可得出这四个点都在圆111,,FED上，充分性得证.0111221FbEaDbaA01A从而该定理得证.224.2线性方程组在高次方程理论中的应用定理4.2.1若为互不相等的四个实数，则四个方程dcba,,,.0cbxaxdx,0baxdxcx,0adxcxbx,0dcxbxax23232323没有公共实根.证明：我们可以用反证法进行证明.假设这四个方程有公共实根，那么关于的齐次线性方程组0xdcba,,,就会有一组非零解，即.0111130020203000203002030xxxxxxxxxxxx可以得到，将带入到原方程组中，就可以得到，这与10x10xcb已知互不相等矛盾，从而可以得出这四个方程没有公共实根.dcba,,,4.3线性方程组在化学中的应用线性方程组在化学中的应用也比较广泛，特别是配平化学方程式，根据守恒规律，可得出相应的方程，从而求解出方程组的解，就可以配出化学方程组的系数了.23例4.2.1在高温条件下，一氧化碳可以还原四氧化三铁生成单质铁和二氧化碳，请配平该化学方程式243COFeOFeCO解：为了配平该化学方程式，我们不妨假设反应物和生成物的量分别为，根据守恒定律，各原子的数目在发生化学反应前后保持不变，4321,,,xxxx则根据铁原子守恒，有,233xx根据碳原子守恒，有,41xx为了平衡氧原子，有,42124xxx将这三个式子综合起来并进行移项，就可以得到.0x2x4x,0xx,0x3x4214123由于方程的个数少于未知量的个数，我们可以知道该线性方程组必定有非零解，且每一个分量必须为正数解，简化其增广矩阵就有.010400013000041020410100100130B24我们将取作自由未知量，就有2x.x4x,x3x,x4x242321特别地，取,1x2此时可求出.4x,3x,4x431此时化学方程式有以下形式：243434COFeOFeCO255总结与展望本文系统介绍了线性方程组相关的知识，主要有两大模块：求解方法和应用.线性方程组求解比较常见的方法有一般消元法、克拉默法则、LU分解法等方法，对这三种方法我们分别进行了举例以便于理解，在其应用方面，我们分为了理论应用和实际应用，并进行了举例.关于线性方程组的求解方法，本文只列举了三种，还有其它的一些方法并没有引用进来，有待进一步了解.关于线性方程组的应用，还有其它方面的应用，比如电路问题、营养搭配问题、维他命的配方等等都可以列举相关的例子.26致谢在论文完成之际，我要特别感谢刘先平老师给予我无私的指导，然后根据老师的指导进行修改.在论文写作的这段时间里，我碰到了各种各样的障碍，包括论文选题、论文写作等方面，在论文写作这方面我还是显得经验不足，有待进一步加强.在这里，我还要感谢各位老师和同学们给我无私的帮助，正是各位老师和同学们兢兢业业的指导个帮助之下，我才能够顺利完成本论文.由于本人学识是有限的，加之时间非常仓促，文中不免有错误和待改进之处，真诚欢迎各位师长、同行提出宝贵意见.27参考文献[1]张科.两类线性方程组的预处理技术及数值求解方法[D].上海大学,2014.[2]郑汉垣.大规模稀疏线性方程组求解的并行GaBP算法研究[D].上海大学,2014.[3]周挺辉,赵文恺,严正,徐得超,江涵.基于图形处理器的电力系统稀疏线性方程组求解方法[J].电力系统自动化,2015,v.39;No.55202:74-80.[4]史文谱,刘迎曦,褚京莲,郭淑红.求解线性方程组的一种新方法[J].计算力学学报,2003,06:715-720.[5]严辉银.求解二乘二分块实线性方程组的块分裂预处理方法[D].兰州大学,2015.[6]邬贵明.FPGA矩阵计算并行算法与结构[D].国防科学技术大学,2011.[7]邢亚斌,卢鹏,李冬媛.线性方程组的求解方法[J].才智,2011,17:103-104.[8]曲婧.一种基于神经网络求解线性方程组的新方法[J].山西电子技术,2010,No.14902:9-10.28[9]徐晓飞,曹祥玉,姚旭,陈盼.一种基于DoolittleLU分解的线性方程组并行求解方法[J].电子与信息学报,2010,v.3208:2019-2022.[10]马成业,杨胜良,黎锁平.求解病态线性方程组的一个正则化方法[J].甘肃科学学报,2010,v.22;No.9404:33-35.[11]苏焱.浮体水动力学边界元分析中线性方程组迭代求解方法研究[D].哈尔滨工程大学,2011.[12]李扬.求解线性方程组及不等式组的ABS方法[D].辽宁师范大学,2007.[13]李文伟.一类线性方程组和矩阵方程的数值求解方法的研究[D].南昌大学,2014.[14]李超.求解多右端线性方程组的块种子投影方法[D].南京航空航天大学,2013.[15]李晓爱,陈玉花,张耘,王新苹.求解大型稀疏线性方程组的Krylov子空间方法的发展[J].科技导报,2013,v.31;No.40111:68-73.[16]王纪林,周钢.求解线性方程组的初参数方法[J].上海交通大学学报,1994,03:100-105.[17]吴丹红.一种加速大规模线性方程组求解的并行方法[J].机电工程,2008,No.15804:58-59+66.29[18]李欣.求解非对称线性方程组的总体拟极小向后扰动方法[J].南京大学学报(自然科学版),2005,04:350-355.[19]张少杰,杨陈东.求解对称正定线性方程组的正交基变换方法[J].河南科学,2016,v.34;No.20803:310-314.[20]李欣.求解线性方程组的总体（拟）极小向后扰动方法[D].南京航空航天大学,2004.[21]王萼芳，石生明.高等代数[M].北京:高等教育出版社,2003.[22]李刚，张显德.高等代数辅导及习题精解[M].吉林：延边大学出版社，2012.[23]闫晓红.高等代数全程导读及习题全解[M].北京：中国时代经济出版社，2006.[24]王向东等.高等代数常用方法[M].北京：科学出版社，1989.[25]StoerJ,BulirschR.IntroductiontoNumericalAnalysis[M].NewYork:Springer-Verlag,1980.30独创性声明本人声明所呈交的论文（设计）是我个人在导师指导下进行的研究工作及取得的研究成果。尽我所知，除文中已经标明引用的内容外，本论文（设计）不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出贡献的个人和集体，均已在文中以明确方式标明.本人完全意识到，本声明的法律结果由本人承担。31学位论文作者签名：日期：年月日学位论文（设计）版权使用授权书本论文（设计）作者完全了解学校有关保留、使用学位论文的规定，即：学校有权保留并向国家有关部门或机构送交论文的复印件和电子版，允32许论文（设计）被查阅和借阅。本人授权湖北民族学院可以将本论文（设计）的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本论文（设计）。保密□，在_____年解密后适用本授权书。不保密□.（请在以上方框内打“√”）学位论文作者签名：指导教师签名：日期：年月日日期：年月日本论文属于