3
SAS IMLを使用して行列の対数を計算するにはどうすればよいですか？SAS IMLを使用して行列の対数を計算する
私はMatlab logm() functionの機能に類似したものを得ようとしています。
A
答えて
1
私はこれを行う事前に書かれた機能を知らないので、あなた自身でそれを書く必要があります。しかし、Nick Highamはこの機能をhis 'mftoolbox' set of functionsに提供しています。 MATLABとSAS/IMLは非常に似ているので、HighamのコードをMATLABからIMLに変換することはそれほど難しくありません。
いくつかのlogm（）アルゴリズムがあります。 Highamのlogm_iss.mファイルに実装されている逆スケーリングと二乗法をお勧めします。この関数はあなたにも実装する必要があるいくつかの他の関数を呼び出します。 sqrtm_dbp（）、 logm_pade_pf（）
背景情報と例が彼のホームページからやから入手できますほとんどがハイアムの論文のいくつかで提供されていますGoogle Scholarまた、Higham（2008）行列の関数：理論と計算を読むことができます。
これを行う方法が組み込まれているかどうかわかりません。 IMLMLIBライブラリ（[EXPMATRIX]（https://support.sas.com/documentation/cdl/en/imlug/68150/HTML/default/viewer.htm#imlug_langref_sect127.htm））に逆の実装がありますが、 LOGMATRIXに相当するものはありません。 Rickはこれを見て、ネイティブにやる方法や、アルゴリズムを実装した他の人の論文などを、より良い方向に伝えることができることを願っています。 – Joe
あなたは実際の行列を持っていると仮定しますか？それは対称的ですか？ – Rick
こんにちはリック、それは現実で対称ではありません。 – jaamor