第15章电路方程的矩阵形式教学课件

结束15-1 割集1.定义 连通图G的一个割集是G的一个支路集合，如果把这些支路移去，将使G(恰好)分离为两个部分，但是少移去其中一条支路，G将仍是连通的。F(a，d，f)这个支路集合就是G的一个割集。adfbceQ1adfbceQ2Q3Q4(a，b，e)(b，c，f)(c，d，e)显然，对右图，汇集于同一结点的支路都是G的一个割集。7/12/20241结束全移，G一分为二；少移一条，G连通。(b,d,e,f)是(a,d,e,f)不是G的割集！Q5adfbceadfbceQ6Q7adfbceadfbceQ8(a,e,c,f)是(a,b,c,d)也是原因：少移去e，G仍为两部分。7/12/20242结束(a,b,c,d,e)不是G的割集！原因：全移，G被分为三部分。2.割集的判断与确定 直观方便的方法是闭合面加定义。adfbceQ9注意：有些割集可能不易用与闭合面相切割的方法表示。abecdf无法作闭合面判断割集(a,b,c,d)。Qabcde与Q相切割的支路集合(a,b,e)不是割集。7/12/20243结束3.独立割集和基本割集KCL适用于任一闭合面。属同一割集的所有支路电流也满足KCL。对于一个连通图 G，总可以列出与割集数量相等的KCL方程。但它们不一定线性独立。(1)独立割集 与一组线性独立的KCL方程相对应的割集，称为独立割集。abecdfQ1Q2Q3Q4当割集的所有支路连接于同一结点时，割集的KCL变为结点的KCL。对较大规模的电路，用观察法选择一组独立割集是困难的。借助于树，就比较方便。7/12/20244结束(2)独立割集的确定选一个树，一条树支与相应的连支可以构成一个割集。由一条树支与相应的连支构成的割集叫单树支割集。对于具有n个结点b条支路的连通图，树支数为(n-1)条。这(n-1)个单树支割集称为基本割集组。btl1l2l3Q独立割集组不一定是单树支割集。就象独立回路不一定是单连支回路一样。而基本割集组是独立割集组。7/12/20245结束树支为2,3,4,6时的基本割集组 树支为5,6,7,8时的基本割集组。12345678Q1Q1(1,2,5,7,8)12345678Q2Q2(1,3,5,8)12345678Q3Q3(1,4,5)Q4Q4(5,6,7,8)12345678Q1Q2Q3Q4同一个图，有许多不同的树，因此能选出许多不同的基本割集组。7/12/20246结束15-2 关联矩阵、回路矩阵、割集矩阵1.关联矩阵的特点 描述结点与支路关联的矩阵。是一个(nb)阶的矩阵。Aa=12341 2 3 4 5 6-1-1+1 0 0 00 0-1-1 0+11i12i23i34i45i5i66+1 0 0+1+1 00+1 0 0-1-1(1)Aa的元素定义ajk=+1，支路k与结点j关联，且方向背离结点；ajk=-1，支路k与结点j关联，且方向指向结点；ajk=0，支路k与结点j无关联。7/12/20247结束(2)降阶关联矩阵A划去Aa中任意一行所得到的(n-1)b阶矩阵。A =12341 2 3 4 5 6-1-1+1 0 0 00 0-1-1 0+1+1 0 0+1+1 00+1 0 0-1-11i12i23i34i45i5i66被划去的行对应的结点可以当作参考结点。a 提示给定A可以确定 Aa，从而画出有向图。若以结点 4 为参考结点，把式中的第 4 行划去，得 A7/12/20248结束(3)用A表示KCL的矩阵形式b(=6)条支路电流可以用列向量表示i=i1,i2,i6 TAi=A=-1-1+1 0 0 0 0 0-1-1 0+1+1 0 0+1+1 0-1-1+1 0 0 0 0 0-1-1 0 +1+1 0 0+1+1 0 i1 i2 i6=-i1 i2 +i3-i3 i4 +i6+i1 +i4 +i5=000Ai=结点1的KCL结点(n-1)的KCL结点2的KCL Ai=01i12i23i34i45i5i667/12/20249结束(4)用A表示KVL的矩阵形式以b(=6)阶列向量表示支路电压：u=u1,u2,u6 T并取某一结点(取)为参考，(n-1=3)个结点电压的列向量：un=un1,un2,un3 T结点电压与支路电压之间的关系为u=ATunu1u2u3u4u5u6=1i12i23i34i45i5i66-un1+un3-un1 un1-un2-un2+un3 un3un2=un1un2un3-1 0 1-1 0 0 1-1 0 0-1 1 0 0 1 0 1 0AT可以认为，这是用A表示KVL的矩阵形式。7/12/202410结束小结 矩阵 A表示有向图结点与支路的关联性质。用 A表示的 KVL 的矩阵形式为 u=ATun 用 A表示的 KCL 的矩阵形式为 Ai=07/12/202411结束2.回路矩阵 描述回路与支路关联的矩阵。是一个(lb)阶的矩阵。(1)B 的元素定义 bjk=+1，支路k与回路j关联，且方向一致；bjk=-1，支路k与回路j关联，且方向相反；bjk=0，支路k与回路j无关联。1231 2 3 4 5 611234560 10-1 10 1 10 0 10 0 01-1 1B=7/12/202412结束(2)基本回路矩阵Bf Bf 反映了一组单连支回路与支路间的关联关系。写Bf时的排列顺序：先连支后树支。Bf=1l Bt(3)用B表示的KVL矩阵形式：u1+u3-u5+u6=0：u2+u3+u6=0：u4-u5+u6=01234561231 2 4 3 5 61 00 1-1 10 10 1 0 10 01 0-1 1Bf =Bu=01 0 1 0 1 10 1 1 0 0 10 0 0 1 1 1u1u2u3u4u5u6=0007/12/202413结束(4)用B表示的KCL矩阵形式 若用列向量表示 l(=3)个独立回路电流：il=il1 il2 ill T 则支路电流与回路电流之间的关系可以表示为i=BTil 可以认为是用B 表示KCL的矩阵形式。123456i1i2i3i4i5i6=il1il2il3 1 0 0 0 1 0 1 1 0 0 0 1-1 0-1 1 1 1 ii1il2il1+il2il3-il1-il3il1+il2+il3=7/12/202414结束3.割集矩阵Q描述割集与支路关联的矩阵。Q是一个(n-1)b阶的矩阵。各元素定义为：qjk=+1，支路k与割集j关联，且方向一致；qjk=-1，支路k与割集j关联，且方向相反；qjk=0，支路k与割集j无关联。若选单树支割集为一组独立割集，则得到基本割集矩阵Qf。排列顺序为先树支后连支。1234561231 2 3 4 5 6Q1-1-1 1 0 001 0 0 1 1 0-1-1 0-10 1Q2Q3Qf=1t Ql Q=7/12/202415结束(1)用割集矩阵Q表示的 KCL的矩阵形式 因属同一割集的所有支路的电流也满足KCL，所以Q i=0123456Q1Q2Q3-1-1 1 0 001 0 0 1 1 0-1-1 0-10 1i1i2i3i4i5i6=-i1 i2 +i3i1 +i4 +i5-i1 -i2-i4+i6=000Qf=1t Ql =100010001-11-1-10-101-13 5 6 1 2 4Q1Q2Q37/12/202416结束(2)用基本割集矩阵Qf表示KVL的矩阵形式式中 ut=ut1 ut2 ut(n-1)T 为树支电压列向量。对右图：ut=ut1 ut2 ut3T u=u3 u5 u6 u1 u2 u4Tu=Qf utT123456Q1Q2Q3u=100-1-10010101001-1-1-1ut1ut2ut3=ut1 ut2 ut3-ut1+ut2-ut3-ut1-ut3 ut2-ut3=u3=u5=u6=u1=u2=u4当选单树支割集为独立割集时，树支电压可视为割集电压。树支电压(割集电压)也是一组完备的独立变量，支路电压可以用树支电压表示。7/12/202417结束*15-3 矩阵A、Bf、Qf之间的关系2.对任一图G，当A、B、Q 的列按相同的支路编号排列时：ABT=0 或 BAT=0 QBT=0 或 BQT=0 3.若A、Bf、Qf 对应同一个树，且支路编号按先树支后连支的相同顺序排列写出。则有：TBtQl=-=-At Al-1TBt=At Al-11.A i=0 Q i=0u=A unTu=Qf utT在形式上相似。所以对某些图G有 Qf =A7/12/202418结束15-4 回路电流方程的矩阵形式一、复合支路 既含阻抗(导纳)，又有电源。(1)支路阻抗Zk是单一的R或L或C，但不是它们的组合；(2)可以缺少某种元件。但不许存在无伴电流源支路。二、支路方程的矩阵形式情况1 电路无互感Zk-+.Usk.Isk.Iek.Ik-+.Uk .Uk=Zk .(Ik+.ISk).-USk式中各量为第 k 条支路的阻抗、独立电流源和独立电压源。无独立源时将其置零。设.I=.I1.I2.IbT .U=.U1 .U2 .UbT.IS=.IS1.IS2.ISbT .US=.US1 .US2 .USbT则 .U=Z .(I+.IS).-US避免造成计算困难。7/12/202419结束Z称为支路阻抗矩阵，Z是对角矩阵，对角元素是各支路阻抗。情况2 电路有互感 设在b条支路中，1g支路之间相互有耦合，则有(g+1)b支路之间无耦合，关系式同情况1。.U=Z .(I+.IS).-USZ1 Z2.Zb00Z=.U1=Z1.Ie1jwM12.Ie2jwM13.Ie3 jwM1g.Ieg .-US1 .U2=jwM21.Ie1+Z2.Ie2jwM23.Ie3 jwM2g.Ieg .-US2 .Ug=jwMg1.Ie1jwMg2.Ie2jwMg3.Ie3 .-USg+Zg.Ieg.Ie1=.I1+.IS1，.Ie2=.I2+.IS2，；M12=M21，。7/12/202420结束有互感和无互感，方程形式相同。有互感时，Z 不再是对角阵。非对角线元素将含互感阻抗，其正负号根据同名端确定。=-.U=Z .(I+.IS).-US .U1 .U2 .Ug .Ug+1 .UbZ1jwM12jwM1g00jwM21Z2 jwM2g00jwMg1jwMg2Zg00000Zg+100000Zb.I1+IS!.I2+IS2.Ig+ISg.Ig+1+IS(g+1).Ib+ISb .US1 .US2 .USg .US(g+1).USb7/12/202421结束情况3 含受控电压源的复合支路受控电压源与无源元件串联，控制量可以是其它支路无源元件的电压或电流。在第十章，我们曾用受控源替代法分析过含有互感的电路。所以当支路含受控电压源时，可以仿照含互感的方法处理。但互感是成对出现的，而受控源可以单个出现。Z的非对角线元素将含有与控制系数有关的元素。其正负号的确定：控制量、被控量与复合支路的参考方向都一致(或都相反)时取“+”，否则取“-”。Zk-+.Usk.Isk.Ik-+.Uk+-.Udk .U=Z .(I+.IS).-US支路方程的矩阵形式仍然是：7/12/202422结束三、回路电流方程的矩阵形式用B表示的KVL：令 Zl=BZBT，则 Zl 称为回路阻抗矩阵。Zl的主对角线元素为自阻抗；非对角线元素为互阻抗。.BU=0 .U=Z .(I+.IS).-US将支路方程代入得：BZ.I .-BUS=0+BZ.IS 用B表示的KCL：.I=BT.Il 代入上式得回路电流方程的矩阵形式为：BZBT.I=.BUS-BZ.IS BZBT是 l 阶方阵。.BUS 和BZ都是l 阶列向量。7/12/202423结束四、回路电流方程的编写步骤 P401例15-1解：(1)作有向图，选树；(2)画基本回路电流，参考方向同连支电流；R1-+.US2.IS1jwL4R2jwL3jwC5112345(3)写基本回路矩阵B、支路阻抗矩阵和电压、电流列向量；B=-1 0 1 0 1 0 1 0 1 -1 Z=diag R1,R2,jwL3,jwL4,jwC51 .US=0 -US2 0 0 0T.IS=IS1 0 0 0 0T7/12/202424结束(4)求回路阻抗矩阵Zl=BZBT(5)并计算整理便得回路电流方程的矩阵形式。Zl=BZBTZl.I=.BUS-BZ.IS 将Zl、.US、.IS 代入式=-1 0 1 0 10 1 0 1 -1R1R2jwL3jwL4jwC51-101010101-1R1+jwL3+jwC51jwC51-jwC51-R2+jwL4+jwC51.Il1.Il2=.R1IS1 .-US2计算得7/12/202425结束15-5 结点电压方程的矩阵形式一、复合支路及其方程的矩阵形式情况1 无受控电流源、无耦合电感Yk-+.Usk.Isk.Iek.Ik-+.Uk.Idk与回路法定义的复合支路相比，增加了受控电流源。但不许存在受控电压源；也不许存在无伴电压源。对第 k 条支路有.Ik=Yk .(Uk+.USk).-ISk对整个电路有.I=Y .(U+.US).-ISY 称为支路导纳矩阵。Y 是一个对角矩阵，对角线元素为各支路导纳。7/12/202426结束情况2 无受控电流源、有耦合电感相当于回路法的情况2：VCR的矩阵形式与情况1相同。差别只是 Y 不再是对角阵。Yk-+.Usk.Isk.Iek.Ik-+.Uk .U=Z .(I+.IS).-US电路有耦合电感时，支路阻抗 Z 不是对角阵，在 Z 的非主对角线元素中将含有互感阻抗。利用上式可求得 .Y U=.I+.IS .-YUS.I=Y .(U .+US).-ISY=Z-1是支路导纳矩阵。7/12/202427结束 情况3 有受控电流源 设：第k条支路受控源受第 j条支路电流(或电压)控制。Yk-+.Usk.Isk.Ik-+.Uk.IdkYj-+.Usj.Isj.Iej.Ij-+.Uj+-.Uej.Idk=gkj Uej或.Idk=bkj Iej因为 .Iej=Yj Uej .Idk=bkj Yj Uej所以无论是流控还是压控，均化成VCCS，且控制系数用Ykj 表示：Ykj=gkjbkj Yj第k条支路方程为：.Ik=Yk .(Uk+.USk).+Idk .-ISk.Ik=Yk .(Uk+.USk)+Ykj .(Uj .+Usj).-ISk.Idk=Ykj .Uej=Ykj(.Uj+.USj)注意它们在复合支路中的方向。7/12/202428结束支路方程的形式同情况1。导纳矩阵Y不是对角阵，非主对角线包含与控制系数有关的元素，可以单个出现。.Ik=Yk .(Uk+.USk).(Uj+Ykj .+Usj).-ISk.Ij.Ik0Ykj.I=Y .(U .+US).-ISk行j列对 b 条支路有.I1.I2.Ib=Y1Y2.Yj.Yk.Yb000000000 .U1+.US1 .U2+.US2 .Uj+.USj .Uk+.USk .Ub+.USb-.ISj.ISk.IS1.IS2.ISb(在情况2中则是成对出现)。7/12/202429结束二、结点方程的矩阵形式描述结点与支路关联的矩阵是A。用A表示的KCL：用A表示的KVL：.A I=0，.U=AT Un支路方程：.I=Y .(U .+US).-IS用结点电压表示支路电流.I=YAT .Un .+YUS .-IS代入用A表示的KCL得 AYAT .Un .+AYUS .-AIS =0结点方程的矩阵形式为 AYAT .Un .=AIS .-AYUS令 Yn=AYAT,.Jn .=AIS .-AYUS则结点方程可以写为 Yn .Un .=JnYn称为结点导纳矩阵。.Jn是由独立源引起的注入结点的电流列向量。7/12/202430结束三、结点电压方程的编写步骤 P405例15-2(1)作有向图，选参考结点；(2)写关联矩阵A、独立电源列相量和支路导纳矩阵；L1R5R4iS4L2R3C6iS3123456A=1 0 1 1 0 0-1 1 0 0 0 1 0-1 0-1 1 0 .Us=0,.Is=0 0 IS3 IS4 0 0 TjwL11jwL21R31R41R51jwC6Y=diag ,AYAT .Un .=AIS .-AYUS(3)求AYAT并代入得到 AYAT .Un .=AIS 7/12/202431结束观察结点导纳矩阵发现 主对角线元素为自导纳，其余为互导纳。相当于第三章所列结点方程等号左边的系数。独立源列向量为注入结点的电流(等号右边的常数)。.Un1 .Un2 .Un3 .IS3+.IS4=0 .-IS4 AYAT .Un .AIS L1R5R4iS4L2R3C6iS3R31+R41+jwL11-jwL11-R41-jwL11jwL11+jwL21+jwC6-jwL21-R41-jwL21R41+R51+jwL217/12/202432结束P406例15-3 设写支路方程的矩阵形式。控制量、被控量与复合支路的参考方向都一致(或都相反)时取“+”，否则取“-”。.Id2=g21U1,.Id4=b46I6R11R21jwC3jwC4jwL51jwL61-g21jwL6b4600000Y=2行1列4行6列解：电路含受控源，但无互感。支路导纳矩阵为：L5R2C3L6R1iS1iS4C4id4-+.US4+-u1+-uS2id2i60+-u65614230注意位置和正负7/12/202433结束可得到独立电源列向量为：将以上所求代入.IS=.US=0.I=Y .(U .+US).-IS.IS1 0 0 .-IS4 0 0 T .-US2 0 .US4 0 0T与复合支路相反取正，否则取负。L5R2C3L6R1iS1iS4C4id4-+.US4+-u1+-uS2id2i60+-u65614230R11R21jwC3jwC4jwL51jwL61-g21jwL6b4600000=-.I2.I3.I1.I4.I5.I6 .U1+0 .U2+.US2 .U3+0 .U4+.US4 .U5+0 .U6+0.IS100 .-IS1007/12/202434结束15-6 割集电压方程的矩阵形式一、关于割集电压 割集电压也是一组完备的独立变量。以割集电压作为电路独立变量的分析方法称为割集电压法。当选单树支(基本)割集作为独立割集时，树支电压就是割集电压。割集电压是一种假想电压，就象假想的回路电流一样。u=Qf utT否则，ut 可理解为一组独立的割集电压。Qadfbce支路a、b、c、d是 G 的一个割集。将它们全部移去，G被分离成两个部分。Qfeut两分离部分之间的电压就是割集电压。7/12/202435结束可以认为，割集电压法是结点电压法的推广。也可以说结点电压法是割集电压法的一个特例。若选一组独立割集 使每一割集都由汇集在一个结点上的支路构成时，割集电压法就成为结点电压法。割集电压法规定的复合支路与结点电压法完全相同，因此支路方程的矩阵形式也相同：.I=Y .(U .+US).-IS12345678Q1Q2Q3Q47/12/202436结束二、割集电压方程描述割集与支路关联的矩阵是 Qf。用Qf 表示的KCL：代入用Qf 表示的KCL方程就得到割集方程的矩阵形式：.Qf I=0，KVL：.U=Qf UtT支路方程：.I=Y .(U .+US).-IS用割集电压表示支路电流.I=Y .Ut .+YUS .-ISTQf Qf YTQf .Ut .=Qf IS .-Qf YUS若令 Yt=Qf Y TQf则Yt 称为割集导纳矩阵。7/12/202437结束三、割集电压方程的编写解：初始条件为零，用运算形式。作有向图，选独立割集组；选支路1、2、3为树支，则单树支割集如图所示。12345Ut1(s)Ut2(s)Ut3(s)树支电压Ut1(s)、Ut2(s)和Ut3(s)就是割集电压，取树支电压的方向为割集方向。写基本割集矩阵 Qf；Qf=1 2 3 4 5Q1Q2Q31 0 00 1 00 0 11 1-1 01 1P408 例15-4L4R1R2L3C5iS1iS2先树支后连支7/12/202438结束写独立源列向量和支路导纳矩阵；US(s)=0 IS(s)=IS1(s)IS2(s)0 0 0T代入L4R1R2L3C5iS1iS2Y(s)=R11,R21,sL31,sL41,sC5 diag Qf Y(s)QfTUt(s)=Qf IS(s)-Qf Y(s)US(s)得割集电压方程Ut1(s)Ut2(s)Ut3(s)=IS1(s)IS2(s)0R11+sL41+sC5-sL41sL41+sC5-sL41R21+sL41-sL41sL41+sC5-sL41sL31+sL41+sC57/12/202439结束Yt 的主对角线元素分别为与割集Q1、Q2、Q3相关联支路的导纳之算术和，称自导纳。其它元素分别是与两相邻割集关联支路的导纳算术和，称互导纳。当两割集方向不一致时，互导纳前加“-”号。L4R1R2L3C5iS1iS2Ut1(s)Ut2(s)Ut3(s)与割集方向相反的割集电流源取“+”，相同取“-”。Ut1(s)Ut2(s)Ut3(s)=IS1(s)IS2(s)0R11+sL41+sC5-sL41sL41+sC5-sL41R21+sL41-sL41sL41+sC5-sL41sL31+sL41+sC57/12/202440结束状态方程的编写 在线性电路中，选独立的电容电压和电感电流作为状态变量列写状态方程和求解最方便。1.直观法的编写步骤 在状态方程中，要包含对状态变量的一次导数，所以：(1)对只含一个C的结点 列KCL方程；C+-R2L2i1uSR1+-uCL1i2iSi2+iS12结点 CdtduC=-i1-i2(2)对只含一个L的回路 列KVL 方程；回路1 L1dtdi1=uC-R1(i1+i2)+uS回路2 L2dtdi2=uC-R1(i1+i2)+uS-R2(i2+iS)(3)列其它方程(如有必要)，消去非状态变量。iR1iR2iR1iR27/12/202441结束(4)整理成矩阵形式直观法适用于不太复杂的电路。对复杂电路采用系统法CdtduC=-i1-i2L1dtdi1=uC-R1(i1+i2)+uSL2dtdi2=uC-R1(i1+i2)+uS-R2(i2+iS)dtduCdtdi1dtdi2=0-C1-C1L11-L1R1-L1R1L21-L2R1-L2R1+R2uCi1i2+0 0L110L21-L2R2uSiS7/12/202442结束2.系统法的编写步骤(1)选特有树；树支包含全部电容和电压源，不包含电感和电流源。(2)对单树支割集列KCL方程；(3)对单连支回路列KVL方程；(4)列其它必要的方程，消去非状态变量；(5)整理并写成矩阵形式。只要电路中不存在仅由电容和电压源构成的回路；也不存在仅由电感和电流源构成的割集。特有树就一定存在。7/12/202443结束例:列出图示电路的状态方程。解：uS1+-L7C2C30C4L8R6G5iS9选特有树：123764589对单树支割集列KCL方程：C2dtdu2=i7C3dtdu3=i6+i7C4dtdu4=i6+i8对单连支回路列KVL方程：L7dtdi7=-u2-u3L8dtdi8=-u4-u5再列两个方程，消去i6和u5：i6=R61u6=R61(-u4-u3+uS1)u5=G51i5=G51(i8+i9)7/12/202444结束整理令u2=x1，u3=x2，u4=x3，i7=x4，i8=x5。写成矩阵形式即可。dtdu2=C21i7dtdu3=C3R61-u3C3R61-u4+1C3i7+C3R61uS1dtdu4=C4R61-u3C4R61-u4+1C4i8C4R61+uS1dtdi7=-L71u2-L71u3dtdi8=-L81u4G5L81-i8G5L81-iS97/12/202445结束.x1.x2.x3.x4.x5=x1x2x3x4x5+uS1iS9dtdx1=C21x4000C21000dtdx2=C3R61-x2C3R61-x3+1C3x4+C3R61uS10C3R61-C3R61-C310C3R610dtdx3=C4R61-x2C4R61-x3+1C4x5C4R61+uS10C4R61-C4R61-0C41C4R610dtdx4=-L71x1-L71x2L71-L71-00000dtdx5=-L81x3G5L81-x5G5L81-iS900L81-0G5L81-0G5L81-uS1+-L7C2C30C4L8R6G5iS9这就是图示电路的状态方程。7/12/202446结束本章结束7/12/202447供娄浪颓蓝辣袄驹靴锯澜互慌仲写绎衰斡染圾明将呆则孰盆瘸砒腥悉漠堑脊髓灰质炎(讲课2019)脊髓灰质炎(讲课2019)供娄浪颓蓝辣袄驹靴锯澜互慌仲写绎衰斡染圾明将呆则孰盆瘸砒腥悉漠堑脊髓灰质炎(讲课2019)脊髓灰质炎(讲课2019)