3回答
我有两个矩阵4x2。如何实现这样的乘法:输出应该是一个矩阵4x1，其中每个元素都是原始矩阵中行中元素的乘积之和。如下所示：
[1 2;
A = 3 4;
5 6;
7 8]
[1 2;
B = 3 4;
5 6;
7 8]
结果C矩阵将为：
[1*1 + 2*2;
C = 3*3 + 4*4;
5*5 + 6*6;
7*7 + 8*8]
浏览 2提问于2012-05-29得票数 3
回答已采纳
我正在试图弄清楚CFG将是什么样的语言，这样的描述：
只有一个a0或更多b0或更多c
我试过这个：S -> a | Sb | Sc
或者类似的事情：
S -> a | B | C
B -> Bb
C -> Cc
但似乎不起作用。还有其他更好的方法来用CFG来描述这种语言吗？
浏览 6提问于2021-03-25得票数 2
我需要在Scilab中创建一个块三对角矩阵。更明确地说，我想设置一个M^2 X M^2矩阵，格式如下：
B C 0 ... 0 0
D B C 0 ... 0
0 D B C ... 0
0 ... ... ... ...
0 0 ... D B C
0 0 0 ... D B
其中B、C和D是M×M矩阵。
我知道如何在"sysdiag“函数的帮助下创建块对角矩阵，但我还没有找到任何参考资料来轻松创建块三对角矩阵。
顺便说一下，给定矩阵B，C和D，我知道如何创建上面的矩阵。
auxA = sysdiag(B,B
浏览 1提问于2014-12-17得票数 0
1回答
我试图用Python创建K MxN matrices，它存储在一个(M,N,K) numpy数组C中，由两个矩阵( A和B )组成，形状分别为(K, M)和(K,N)。第一个矩阵被计算为C0 = a0.T x b0，其中a0是A的第一行，b1是第一行B，第二个矩阵是C1 = a1.T x b0等等。
现在我用一个for循环来计算矩阵。
import numpy as np
A = np.random.random((10,800))
B = np.random.random((10,500))
C = np.zeros((800,500,10))
for k in range(10):
浏览 1提问于2019-04-21得票数 2
回答已采纳
2回答
C语言中的稀疏矩阵乘法
、、、
我有两个矩阵市场格式的稀疏矩阵文件：
row col val
1 1 3.0
1 2 1.0
2 3 2.0
etc...
目前，我已经将文件分成了6个数组：
row_A[], col_A[], val_A[], row_B[] …
它们分别包含行索引、列索引和值。
我想轻松地将这两个矩阵相乘，而不必首先将它们转换为密集矩阵格式。有没有这样做的算法？
我在Quora上找到了这个伪代码，但我不确定它是不是最好的实现，或者它将如何在C：中实现
multiply(A,B):
for r in A.rows:
for c in A.rows[r]:
fo
浏览 1提问于2018-10-26得票数 2
1回答
我的工作-我得到了相同的答案，除了最后一个问题#什么是A^(-1)*A？
将numpy导入为np
初始化矩阵A
A = np.array([[1, 2, 0], [0, 5, 6], [7, 0, 9]])
转置A
A_trans = np.transpose(A)
取A的倒数
A_inv = np.linalg.inv(A)
什么是A^(-1)*A？
A_invA = (np.linalg.inv(A))*A
print(A)
print(A_trans)
print(A_inv)
print(A_invA)
[[1 2 0]
[0 5 6]
[7 0 9]]
[[1 0 7]
浏览 2提问于2020-12-24得票数 0
1回答
我有一个网站，教学生有限数学。
在这个网站上，我有一个计算器，可以进行组合和排列，现在我正在尝试包括一个矩阵计算器，它将对矩阵进行加、减、乘和逆矩阵运算。
我使用Javascript和西尔维斯特的库来做计算。我成功地创建了一个程序，它将接受用户输入到表单中的值并进行计算，但这只适用于特定大小的矩阵(4x4)。
我搞不懂的是，如何从表单中只获取非空的值，并从中创建一个矩阵，然后将这些值输出到结果表单中的相应字段中。
我正在使用的一些Sylvester方法
// create matrix from embedded array and assign to var A
var A = $M([
浏览 2提问于2011-10-20得票数 3
2回答
我有一个很大(500k X 500k)的稀疏矩阵。我想知道它的主要组成部分(实际上，即使只计算最大的PC也可以)。随机化的PCA效果很好，除了它本质上是寻找协方差矩阵的特征向量，而不是相关矩阵。使用大型稀疏矩阵的协方差矩阵找到PCA的软件包有什么想法吗？最好是在python中，尽管matlab和R也可以工作。
(作为参考，向提出了类似的问题，但方法涉及协方差矩阵)。
浏览 1提问于2012-11-28得票数 3
我即将用C#、MVC和AngularJS编写一个项目。在这个项目中，我有一个大约100-200个地点的列表，我需要分成组，每组3-4个目的地最多。组的划分方式是其中的目的地最少。然后必须建立从单个来源到每个组的路由(为每个组单独建立路由)。我对谷歌地图api有点陌生，迷失在其中，有没有一种服务可以获取多个目的地，并检索它们之间的最小距离？或者我必须在我的算法中通过循环遍历所有目的地来找到它？我所找到的只是起点和单一目的地之间的距离矩阵。
浏览 0提问于2016-12-05得票数 0
2回答
我想要计算向量，
S=A B u，
其中s和u是N维复向量，A是N×M复矩阵，B是M×N复矩阵。当A、B和u的元素表示为浮点数时，以下哪两种方法具有更好的准确性(更有效的数字)？
(1)先计算B。
首先做矩阵向量乘法，
Y=B u
然后，另一个矩阵向量乘法
S=A y
(2)先计算A、B。
首先做矩阵-矩阵乘法，
C=A B
然后，矩阵向量乘法
S=C U
有什么已知的一般规则吗？
顺便说一下，我知道方法(1)比方法(2)效率高得多。
浏览 2提问于2011-09-26得票数 3
回答已采纳
我有一个在PHP中表示为数组的矩阵：
array(
array('a','b'), // a | b
array('c','d') // c | d
)
我有一个用类似方式表示的第二个矩阵(这里没有内部数组，但我不介意添加它们)：
array(
'e', // e
'f' // f
)
我想要这样的组合：
array(
array('a','b',&
浏览 0提问于2010-11-14得票数 2
回答已采纳
在Maxima中，如何得到两个矩阵的元素积(Hadamard)？
e.g
a b
c d
*
e f
g h
->
a*e b*f
c*e d*h
浏览 1提问于2019-01-25得票数 1
回答已采纳
2回答
我真的很困惑。我找了很多教程，但找不到明确的答案。
A B B D
1 X x 5
2 x y 6
x 4
我想跨越这两个表。A，B，B，d是属性。
A B B D
1 X x 5
2 x x 5
1 X y 6
2 x y 6
1 X x 4
2 x x 4
按照笛卡尔的规则，这应该是正常的答案。横过所有的行。但我对同一栏B感到困惑。同一栏似乎是两次？
浏览 3提问于2015-11-26得票数 2
回答已采纳
2回答
我必须有两个大小均匀的非常大的向量(列)A和B。我想把向量A除以向量B。这将给我一个包含零的大矩阵AxB，除了最后一列。本专栏包含我感兴趣的值。当我简单地在Matlab脚本中对向量进行除法时，内存就用完了。可能是因为矩阵AxB变得非常大。也许我可以通过重复以下步骤来防止这种情况发生：
计算矩阵AxB的第一行
过滤最后一个值，并将其放入另一个向量C中。
删除使用的矩阵AxB行
对向量A中的所有行重做步骤1-4
我怎么才能做一个循环来做这个呢？
浏览 3提问于2013-02-12得票数 0
3回答
为什么BLAS有一个用于矩阵-矩阵乘法的gemm函数和一个用于矩阵-向量乘法的单独的gemv函数？矩阵-向量乘法不只是矩阵-矩阵乘法的特例，其中一个矩阵只有一行/列吗？
浏览 3提问于2011-08-16得票数 24
回答已采纳
我正在编写一个简单的程序来解决问题。
该程序创建了一个900 x 900矩阵。在程序创建此矩阵c之前，我按以下方式预先分配它：
c = zeros(900,900);
这将产生一个橙色错误消息：“分配给变量'c‘的值可能未使用”。
在以后的程序中，矩阵c中填充了数字。那么为什么会出现错误消息呢？
浏览 5提问于2010-02-23得票数 2
回答已采纳
1回答
矩阵计算
、、、、
我试图用numPy广播计算Rij = Aij /Cij。如果矩阵大小不相同(n×n)，也会引发异常。
我不太确定这是正确的，还是我应该做元素智慧或矩阵智慧。有人能告诉我怎么做吗？
A = [[(i+j)/2000 for i in range(500)] for j in range(500)]
B = [[(i-j)/2000 for i in range(500)] for j in range(500)]
C = [[((i+1)/(j+1))/2000 for i in range(500)] for j in range(500)]
def matrix_R(A,B,C):
浏览 2提问于2022-11-17得票数 0
1回答
我想对包含大量零的大型矩阵进行奇异值分解。特别是从对称矩阵A的对角化得到的U和S，这意味着A=U*S*转置(U^*)，其中S是对角矩阵，U包含所有特征向量作为列。
我在网络上搜索将奇异值分解和稀疏矩阵结合在一起的c++库，但是只能找到少数的库，而不是所有的特征向量。有没有人知道是否有这样的图书馆？
另外，在得到U和S之后，我需要把它们乘成一些稠密的向量。
浏览 5提问于2015-11-22得票数 0
问题如下：
我有一套20个方程式：
r1 = s1*h1_1 + s2*h1_2 + ... s20*h1_20
r2 = s1*h2_1 + ...
...
r20 = s1*h20_1 + ...
其中r，s和h是矩阵，*表示点态乘积。它可以用矩阵形式R= H*S重写，我想为S解这个方程，所以我需要计算inv(H)*R，但是如果H的每个元素都是矩阵，我如何计算inv(H)？我不能简单地将这些较小的矩阵串联成一个更大的矩阵H，然后再将其反转--因为这将给出一个与之不同的结果，例如用符号值反演矩阵H，然后用较小的矩阵代替这些符号值(因为方程组中存在点态积)。
到目前为止，我想出了一个解决
浏览 0提问于2015-07-23得票数 4
回答已采纳
1回答
我正确地将一个名为A的20x60矩阵除以一个名为B的1x60向量。因此，C = A/B，其中C是一个20x1向量。
MATLAB在A/B中做什么？我在上找到了答案
如果A是带m ~= n的矩形m×n矩阵，B是带有n列的矩阵，则x = B/A返回方程组x*A = B的最小二乘解。
但是，当我尝试C*B或B*C时，我会得到一个维度错误，为什么？
浏览 0提问于2018-09-26得票数 3
回答已采纳
1回答
将矩阵扩展为块矩阵-索引问题
、、、、
我希望用两种方式展开矩阵，但我无法计算出它的代码；我不需要精确的代码，而只是在文献中寻找一个指针。下面的问题涉及量子力学中换位方程的线性化，我从物理学的角度了解理论，但我不知道我的问题在编程世界中的名字。
我有矩阵A，它是NxN矩阵，在这个例子中，设为N=3；所以A是：
A=[a11 a12 a13
a21 a22 a23
a31 a32 a33]
我想做一个矩阵B，它是N^2xN^2大小。B有块状：
B=[A11 A12 A13
A21 A22 A23
A31 A32 A33]
A11=[a11 0 0
0 a11 0 ...
0
浏览 1提问于2016-07-28得票数 2
回答已采纳
1回答
float(multiply(colVec1,colVec2).T * (matrix*matrix[i,:].T)) + c
我是Python和numpy的新手，我正在努力理解上面的代码片段所做的事情。
multiply().T部分执行逐元素的乘法，然后执行转置，这样结果就变成了行向量。
我正在努力理解matrix[i,:]的工作。它是通过选择第一行向量来创建子矩阵，还是从第一行向量一直到矩阵的末尾创建子矩阵？
*执行一个点积，然后使用float()将其转换为浮点。
浏览 2提问于2016-09-02得票数 0
回答已采纳
1回答
用迭代法计算矩阵的行列式
、、、、
我编写了C++代码，使用递归方法计算矩阵的行列式。
现在，我想用迭代方法重写这个递归方法。但我想不出如何做到这一点。
问题是：如何将特定的递归方法(int行列式(int *& const int))重写为迭代方法？
以下是我的C++代码：
// Determinant C++
#include <iostream>
#include <ctime>
using namespace std;
void RandInitArray(int *, const int, int = 0, int = 20);
void Display(int *, const int
浏览 0提问于2018-06-09得票数 2
1回答
我需要为我的工作生成大量的随机均值不变正交矩阵。均值不变矩阵具有A*1_n=1_n性质，其中1_n是标量1基本np.ones(n)的大小为n的向量.我使用Python，特别是Numpy来创建我的矩阵，但是我想确保我的方法是正确的，也是最有效的。另外，我还想介绍我尝试过的三种不同的正交化方法的结果，并希望能解释为什么一种方法比另一种方法更快。在这篇文章的结尾，我问了四个关于我的发现的问题。
一般情况下，为了创建一个均不变的随机正交矩阵A，需要创建一个随机平方矩阵M1，将其第一列替换为1列，并将该矩阵正交化。然后，再用另一个矩阵M2这样做，最终的均值不变随机正交矩阵A= M1*(M2.T)。这一过
浏览 20提问于2021-02-26得票数 0
1回答
问题是，对于矩阵A、B、C和D --它们是长度n的向量n*n和x --在Matlab上以最有效的方式和最不有效的方式找到E = DCBAx。
计算E最明显的方法就是把它们直接乘以。
，这是最有效的方法吗？最不有效的方法是什么？
浏览 0提问于2014-09-25得票数 3
回答已采纳
1回答
虽然我使用的是C语言，但我认为编程语言和代码是不相关的。
我想用CSR格式表示一个对称矩阵，用于快速矩阵-向量积。
我读取的输入流只给出了矩阵的上三角部分，一个数字接一个数字，从左到右，行到行。
示例：
0 2 3 5
2 0 1 4
3 1 0 9
5 4 9 0
输入流：
0 2 3 5 0 1 4 0 9 0
比方说，如果我要从二维数组中构建一个常规矩阵，那就很容易了。我需要的只是我得到的数据。但CSR代表的情况并非如此。
浏览 9提问于2017-05-05得票数 0
1回答
gpu上的大矩阵乘法运算
、、、
我需要在GPU上使用CUDA实现大矩阵的矩阵乘法。每个矩阵的大小都比GPU内存大。因此，我认为我需要一个算法来有效地实现这一点。我在网上找了找，但什么也找不到。谁能给我这样的算法的名称或链接。
谢谢
浏览 0提问于2013-01-28得票数 8
回答已采纳
1回答
需要一些关于艾根的建议。
我想用这种方式从向量系数(一列矩阵) 'a‘和'b’创建矩阵：
| a0*b0 | a0*b1 | a0*b2 | ... | a0*bn |
| a1*b0 | a1*b1 | a1*b2 | ... | a1*bn |
| a2*b0 | a2*b1 | a2*b2 | ... | a2*bn |
| ... | ... | ... | ... | ... |
| am*b0 | am*b1 | am*b2 | ... | am*bn |
如果简单的话，我想创建由两个向量组成的乘法表，并将它添加到另一个矩阵中。
我怎样才能以最有效率
浏览 6提问于2020-02-06得票数 1
回答已采纳
我实现了一个算法，它以R行C列的矩阵作为输入。我说算法的最坏情况时间复杂度是O(C√C * R^3)或O(C^1.5 * R^3) 现在有人问我，它不能简单地表示为最坏的情况的O(R^3)吗？我想说的是，由于有两个输入(不是一个)，有时C可能很大，有时R可能很大，所以我们不能将其简化为简单的O(R^3)，应该同时考虑C和R。 我的答案正确吗？若否，原因为何？
浏览 40提问于2020-08-05得票数 1
回答已采纳
1回答
特征稀疏矩阵行列式为零
、、、
我试图计算我所操作的稀疏矩阵是否是正定的。为此，我试图使用西尔维斯特标准，这意味着领先的未成年人是积极的。
为了计算矩阵的行列式，我构造了一个矩阵的每个块的sparseLU求解器，它可以给出矩阵的行列式。但从某个维度(约130*130)开始，我得到的结果是，所有决定因素都是0。这不是我问题中的一些特殊维数(矩阵的块为32*32)，因此我认为这个问题与本征元所应用的某些截断算法有关，其行列式只是低于某个阈值。
我对这样一个机制的探索没有取得好的结果。我的矩阵的维数约为16k*16k，所有非零元素都在对角线附近的96个元素上。
是否有任何截断机制在本征中实现，我能以某种方式控制它的阈值吗？
浏览 1提问于2018-12-20得票数 4
回答已采纳
1回答
我试图在Julia中建立一个随机稀疏对称正定矩阵。
我尝试过的一种方法是：
n = 10
density = 0.1
X = sprand(n,n,density)
X = 0.5*(X+X') + n*eye(n)
@test isposdef(X) == true
这是可行的，但它给你一个对角占优矩阵，我不想要。所有条目都应该来自相同的随机分布。
浏览 8提问于2017-12-11得票数 0
回答已采纳
2回答
我创建了一个向量和矩阵类，我试图执行一些操作，例如矩阵和向量的乘法，矩阵和矩阵的乘法，矩阵和浮点数的乘法(标量)。我似乎很难得到矩阵*向量和矩阵*矩阵的正确乘积。下面是用于处理这些操作的矩阵类的一部分：
// Matrix * vector, result vector
Vector Matrix::operator*(const Vector & other) const
{
if (other.getDimensions() == 4)
{
float floats[4];
const float* temp = other.getData();
for (
浏览 6提问于2017-03-03得票数 0
回答已采纳
1回答
我想把sympy中的两个矩阵相乘。但是，.dot()函数似乎不能正常工作，因为当使用转置矩阵以便尺寸匹配时，出现的错误与不使用它时相同。
复制代码：
from sympy import *
x1, x2, x3 = symbols('x1 x2 x3')
x = Matrix([x1, x2, x3])
m = Matrix([1, 2, 3])
xm = (x-m)
print(xm.shape)
xmt = xm.T
print(xmt.shape)
c = Matrix([[1, 2, 3], [2, 1, 4], [3, 4, 5]])
print(c.shape)
浏览 1提问于2022-11-05得票数 2
回答已采纳
我有一个矩阵m，它有一些float数。当我把它与恒等矩阵相乘时，当得到的矩阵变成零的时候，一些元素。实际上，矩阵不应该有变化。有谁可以帮我？
在上面的图像中，在矩阵中没有0.013。
浏览 1提问于2021-01-03得票数 0
1回答
我正在为稀疏线性代数计算编写一个库，作为我论文工作的后端，我已经到了一个十字路口。我正在使用现代的Fortran (不要呻吟，它已经继承和多态以及所有的爵士乐已经10年了)。
从软件设计的角度来看，我的主要问题是使迭代求解器能够使用稀疏矩阵，而不知道它们的存储格式。迭代求解者需要知道的唯一功能是如何将矩阵乘以向量。
为此，我使用了一个抽象类sparse_matrix和用于矩阵向量乘法的虚拟方法matvec；然后有几个子类，表示每种存储格式，它们用自己的实现覆盖父类matvec。我相信这就是所谓的“模板”模式，是吗？
我正在考虑重构我的代码以使用组合而不是继承。为了达到这个目的，稀疏矩阵由一个有
浏览 0提问于2013-07-29得票数 3
2回答
我试图为矩阵中的每个单元格条目编写最小的值，无论是它的行和值还是同维新矩阵中的colsum值。
例如：
假设我有矩阵c，如下所示：
x <- matrix(seq(1:6),2)
x
[,1] [,2] [,3]
[1,] 1 3 5
[2,] 2 4 6
它的龙和钶钽铁矿石是：
rowSums(x)
[1] 9 12
colSums(x)
[1] 3 7 11
基于这一信息，新的矩阵应该如下所示：
[,1] [,2] [,3]
[1,] 3 7 9
[2,] 3 7 11
我一直在考虑
浏览 3提问于2013-12-28得票数 3
回答已采纳
在下面的代码中，为什么相同的符号*对矩阵和数组显示不同的行为？为什么它在矩阵的情况下不显示元素乘法？
a= np.arange(2*3).reshape(2,3)
print(a)
b= np.array([[0,0,0],[1,1,1]])
print(b)
c= a*b
print(c)
d=np.matrix([[1,2],[3,4]])
e= np.matrix([[0,0],[1,1]])
f= d*e
print("d: ",d,"e: ",e)
print(f)
浏览 0提问于2018-08-23得票数 3
3回答
用本征C++库使用GDB
、、、、
我正在使用艾根C++库，可从下载。这是一个更容易处理矩阵和数组的C++库。我使用g++编译器和gdb进行调试。但是，我发现在使用gdb时无法打印矩阵的内容(由特征提供)。
浏览 7提问于2014-08-01得票数 20
回答已采纳
我有一个df，有几个列的值为0或1。
a b c d e
1 0 0 0 0
0 1 0 1 0
0 1 0 1 0
1 0 1 0 1
我想要创建一个5乘5的矩阵，显示如果列中有1在同一行的总计数。我只想考虑1，在对角线情况下，它会自动地反映该列中的一行总数1。输出如下：
a b c d e
a 2 0 1 0 1
b 0 2 0 2 0
c 1 0 1 0 1
d 0 2 0 2 0
e 1 0 1 0 1
谢谢。
苏迪尔
浏览 1提问于2016-05-12得票数 1
回答已采纳
2回答
我正在做一个项目，它需要计算一个矩形矩阵的正交正交基，这个矩阵可能是秩不足的，也可能不是。在matlab中，我们可以调用基于svd的orth()函数来处理这个问题，但我需要用C实现它。我尝试使用Gram-Schmidt，当矩阵是满秩时，它工作得很好。C语言中有没有可以解决这个问题的库？或者关于用C实现它的一些提示？非常感谢。
浏览 0提问于2013-07-01得票数 0
回答已采纳
1回答
我有一个PySpark数据
simpleData = [("person0",1, 1, 0), \
("person1",1, 1, 1), \
("person2",1, 0, 0), \
("person3",0 ,0, 0 ), \
]
columns= ['persons_name','A', 'B', 'C']
exp = spark.createDataFrame(data = simpleData, schema = col
浏览 1提问于2020-12-02得票数 0
回答已采纳
1回答
我写了一个函数来计算两个动态数组的乘积。我的函数给出的是随机数，而不是预期的结果。我有两个初始化动态矩阵A和B的函数。现在，当我试图计算动态矩阵C，使得C=A*B时，我的函数没有显示正确的结果。我已经检查了我的代码，似乎我在某些地方做错了，但我不能准确地找出问题所在。下面，我粘贴了计算C=A*B的函数。我用Dev-C++用C++编写了它。您的帮助将不胜感激。
void CalculateDynamicMatC( double *DynMatA, double *DynMatB, int n)
{
int i, j, k, m=10;
double *DynMatC=new dou
浏览 4提问于2012-04-29得票数 0
1回答
矩阵操作的问题
、、
我有两个矩阵A和B和一个矩阵C，其中
C=A/B
现在我想取C的导数，它是
(B*dA/dx - A*dB/dx)/B^2
我不确定上面的内容会是什么样子。我发现在某个地方它看起来像
(B/B)*((dA/dx)/B) - (A/B)*(dB/dx)/B
我很困惑我能执行什么类型的操作。我可以这样做吗？
(B/B^2) * (dA/dx) - (A/B^2)*(dB/dx)
我尝试了上面两种方法，他们给了我不同的结果。我不确定什么运算是允许的，在矩阵的情况下，特别是乘法和除法。所以有谁能澄清一下吗？
浏览 3提问于2012-04-30得票数 0
回答已采纳
1回答
假设我有两个矩阵：M是2x3，C是3x4…
import numpy as np
M = np.matrix([ [1,2,3], [4,5,6]])
C = np.matrix([ [1,1,2,2], [1,2,1,2], [2,1,1,1]])
我需要将这些矩阵相乘，得到大小为2x4的结果。使用np.dot(M, C)很容易做到这一点，但我希望将向量积相乘，而不是求和。
例如，一个常规的点积应该是：
result[0,0] = (M[0,0]*C[0,0]) + (M[0,1]*C[1,0]) + (M[0,2]*C[2,0])
# result[0,0] = (1*1) + (2*1)
浏览 3提问于2016-01-08得票数 3
4回答
复矩阵指数在C++中的应用
、、、、
在c/ c++中真的可以计算复矩阵的矩阵指数吗？
我已经设法使用GNU Science Library中的blas函数求出了两个复矩阵的乘积。对于matC = matA * matB：
gsl_blas_zgemm (CblasNoTrans, CblasNoTrans, GSL_COMPLEX_ONE, matA, matB, GSL_COMPLEX_ZERO, matC);
我已经设法得到了一个矩阵的矩阵指数通过使用未记录的
gsl_linalg_exponential_ss(&m.matrix, &em.matrix, .01);
但这似乎不接受复杂的论点。
有没有办法做
浏览 0提问于2012-04-07得票数 4
2回答
矩阵乘法
、
例如，我有一个矩阵，A = np.array([[1,0],[0,1]])和B = np.array([[1,2],[3,4],[5,6]])形式的ndarray。我想要一个矩阵乘法的每个数组的数组B由矩阵A。
我做了一个循环，就像这样
C = []
for b in B:
C.append(np.matmul(A,b))
C=np.array(C)
out[]：
array([[1, 2],
[3, 4],
[5, 6]])
但我知道这对循环来说是很费时的。有更好的方法吗？
浏览 4提问于2021-07-09得票数 1
回答已采纳
1回答
假设我有4个不变的矩阵和1个可变的矩阵：
Eigen::MatrixXf M,N,O,P; // Fixed
Eigen::MatrixXf Q; // Value changes
我有一个想要计算的表达式：
auto myExpression = M*N*O*P*Q;
如果我执行以下操作：
for (int i = 0; i<5; ++i) {
Q = MatrixXf::Random(SomeSize,SomeSize);
MatrixXf result = myExpression.eval();
}
这会像预期的那样工作吗？这意味着我从Q的新值中得到了正确的结果？
浏览 10提问于2016-09-25得票数 0
1回答
我希望这不是一个重复的问题。我想用矩阵的形式计算向量之和的平方。例如，我有一个向量：
v = c(a, b, c, d)
我想要实现的是方程展开(a+b+c+d)^2 = a^2 + b^2 + c^2 + d^2 + 2ab + 2ac + 2ad + 2bc + 2bd + 2cd的矩阵形式：
a^2 ab ac ad
ab b^2 bc bd
ac bc c^2 cd
ad bd cd d^2
在r中是否有一个函数可以这样做，而不用手工写出它呢？
非常感谢你的帮助。
浏览 4提问于2021-10-04得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
