• 沒有找到結果。

第10讲线性方程组

N/A
N/A
Protected

Academic year: 2021

Share "第10讲线性方程组"

Copied!
51
0
0

加載中.... (立即查看全文)

全文

(1)

传媒与信息工程学院 欧 新 宇

第10讲 线性方程组

第5章 线性方程组

(2)

线性方程概述

线性方程组解的几何意义

高斯消元法与行阶梯方程组

矩阵的初等变换

线性方程组应用实例(Python)

(3)

线性方程概述

(4)

【例5.1】食品配方的应用

某食品厂收到某种食品的订单,要求这种食品由甲、乙、丙、丁四种原料 做成,且该食品中含蛋白质、脂肪和碳水化合物的比例分别为15%、5%和12

%。而甲、乙、丙、丁原料中含蛋白质、脂肪和碳水化合物的百分比由下表给 出。那么,如何用这四种原料配置出满足要求的食品呢?

𝑥 1 + 𝑥 2 + 𝑥 3 + 𝑥 4 = 1

20%𝑥 1 + 16%𝑥 2 + 10%𝑥 3 + 16%𝑥 4 = 15%

3%𝑥 1 + 8%𝑥 2 + 2%𝑥 3 + 5%𝑥 4 = 5%

10%𝑥 1 + 25%𝑥 2 + 20%𝑥 3 + 5%𝑥 4 = 12%

𝑥 1 + 𝑥 2 + 𝑥 3 + 𝑥 4 = 1

20𝑥 1 + 16𝑥 2 + 10𝑥 3 + 15𝑥 4 = 15 3𝑥 1 + 8𝑥 2 + 2𝑥 3 + 5𝑥 4 = 5 10𝑥 1 + 25𝑥 2 + 20𝑥 3 + 5𝑥 4 = 12

解出上述方程组的解,即可获得四种原料的配比。

1. 线性方程概述

某食品

蛋白质(%) 20 16 10 15 15

脂肪(%) 3 8 2 5 5

碳水化合物(%) 10 25 20 5 12

(5)

线性代数的任务之一——解方程 1. 线性方程概述

【定义】若线性方程组无解 ,则称该方程组是不相容的;如果

线性方程组至少存在一个解,则称该方程是相容的。

(6)

线性方程组解的求解和几何意义

(7)

适定二元线性方程组的求解 2. 线性方程组解的几何意义

(𝑎) ൜𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 3𝑥

2

= 3 ⇒ 𝑟

2

= 𝑟

2

+ 𝑟

1

⇒ ቊ𝑥

1

− 2𝑥

2

= −1 𝑥

2

= 2

(𝑏) ቊ𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 2𝑥

2

= 3 ⇒ 𝑟

2

= −𝑟

2

⇒ ቊ𝑥

1

− 2𝑥

2

= −1 𝑥

1

− 2𝑥

2

= −3

(𝑐) ቊ𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 2𝑥

2

= 1 ⇒ 𝑟

2

= 𝑟

2

+ 𝑟

1

⇒ ቊ𝑥

1

− 2𝑥

2

= −1 0𝑥

1

+ 0𝑥

2

= 0

矛盾方程,

无解

无穷解

可以由下而上

地回 代 解 出 𝒙 𝟐 =2 ,

𝒙 𝟏 =3 , 这 是 解 线

性 方 程 组 的 规 范

方法。

(8)

适定二元线性方程组的求解(几何形态)

2. 线性方程组解的几何意义

ቊ 𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 2𝑥

2

= 1 (c) 有无穷个解 ቊ 𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 3𝑥

2

= 3 (a) 有唯一解

ቊ 𝑥

1

− 2𝑥

2

= −1

−𝑥

1

+ 2𝑥

2

= 3

(b) 无解

(9)

超定二元方程组的近似解

三个方程,只有两个变量。它们所对应的三根直线并不共点,

即方程组 不相容 ,称为超定方程组。它没有精确解,但有近似 解——最小二乘解。

2. 线性方程组解的几何意义

𝑥 1 + 𝑥 2 = 1 𝑥 1 − 𝑥 2 = 3

−𝑥 1 + 2𝑥 2 = −3

(10)

三元方程组的求解

从第2个方程中减去第1个方程的2倍,得到 x 系数为零的新的 第2个方程,再从第3个方程中减去第1个方程的-5倍,得到 x 系数为零的新的第3个方程。

再将第2个方程乘7/5与第3个方程相加,在第3个方程中消去 y。

于是形成了阶梯形的结构。可以由下而上地回代解出z, y, x,

这是解线性方程组的规范方法。

2. 线性方程组解的几何意义

(11)

三元方程组解的几何意义

这三个方程在笛卡尔坐标系中的图形是三个平面,方程组的 解就是它们的交点坐标。

若将第三个方程改一下,消元后剩了两个方程,交点就成了 交线,说明有无数个解,构成了一根直线。

2. 线性方程组解的几何意义

3 2

1

2 1

3

2 4

4 4 4

2 3 3 5 3 5 5 3 5

4 11 5 3 5 0 0

r r

r r

r r

x y z x y z x y z

x y z y z y z

x y z y z

− −

+ − = + − = + − =

  

 − + = ⎯⎯⎯⎯ →  − + = − ⎯⎯⎯→  − + = −

  

 − − =  − + = −  =

  

(12)

三元方程组解的几何意义

对于更多元的线性方程组,不可能想象出其空间的几何图形,

但关于欠定、适定和不相容方程(超定)的基本概念是一脉相承的,

它们的解的特性也都可以推广到高维空间。

2. 线性方程组解的几何意义

(13)

课堂互动一 Link

(14)

高斯消元法与行阶梯方程组

(15)

将求解方法推广到高阶系统,需要借助于矩阵,并用计算机进行求 解,这是线性代数与初等代数的区别。设m为方程的个数,n为未知数

的个数,则n元线性方程组可以表示为:

𝑎 11 𝑥 1 + 𝑎 12 𝑥 2 + ⋯ + 𝑎 1𝑛 𝑥 𝑛 = 𝑏 1 𝑎 21 𝑥 1 + 𝑎 22 𝑥 2 + ⋯ + 𝑎 2𝑛 𝑥 𝑛 = 𝑏 2

𝑎 𝑚1 𝑥 1 + 𝑎 𝑚2 𝑥 2 + ⋯ + 𝑎 𝑚𝑛 𝑥 𝑛 = 𝑏 𝑚

线性方程组中解的全体称为方程组的解集(solution set)。解方程组 就是求其全部解,即解集。一般情况下,变量的个数n与方程的个数m 不一定相等。

适定方程组

(m=n) :存在着唯一的一组解;

欠定方程组

(m<n) :其解存在,但不唯一;

超定方程组

(m>n) :不存在精确解,可以求出其近似解。

3. 高斯消元法与行阶梯方程组

(16)

等价方程组

【定义】若两个含有相同变量的方程组具有相同的解集,则称它们是等 价的,这两个方程被称为同解方程组。

考虑以下两个方程组:

(𝑎) ቐ

3𝑥

1

+ 2𝑥

2

− 𝑥

3

= −2 𝑥

2

= 3

2𝑥

3

= 4

(𝑏) ቐ

3𝑥

1

+ 2𝑥

2

− 𝑥

3

= −2

−3𝑥

1

− 2𝑥

2

+ 𝑥

3

= 5 3𝑥

1

+ 2𝑥

2

+ 𝑥

3

= 2

不难得到它们的解都是(-2, 3, 2),因此这两个方程称为同解 方程组。

【定义】消元法: 通过消元变换把方程组化为容易求解的阶梯形结

构的同解方程组。

(17)

【例5.2】 利用消元法求解方程组

求解方程组 ቐ3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 3𝑥

1

+ 3𝑥

2

− 𝑥

3

= −5 2𝑥

1

+ 2𝑥

2

− 𝑥

3

= 4

(I)

解:对(I)式进行消元变化可以得到(II), (III)式。

(II)ቐ

3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 𝑥

2

+ 𝑥

3

= −1 2/3𝑥

2

+ 1/3𝑥

3

= 20/3​​​​

(III) ቐ 3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 𝑥

2

+ 𝑥

3

= −1

− 1/3𝑥

3

= 22/3​​​​

1. 将(I)式中的第一个方程分别乘以-3/3及-2/3,加到第二、三方程上,可以 消去后两个方程中的变量𝑥

1

;

2. 将(II)式中的第二个方程乘-2/3,加到第三个方程中,消去其中的𝑥

2

,得 到第(III)个方程组;

3. 第(III)个方程组称为行阶梯形方程组。这样的阶梯形方程组可以用回代法 方便地逐个求出它的解。

3. 高斯消元法与行阶梯方程组

(18)

【例5.2】 利用消元法求解方程组 3. 高斯消元法与行阶梯方程组

回代过程如下:

1. 由最后一个方程解出x 3 = -22 2. 代入第二个方程,解得:

x 2 =21

3. 再将x 2 ,x 3 代入第一个方程,

得到:x 1 = −30。

最后的行阶梯形方程组(VI)只保留了系数均为1的对角项。

得到它就等于求出了方程组的解。

(IV)ቐ

3𝑥

1

+ 2𝑥

2 − 1/3𝑥 2

= −48/3

3𝑥 1 + 2𝑥 2 − 1/3𝑥 2

= −12/3

3𝑥 1 + 2𝑥 2

− 1/3𝑥

3

=

−22/3

(V) ቐ

3𝑥

1 + 2𝑥 2 − 1/3𝑥 2

= −90/3

3𝑥 1 + 2𝑥 2 − 1/3𝑥 2

= 212/3

3𝑥 1 + 2𝑥 2

− 1/3𝑥

3

= 22/30 (VI) ቐ

3𝑥 1 + 2𝑥 2 − 1/3𝑥 2

= −30/3

3𝑥 1 + 2𝑥 2 − 1/3𝑥 2

= 212/3

3𝑥 1 + 2𝑥 2 − 1/3𝑥 3

= −220

(19)

行阶梯矩阵

【定义】若一个矩阵满足:

(1)每一个非零行中的第一个非零元为1;

(2)第k行的元不全为零时,第k+1行首变量之前零的个数多于第k行首变量 之前零的个数;

(3)所有元素均为零的行必在不全为零的行之后;

则称其为行阶梯形矩阵。

【例5.3】下列矩阵为行阶梯形矩阵 1 4 2

0 1 3 0 0 1

, 1 2 3 0 0 1 0 0 0

, 1 3 1 0 0 0 1 3 0 0 0 0

, 1 0

0 1 , 0 1 0 0 0 1

【例5.4】下列矩阵不是行阶梯形矩阵 2 4 2

0 3 3 0 0 4

, 1 2 3 0 0 2 0 0 0

, 0 0 0 0 0 0 1 3 0 1 0 0

, 0 1 1 0

3. 高斯消元法与行阶梯方程组

(20)

消元法的三种同解变换

在消元过程中,主要对原方程组进行了三种变换:

① 互换两个方程的位置;称为位置变换。

② 用一个非零数 k 乘某个方程;称为数乘变换。

③ 把一个方程的 k 倍加到另一个方程上。称为消元变换。

这三种变换称为线性方程组的同解变换。因为对方程组而言,

这些变换不会改变方程组的解。

消元法规则刻板,容易程序化,其计算量又是最少的,所以 可以利用计算机程序实现。消元法是一切线性方程组求解的基础。

3. 高斯消元法与行阶梯方程组

(21)

课堂互动二 Link

(22)

矩阵的初等变换

(23)

重新审视【例5.2】可以发现,从方程组(I)变换到方程组(VI)的全过程中,

方程组中的变量𝒙

𝒏

并没有参与任何运算,参与运算的只是方程组的系数和常 数。于是可省略变量,只提取方程组等式左端的系数和右端常数,分别写成 如下系数矩阵A和常数向量b。

3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 3𝑥

1

+ 3𝑥

2

− 𝑥

3

= −5 2𝑥

1

+ 2𝑥

2

− 𝑥

3

= 4

⇒ 𝑨 =

3 2 −2 3 3 −1 2 2 −1

, 𝒃 =

−4

−5 4

其中,矩阵A的行号表示方程的序号,列号表示变量 𝒙 的序号。由矩阵A,

b很容易恢复出方程组原型。所以,可以通过这两个矩阵来研究线性方程组。

学习线性代数的主要目标就是,要学会利用矩阵来描述系统,并用矩阵软件 工具去解决各种问题。

4. 矩阵的初等变换

(24)

基于增广矩阵的方程组

【定义】由线性方程组所有系数所构成的矩阵,称为线性方程组的系数 矩阵。系数矩阵为 n 阶方阵的方程组也称为

n阶方程组。

【例5.2】

中, A, b是三阶线性方程组的系数矩阵。把A,b并排起来 组成的矩阵C 称为方程组的增广矩阵,知道了C也就知道了线性方程组 的全部参数。所有针对方程组所做的消元变换,都可以表示为对系数增 广矩阵的行变换。

4. 矩阵的初等变换

𝑨 =

3 2 −2 3 3 −1 2 2 −1

, 𝒃 =

−4

−5 4

, 𝑪 = [𝑨, 𝒃] =

3 2 −2 3 3 −1 2 2 −1

−4

−5

4

(25)

矩阵的初等行变换 4. 矩阵的初等变换

【定义】下面三种变换称为矩阵的初等行变换:以下的

r

表示行

(row)

(1)交换两行的位置(交换第i, j行,记作r

i ←→r j

);

(2)以非零数k乘某行(以k乘第i行,记作kr

i

);

(3)把某一行的k倍加到另一行上(把第j行的k倍加到第i行上,记作r

i +kr j

) 矩阵的这三种初等行变换就对应于方程组的三种初等变换(位置、数乘 和消元)。它们都是可逆的,且其逆变换是同一类型的初等变换,由此可知,

初等行变换是同解变换。

如果矩阵

A经有限次初等行变换变成

矩阵

B,就称

矩阵

A与

矩阵

B等价。

(26)

方程消元等价于行变换 4. 矩阵的初等变换

线性方程组和它的增广矩阵是一一对应的,对线性方程组进行初等行变换

就是对其增广矩阵进行初等行变换。于是解方程组的过程可用矩阵的初等行变 换来一一对照。

前向消元过程

3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 3𝑥

1

+ 3𝑥

2

− 𝑥

3

= − 5 2𝑥

1

+ 2𝑥

2

− 𝑥

3

= 4

⇒ ቐ

3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 𝑥

2

+ 𝑥

3

= −1 2/3𝑥

2

+ 1/3𝑥

3

= 20/3

⇒ ቐ

3𝑥

1

+ 2𝑥

2

− 2𝑥

3

= −4 𝑥

2

+ 𝑥

3

= −1

− 1/3𝑥

3

= 22/3

系数矩阵

(27)

方程消元等价于行变换 4. 矩阵的初等变换

对最后一个矩阵按照消元回代的思想继续进行初等行变换,

最终可以获得一个对角矩阵;之后再将各行除以对角元素,最终 可以获得一个单位矩阵,此时最右一列就是方程组的解。

3 2 −2 − 4 0 1 1 −1 0 0 −1/3 22/3

3 2 0 −48 0 1 0 21 0 0 −1/3 22/3

3 0 0 −90 0 1 0 21 0 0 −1/3 22/3

1 0 0 −30 0 1 0 21 0 0 1 −22

⇒ ቐ

𝑥 1 = −30

𝑥 2 = 21

𝑥 3 = −22

(28)

行最简形矩阵 4. 矩阵的初等变换

事实上,在进行消元变化的时候,如果给出的原始矩阵不是满秩的方阵,

则所获得的矩阵将不再是一个单位矩阵,而是行最简形矩阵。这意味着对于系 数矩阵为非适定矩阵或奇异矩阵的方式组来说,方程组将无法获得唯一解。

【定义】若一个矩阵满足

(1)矩阵式行阶梯型的;

(2)每一行的第一个非零元是该列唯一的非零元,

则称该矩阵为行最简形矩阵(reduced row echelon form)。

下列矩阵为行最简形:

1 0 0 1

1 0 2 0 1 0 0 0 0

1 0 0 0 1 2

0 1 2 0 0 1 0 0 0 0 0 0

1 3 0 0 0 0

0 1 1 5 0 0

(29)

课堂互动三 Link

(30)

线性方程组应用实例(Python)

(31)

Numpy.linalg库

numpy.linalg是Numpy库中提供的线性代数库,主要包括 矩阵的点积、内积、矩阵积、行列式、逆矩阵和求解矩阵方程的 函数。其中np.linalg.solve(A,b)可以实现线性方程的求解。但对 于非适定方程或系数矩阵为奇异矩阵时,该库无法进行运算。

5. 线性方程组应用实例(Python)

3𝑥 1 + 2𝑥 2 − 2𝑥 3 = −4 3𝑥 1 + 3𝑥 2 − 𝑥 3 = −5 2𝑥 1 + 2𝑥 2 − 𝑥 3 = 4

𝑥 1 = −30

𝑥 2 = 21

𝑥 3 = −22

(32)

5. 线性方程组应用实例(Python)

如何解决?

使用系数增广矩阵的行变换法求解方程组

将系数矩阵化简为行最简形矩阵 1 0 0

0 1 0 0 0 1

0 0 0 0 0 0 0 0 0

0 0 0

1 0

行最简形矩阵 0 1

(33)

Sympy 库

SymPy是一个符号计算的Python库。它的目标是成为一个全 功能的计算机代数系统,同时保持代码简洁、易于理解和扩展。它 完全由Python写成,不依赖于外部库。

此处,我们使用sympy库中最基本的运算函数:

解线性方程组(求阶梯矩阵):Sympy.Matrix.rref()

求矩阵的秩: Sympy.Matrix.rank()

5. 线性方程组应用实例(Python)

(34)

【计算题】用计算机求解线性方程组

【例5.5】求下列方程组的解

2𝑥 1 − 2𝑥 2 + 2𝑥 3 + 6𝑥 4 = −16 2𝑥 1 − 𝑥 2 + 2𝑥 3 + 4𝑥 4 = −10 3𝑥 1 − 𝑥 2 + 4𝑥 3 + 4𝑥 4 = −11 𝑥 1 + 𝑥 2 − 𝑥 3 + 3𝑥 4 = −12

5. 线性方程组应用实例(Python)

(35)

【计算题】用计算机求解线性方程组

【例5.6】设方程组的系数矩阵𝑨,𝒃如下,判断它的解的性质及 𝑨的秩。

𝐴 =

−2 −2 2 1 −5 1

2 −2

−3 −1

−1 2 −5

−1 2 1

6 5 0 −1

, 𝑏 =

−2

−1 2 0

5. 线性方程组应用实例(Python)

(36)

【计算题】用计算机求解线性方程组

【结果分析】

在本例中,方程组是

欠定矩阵

,无法用

numpy.linalg.solve()方法进行求解。在后续的

例题中,我们都将使用基于行最简形矩阵化简 的方法rref()完成。

5. 线性方程组应用实例(Python)

函数 sympy.matrix.rref()方法的输出是一个2维元组,其中:

⚫ 第一项,为化简后的 行最简形矩阵;

⚫ 第二项,为行最简形中 “ 1”所在列的序号 。此例中,该项的的 长度为4,说明原矩阵4个主元和主元行,即矩阵的秩 𝑟=4 。 该结论也可以通过M.rank()方法进行验证。

此处, 秩<方程未知数个数 ,因此方程为欠定矩阵。

(37)

【计算题】用计算机求解线性方程组

【结果分析】

既然方程组是欠定方程,那 么不难发现未知数𝑥 3 ,𝑥 4 不存在唯 一解,即行最简形第三行的结果:

𝑥 3 − 𝑥 4 = −2/3

5. 线性方程组应用实例(Python)

设 𝑥 4 = 𝑐 , 𝑐 是任意常数,则可以得到该方程组的解为:

𝑥 1 = −2/9 𝑥 2 = 2/9

𝑥 3 = c − 2/3

𝑥 4 = −1/3

(38)

【应用题】计算插值多项式

【例5.7】求插值多项式𝑝(𝑡) = 𝑎 0 + 𝑎 1 𝑡 + 𝑎 2 𝑡 2 + 𝑎 3 𝑡 3 的各系数,

使它能通过下表中各点。并求多项式在 𝑡 =1.5 时的值。

解:根据题目要求,列出方程组:

𝑎 0 𝑎 0 𝑎 0 𝑎 0

+ + + +

0 𝑎 1 2𝑎 1 3𝑎 1

+ + + +

0 𝑎 2 4𝑎 2 9𝑎 2

+ + + +

0 𝑎 3 8𝑎 3 27𝑎 3

=

= =

=

3 0

−1 6

⇒ 𝐴𝒙 = 𝑏 ⇒ 𝒙 = 𝐴\𝑏

5. 线性方程组应用实例(Python)

𝒕 𝒊

0 1 2 3

𝑓(𝑡 𝑖 ) 3 0 −1 6

(39)

【应用题】 计算插值多项式

由输出结果可以得到方程组的解: 𝑎

0

= 3, 𝑎

1

= −2, 𝑎

2

= −2, 𝑎

3

= 1。

将结果代入多项式,可以得到插值多项式:

𝑝 𝑡 = 3 − 2𝑡 − 2𝑡

2

+ 𝑡

3

⚫ 当t=1.5时,可以求得多项式的值:

𝑝 1.5 = 3 − 2 1.5 − 2 1.5

2

+ 1.5

3

= −1.125

5. 线性方程组应用实例(Python)

-5

0

5

10

15

(40)

【应用题】平板稳态温度的计算

【例5.8】如图所示,假设存在一个密度均匀 的金属元器件,在其四周我们获得了8个点的 瞬时温度。已知元器件各点的温度等于其四周 温度的平均值,求abcd四个点的温度。

解:按要求可以得到abcd四点的温度方程式,整理以后可以得 到其系数矩阵 Ax = b:

𝑥

𝑎

= (10 + 20 + 𝑥

𝑏

+ 𝑥

𝑐

)/4 𝑥

𝑏

= (20 + 40 + 𝑥

𝑎

+ 𝑥

𝑑

)/4 𝑥

𝑐

= (10 + 30 + 𝑥

𝑎

+ 𝑥

𝑑

)/4 𝑥

𝑑

= (40 + 30 + 𝑥

𝑏

+ 𝑥

𝑐

)/4

,𝐴𝑥 =

4 −1

−1 4

−1 0 0 −1

−1 0 0 −1

4 −1

−1 4

𝑥

𝑎

𝑥

𝑏

𝑥

𝑐

𝑥

𝑑

=

30 60 40 70

5. 线性方程组应用实例(Python)

(41)

【应用题】平板稳态温度的计算 5. 线性方程组应用实例(Python)

从程序输出结果可以得到abcd四个 点的温度分别是:a=20度, b=55/2度, c=45/2度, d=30度。

(42)

【应用题】交通流量分析

【例5.9】如图所示,某城市市区的交叉路 口由两条单向车道组成。图中给出了在交通 高峰时段每小时进入和离开路口的车辆数。

计算在四个交叉路口间车辆的数量。

解:在每一个路口,必然存在进入的车辆与 离开的车辆数相等,于是有:

节点A : x 1 + 450 = x 2 + 610

节点B : x 2 + 520 = x 3 + 480

节点C: x 3 + 390 = x 4 + 600

节点D: x 4 + 640 = x 1 + 310

5. 线性方程组应用实例(Python)

(43)

【应用题】交通流量分析

将节点信息转换为方程式,并得到增广矩阵。

𝑥 1 𝑥 3 4𝑎 2

−𝑥 1

− 4𝑎 2

0

𝑥 2 𝑥 2 4𝑎 2

0

𝑥 2

− 4𝑎 2

0

𝑥 2 𝑥 3 𝑥 3 0

𝑥 2 𝑥 2

− +

𝑥 2 𝑥 2 𝑥 4 𝑥 4

=

= =

=

160

−40 210

−330

1 − 1 0 0 0 1 − 1 0 0 0 1 − 1

−1 0 0 1

160

−40 210

−330

5. 线性方程组应用实例(Python)

(44)

【应用题】交通流量分析

【结果分析】

从结果不难看出,该方程是一个欠定方程组,其秩=3<未知 数个数4。即得不出唯一解。如果设 𝑥4 为自由变量,并将它移到 增广项位置,就可以求出解。其物理意义是,给定任意 𝑥4 ,均 可得到每个路口的交通流量。例如,设 𝑥4=150 ,则:

5. 线性方程组应用实例(Python)

(45)

【应用题】化学方程式的配平

【例5.10】给定化学方程式,试将其进行配平。即通过设置系 数,实现每个原子左右数量都相等。

(𝑥 1 )𝐶 3 𝐻 8 + (𝑥 2 )𝑂 2 → (𝑥 3 )𝐶𝑂 2 + (𝑥 4 )𝐻 2 𝑂

⚫ 四种物质的成分列向量:𝐶 3 𝐻 8 3 8 0

, 𝑂 2 0 0 2

, 𝐶𝑂 2 1 0 2

, 𝐻 2 𝑂 0 2 1

⚫ 配平方程为:𝑥 1 ⋅ 3 8 0

+ 𝑥 2 ⋅ 0 0 2

= 𝑥 3 ⋅ 1 0 2

+ 𝑥 4 ⋅ 0 2 1

⚫ 移项化后可知:𝐀 = 3 0 − 1 0 8 0 0 − 2 0 2 − 2 − 1

, 𝐛 = 0 0 0

5. 线性方程组应用实例(Python)

(46)

【应用题】化学方程式的配平

从结果不难看出,该方程是一个欠定方程组,其秩=3<未知数个数4。即得 不出唯一解。为了实现化学方程式的配平,我们可以设 𝑥4=4 (通常取分母的 最小公倍数),则方程的解可以变为唯一解。

由此,可以得到配平后的方程:𝑪

𝟑

𝑯

𝟖

+ 𝟓𝑶

𝟐

= 𝟑𝑪𝑶

𝟐

+ 𝟒𝑯

𝟐

𝑶

5. 线性方程组应用实例(Python)

(47)

【应用题】电路图分析

【例5.11】在一个电路中,可根据电阻大小和电源电压来确定电路中各分支

的电流。如图所示,根据基尔霍夫定律,

1. 任一节点上流出电流的量等于流入电流的量

2. 任一回路上电压的代数和等于各元件压降的代数和 3. 压降E满足欧姆定律𝐸 = 𝑖𝑅

试求三条支路上的电流𝑖

1

, 𝑖

2

, 𝑖

3

。 根据以上定律可以得到:

⚫ 节点A:𝑖

1

− 𝑖

2

+ 𝑖

3

= 0

⚫ 节点B: − 𝑖

1

+ 𝑖

2

− 𝑖

3

= 0

⚫ 上层回路:4𝑖

1

+ 2𝑖

2

= 8

⚫ 下层回路:2𝑖

2

+ 5𝑖

3

= 9

5. 线性方程组应用实例(Python)

1 −1 1

−1 1 −1 4 2 0 0 2 5

0 0 8 9 增广矩阵

(48)

【应用题】电路图分析

根据以上定律可以得到:

⚫ 节点A:𝑖

1

− 𝑖

2

+ 𝑖

3

= 0

⚫ 节点B: − 𝑖

1

+ 𝑖

2

− 𝑖

3

= 0

⚫ 上层回路:4𝑖

1

+ 2𝑖

2

= 8

⚫ 下层回路:2𝑖

2

+ 5𝑖

3

= 9

5. 线性方程组应用实例(Python)

(49)

⚫ 了解二阶和三阶线性方程组在笛卡尔坐标系中的几何意义

⚫ 掌握利用初等变换把增广矩阵化为行最简形的方法

⚫ 秩表明独立方程的个数。系数矩阵与增广矩阵的秩相等是线性 方程组有解的充分必要条件(适定方程)

⚫ MATLAB实践

(1)能够利用Python构造矩阵

(2)能够利用Python求解线性方程组的解

(3)能够利用Python求矩阵的秩

本章小节

(50)

课后作业 Link

(51)

读万卷书 行万里路 只为最好的修炼

QQ:14777591 (宇宙骑士)

Email:ouxinyu@alumni.hust.edu.cn Tel:18687840023

參考文獻

相關文件

明龍計算一題兩個數相加的數學題目,不小心算成了相減,所得到的答

甲型禽流感 H7N9 H7N9 H7N9 H7N9 H7N9 H7N9 H7N9 H7N9 - - 疾病的三角模式 疾病的三角模式 疾病的三角模式 疾病的三角模式 疾病的三角模式

一階隨機差分方程式.

1871—1946)和稍晚的塞梵利(F.Severi,1879—1961),他们主要的结果是代数曲面的分类.头一个

笛卡儿企图通过坐标系给几何引进新方 法, 他的成就远远超出他的期望. 坐标系是数 学中的双刃剑, 使得几何的目的可以通过代 数达到, 反过来,

一般说来, 对于一个区间上单调的函数的 图形都存在一个需要判别弧段位于相应的弦线

4、任给空间中某一点,及某一方向,过该定点且垂

• 孤波解的性质:1- x 趋向无穷时 u 趋于0;2- 波是单向传播的;3- 波幅与波速成正比;4- KdV方程的色散与非线性刚好抵消,导致