c - C中的特征向量、特征值不动点计算

Question

** 已编辑 **

我尝试使用 libfixmath 将提到的 Jacobi 算法更改为定点算法但我没有得到正确的结果。我错过了什么？？

edited Jacobi code

makefile

---

这里是 C 新手。我不知何故陷入了深渊，找不到出路。如果您能提供帮助，那就太棒了!

情况:我正在尝试在 C 中实现 ICA 算法。我是使用浮点运算( double 、浮点)来实现的。现在我希望此代码将其转换为定点，以便我可以将其导入 32 位的 ARM 微 Controller (这就是我不能使用 double 、浮点等的原因)。

我找到了四个我认为可以帮助我的库:

1. http://sourceforge.net/projects/avrfix/files/
2. http://www.dsprelated.com/showcode/40.php
3. http://sourceforge.net/p/fixedptc/code/ci/default/tree/
4. 库修复数学

我没有使用 1. 2. 或 3。我目前正在尝试 libfixmath，因为几乎所有的计算都是用矩阵完成的。

我的问题是试图找到协方差矩阵(正对称 3x3 矩阵)的特征值和特征向量。我四处搜索了执行特征分解或 SVD 等的库或函数，但我没有找到任何东西。

你如何在定点中进行那种计算？？有没有我没发现的功能/库？我是否必须更改我在浮点中的特征函数(逐行转换为定点 - 即 fix16_from_dbl() )？

我当前的特征函数(当然不是我的，我认为它来自 Numerical Recipes)

jacobi.h

其他相关问题:这是 StackOverflow fixed SVD

****这是我的第一个问题，如果我的问题有任何需要更正的地方，请说出来....不要活着吃我! :)

** 已编辑 **

libfixmath 做 Cholesky 1 , 维基 2和QR分解。

我尝试玩数学游戏并使用此数据(来自上述函数的数据)生成矩阵的特征向量或特征值，但我失败了。

如果有人知道该怎么做，那么问题就解决了。

*我应该将其发布到数学堆栈网站吗？

1 http://rosettacode.org/wiki/Cholesky_decomposition

2 http://en.wikipedia.org/wiki/Cholesky_decomposition#The_Cholesky.E2.80.93Banachiewicz_and_Cholesky.E2.80.93Crout_algorithms

Answer 1

GNU Scientific Library怎么样？ ？它有一个 number of functions related to eigenvalue decomposition

你可以使用类似 IQMath 的东西为了从 float 转换为定点数。类似问题也有很多答案here