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)で計算。 »
「学問・資格」カテゴリの記事
- 熱力学の関数の関係を覚えるBorn Square （まあ熱力学のはじき（きはじ)?) もちろんMax Bornはこんなの使わずに導出できただろうが、、、(2021.04.29)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年4月19日 IEEE Microwave Magazineの記事で細胞の特徴づけにマイクロ波、Microwave Journalの特集は機械学習を使ったアンプ設計、QorvoがBAW技術を新型コロナウイルスの抗原検査に、など。(2021.04.19)
- 2つのSパラメータのカスケード接続による合成をWolframAlphaに手伝ってもらって計算した。(2021.04.14)
- iPhoneで使える数学ソフトMaple Calculatorで手書き入力でMathematicaの常微分方程式の例を試す（その１）(2021.03.29)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年3月16日 IEEE Microwave magazineの特集は中国のマイクロ波技術、Microwave Journalの特集はCATRでミリ波ビームスキャンを測定、SIJの特集はなんと450ドルのネットアナ、フリーの高周波シミュレータ（Qucsなど）。(2021.03.16)
「日記・コラム・つぶやき」カテゴリの記事
- 地球と月の間の万有引力より太陽と月の間の万有引力の方が強いのか！知らなかった、、、(2021.05.05)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（5/2更新） 日本は第4波が進行中、他の国もイギリスを除いて増加している。 (2021.05.03)
- AppleのUWB利用を忘れものタグ（位置タグ) AirTagの分解（iFixitの写真を解説、5/4追記）(2021.05.02)
- 万有引力と重力、そしてそのエネルギーをちゃんと地球のスケールでグラフ化したらわかるのでは（ひろゆきが位置エネルギーが存在しないといったことに対して）(2021.04.30)
- 新型コロナウイルス、中国、日本、韓国、アメリカ、ドイツ、フランス、イギリスでの感染者数を指数関数＆ロジスティック関数＆Log-Logプロットでべき関数フィッティングした（4/25更新） 日本は完全に第4波なので色を変えて表示。フランスもワクチン接種進んでいるはずなのに全然減らない。(2021.04.26)
« Python+SymPyで2*NポートのSパラメータをTパラメータに変換（1） 2ポートの場合 | トップページ | 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 »
コメント