我有一个矩阵,例如这个例子(我的实际矩阵可以大得多)
A = [-1 -2 -0.5;
0 0.5 0;
0 0 -1];
只有两个线性无关的特征值(重复特征值-1).我想用广义特征向量获得完整的基础.我知道如何做到这一点的一种方法是使用Matlab jordan在Symbolic Math工具箱中的功能,但我更喜欢为数字输入设计的东西(实际上,有两个输出,jordan大型矩阵失败:"MuPAD命令出错:相似矩阵太大").我不需要Jordan规范形式,这在数字上下文中是非常不稳定的,只是一个广义特征向量的矩阵.是否有函数或函数组合以数字稳定的方式自动执行此操作,或者必须使用通用的手动方法(这样的过程有多稳定)？
注意:通过"广义特征向量",我指的是一个非零向量,可用于增加所谓的有缺陷矩阵的不完全基础.我并不是指使用或通过解决广义特征值问题而得到的特征值的特征向量(尽管后一种用法很常见,但我认为最好避免使用).除非有人能纠正我,否则我不相信这两者是一样的.eigqz
更新1 - 五个月后:
请参阅我的答案,了解如何为大于82乘82的矩阵(本问题中的测试矩阵的极限)符号化地获取广义特征向量.
我仍然对数字方案感兴趣(或者如果这些方案都与计算Jordan形式有关,那么这些方案可能会如何不稳定).我不希望盲目地实现已被标记为此问题的副本的线性代数101方法,因为它不是数值算法,而是用于评估学生的铅笔和纸张方法(我认为它可以实现但象征性地).如果有人能指出我对该方案的实施或对其进行数值分析,我会对此感兴趣.
更新2 - 2015年2月:在R2014b中测试的所有上述内容仍然是正确的.