问题标签 [symmetric]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
cuda - 用 CUDA 计算 (A - B(D^-1)B^T )^-1
使用CUDA C计算以下表达式的最有效方法可能是什么?
(A - B(D^-1)B^T)^-1
其中D是一个非常大的对称矩阵,A是一个小的对称矩阵,这使得B和B^T是中等大小的矩形非对称矩阵。当然 (^-1) 和 (^T) 分别是逆运算和转置运算。
python - python中的复杂对称矩阵
我试图在 python中对角化一个复杂的对称矩阵。
我查看了 numpy 和 scipy linalg 例程,但它们似乎都处理厄米特矩阵或实对称矩阵。
我正在寻找的是某种方法来获得我的起始复杂和对称矩阵的高木分解。这基本上是标准的特征分解 S = VDV^-1 但是,由于起始矩阵 S 是对称的,因此生成的 V 矩阵应该自动正交,即 VT = V^-1。
有什么帮助吗?
谢谢
matlab - matlab中的对称FFT
我将在matlab中实现OFDM系统。我需要从数据中获取 IFFT 对称,然后再从结果中获取 FFT。IFFT 正确,但 FFT 不正确,结果数字的前一半类似于 IFFT 之前的数据,但后半部分是错误的。当我知道 IFFT 采用“对称”时,我只是不知道是否应该使用 FFT 函数。
这是我使用的功能:
谢谢
linear-algebra - 对于一个只有 3x3 的对称正定线性系统,Cholesky 是否仍然比 Householder 快?
我正在尝试解决对称正定的线性系统Ax=b
。A
3x3
尽管它的规模很小,但我将不得不重复它A
数百万次。所以效率还是很重要的。
线性系统有很多求解器(C++,通过 Eigen)。我个人更喜欢:HouseholderQr().solve()
, 和llt().solve()
, ldlt().solve()
。
我知道当n
非常大时,基于 Cholesky 分解的求解器比 Householder 的求解器更快。但是对于我的情况何时n
只有 3,我如何比较它们的相对效率?有没有精确float operation
分析的公式?
谢谢
r - R在对称矩阵上应用函数
我想利用一apply
组函数中的一个来进行一些计算。
首先,我有两个矩阵,mat1
并且mat2
:
mat1
:
mat2
:
mat2
是使用与此处无关的函数计算mat1
的,本质上我想应用一个加权函数,当数据较少(因此不太准确)时mat1
对结果进行惩罚。mat2
因此,为了实现这一点,我想对于 中的某个坐标x,y
,mat2
计算 的两列的成对完整性mat1
。
例如:mat2["a","b"]
或mat2["b","a"]
(应该相同)将成为原始值 *(和的完整行和mat1
/a
的总行b
)。mat1
a
b
所以真正的问题是如何将一个函数应用于一个矩阵,该矩阵为每一列循环每一列(双循环)并将其存储在权重矩阵中以与另一个矩阵相乘?
我已经可以使用rollapply
fromzoo
包比较两行,如下所示:
我得到:
如,比较 a 和 b,2 行是完整的,比较 b 和 c,1 行是完整的。那么如何比较 a 与 b、a 与 c 以及 b 与 c 的比较呢?
谢谢。
r - R在对称矩阵中获得最高的x单元格及其行名/列名?
我有一个对称矩阵mat
:
我想计算它的两个最高元素。现在因为它是一个对称矩阵,我想像upper.tri
这样使用:
所以上面是这样一种迂回的做事方式,即使这样,输出也有“重复”的行,因为它们只是被转置了..
有人有更直观的方法吗?
谢谢。
r - 在 R 上整数值有限的对称矩阵中找到最大行列式
我正在使用 R 来尝试找到对称 3x3 矩阵的最大可能行列式,其所有条目都位于整数集合 {1,2,3,4,5,6} 可以具有?
有什么想法我会怎么做?
r - R获取列的所有组合并对它们执行返回对称矩阵的函数?
我有一个mat
看起来像这样的矩阵:
我想通过应用函数来构建一个对称矩阵Sum(Col1-Col2)
。最终结果将如下所示:
这样空格代表差异的总和。IE[1,2] = Sum(A-B)
。
我研究了以下方法:
但它并没有给我一个很好的对称矩阵。
有人有想法么?
谢谢。
编辑 - 感谢特洛伊上述作品。但是如果我想计算呢Sum((Col1-Col2)^2)
(Sum(((A_1,A_2,..,A_n)-(B_1,B_2,..,B_n))^2)
所以不能先求和A
然后B
再减去,否则答案将被关闭)。
r - R加速方阵的矢量化
任何人都可以帮助我加快一些代码:
哪里mat
是观察和对象的NxM
矩阵,例如:N
M
LinReg
定义为:
基本上,我正在对 中的每个对象组合(即Obj1, Obj2
和Obj2,Obj3
和Obj1, Obj3
)执行回归分析(OLS) mat
,然后使用包中的adf.test
函数tseries
并存储p-value
. 最终结果sym.pr
是所有的对称矩阵p-values
(但实际上它不是 100% 对称的,请参阅此处了解更多信息),但它就足够了。
使用上面的代码,在一个600x300
矩阵(600 个观察和 300 个对象)上,大约需要 15 分钟。
我想可能只计算对称矩阵的上三角形,但不知道如何去做。
有任何想法吗?
谢谢。
r - R - 列表与重复项的组合?
如何获得具有重复项的列表的所有组合。重复我的意思是一个元素本身。我正在构建一个对称矩阵。
结果是:
当我想要:
甚至:
但是我的矩阵很大,所以我想将组合保持在最低限度(所以最好是第二个结果),因为更多的组合 = 更多的计算 = 更长的运行时间..
谢谢。