多元一次方程组

早在初中时,我们就学习过多元一次方程组的求程。

所谓多元一次组指的是含有有个未知数,并且含有未知数的项的次数都是1方程组。

一个简单的例子:
学校的篮球数比排球数的2倍多3个,足球数与排球数的比是2:3,三种球共48个,求三种球各有多少?
假设篮球个数为x1,排球个数为x2,足球个数为x3,
根据题目有以下的关系式成立:
cd5d438fadc34af9afd4a5244e61619a?from=pc.jpg
图1. 多元一次方程组

对于多元一次方程组,最常用消元法来进行求解。

思路是经过各个方程组之间的线性运算,消去一些未知数,最终得到一个未知数的方程。

根据该方程解出未知数,再代回其它方程得到其它未知数的值。

整个过程非常繁琐,一不小心就会出错;

矩阵运算

在学习了矩阵运算之后,这一求解过程可以化为矩阵运算而大大简化;

跟初等数学里面学习的加、减、乘、除运算规则一样,伟大的数学家们也为矩阵运算制定了一些规则。

对于矩阵的加、减,限于文章篇幅,在此不作赘述。

对于矩阵的乘法,有以下的规则:
设有一个m行s列的矩阵表示为 050b88dfb3d94bef96bb8296f9c35c7b.jpg
一个s行n列的矩阵表示为 f590d295ed034ef1bea58ba4d28c9377.jpg
则矩阵A左乘于矩阵B所得到的矩阵C是这样一个矩阵:
  (1) 行数与(左矩阵)A相同,列数与(右矩阵)B相同,即C= afa450a31f524caab2b3ecf7939cadb1.jpg
  (2) C的第i行第j列的元素由A的第i行元素与B的第j列元素对应相乘,再取乘积之和。


比如,两个矩阵:
7f77cd8c3eea47eab1150099928966cc?from=pc.jpg
A、B矩阵

计算 8d0bfcc18a0b4633aecabbecf6cbef06.jpg ,从规则(1),可以知道C是2行3列的矩阵。

每个元素的数值由规则(2)计算得到,比如:

第0行0列的数值为:
c00=1*1+2*-1=-1;

第0行1列的数值为:
c01=1*2+2*1=4;

依次算出各元素的数值,并表示为矩阵,得到:
cf134a6c693b4b0eb7f1b68ca25fadf3?from=pc.jpg
C=AB

对图1的多元一次方程组,我们做一些变形,

并把未知数表示为3行1列的矩阵,可以得到:
c774a8481ba44ed3a7a1c6610b6e357d?from=pc.jpg
矩阵表达

表示为: 87661d22823f447f98573931ae1addfd.jpg ,然后就可以经过矩阵运算求解方程组。
左右两边同时左乘A的逆矩阵,得到: e9970b8e4bef43a89bd0c381cb71522d.jpg

N年前,matlab还能随意使用,做矩阵运算非常只是几条命令的事情,非常简单。

在被matlab抛弃的日子里,虽然对matlab还怀有深深的眷恋,无奈落花有意、流水无情;

我只能退而求其次用python来排遣寂寞。

利用,以下python脚本,
import numpy.matlib
  • import numpy as np
  • A=np.array([[1,-2,0],[0,3,-2],[1,1,1]])
  • Y=np.array([[3],[0],[48]])
  • X=np.matmul(np.linalg.inv(A),Y);
  • 复制代码
    求得,[x1,x2,x3]=[23, 10, 15];

    矩阵运算在电路计算、分析中的应用

    在线性电路中,首先通过电路理论得到电路的线性方程组。

    把线性方程组的求解转换为矩阵运算。

    再用python编写脚本进行矩阵运算,就可以算出电路中的各电气参数;

    如下图的电路分析题目,求解i1和i2;
    22cd31df55cd4731b16ab6623ca8401a?from=pc.jpg
    电路分析实例

    根据网孔电流法,有:
    6694232ce6d04cb685c0957b289a5671?from=pc.jpg
    网孔电流法列出方程组

    化简得到:
    918a0c9e4ddf4adaa850a2c33311463e?from=pc.jpg
    化简后的方程组

    编写python脚本,
    import numpy.matlib
  • import numpy as np
  • A=np.array([[30,-20],[-20, 80]])
  • Y=np.array([[-10],[20]])
  • X=np.matmul(np.linalg.inv(A),Y)
  • 复制代码
    计算得到[i1,i2]=[-0.2, 0.2];

    总结

    本文旨在介绍矩阵的概念、矩阵运算的规则,以及与多元一次方程组的关系。

    并以一个简单的电路为例,介绍电路计算过程;

    矩阵运算以及电路分析远不限于此,

    另外通过符号的矩阵运算,可以得到电参数的解析表达式并对其化简,可以大大简化滤波电路的分析计算。

    来源:物联网全栈开发