Python+SymPy（ブロック行列）で2*NポートのSパラメータをTパラメータに変換（2） 4ポートの場合
先日は2ポートの計算でしたが、こちらが本題。
4ポートの場合のSパラメータ・Tパラメータ変換を以下のように考える。
（ポート1とポート2、ポート3とポート4が組の場合）
これはブロック行列の計算が必要だが、なんとSymPyはその計算ができる。
4x4行列を2x2行列4個に分解して計算を進める。
MatrixSymbok, BlockMatrix, block_collapseを使う。こんな感じ。
from sympy import *
init_printing()
SA=MatrixSymbol('SA',2,2)
SB=MatrixSymbol('SB',2,2)
SC=MatrixSymbol('SC',2,2)
SD=MatrixSymbol('SD',2,2)
I=Identity(2)
O=ZeroMatrix(2,2)
A=BlockMatrix([[I,-SA],[O,-SC]])
B=BlockMatrix([[SB,O],[SD,-I]])
block_collapse(Inverse(A)*B)
すると、結果はこうなる。
おお、完全にあってる。
2*Nポートも全く同様にできます。
案外、この計算している文献が少ない、、、
« Python+SymPyで2*NポートのSパラメータをTパラメータに変換（1） 2ポートの場合 | トップページ | 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス）を生むThomas-Rössler方程式のパラメータbを色々変えて、Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP！ インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
「日記・コラム・つぶやき」カテゴリの記事
- ExcelのLET関数＋SEQUENCE関数で数値計算シリーズ（その5）ワンライナーでラマヌジャンの円周率公式を使って計算する。n=2で収束した。(2021.03.01)
- mRNAワクチンやCRISPR cas9のことが全然わからんので昔買ったブルーバックスの「アメリカ版 大学生物学の教科書」を3巻まで読む、、、と思ったら11年ぶりに完全改訂されたものが出てた！(2021.03.02)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（2/28更新）さすがにどの国も伸びは鈍ったか。(2021.03.01)
« Python+SymPyで2*NポートのSパラメータをTパラメータに変換（1） 2ポートの場合 | トップページ | 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 »
コメント