- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
以前在这里问过类似的问题,但似乎没有人回答我的例子。我使用 Mathematica 和 SciPy 计算矩阵 A 的特征值和特征向量;特征值一致,但特征向量并非如此:
(1) 最低(特征值)特征向量一致
(2) Mathematica 和SciPy 的其余对应特征向量不存在乘法关系
(3) 我可以使用外积计算将 SciPy 的特征向量发送到 Mathematica 的对应特征向量的变换矩阵 T
T = numpy.outer(MathematicaEigenvector, SciPyEigenvector)
这样
MathematicaEigenvector = numpy.dot(T, SciPyEigenvector)
我希望所有 SciPy-Mathematica 特征向量对的变换矩阵 T 都相同,因为T 只是将矩阵 inv(T).A.T 的特征向量与原始矩阵 A 的特征向量相关联的矩阵。然而,对每个特征向量对执行步骤 (2) 会给出不同的 T 矩阵。
谁能解释一下?如果需要,我可以发布矩阵。
更新:python代码和矩阵如下:
S = [[0., -1, -1, -1, 0, 0, -1, 0, 0],
[-1, 0., -1, 0, -1, 0, 0, -1, 0],
[-1, -1, 0., 0, 0, -1, 0, 0, -1],
[-1, 0, 0, 0., -1, -1, -1, 0, 0],
[0, -1, 0, -1, 0., -1, 0, -1, 0],
[0, 0, -1, -1, -1, 0, 0, 0, -1],
[-1, 0, 0, -1, 0, 0, 0., -1, -1],
[0, -1, 0, 0, -1, 0, -1, 0., -1],
[0, 0, -1, 0, 0, -1, -1, -1, 0.]];
eig_val,eig_vec = scipy.linalg.eig(S)
idx = eig_val.argsort()
eig_val = np.array(eig_val[idx])
eig_vec = np.array(eig_vec[:,idx])
Mathematica 特征向量是:
[-0.333333, -0.333333, -0.333333, -0.333333, -0.333333, -0.333333, -0.333333, -0.333333, -0.333333],
[0.0385464, 0.570914, 0.371276, -0.570914, -0.0385464, -0.238184, -0.33273, 0.199638, 0.],
[0.570246, -0.0269007, 0.197029, 0.0269007, -0.570246, -0.346316, 0.373217, -0.22393, 0.],
[-0.0816497, 0.0816497, -0.489898, -0.0816497, 0.0816497, -0.489898, 0.408248, 0.571548, 0.],
[-0.333333, -0.333333, 0.166667, -0.333333, -0.333333, 0.166667, 0.166667, 0.166667, 0.666667],
[-0.288675, 0.288675, 2.498e-16, -0.288675, 0.288675, -1.94289e-16, 0.57735, -0.57735, 0.],
[-0.5, 0.5, -2.04678e-16, 0.5, -0.5, 2.41686e-16, -9.25186e-17, 5.55112e-17, 0.],
[0.166667, 0.166667, -0.333333, 0.166667, 0.166667, -0.333333, -0.333333, -0.333333, 0.666667],
[0.288675, 0.288675, -0.57735, -0.288675, -0.288675, 0.57735, 4.02456e-16, -2.08167e-16, 0.]
而 SciPy 特征向量是:
[-0.33333333 -0.33333333 -0.33333333 -0.33333333 -0.33333333 -0.33333333 -0.33333333 -0.33333333 -0.33333333]
[ 0.12054181 -0.17813781 0.50013951 0.08577902 -0.21290061 0.4653767 -0.2872389 -0.58591853 0.0923588 ]
[ 0.12191583 -0.21327897 0.26215377 -0.28683603 -0.62203084 -0.1465981 0.35987707 0.02468226 0.500115 ]
[ 0.66666667 0.16666667 0.16666667 0.16666667 -0.33333333 -0.33333333 0.16666667 -0.33333333 -0.33333333]
[-0.16604424 -0.59504716 -0.43689399 0.43294845 0.00394553 0.16209871 0.43294845 0.00394553 0.16209871]
[-0.01305419 0.07446538 -0.0614112 -0.54881726 0.36347168 0.18534558 0.56187145 -0.43793706 -0.12393438]
[-0.66666667 0.33333333 0.33333333 0.33333333 -0.16666667 -0.16666667 0.33333333 -0.16666667 -0.16666667]
[-0.21052033 0.65306873 -0.4425484 0.10526016 -0.32653437 0.2212742 0.10526016 -0.32653437 0.2212742 ]
[-0.02303417 0.0714558 -0.04842162 0.09679298 0.41311466 -0.50990763 -0.0737588 -0.48457045 0.55832926]
[ 4.67737437 0.12612917 0.75157798 -0.09378424 0.91674876 2.36234989 1.03706802 -9.0725069 0. ]
以上均按特征值排序 [-4.+0.j, -1.+0.j, -1.+0.j, -1.+0.j, -1.+0.j, 2.+0.j, 2. +0.j, 2.+0.j, 2.+0.j]
最佳答案
我认为原因如下:因为存在重复的特征值,所以变换矩阵 T 必须作用于该子空间中特征向量的线性组合,而不是单个特征值。也就是说,我的第一个代码片段应该修改为:
T = numpy.outer(MathematicaEigenvectorSubspace, SciPyEigenvectorSubspace)
我还没有通过找到使两个子空间等价的线性组合来检查这是否明确有效。
关于python - 特征向量变换的差异 : Mathematica vs. SciPy,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/25723478/
设置 我希望能够定义一个特征,使得任何实现该特征的结构不仅必须实现函数,而且还必须为某些常量指定值。所以也许是这样的: trait MyTrait { const MY_CONST: u8;
在我的 Web 应用程序中,授权用户至少有 4 个“方面”:http session 相关数据、持久数据、facebook 数据、运行时业务数据。 我决定使用案例类组合而不是特征至少有两个原因: 性状
我正在尝试使用以下代码从类中获取完整数据成员的列表: import std.stdio; import std.traits; class D { static string[] integr
我正在尝试实现 From对于我的一种类型。它应该消耗任意长度的行(仅在运行时已知)并从行中获取数据。编译器提示 &[&str; 2]不是 &[&str] ,即它不能将固定大小的切片转换为任意长度的切片
有人可以请你这么好心,并指出一种提取拟合树中使用的列/特征的方法,使用如下代码: library(dplyr) library(caret) library(rpart) df % dplyr
假设我定义了一个 Group所有组操作的特征。是否可以创建一个包装器AGroup超过 Group无需手动派生所有操作? 基本上,我想要这个: #[derive (Copy, Debug, Clone,
最近浏览了Markus Stocker的博客他很好地解释了如何在使用 observation 时表示传感器观察结果。 SSN 的模块本体论。我完全理解他的解释,但我发现有一件事多余地代表了一个的两个特
我有以下情况/代码; trait Model { def myField: String } case class MyModel(myField: String) extends Model
我想让一个案例类扩展一个特征 以下是我的要求: 我需要为 child 使用案例类。这是一个硬性要求,因为 scopt ( https://github.com/scopt/scopt ) parent
最近浏览了Markus Stocker的博客他很好地解释了如何在使用 observation 时表示传感器观察结果。 SSN 的模块本体论。我完全理解他的解释,但我发现有一件事多余地代表了一个的两个特
我有以下情况/代码; trait Model { def myField: String } case class MyModel(myField: String) extends Model
不确定标题是否完全有意义,对此感到抱歉。我是机器学习新手,正在使用 Scikit 和决策树。 这就是我想做的;我想获取所有输入并包含一个独特的功能,即客户端 ID。现在,客户端 ID 是唯一的,无法以
我想读取具有 Eigen 的 MNIST 数据集,每个文件都由一个矩阵表示。我希望在运行时确定矩阵大小,因为训练集和测试集的大小不同。 Map> MNIST_dataset((uchar*)*_dat
在 MATLAB 中,我可以选择一个分散的子矩阵,例如: A = [1 ,2 ,3;4,5,6;7,8,9] A([1,3],[1,3]) = [1,3;7,9] 有没有用 Eigen 做到这一点的聪
我在执行 Into 时遇到问题Rust 中通用结构的特征。下面是我正在尝试做的简化版本: struct Wrapper { value: T } impl Into for Wrapper {
我有这段 matlab 代码,我想用 Eigen 编写: [V_K,D_K] = eig(K); d_k = diag(D_K); ind_k = find(d_k > 1e-8); d_k(ind_
我正在使用 Eigen C++ 矩阵库,我想获取对矩阵列的引用。文档说要使用 matrix_object.col(index),但这似乎返回了一个表示列的对象,而不是简单地引用原始矩阵对象中的列。我担
在乘以很多旋转矩阵之后,由于舍入问题(去正交化),最终结果可能不再是有效的旋转矩阵 重新正交化的一种方法是遵循以下步骤: 将旋转矩阵转换为轴角表示法 ( link ) 将轴角转换回旋转矩阵 ( lin
定义可由命名空间中的多个类使用的常量的最佳方法是什么?我试图避免太多的继承,所以扩展基类不是一个理想的解决方案,我正在努力寻找一个使用特征的好的解决方案。这在 PHP 5.4 中是可行的还是应该采用不
定义可由命名空间中的多个类使用的常量的最佳方法是什么?我试图避免太多的继承,所以扩展基类不是一个理想的解决方案,我正在努力寻找一个使用特征的好的解决方案。这在 PHP 5.4 中是可行的还是应该采用不
我是一名优秀的程序员,十分优秀!