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駆け出しエンジニアの豆知識: ⑩パワーアンプの近くにある電源系のICはエンベロープトラッカー。ETモードとAPTモードがある。(2021.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年9月15日 iPhone13の内部構造、IEEE Microwave magazineの特集はマイクロ波と機械学習、Microwave Journalは車載へのUWB応用、Huaweiは2030年に6Gを狙う、など。(2021.09.16)
- 高周波回路シミュレータQucsStudioを使ってみる（その４） 電磁界シミュレーション（FDTD）機能を使ってKeysight ADSのMomentumの例題のフィルタ（スタブLPF)を計算。(2021.09.13)
- コラッツ予想を複素数に拡張したものをPythonで図示してみる。(2021.09.07)
- コラッツ予想(Collatz conjecture)の図示をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! 証明出来たら1億円だそうですが、過去の研究結果とかタオさんの結果とか知ってからやった方がいいと思います。(2021.09.05)
「日記・コラム・つぶやき」カテゴリの記事
- スマホRF駆け出しエンジニアの豆知識: ⑩パワーアンプの近くにある電源系のICはエンベロープトラッカー。ETモードとAPTモードがある。(2021.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年9月15日 iPhone13の内部構造、IEEE Microwave magazineの特集はマイクロ波と機械学習、Microwave Journalは車載へのUWB応用、Huaweiは2030年に6Gを狙う、など。(2021.09.16)
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（9/12更新）。先週より減ったといっても先週がひどすぎるので全然減速して見えない。韓国もコンスタントに上がってる。(2021.09.12)
- 高周波回路シミュレータQucsStudioを使ってみる（その４） 電磁界シミュレーション（FDTD）機能を使ってKeysight ADSのMomentumの例題のフィルタ（スタブLPF)を計算。(2021.09.13)
- タカラトミーの電池なし、手で転がしてキャパシタを充電して走るテコロチャージのプラレールが発売される。過去の特許(チョロQ、Q-eyes)を見てみる。電気二重層キャパシタを使ってるっぽい？(2021.09.08)
« Python+SymPyで2*NポートのSパラメータをTパラメータに変換（1） 2ポートの場合 | トップページ | 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 »
コメント