如何使用 SAS IML 计算矩阵的对数？
我正在尝试获得类似于 Matlab logm() 函数的功能。
问问题
148 次
1 回答
1
我不知道会执行此操作的预先编写的函数，因此您需要自己编写它。但是，Nick Higham 在他的“mftoolbox”函数集中提供了这个函数。MATLAB 和 SAS/IML 非常相似，因此将 Higham 的代码从 MATLAB 转换为 IML 并不太难。
有几种 logm() 算法。我推荐在 Higham 的 logm_iss.m 文件中实现的逆缩放和平方方法。该函数调用您还需要实现的其他几个函数：sqrtm_dbp()、logm_pade_pf()
背景信息和示例在 Higham 的几篇论文中提供，其中大部分可从他的主页或 Google Scholar 获得。您还可以阅读 Higham (2008)矩阵函数：理论和计算。
于 2016-05-09T13:50:59.143 回答