本文介绍了如何在Windows环境下使用UMFPACK库解决大规模稀疏矩阵方程。UMFPACK是一个经过验证的库，通过编译提供的源代码或使用预编译的库文件，可以将其集成到工程中。文章详细阐述了矩阵的压缩表示方法，包括Ap、Ai和Ax数组的含义，并展示了简单的使用示例。在实际应用中，UMFPACK解算4611x4611矩阵方程仅需0.025ms，显示出良好的性能。
最近在对模型进行PRT处理时需要用到最小二乘法来做采样点的优化拟合，对最小二乘法有了解的可能就会知道（
http://en.wikipedia.org/wiki/Least_squares
），一般情况下为了便于目标点的计算通常会引入矩阵并得到一个如下的矩阵方程：
Ax = b
其中的x就是一个关于目标点的参数列向量，b是对应的观测值向量，而A就是根据应该背景得到的数据矩阵。求解这个矩阵形式的线性方程组就可以得到最小二乘法下的x优化估值。
Ax = b
其中的x就是一个关于目标点的参数列向量，b是对应的观测值向量，而A就是根据应该背景得到的数据矩阵。求解这个矩阵形式的线性方程组就可以得到最小二乘法下的x优化估值。
可能在大多数情况下得到的A会是一个规模较大且较为稀疏的矩阵。当然，这里可以自己写程序来求解，但能用经过验证的库岂不更好？UMPACK就是求解类似于Ax=b这样问题的一个库，来自佛罗里达州立大学。可以直接到http://www.cise.ufl.edu/research/sparse/umfpack/去下载对
最低0.47元/天 解锁文章
1499
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
