现代电路第1章.docVIP

第1章 矩阵运算的计算机方法及稀疏矩阵 利用基尔霍夫电流定律（KCL）、基尔霍夫电压定律（KVL）与元件约束方程结合起来建立的方程组，大部分均为线性代数方程组。本章主要讨论对线性方程组及相应矩阵的运算方法。一般而言，如果网络是线性网络，那么方程式也是线性的；非线性网络则导致非线性方程组，但其解也是某些工作点附近、通过使这些方程式线性化来求得，或者利用分段线性化来求解。所以，解线性方程组的方法是解所有问题的基础。解线性方程组的方法一般分为直接法、迭代法，在本章中仅研究直接法。对于大型的矩阵，如果其元素大部分为零，则求解要用稀疏矩阵的方法。 1.1 高斯消元法解线性方程组 高斯消元法是线性方程组直接解法中优秀的解法之一，它是建立在这样一个事实基础上的，即将一个方程乘以一个常数加到另一个方程上，方程组的解不变。下面以一实例来说明该方法。 1.1.1 例题分析 [例1-1] 试解线性方程组 (1-1a) 解：将式(1-1a)中的1式与2式对调，得到 (1-1b) 将式(1-1b)中的2式减去1式的2倍，3式减去1式的5倍，这样后两个方程中的都被消去，得到 (1-1c) 将式(1-1c)中的2式乘以1/6，得到 (1-1d) 将式(1-1d)中的3式减去2式的17倍，得到 (1-1e) 式（1-1e）是一个阶梯形线性方程组，从中可以求出，，。过程如下： 由式（1-1e）中的3式解得，将代入2式得到，最后将，代入1式，得到。 显然，方程组（1-1a）～方程组（1-1e）都是同解方程组，因此方程组（1-1e）的解也就是方程组（1-1a）的解。这就是高斯消元法。由式（1-1a）～式（1-1e）叫做消去过程，而由方程组（1-1e）的3式、2式、1式自下而上逐个求出、、的过程叫做回代过程。 1.1.2 算法说明 线性方程组 (1-2) 式中，A是n×n维常数矩阵，B是n维已知向量，x为n维未知向量。将式（1-2）展开，得到 (1-3) 从形式上看，对矩阵A求逆可解方程组（1-3），得到 (1-4) 如果只需要一个输出变量，这时则用称之为克莱姆法则的方法来求得其解。 克莱姆法则可以这样来叙述，对式(1-3)表示的方程组，其解向量x的第k个分量是两个行列式的值之比，其分母是|A|，分子是用B代替A的第k列而得到的行列式的值。 当解小规模问题时，克莱姆法则是可以应用的，而且在理论研究上也是需要的；然而它的计算量很大，所以在计算机上很少使用。这些也同样适用于矩阵求逆。 高斯消元法是最好的算法之一，它是建立在这样一个事实基础上的，即将一个方程乘以常数加到另一个方程上，方程组的解不变。 考虑式(1-3)，把它改写成方程组，并且用表示右侧向量的元素，这样简化了表示，方程组变成 (1-5) 将式(1-5)中的第一个方程除以，并把它写成为 式中，=,。将其乘以，并且与 式(1-5)中的第二个方程式相加，则新的第二个方程的系数成为 （） 这样选择的相乘因子使得变为零。相似地对其他方程式，使 （） 则使第一列的系数除去为１之外，其他皆为零。实际上，不必要计算变为零的元素，这样在计算机存储中将删去，并且从开始计算。按照上面方式处理后，方程组(1-5)变成下列形式 (1-6) 在括号内上标表明其值已变换了一次。 下一步排除第一行和第一列，对式（1-6）第二个方程至第个方程施以同样的处理，其公式变为 为了实效，省略计算＝0，对所有行重复这个过程。由于起始状态为，则高斯消元法的一般公式成为 (1-7) 那么最后所得的方程组成为下列形式 (1-8) 求解方程组（1-8）所需的以下步骤称为回代过程。将最后一个变量xn代入第个方程式得到xn-1，并且照此做下去。一般讲，回代过程为 （） (1-9) 这里把上标去掉了。 在计算机的术语中，一种算法的代价用计算量来度量，每个计算量是一个乘法和一个减法的组合。可以证明，高斯消元法需要大约个计算量，而回代过程则用大约个计算量就可以完成（为矩阵的规模）。 1.1.3 选主元素 列主元素或者全主元素法可以改善线性方程组解的精度。考虑式(1-3)并假设，则无论高斯消元法还是LU分解法都不能应用，因为