MATLAB的答案

最小非零特征值的广义特征值问题

18视图(30天)
Zoltan Csati
Zoltan Csati 2018年4月27日
回答: 安德鲁·克尼亚杰夫 2018年8月12日
我有两个矩阵, 一个 和 B ,我想解决广义特征值问题 斧头 =λ* Bx 。事实上我只需要 最小非零特征值 。矩阵的性质:
一个 是对称的,奇异与已知无效(但没有先验已知的内核),稀疏的吗
B 是对称的,奇异,半正定与已知的内核,稀疏,甚至线性无关的部分是坏脾气的
最小的非零特征值(由于病态)将数值结果类似1 e15汽油。然而,我知道我感兴趣的特征值不舍入高原附近。如果我知道这个值约,我可以使用
λ= eigs (A, B, k,猜);
k是特征值的数量我请求和猜接近最小的非零特征值我正在寻找。
因为我没有猜的信息,目前我转换 一个 和 B 完整的矩阵和调用eig:
λ= eig(完整(A)、(B));
然而,这是非常缓慢的。什么好主意吗?

在回答这个问题。

接受的答案

安德鲁·克尼亚杰夫
安德鲁·克尼亚杰夫 2018年8月12日
从矩阵A和B都是单一的,它不是一个简单的数值问题。甚至eig(完整(A)、(B))可能会给你错误的答案。叫EIGS EIGS (A, B, k,猜)不是可能的工作,因为它依赖于计算矩阵的逆,猜B,这是单一的在你的情况中。
如果B的内核是内核的一个子集,和维度的差异不是太大,你可能想尝试//www.tatmou.com/matlabcentral/fileexchange/48-lobpcg-m约束blockVectorY是向量的集合生成已知的k B . LOBPCG给你最小的内核,无论0,不信,不过,所以您可能需要选择k足够大到你需要的地方。LOBPCG k的值是线性无关的初始近似特征向量的数量。如果k太大,相对于矩阵的大小,它不会工作得很好…

更多的答案(0)

类别

找到更多的在线性代数在帮助中心和文件交换

标签

下载188bet金宝搏

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!