使用SAS IML计算矩阵的对数

Question

如何使用SAS IML计算矩阵的对数？

我想尝试类似于Matlab logm() function的功能。

Answer 1

我不知道会有这样做的预先编写的功能，因此您需要自己编写。但是，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学术搜索获得。您还可以阅读Higham（2008）矩阵函数：理论与计算。