行列の平方根（X^2 =AとなるX）をPython+SciPy(linalg.sqrtm)で計算。
ある行列Aがあって、A=X^2となるようなX,つまり行列の平方根の計算が必要となって調べてみた。
英語版のWikipediaに結構詳しく書いてある。
単純に成分の平方根を取ったらいいとかそんな簡単じゃなくて、
無限にあったり複数あったりなかったり、半正定値なら決まるとか、かなり面倒。
こういうのは既存のライブラリを使う方が簡単、ということでまずはmatlab。
sqrtmというのを使うようだ。
さらに同じsqrtmがSciPyでも使える。
どちらも同じアルゴリズム”Blocked Schur algorithms”を使っている。
こちらに説明があります。
ということで早速Pythonで計算してみよう。説明資料とmatlabの例を計算してみた。
scipy.linalg.sqrtmを使ってこんな感じで計算できる。
« Python+SymPy（ブロック行列）で2*NポートのSパラメータをTパラメータに変換（2） 4ポートの場合 | トップページ | Excel VBAで複素N×N行列（動的配列でNを決める）のさらに動的配列を作る方法（NポートSパラメータの周波数特性みたいな） »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年11月27日 Ericssonのレポートで5Gは2025年に26億台に、インテルとメディアテックが5Gで協業、iPhone11のスマートバッテリーケースのX線写真、Qorvoのミリ波フロントエンド、など。(2019.11.28)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年11月11日 Huawei Mate 30 Pro 5G分解、AppleのUWBチップU1解析、Qualcommの5G RFモジュール解説など。(2019.11.11)
- 又吉直樹のヘウレーカ 「世界最速！？量子コンピューターって何？」を見てました。リアルタイムツイートのまとめ。(2019.11.07)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年10月28日 AppleのUWBモジュール（U1チップ内蔵）解析、GSAのミリ波レポート、Infineonの60GHzレーダセンサの配線とか。(2019.10.28)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2019年10月24日 Qualcommが5G関連企業に200Mドル出資＆5Gの利点を プレゼン、放熱機構がすごいMicrosoft Surface Laptopや60GHzセンサが搭載のGoogle Pixel 4 分解(2019.10.24)
「日記・コラム・つぶやき」カテゴリの記事
- エボラ出血熱治療に関するクラウドファンディングに出資したら、エボラウイルスとインフルエンザウイルスのストラップいただいた。(2019.11.26)
- eetimesが選ぶ、ハロウィンに向けた面白電子工作プロジェクト10選。(2019.10.26)
- 道の真ん中を亀が歩いていた、、、(2019.08.13)
- ようやく探していたストロガッツ先生の「Infinite powers」買った！微分・積分の歴史から、ウサイン・ボルトやHIVとの闘い、など興味深い題材がたくさん。(2019.06.15)
- 「量子革命」を読んだ。ボーアがなぜあんなにアインシュタインとの論争で有名か、とか前期量子論だけじゃなく偉大な貢献をしたというのを知らなかったので勉強になる。ゾンマーフェルトが方位量子数と磁気量子数を導入したのも。(2019.06.02)
« Python+SymPy（ブロック行列）で2*NポートのSパラメータをTパラメータに変換（2） 4ポートの場合 | トップページ | Excel VBAで複素N×N行列（動的配列でNを決める）のさらに動的配列を作る方法（NポートSパラメータの周波数特性みたいな） »
コメント