行列の平方根（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パラメータの周波数特性みたいな） »
「学問・資格」カテゴリの記事
- 熱力学の関数の関係を覚える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パラメータに変換（2） 4ポートの場合 | トップページ | Excel VBAで複素N×N行列（動的配列でNを決める）のさらに動的配列を作る方法（NポートSパラメータの周波数特性みたいな） »
コメント