eigenvector相关内容

如何计算 2^n 模 1000000007 , n = 10^9

计算这个最快的方法是什么,我看到有些人使用矩阵,当我在互联网上搜索时,他们谈到了特征值和特征向量(不知道这个东西)......有一个问题减少了到递归方程f(n) = (2*f(n-1)) + 2 ,并且 f(1) = 1,n 可能高达 10^9....我已经尝试过使用 DP,存储多达 1000000 个值并使用常见的快速求幂方法,但都超时了我在这些需要计算大值的模数问题上通常很弱 解决方案 ..
发布时间:2021-12-19 08:19:43 其他开发

MATLAB eig 有时会返回反转符号

我正在尝试编写一个程序来获取任意大小的矩阵 A,然后 SVD 对其进行分解: A = U * S * V' 其中A是用户输入的矩阵,U是A * A'的特征向量组成的正交矩阵,S是奇异值的对角矩阵,V是A' * A的特征向量的正交矩阵. 问题是:MATLAB 函数eig 有时会返回错误的特征向量. 这是我的代码: 函数 [U,S,V]=badsvd(A)W=A*A';[U,S]=e ..
发布时间:2021-12-08 15:05:57 其他开发

Matlab eig 是否总是返回排序值?

我在 Matlab 中使用了一个函数: [V,D] = eig(C); 我看到 V 和 D 总是按升序排序.它总是这样还是应该在我获得 V 和 D 值后对它们进行排序? 解决方案 V 不按任何顺序排序,除了对应于关联特征值的顺序.但也许你不是那个意思. 特征值倾向于按降序排列,但这根本不能保证.它们往往是有序的,因为最大的往往会从最上面的算法中流出.Eig 最后没有排序来确保这一 ..
发布时间:2021-12-08 15:04:52 其他开发

在 scipy 中获取非标准化特征向量

Scipy 和 Numpy 返回归一化的特征向量.我正在尝试将向量用于物理应用程序,我需要它们不被标准化. 例如a = np.matrix('-3, 2; -1, 0')W,V = spl.eig(a) scipy 返回 [-2,-1] 的特征值 (W) 和模态矩阵 (V)(特征值作为列)[[ 0.89442719 0.70710678][ 0.4472136 0.70710678]] ..
发布时间:2021-07-16 20:59:35 其他开发

跟踪 1 参数矩阵族的特征向量

我的问题是:我正在尝试通过(截断的)Karhunen-Loeve 变换对随机过程进行谱分解,但我的协方差矩阵实际上是一个 1 参数矩阵族,我需要一种方法来估计/可视化我的随机过程如何取决于此参数.为此,我需要一种方法来跟踪 numpy.linalg.eigh() 生成的特征向量. 为了让您了解我的问题,这里有一个示例玩具问题:假设我有一组点 {xs},以及一个具有协方差 C(x,y) = 1 ..
发布时间:2021-06-10 19:36:28 Python

特征向量中心性numpy的不同结果

以下示例给出了使用 eigenvector_centrality 和 eigenvector_centrality_numpy 获得的不同结果.有没有办法让这种计算更加健壮?我正在使用 networkx 2.4、numpy 1.18.5 和 scipy 1.5.0. 将 numpy 导入为 np将 networkx 导入为 nx邻接矩阵 = {0:{1:0.6,},1:{2:0,3:0,},2: ..
发布时间:2021-06-07 19:54:03 其他开发

NumPy与MATLAB之间的特征向量存在细微差异

我有刚度矩阵和质量矩阵.我想计算我的结构振动形状和周期(特征值/矢量),所以我为此使用了NumPy.特征值与MATLAB给出的特征值相同,但是当我将特征向量与MATLAB给出的特征值进行比较时,发现一些小的差异(小于1E-10). 这是为什么,如何使两个结果相等? 我试图提高NumPy的精度,但是没有用. 将numpy导入为np#S刚度矩阵#M质量矩阵w,f = np.linalg. ..
发布时间:2021-05-31 21:01:20 Python

如何找到具有不同特征值的两个矩阵的公共特征向量

我正在寻找或更确切地说在两个矩阵 A 和 B 之间建立公共特征向量矩阵X,例如: AX = aX带有"a"特征值对应的对角矩阵BX = bX,其中"b"为0.特征值对应的对角矩阵 其中 A 和 B 是正方形和可对角线化的矩阵. 我看过类似的帖子,但没有设法得出结论,即当我构建由以下内容定义的最终想要的同构异形 F 时具有有效的结果: F = PDP ^ -1 我还阅读了维基百科 ..
发布时间:2021-05-29 20:53:54 其他开发

将Eigen :: MatrixXf转换为2D std :: vector

有没有比点对点复制值更优雅的解决方案了? 像这样的东西适用于一维矢量... vectorvec(mat.data(),mat.data()+ mat.rows()* mat.cols()); 我尝试了GCC编译器针对vector,但没有解决... 解决方案 Eigen :: MatrixXf 使用有效的线性内存,而 vector 为 vector ..
发布时间:2021-05-03 20:09:50 其他开发

按特征向量对特征向量排序(关联排序)

我有一个未分类的特征值向量和一个相关的特征向量矩阵.我想针对特征值的排序集对矩阵的列进行排序.(例如,如果eigenvalue [3]移至eigenvalue [2],我希望特征向量矩阵的第3列移至第2列.) 我知道我可以通过 std :: sort 在 O(N log N)中对特征值进行排序.不用滚动我自己的排序算法,如何确定矩阵的列(相关的特征向量)以及它们的特征值在排序时就遵循了? ..
发布时间:2021-04-19 19:31:01 C/C++开发

Zheev特征值和向量不正确

我在Fortran科学代码中使用LAPACK zheev例程来计算不太大(可能永远不会超过1000大小)的矩阵的特征值和向量. 由于此步骤发生在计算的开始,所以我必须获得很高的精度,以避免重要的错误传播.问题在于,在我的测试用例中(仅使用12x12矩阵),计算精度仅为1e-9左右,这根本不够. 我与numpy.linalg.eigh进行了比较,后者给出了可笑的更好结果,我想知道如何使用 ..
发布时间:2020-11-10 06:54:34 其他开发

R中的特征向量复数非对称矩阵不同于Matlab:如何解决这个问题?

我想计算R中复杂矩阵的特征值. 比较在MATLAB中获得的值,我没有从R中从完全相同的矩阵进行计算而获得的本征值相同. 在R中,我使用了eigen().在MATLAB中,使用eig()或eigs()(两个函数都给出相同的特征值,但不同于R的特征值). 对于真实矩阵,R和MATLAB是一致的. 如何在R中获得与MATLAB相同的结果? Matlab 实矩阵中的示例: sque ..
发布时间:2020-08-13 19:42:29 其他开发

计算R中的稀疏矩阵的特征向量

我正在尝试计算R中一个大型稀疏矩阵的m个第一特征向量.使用eigen()是不现实的,因为在这里,大型表示N> 10 6 . 到目前为止,我已经确定应该使用igraph软件包中的ARPACK,该软件包可以处理稀疏矩阵.但是我无法在一个非常简单的(3x3)矩阵上工作: library(Matrix) library(igraph) TestDiag ..
发布时间:2020-08-06 02:16:13 其他开发

当协方差矩阵为零时,如何在R中使用princomp()函数?

在R中使用princomp()函数时,遇到以下错误:"covariance matrix is not non-negative definite". 我认为,这是由于协方差矩阵中的某些值为零(实际上接近零,但在舍入过程中变为零)所致. 当协方差矩阵包含零时,是否可以进行PCA处理? [FYI:获得协方差矩阵是princomp()调用内的中间步骤.可以从此处下载用于重现此错误的数 ..
发布时间:2020-07-31 04:10:21 其他开发

计算R中前两个主成分的最快方法是什么?

我正在R中使用princomp来执行PCA.我的数据矩阵很大(10K x 10K,每个值最多4个小数点).在Xeon 2.27 GHz处理器上需要约3.5个小时和约6.5 GB的物理内存. 由于我只想要前两个组件,有没有更快的方法呢? 更新: 除了速度之外,还有一种高效的内存存储方式吗? 使用svd(,2,)计算头两个组件需要花费大约2个小时和大约6.3 GB的物理内存. ..
发布时间:2020-07-31 04:09:49 其他开发