番茄的随笔12：从矩阵变换看电力系统模变换与序变换（上）

1. 概述

        蛇年十月，凌晨三点，喂完崽后，辗转反侧。崽半岁有余，憨态可掬，妻鼾声如雷，不绝于耳，觉吾之所幸，兴起，遂不能寐。

        哈哈，番茄文笔还是不错的，早先说过，鄙人高中可是语文课代表，深以为幸。既然睡不着觉，脑子就干点活吧。

        想想一直以来，困扰自己的电力系统故障分析中一些概念，结合AI搜索的资料，想明白一些事情，这也不能算是失眠了，这叫思考。

2. 矩阵对角化与解耦

      让我们来回想一下矩阵的相关概念：特征值、特征相量和矩阵对角化。

      对于一个n×n的方阵A，我们在看到这n×n个数的时候，难免不知所云。这玩意有个啥用，如果Y=AX，而A是个方阵，你看到的是个这玩意：

        但凡是个地球人，想一眼看懂y1和x1的关系，都得直摇头，但是，如果B是个对角矩阵，而Y* = BX*，你看到的就是这个：

       再问你y1*和x1*的关系，就很容易回答了。

       所以，看到了，这就是矩阵对角化的魅力。搁在线性代数里你管它叫对角化，不挑你毛病，但是，放在故障分析里面，你得叫人家什么，叫人家什么，对喽，解耦。

       回到矩阵对角化，方针A可对角化（对角矩阵B），也就是存在一个可逆矩阵P，满足：

      对角化后，对角矩阵B对角线上的每个元素就是A的特征值，矩阵P的每列的列向量P1、P2…称为特征向量。

      什么样的矩阵可以对角化，别的不用管：实对称矩阵一定可以对角化。这个结论就够了。

      考研线性代数的朋友，应该回忆起了求解矩阵的特征值和特征向量的惨痛记忆，这是考研数学的大题里的一个步骤，属于圈起来要考的考点。但是这个求解的过程中，我们也知道了，如果特征值还能说是唯一的，那特征向量绝对不能说是唯一的，这也可以理解，条条大路通罗马，罗马就一个，但是去罗马的路多得很，高兴了你先跑到太阳上，再去罗马都没问题。也就是说，方法有很多，但目的就一个，电力系统里面的变换，目的就一个，解耦。

        我们知道三相电力系统里面电压U和电流I的关系有U=ZI：

      让地球人直接写三相电压和三相电流的关系，这很难弄明白，于是，就用到了我们前面说的“解耦”，把这个阻抗矩阵对角化以后，就容易多了。

3. 阻抗矩阵的罗马在哪里

        我们来看看这个阻抗矩阵Z，第一眼什么鬼东西，再看看，好像有规律，这玩意有点意思啊，看起来是9个元素，但是聪明的你肯定首先会猜：

1. Zab和Zba是不是相等啊；
2. Zaa、Zbb、Zcc是不是也相等啊；
3. Zab、Zbc、Zca是不是也相等啊；
4. 如果这些都相等，该多好啊。

     恭喜你，兄弟，三相电力系统里面，这些真的能看成相等的（近似），你有了自阻抗和互阻抗的概念。

      自阻抗和互阻抗来源于多端口网络，配电网中任意一条线，都可以看成多端口网络；对于多端口网络，自阻抗描述的是一个端口自身的电压与电流关系，一个端口的自阻抗是当其他所有的端口都开路时，该端口的电压和电流的比值；互阻抗描述的是一个端口电压与另一个端口电流的关系，反映的是两个端口网络间的耦合强度。之前番茄的文章（随笔5）中也尝试理解过自阻抗和互阻抗，可以温故一下而知心-新。

   查阅了一下电力系统线路的自阻抗和互阻抗计算公式：    

        可见，假设理想的配电网系统，三相导线认为完全相同，呈现等角三角形排列，那么猜想2）和3）也是成立的。我们假定三相导线各自的自阻抗为Zs，互阻抗为Zm，那么阻抗矩阵转换为：

          利用卡尔丹公式求解这个方程组：

      是不是很熟悉，这三个特征值分别就是对应了零序阻抗，正序阻抗和负序阻抗。

      于是我们费劲千辛万苦终于找到了阻抗矩阵Z的对角化矩阵，满足

，对角矩阵为：