使用 SAS IML 计算矩阵的对数
compute the logarithm of a matrix using SAS IML
如何使用 SAS IML 计算矩阵的对数？
我正在尝试获得类似于 Matlab 功能的东西 logm() function。
我不知道可以执行此操作的预写函数，因此您需要自己编写。然而，Nick Higham 在 his 'mftoolbox' set of functions 中提供了这个功能。 MATLAB 和 SAS/IML 非常相似，因此将 Higham 的代码从 MATLAB 转换为 IML 并不太困难。
有几种 logm() 算法。我推荐逆比例缩放和平方方法，它在 Higham 的 logm_iss.m 文件中实现。该函数调用您还需要实现的其他几个函数：
sqrtm_dbp(),
logm_pade_pf()
Higham 的几篇论文中提供了背景信息和示例，其中大部分可以从他的主页或 Google Scholar 获得。您还可以阅读 Higham (2008) 矩阵函数：理论与计算。
如何使用 SAS IML 计算矩阵的对数？
我正在尝试获得类似于 Matlab 功能的东西 logm() function。
我不知道可以执行此操作的预写函数，因此您需要自己编写。然而，Nick Higham 在 his 'mftoolbox' set of functions 中提供了这个功能。 MATLAB 和 SAS/IML 非常相似，因此将 Higham 的代码从 MATLAB 转换为 IML 并不太困难。
有几种 logm() 算法。我推荐逆比例缩放和平方方法，它在 Higham 的 logm_iss.m 文件中实现。该函数调用您还需要实现的其他几个函数： sqrtm_dbp(), logm_pade_pf()
Higham 的几篇论文中提供了背景信息和示例，其中大部分可以从他的主页或 Google Scholar 获得。您还可以阅读 Higham (2008) 矩阵函数：理论与计算。