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ニュース 2025年7月17日 IEEE Microwave Magazineで地上監視レーダ解説など、Microwave Journalでミリ波AiPのOTA測定解説など、5G Americasが5G Advancedホワイトペーパー発行、TDKのAI 微細欠陥検出edgeRX Vision、IEEE Trans. on THz の今月号発行(2025.07.17)
- Xのタイムラインが何かランダムウォークのガウス分布じゃなくて逆正弦法則のように見えてきている…仕方ないのでJulia言語でモンテカルロシミュレーションで試行回数を変えて逆正弦法則とガウス分布を見る。(2025.07.16)
- 高周波・RFニュース 2025年7月13日 Pythonの高周波ライブラリscikit-rfがv1.8.0に、SamsungがGalaxy Z Fold7など発表→QualcommがSnapdragon 8 Eliteが使われていると発表、NGMNが基地局アンテナの推奨事項をまとめる、STMicroとMetalenzがメタサーフェス光学のライセンス締結(2025.07.14)
- Interface2025年8月号Pythonで体験！はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10／10、RCR wireless newsのウェビナー2件（6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2025年7月17日 IEEE Microwave Magazineで地上監視レーダ解説など、Microwave Journalでミリ波AiPのOTA測定解説など、5G Americasが5G Advancedホワイトペーパー発行、TDKのAI 微細欠陥検出edgeRX Vision、IEEE Trans. on THz の今月号発行(2025.07.17)
- Xのタイムラインが何かランダムウォークのガウス分布じゃなくて逆正弦法則のように見えてきている…仕方ないのでJulia言語でモンテカルロシミュレーションで試行回数を変えて逆正弦法則とガウス分布を見る。(2025.07.16)
- 高周波・RFニュース 2025年7月13日 Pythonの高周波ライブラリscikit-rfがv1.8.0に、SamsungがGalaxy Z Fold7など発表→QualcommがSnapdragon 8 Eliteが使われていると発表、NGMNが基地局アンテナの推奨事項をまとめる、STMicroとMetalenzがメタサーフェス光学のライセンス締結(2025.07.14)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10／10、RCR wireless newsのウェビナー2件（6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
- 高周波・RFニュース 2025年7月2日 5G Americasが6Gに向けセンシングと通信ホワイトペーパー発行、KYOCERA AVXが3dBハイブリッドカプラ発表、TDKが車載薄膜インダクタ発表、Nordicが1次電池向けPMIC発表、ローデ・シュワルツの6GとAI／ML解説記事(2025.07.02)
コメント
« Python+SymPyで2*NポートのSパラメータをTパラメータに変換（1） 2ポートの場合 | トップページ | 行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。 »
参考文献を教えていただけますか？
S-->Tの変換は２端子対回路の説明ばかりで、それ以上だとどうやるのか知りたいです。
行列の変数(a1,a2,a3,a4,b1,b2,b3,b4)が変換後は入れ替わっていますが、ここを任意に入れ替えた場合も同じ計算処理でしょうか。
投稿: ひろつぐ | 2022年4月 8日 (金) 22時19分