矩阵论在电路网络分析中的应用
摘要:电路网络分析中,运用矩阵论的相关知识可以直观的解决一些复杂问题,
比如所在支路存在无伴电压源的情况,而且矩阵运算方便进行计算机算法,在解决含大量节点的电路时是人工计算无法比拟的。若电路中存在无伴电压源支路时,由于该支路的导纳为无穷大,这给节点电压方程和割集方程的建立带来困难。解决这一问题的方法之一是将无伴电压源的支路电流也作为网络变量。因此,在改进的节点方程中是以节点电压和某些支路电流作为未知量。所述的支路电流包括无伴电压源支路电流和直接求解的支路电流。
改进节点法
将网络的支路划分为三类,一类是一般支路,另两类是无伴电压源支路和直接求电流的支路。后两类支路都可以以二端元件作为一条支路,支路电压和支路电流选择关联参考方向。
网络中的支路编号按照一般支路、无伴电压源支路和直接求电流支路,可将网络的关联矩阵A写成如下分块矩阵形式:
A A0
式中
AE
Ax
A
是反映一般支路与节点之间的关联关系的子阵。E是反映无伴电压源支
A
路与节点之间的关联关系子阵。x是反映直接求电流支路与节点之间关联关系子阵。
将支路电流向量和支路电压向量也按同样的顺序分块:
A0
Ib(s) I0(s)Ub(s) U0(s)
根据基尔霍夫电流定律,有
IE(s)UE(s)
Ix(s) Ux(s)
T
T
A0
AE
I0(s)
0Ax I(s) E
Ix(s)
A0I0(s) AEIE(s) AxIx(s) 0
根据基尔霍夫电压定律,有
U0(s)
U(s) A E 0 U0(s)
AE
Ax Un(s)
T
一般支路、无伴电压源支路和直接求电流支路的电流电压关系方程分别为
I0(s) Y0(s)U0(s) Yo(s)Us(s) Is(s)UE(s) USE(s)Ix(s) Yx(s)Ux(s)
将基尔霍夫电压方程带入得
I0(s) Y0(s)A0TUn(s) Y0Us(s) Is(s)
Ix(s) Yx(s)AxTUn(s)Yn0(s) A0Y0(s)A0T
将以上方程式列写为矩阵形式为
Yn0(s) AT
E
T
Y(s)Ax x
AE00
Ax Un(s) In0(s)
I(s) U(s) 0 E SE E Ix(s) 0
该式即为改进节点方程的一般形式,改进的节点法是以增加网络变量数为代价,避开了写无伴电压源支路的支路导纳。设网络有N+1个节点、p个无伴电压源支路和r个直接求电流源支路,则改进的节点方程的网络变量数为N+p+r个,系数矩阵为N+p+r阶方阵。虽然系数矩阵的维数增加了,但矩阵是稀疏的,利用稀疏矩阵(非零元素占全部元素的百分比很小(例如5%以下)的矩阵。有的矩阵非零元素占全部元素的百分比较大(例如近50%),但它们的分布很有规律,利用这一特点可以避免存放零元素或避免对这些零元素进行运算,这种矩阵仍可称为稀疏矩阵)技术计算仍很方便。
这种改进的节点电压矩阵可以应用到实际电路中,例如一个三阶电路含有两个电容一个电感,b1支路含有有伴电压源电导G1,b6、b7为直接求电容C1、C2支路,b2是电感L支路,b3是一个单纯电导G2支路,b4是电流源和电导G3的并联支路,b5是无伴电压源支路,以节点⑤为参考节点它的拓扑结构图可以如图所示:
1
5
列出关联矩阵
1 0A
0 0
0100
001 1
0001
00 10
1 100
0 1 1 0
前四列为A0阵,第五列为AE阵,后两列为直接求解阵Ax。 一般支路的导纳矩阵Y0(s)为对称阵即
G1 0
Y0(s)
0 0
01
SL00
00G20
0 0 0 G3
直接求电流支路的支路导纳矩阵Yx(s)为
SC1
Yx(s)
0
因为Yn0(s) A0Y0(s)A0T
0
SC2
G1
0T
Yn0(s) A0Y0(s)A0
0 0
01
sL00
00G2 G2
0
G2
G2 G3
00 sC1 sC1
Yx(s)Ax
0sC sC0 22
T
节点电流注入为
In0(s) GUIs(s) 1s(s)00
改进后的节点矩阵可写为
T
0000 G1
1 0000
sL 00G2 G2 1
0 G2G2 G30 0 00100
000 sC1 sC1
00 0sC2 sC2
求出所求支路的电流量。
0
1s(s) Un1(s) GU
0 11 Un2(s)
Un3(s) 0 0 1
Un4(s) Is(s)
00
IE(s) USE(s) 00
I(s) 0
10 c1
0 Ic2(s) 0 1
1
此方法可以有效地避开无伴电压源,支路导纳无穷大的病态情况,最终可以对于此电路的拓扑结构可以利用不定导纳矩阵对端子③,④进行隐藏。定义e为无论多么小的正数, e>0。 其不定导纳矩阵
G1 sC1
sC1 Yi 0
0
G
1
sC1
sC1 sC2
sC20 1
sL
1sL
0 sC2sC2 G2
G2 1
1
00
1
sL
1
G3
11
G1 G3
sL
G2G2 G3 G3
由于端子⑤接地,将第五行、第五列划去。得到Yid
G1 sC1
sC1 Yid
0 0
sC1
sC1 sC2
sC20
1sL
0 sC2sC2 G2
G2
1
00
G2G2 G3
由节点电流方程
Yid(s)Un(s) In(s)
Y11Y12
将Yid矩阵分块
YY 2122
Y11Y12 U I
YY U I 2122 2 2
I2 0
Y21U Y22U2 0,U2 Y22 1Y21U (Y11 Y12Y22 1Y21)U I
1
所以:I (Y11 Y12Y22Y21)U
故消去了节点③、④,得到只含节点①、 …… 此处隐藏:489字,全部文档内容请下载后查看。喜欢就下载吧 ……