電気回路理論/回路の行列表示と節点解析のソースを表示

この節では[[../グラフ理論の基礎|グラフ理論の基礎]]で学んだ事柄を用いて、電気回路の行列表示について説明し、また[[../節点解析法|節点解析法]]との関連について述べる。

== 枝電流と接続行列 ==
回路のある節点に接続する枝の集合はひとつのカットセットをなす。なぜならば、回路をその節点とその他の部分とに分離することができるからである。さらに、KCLよりこのカットセットの枝に流れる電流(枝電流)の和は0となる。

''回路図が無いため、しばらく記号で代用します''
<pre>
　　　　　A 　　　　→i1       B
  ┌────┬────R1────┐
 +│　　   +│  　　+ v1 -　　 +│ 
v ↑i    v2 R2↓i2  　　　　 v3 R3↓i3
 -│  　　 -│ 　　 　　　　　 -│
  └────┴─────────┘
          C
</pre>
たとえば上図のπ形回路で考えてみよう。この回路には3つの節点A、B、Cがある。ここで各節点に接続する枝は1つのカットセットをなし、各カットセットに流れる電流の和は0である。すなわち各節点にKCLを適用すれば、節点から流れ出る電流を正、節点へ流れ込む電流を負として、
:節点A: <math>i_1 + i_2 - i = 0</math>
:節点B: <math>i_3 - i_1 = 0</math>
:節点C: <math>i - i_2 - i_3 = 0</math>
これを行列を用いて書き直せば、次のように書くことができる。
:<math>\begin{pmatrix}
-1 &  1 &  1 &  0 \\
 0 & -1 &  0 &  1 \\
 1 &  0 & -1 & -1 
\end{pmatrix}\begin{pmatrix}
i \\ i_1 \\ i_2 \\ i_3
\end{pmatrix} = O</math>
左辺の行列
:<math>A = \begin{pmatrix}
-1 &  1 &  1 &  0 \\
 0 & -1 &  0 &  1 \\
 1 &  0 & -1 & -1 
\end{pmatrix}</math>
を'''接続行列'''(incidence matrix)という。接続行列の各行が節点に、各列が枝に対応し、''i''行''j''列の成分は''j''番目の枝が''i''番目の節点に接続している様子を表す。0であれば接続しておらず、1ならば電流が流れ出る向きに、-1ならば電流が流れ込む向きに接続していることを示している。また、ベクトル
:<math>\mathbf{i}_b = \begin{pmatrix}
i \\ i_1 \\ i_2 \\ i_3
\end{pmatrix}</math>
を枝電流ベクトルという。

ここで、接続行列の各行を加えると0になることに注意する。これは、この行列の[[線型代数学/行列の基本変形#階数|階数]]が3でない、すなわち先に立式した3つの式が独立ではないことを示している。しかし、Aの第3行を消去した行列
:<math>A' = \begin{pmatrix}
-1 &  1 &  1 &  0 \\
 0 & -1 &  0 &  1
\end{pmatrix}</math>
はこれ以上基本変形を施しても第2行あるいは第1行を消去できないので、この行列の階数は2である。すなわち、3本の式のうち2本が独立であり、2本の方程式のみを考えればよいことが分かる。このように接続行列から1行を消去した行列を'''既約接続行列'''(reduced incidence matrix)という。

一般に節点が''n''個ある場合、''n''個のカットセットを考えてその全てに対してKCLを適用することができる。しかし、それによって立式される''n''個の式は独立ではなく、任意の1つの式を省いた''n''-1個の式が独立となる。また、接続行列の階数は''n''-1となる。

そこで1つの節点を'''基準節点'''(datum node)とし、残りの''n''-1個の節点に対してKCLを適用させればよい。例えば節点Cを基準にすれば、節点Aと節点Bの2式のみ考えればよい。これは節点解析法においても重要な点である。

== 枝電圧と接続行列 ==
節点電位についてはどうだろうか。各節点の電位を<math>v_a, v_b, v_c</math>とすると、節点電位と枝電圧の関係として、
:<math>v = -v_a + v_c</math>
:<math>v_1 = v_a - v_b</math>
:<math>v_2 = v_a - v_c</math>
:<math>v_3 = v_b - v_c</math>
となる。これを行列を用いて書き表せば、
:<math>\begin{pmatrix}
v \\ v_1 \\ v_2 \\ v_3
\end{pmatrix} = \begin{pmatrix}
-1 &  0 &  1 \\
 1 & -1 &  0 \\
 1 &  0 & -1 \\
 0 &  1 & -1
\end{pmatrix}\begin{pmatrix}
v_a \\ v_b \\ v_c
\end{pmatrix}</math>
と書くことができる。右辺の係数行列は先に見た接続行列Aの転置行列となっている。ただし先に見たとおり、Aは1つ余計な行を含んでいる。そこでこの場合も節点Cを基準節点にとり、既約行列を用いて、
:<math>\begin{pmatrix}
v \\ v_1 \\ v_2 \\ v_3
\end{pmatrix} = \begin{pmatrix}
-1 &  0 & \\
 1 & -1 & \\
 1 &  0 & \\
 0 &  1 & 
\end{pmatrix}\begin{pmatrix}
v_a \\ v_b
\end{pmatrix}</math>
と書いてよく、このときこの式は
:<math>\mathbf{v}_b = A'^T \mathbf{v}_n</math>
と書くことができる。行列<math>\mathbf{v}_b</math>を枝電圧ベクトルといい、<math>\mathbf{v}_n</math>は節点電圧ベクトルである。

== 節点解析法 ==
電気回路の各枝は、次のように一般化して書くことができる。

''回路図が無いため、しばらく記号で代用します''
<pre>
　　　　　isk
ik┌───←───┐
→│　　　　+　-　│
─┴─Rk─-─vsk─┴─
+         vk         -
</pre>
''k''番目の枝に抵抗<math>R_k</math>、電流源<math>i_{sk}</math>、電圧源<math>v_{sk}</math>があるとし、このときの枝電流を<math>i_k</math>、枝電圧を<math>v_k</math>とする。もし枝に電源が含まれていなければこれらの値を0にとればよいし、あるいはたとえば電圧源のみで構成されているのであれば抵抗を0として考えればよい。

さて、この枝は[[../ノートンの定理|ノートンの定理]]によって次の等価電流源に置き換えることができる。

''回路図が無いため、しばらく記号で代用します''
<pre>
　　　　　jsk
ik┌───←───┐
→│　　　　　　　│
─┴───Rk───┴─
+         vk         -
</pre>
電源を除去した時の内部抵抗は<math>R_k</math>のみであるから、ノートン等価電流源の内部抵抗は<math>R_k</math>である。また等価電流源の電流を<math>j_{sk}</math>とすると、これは元の枝の両端を短絡した時に流れる電流に等しいから、
:<math>j_{sk} = i_{sk} - \frac{v_{sk}}{R_k}</math>
となる。

このモデルによって、枝電流を計算することができる。枝電流<math>i_k</math>は
:<math>i_k = \frac{v_k}{R_k} + j_{sk} = \frac{v_k}{R_k} + i_{sk} - \frac{v_{sk}}{R_k}</math>
となる。

さて、これを用いると、[[#枝電流と接続行列|前述]]の枝電流ベクトル<math>\mathbf{i}_b</math>は
:<math>\mathbf{i}_b = G\mathbf{v}_b + \mathbf{i}_{sb} - G\mathbf{v}_{sb}</math>　　(1)
と書くことができる。ここで''G''は
:<math>G = \begin{pmatrix}
G_1 & & & \\
& G_2 & & \\
& & \ddots & \\
& & & G_n 
\end{pmatrix} = \begin{pmatrix}
\frac{1}{R_1} & & & \\
& \frac{1}{R_2} & & \\
& & \ddots & \\
& & & \frac{1}{R_n}
\end{pmatrix}</math>
である。また、<math>\mathbf{i}_{sb}</math>と<math>\mathbf{v}_{sb}</math>それぞれ枝にある電流源と電圧源を縦に並べたベクトルである。

[[#枝電流と接続行列|枝電流と接続行列]]で見たように、枝電流と接続行列との間には
:<math>A\mathbf{i}_b = O</math>
の関係がある。したがって、(1)式より
:<math>AG\mathbf{v}_b + A\mathbf{i}_{sb} - AG\mathbf{v}_{sb} = O</math>　　(2)
が導かれる。

[[#枝電圧と接続行列|枝電圧と接続行列]]で見たように、枝電圧と節点電位とは
:<math>\mathbf{v}_b = A^T\mathbf{v}_n</math>
の関係があるから、この式の両辺を入れ替えて左から''AG''をかけると、
:<math>AGA^T\mathbf{v}_n = AG\mathbf{v}_b</math>
となり、(2)式とあわせると、
:<math>AGA^T\mathbf{v}_n = AG\mathbf{v}_{sb} - A\mathbf{i}_{sb}</math>
となる。ふつうは<math>AGA^T = Y_n</math>とおいて、
:<math>Y_n\mathbf{v}_n = AG\mathbf{v}_{sb} - A\mathbf{i}_{sb}</math>
とする。さらに
:<math>\mathbf{j}_{sb} = AG\mathbf{v}_{sb} - A\mathbf{i}_{sb}</math>
と書くことにすれば、この方程式は
:<math>Y_n\mathbf{v}_n = \mathbf{j}_{sb}</math>
と書くことができる。これを'''節点方程式'''(node equations)という。いま''A''は既約接続行列を用いている。<math>Y_n</math>はコンダクタンス行列と呼ばれる''n''-1次の正方行列であり、また<math>\mathbf{j}_{sb}</math>は節点電流源ベクトルと呼ばれる''n''-1次のベクトルである。

節点電流源ベクトルは各節点へ流れ込む電流源を表している。ただし、節点へ流れ込む方向を正、流れ出る方向が負であることに注意する。また、コンダクタンス行列は次のように書かれる。
* 対角成分は''k''番目の節点に接続するすべての枝のコンダクタンスの和となる。これを'''自己コンダクタンス'''(self-conductance)という。
* それ以外の成分は''k''番目の節点と''l''番目の節点を結ぶ枝にあるコンダクタンスの和に負の符号をつけたものになる。これを'''相互コンダクタンス'''(mutual-conductance)という。

したがって[[../節点解析法|前節]]で見たようにして機械的に節点方程式を立てることができる。

[[カテゴリ:電気工学|かいろのきようれつひようしとせつてんかいせき]]