行列の平方根（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駆け出しエンジニアの豆知識: ⑩パワーアンプの近くにある電源系の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パラメータに変換（2） 4ポートの場合 | トップページ | Excel VBAで複素N×N行列（動的配列でNを決める）のさらに動的配列を作る方法（NポートSパラメータの周波数特性みたいな） »
コメント