浅谈线性方程组的求解及其应用

浅谈线性方程组的求解及其应用

【摘要】 线性代数是代数学的一个重要组成部分,广泛应用于现代科学的许多

分支。其核心问题之一就是线性方程组的求解问题。本文先简要介绍了线性方程

组求解的历史,然后给出线性方程组解的结构。重点介绍了解线性方程组的几种

方法:消元法,克拉默法则求解线性方程组的方法。最后介绍了如何利用Matlab

常用电脑软件解线性方程。

关键词:线性方程组 克拉默法则 Matlab

1. 线性方程组求解的历史

线性方程组的解法,早在中国古代的数学著作《九章算术》方程章中已作

了比较完整的论述。其中所述方法实质上相当于现代的对方程组的增广矩阵施行

初等行变换从而消去未知量的方法,即高斯消元法。在西方,线性方程组的研究

是在17世纪后期由莱布尼茨开创的。他曾研究含两个未知量的三个线性方程组组

成的方程组。麦克劳林在18世纪上半叶研究了具有二、三、四个未知量的线性方

程组,得到了现在称为克莱姆法则的结果。克莱姆不久也发表了这个法则。18

世纪下半叶,法国数学家贝祖对线性方程组理论进行了一系列研究,证明了一元

齐次线性方程组有非零解的条件是系数行列式等于零。

2. 线性方程组解的结构

n元线性方程组的一个解(c1,c2,……cn)是一个,维向量,当方程组有无穷多

个解时,需要研究这些解向量之间的关系,以便更透彻地把握住它们。

关于齐次线性方程组的解的结构有以下结论:

1) 定义1齐次线性方程组的一组解η1,η2……ηt称为该方程组的一个基础解

系,如果

a) 该方程组的任一解都能表成η1,η2……ηt的线性组合。

b) η1η2……ηt线性无关。

2) 齐次线性方程组的两个解的和还是解,一个解的倍数还是解。

3) 齐次线性方程组有非零解时必定存在基础解系,并且一个基础解系里有n-r

个解,其中n是未知量的个数,r是系数矩阵的秩。如果系数在数域P中的齐

次线性方程组

的一个基础解系是: η1,η2……ηn-r,则①的全部解为

k1η1+k2η2+……+kn-rηn-r

其中k1,k2,……kn-r取遍取遍数域P中全部数。

3. 线性方程组的求解方法

3.1消元法

解线性方程组的最基本最有效的方法是消元法。它的做法是:先把线性方程

组的增广矩阵经过矩阵的初等行变换化成阶梯形,然后去解相应的阶梯形方程

组。或者把线性方程组的增广矩阵经过初等行变换化成行简化阶梯形,从而可立

即写出方程组的解。

消元法是求解低阶多元线性方程组的方法,此时线性方程组必须是适定方程

组,一般是用于二元一次或三元一次方程组,当未知元增多时,计算效率低甚至

无法求解。

3.2克拉默法则

当系数行列式小为零时,适定方程组有惟一解,其解为:

(i=1,2,……,n)

其中D是系数行列式,D是在系数行列式基础之上结合方程组右边常数形成的

新行列式)在此法则中,行列式的计算显得非常重要)利用行列式的性质计算行列

式最为有效,对于二、三阶行列式可以利用对角线法则计算。

克拉默法则克服了消元法计算效率低甚至无法计算多元一次方程组的缺点,

但是对于系数行列式等于零以及欠定或者超定方程组的情况,它是无能为力的)

事实上,当未知元数过多时,克拉默法则的计算效率就很低。

4. 线性方程组的解法在MATLAB中的实践

MATLA B语言是一种以矩阵运算为基础的计算语言,对于实现线性方程组的求

解非常方便、对一个四兀一次方程组的求解,可以用克拉默法则和逆阵乘积法来

实现,程序如下:

tic;

D=[1 1 1 1 ; 1 2 -1 4 ; 2 -3 -1 -5;3 1 2 11 ];

det(D)

b=[5 -2 -2 0];

D1=[5 1 1 1 ; -2 2-1 4 ; -2 -3 -1 -5 ; 0 1 2 11];

D2=[1 5 1 1 ; 1 -2-1 4 ; 2 -2 -1 -5 ; 3 0 2 11];

D3=[1 1 5 1 ; 1 2 -2 4 ; 2 -3 -2 -5 ; 3 1 0 11];

D4=[1 1 1 5 ; 1 2 -1 -2 ; 2 -3 -1 -2 ; 3 1 2 0];

X1=det(D1)/det(D);

X2==det(D2)/det(D);

X3=det(D3)/det(D);

X4=det(D4)/det(D);

X5=inv(D)*b;

toc

其中克拉默法则用行列式除法Xi=det(Di)/det(D)来实现;逆阵乘积法用

X=inv(D)*b来实现; det(D)是系数矩阵D的行列式运算; inv(D)是D的逆阵运算。 上例中,系数矩阵D不为零,可以用克拉默法则和逆阵乘积法来求解。当系数行

列式为零时,只能用初等变换来求解。对于初等变换,利用阶梯生成函数命令rref

也可以轻松地实现.

可见,MATLAB语言实现线性方程组的求解具有程序简单、直观的特点,同时

还具有计算效率高的优点,在实际计算巾摆脱了系数矩阵阶数未知元数等的限

制。

5. 总结

本文首先介绍了一些线性组解法的历史,然后主要介绍了线性方程组的几种

解法。消元法,克拉默法则求解线性方程组。线性方程组是线性代数中一个最重

要的内容,它除了本文介绍的几种解题的方法以及应用外,还有很多其它的解题

方法以及更多广泛的应用.它是数学以及其它理工科解题时必不可少的知识之

一。对线性方程组的解题方法以及它的应用,也会一直不断的研究下去。

参考文献

【1】 李尚志,线性代数.[M].高等教育出版社

【2】 李排昌.矩阵与解线性方程组[J].中国人民公安大学学报.2011年

【3】 李桂荣、韩忠月、梁超.线性方程组的简便解法[J].德州学院学报.

2007年8月.第23卷(第4期)

【4】 姬五胜.线性方程组解法及其MATLAB实践[J]天水师范学院学报.2009年3

月, 第29卷(第2期)

【5】 陈建莉.线性方程组解法新探[J]纺织高校基础科学学报.2008年6月.第21

卷(第2期)

浅谈线性方程组的求解及其应用

【摘要】 线性代数是代数学的一个重要组成部分,广泛应用于现代科学的许多

分支。其核心问题之一就是线性方程组的求解问题。本文先简要介绍了线性方程

组求解的历史,然后给出线性方程组解的结构。重点介绍了解线性方程组的几种

方法:消元法,克拉默法则求解线性方程组的方法。最后介绍了如何利用Matlab

常用电脑软件解线性方程。

关键词:线性方程组 克拉默法则 Matlab

1. 线性方程组求解的历史

线性方程组的解法,早在中国古代的数学著作《九章算术》方程章中已作

了比较完整的论述。其中所述方法实质上相当于现代的对方程组的增广矩阵施行

初等行变换从而消去未知量的方法,即高斯消元法。在西方,线性方程组的研究

是在17世纪后期由莱布尼茨开创的。他曾研究含两个未知量的三个线性方程组组

成的方程组。麦克劳林在18世纪上半叶研究了具有二、三、四个未知量的线性方

程组,得到了现在称为克莱姆法则的结果。克莱姆不久也发表了这个法则。18

世纪下半叶,法国数学家贝祖对线性方程组理论进行了一系列研究,证明了一元

齐次线性方程组有非零解的条件是系数行列式等于零。

2. 线性方程组解的结构

n元线性方程组的一个解(c1,c2,……cn)是一个,维向量,当方程组有无穷多

个解时,需要研究这些解向量之间的关系,以便更透彻地把握住它们。

关于齐次线性方程组的解的结构有以下结论:

1) 定义1齐次线性方程组的一组解η1,η2……ηt称为该方程组的一个基础解

系,如果

a) 该方程组的任一解都能表成η1,η2……ηt的线性组合。

b) η1η2……ηt线性无关。

2) 齐次线性方程组的两个解的和还是解,一个解的倍数还是解。

3) 齐次线性方程组有非零解时必定存在基础解系,并且一个基础解系里有n-r

个解,其中n是未知量的个数,r是系数矩阵的秩。如果系数在数域P中的齐

次线性方程组

的一个基础解系是: η1,η2……ηn-r,则①的全部解为

k1η1+k2η2+……+kn-rηn-r

其中k1,k2,……kn-r取遍取遍数域P中全部数。

3. 线性方程组的求解方法

3.1消元法

解线性方程组的最基本最有效的方法是消元法。它的做法是:先把线性方程

组的增广矩阵经过矩阵的初等行变换化成阶梯形,然后去解相应的阶梯形方程

组。或者把线性方程组的增广矩阵经过初等行变换化成行简化阶梯形,从而可立

即写出方程组的解。

消元法是求解低阶多元线性方程组的方法,此时线性方程组必须是适定方程

组,一般是用于二元一次或三元一次方程组,当未知元增多时,计算效率低甚至

无法求解。

3.2克拉默法则

当系数行列式小为零时,适定方程组有惟一解,其解为:

(i=1,2,……,n)

其中D是系数行列式,D是在系数行列式基础之上结合方程组右边常数形成的

新行列式)在此法则中,行列式的计算显得非常重要)利用行列式的性质计算行列

式最为有效,对于二、三阶行列式可以利用对角线法则计算。

克拉默法则克服了消元法计算效率低甚至无法计算多元一次方程组的缺点,

但是对于系数行列式等于零以及欠定或者超定方程组的情况,它是无能为力的)

事实上,当未知元数过多时,克拉默法则的计算效率就很低。

4. 线性方程组的解法在MATLAB中的实践

MATLA B语言是一种以矩阵运算为基础的计算语言,对于实现线性方程组的求

解非常方便、对一个四兀一次方程组的求解,可以用克拉默法则和逆阵乘积法来

实现,程序如下:

tic;

D=[1 1 1 1 ; 1 2 -1 4 ; 2 -3 -1 -5;3 1 2 11 ];

det(D)

b=[5 -2 -2 0];

D1=[5 1 1 1 ; -2 2-1 4 ; -2 -3 -1 -5 ; 0 1 2 11];

D2=[1 5 1 1 ; 1 -2-1 4 ; 2 -2 -1 -5 ; 3 0 2 11];

D3=[1 1 5 1 ; 1 2 -2 4 ; 2 -3 -2 -5 ; 3 1 0 11];

D4=[1 1 1 5 ; 1 2 -1 -2 ; 2 -3 -1 -2 ; 3 1 2 0];

X1=det(D1)/det(D);

X2==det(D2)/det(D);

X3=det(D3)/det(D);

X4=det(D4)/det(D);

X5=inv(D)*b;

toc

其中克拉默法则用行列式除法Xi=det(Di)/det(D)来实现;逆阵乘积法用

X=inv(D)*b来实现; det(D)是系数矩阵D的行列式运算; inv(D)是D的逆阵运算。 上例中,系数矩阵D不为零,可以用克拉默法则和逆阵乘积法来求解。当系数行

列式为零时,只能用初等变换来求解。对于初等变换,利用阶梯生成函数命令rref

也可以轻松地实现.

可见,MATLAB语言实现线性方程组的求解具有程序简单、直观的特点,同时

还具有计算效率高的优点,在实际计算巾摆脱了系数矩阵阶数未知元数等的限

制。

5. 总结

本文首先介绍了一些线性组解法的历史,然后主要介绍了线性方程组的几种

解法。消元法,克拉默法则求解线性方程组。线性方程组是线性代数中一个最重

要的内容,它除了本文介绍的几种解题的方法以及应用外,还有很多其它的解题

方法以及更多广泛的应用.它是数学以及其它理工科解题时必不可少的知识之

一。对线性方程组的解题方法以及它的应用,也会一直不断的研究下去。

参考文献

【1】 李尚志,线性代数.[M].高等教育出版社

【2】 李排昌.矩阵与解线性方程组[J].中国人民公安大学学报.2011年

【3】 李桂荣、韩忠月、梁超.线性方程组的简便解法[J].德州学院学报.

2007年8月.第23卷(第4期)

【4】 姬五胜.线性方程组解法及其MATLAB实践[J]天水师范学院学报.2009年3

月, 第29卷(第2期)

【5】 陈建莉.线性方程组解法新探[J]纺织高校基础科学学报.2008年6月.第21

卷(第2期)


相关内容

  • 关于常微分方程求解公式的注记
  • 2008年12月 阴山学刊 Dee.2008第22卷第4期 YINSHANACADEMICJOURNAL V01.22 No.4 关于常微分方程求解公式的注记 穆 勇 (长江大学信息与数学学院.湖北荆州434023) 摘 要:要对常微分方程中的一阶线性齐次方程以及一阶线性非齐次方程的求解公式进行一下 ...

  • 变系数二阶线性微分方程的求解
  • 第22卷第3期2009年6月 高等函授学报(自然科学版) JournaI of VoI.22No.3 2009 HigherCorrespondenceEducation(NaturalSciences) ・高职高专教学・ 变系数二阶线性微分方程的求解 范小勤1,李金洋2 (1.广州番禺职业技术学院 ...

  • 线性微分方程组初值问题的求解公式
  • 第! ! 卷第" 期宝鸡文理学院学报#自然科学版$ #$, &-. /0'&120&345&''676&18. 9:0/; (! ! ) &(" %&' (! **! >0. ! **! 年+月 ? ? ? ? ? ? ...

  • 考研数学一常考题型及重点汇总
  • 第二篇 高等数学 第一章 函数.极限.连续 思考的鱼点拨 "函数.极限.连续"这一部分的概念及运算是高等数学的基础,它们是每年必考的内容之一,数学一中本部分分数平均每年约占高等数学部分的10%. 本章的考题类型及知识点大致有: 1.求函数的表达式: (1)给出函数在某一区间上的表 ...

  • 考研数学三考试大纲 全国统考
  • 考研数学三大纲 编辑 考试科目 微积分.线性代数.概率论与数理统计 1试题结构 考试形式 试卷内容结构 试卷题型结构 2考试内容 微积分 线性代数 概率统计 试题结构编辑 考试形式 1.试卷满分及考试时间 试卷满分为150分,考试时间为180分钟. 2.答题方式 答题方式为闭卷.笔试. 试卷内容结构 ...

  • 数学三考研大纲
  • 考研数学三考试大纲 姓名:曹辉 手机号:[1**********] 考试科目:微积分.线性代数.概率论与数理统计 考试形式和试卷结构 一.试卷满分及考试时间 试卷满分为150分,考试时间为180分钟. 二.答题方式 答题方式为闭卷.笔试. 三.试卷内容结构 微积分 约56% 线性代数 约22% 概率 ...

  • 有限元法在结构力学领域的发展及应用
  • 科 研 训 练 论文题目 院 (系) 专业名称 学生姓名 学生学号 指导教师 年 月 日 有限元在结构力学领域的发展和应用 摘要:有限元法是伴随着电子计算机技术的进步而发展起来的一种新兴数值分析方法,是力学.应用数学与现代计算技术相结合的产物.有限元法是一种高效能.常用的计算方法.有限元法最早应用于 ...

  • 华南科技大学硕士研究生入学考试
  • 华南科技大学硕士研究生入学考试<数学>(含高等数学.线性代数) 考试大纲 一.函数.极限.连续 考试内容 函数的概念及表示法 函数的有界性.单调性.周期性和奇偶性 复合函数.反函数.分段函数和隐函数 基本初等函数的性质及其图形 初等函数 简单应用问题的函数关系的建立. 数列极限与函数极限 ...

  • 类二阶常微分方程组特解形式的探讨
  • 一类二阶常微分方程组特解形式的探讨 作者: 作者单位:杜增吉徐州师范大学,数学科学学院,江苏,徐州,221116 相似文献(10条) 1.期刊论文 樊自安.FAN Zi-an 二元常系数微分方程组解的表达式 -石家庄学院学报2009,11(3) 对于二元常系数微分方程组给出了解的表达式,利用解的表达 ...