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
回答已采纳
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
在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
回答已采纳
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回答
我想要计算向量，
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
回答已采纳
2回答
我必须有两个大小均匀的非常大的向量(列)A和B。我想把向量A除以向量B。这将给我一个包含零的大矩阵AxB，除了最后一列。本专栏包含我感兴趣的值。当我简单地在Matlab脚本中对向量进行除法时，内存就用完了。可能是因为矩阵AxB变得非常大。也许我可以通过重复以下步骤来防止这种情况发生：
计算矩阵AxB的第一行
过滤最后一个值，并将其放入另一个向量C中。
删除使用的矩阵AxB行
对向量A中的所有行重做步骤1-4
我怎么才能做一个循环来做这个呢？
浏览 3提问于2013-02-12得票数 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
回答已采纳
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
回答已采纳
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
回答已采纳
我需要在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回答
我写了一个函数来计算两个动态数组的乘积。我的函数给出的是随机数，而不是预期的结果。我有两个初始化动态矩阵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 (mxn)和b (nxp)。对于每个n，我想用b的第n行乘以a的第n列，给出n (mxp)矩阵。然后，我想通过取矩阵中每个元素的平均值，将这些矩阵“折叠”成一个(mxp)矩阵。在numpy中，是否有一种合理有效的方法(就像在无循环中那样)来做到这一点？
浏览 2提问于2015-01-31得票数 0
回答已采纳
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
回答已采纳
2回答
我有两个矩阵，A和B，我想为它们的每个列做一个乘法，以生成一个新的矩阵。我想到的第一件事是
A = rand(4,3);
B = rand(4,3);
for J=1:SIZE(A,2)
for jj=1:size(B,2)
C(:,:,m) = A(:,j)*B(:,jj)' ;
m = m+1 ;
end
end
但我不想使用for循环，这会使它变慢。有办法吗？
我将使用C的三维矩阵，即由A和B的列相乘而成的矩阵，首先构建C，然后在每个循环中使用它的三维矩阵，还是只在每个循环中进行乘法？
浏览 0提问于2015-04-16得票数 3
回答已采纳
1回答
我在Matlab中有两个矩阵：A的GxM和B的size DxM。我想要创建一个大小为C的矩阵GxD，这样C(g,d)=sum(A(g,:).*B(d,:))/M就可以了。
让我用一个例子来更好地解释
clear
A=[1 2 3;
4 5 6;
7 8 9;
10 11 12]; %GxM
B=[14 15 16;
17 18 19]; %DxM
G=size(A,1);
D=size(B,1);
M=size(A,2);
矩阵C如下所示：
C=zeros(G,D);
for g=1:G
for d=1:D
C(g,d)=sum(
浏览 0提问于2018-09-01得票数 2
回答已采纳
我有一个在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
回答已采纳
我在Matlab中有两个稀疏矩阵，A和B，我想要计算一个三维矩阵C，这样
C(i,j,k) = A(i,j) * B(j,k)
我不需要循环就能这么做吗？
(附带问题:是否有此操作的名称？)
编辑:看来我的问题已经被问到了(只是为了完整的矩阵)：
浏览 4提问于2013-11-30得票数 1
回答已采纳
我有两个n乘m的矩阵，A和B。我想创建一个新的矩阵C，它类似于：
for i = 1:n
C = C + outerProduct(A(i,:), B(i,:));
end
即C是一个大小为m×m的矩阵，是A和B行的所有外积的和。
有没有一种不使用for循环的快速方法(考虑到for循环在Matlab中的速度非常慢)？
浏览 2提问于2012-01-10得票数 2
回答已采纳
1回答
矩阵积分乘法
、、
我有两个矩阵，我分别计算。
其中一个矩阵是Y，它的大小是2x25881，第二个是U，它由两个独立的向量组成，如下所示：
U = [ArrayAT ; ArrayAD]
U的尺寸为2x25881。
现在，我要计算(我的dt =0.001)的积分：
我的做法是：
Y_Int_Sum = 0;
for i = 1 : length(T)
Y_Int_Sum = Y_Int_Sum + Y_Int(:,i)'*[ArrayAT(i) ; ArrayAD(i)]*dt;
end
有什么更好的方法吗？
谢谢。
浏览 2提问于2017-09-09得票数 0
回答已采纳
我正在尝试在CUPY中使用多个GPU并行化多个矩阵乘法。 Cupy加速矩阵乘法(例如$A\times B$)。我想知道我是否有四个方阵A，B，C，D。我想在两个不同的本地GPU上计算AB和CD。如何在CUPY中完成此操作？ 例如，在tensorflow中， for i in xrange(FLAGS.num_gpus):
with tf.device('/gpu:%d' % i): 在CUPY中也有类似的方法吗？关于Cupy的事情是，它直接执行代码，所以它不能运行下一行(例如$C\times D$)，直到当前行结束(例如$A\times B$)。 谢谢Tos的帮助。现在新
浏览 379提问于2019-09-19得票数 1
1回答
我想将矩阵的下面两个矩阵相乘，并设置为我的模型的目标：
m = gp.Model("matrix")
x = m.addMVar((9, 9), vtype=GRB.BINARY, name="x")
c = np.random.rand(9,9)
m.setObjective(x @ c, GRB.MINIMIZE)
以下是我想要达到的目标
这给我带来了以下错误：
错误代码-1:变量不是一维MVar对象
我该怎么解决呢？我想Gurobi不接受2D Mvar对象乘法
浏览 8提问于2021-09-17得票数 0
回答已采纳
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回答
假设我有两个矩阵：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
问题如下：
我有一套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回答
我得到了一个用C++编写Strassen Vinograd算法的任务。我已经写了两次，但我的代码的第一个版本甚至不能工作。结果矩阵左下角的结果是正确的。我的第二个版本运行速度比朴素算法慢，即使N= 64+也是如此。所以我需要帮助，我做错了什么？重要提示:我不允许在递归和结构中使用动态矩阵。此外，最好在不复制的情况下，使用子矩阵的角元素的坐标进行乘法。不正确：
<code>A0</code>
慢，截止日期默认值1：
<code>A1</code>
浏览 18提问于2018-12-16得票数 3
1回答
两个Toeplitz矩阵的乘积？
、、、、
Toeplitz矩阵与长度正确的向量的乘积的O(n log n)算法是众所周知的:将其放入循环矩阵中，将其乘以向量(以及随后的零)，然后返回乘积的顶部n元素。
我发现很难找到最好的(时间)算法来乘以两个相同大小的Toeplitz矩阵。
谁能给我一个解决这个问题的算法？
浏览 1提问于2013-04-09得票数 6
回答已采纳
我有两个矩阵:矩阵a (6*6主矩阵)和矩阵b (14*2包含矩阵a中元素之间的关系)：
a = [
0 1 0 0 1 1
1 0 1 1 0 0
0 1 0 1 0 0
0 1 1 0 0 1
1 0 0 0 0 0
1 0 0 1 0 0 ]
b= [
1 2
1 5
1 6
2 1
2 3
2 4
3 2
3 4
4 2
4 3
4 6
5 1
6 1
6 4 ]
我需要计算以下内容：
TS(i,j)=(a(i,j)+1)/(the number of elements related to i(from matrix b)
* the
浏览 0提问于2012-06-26得票数 0
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回答
当我将numpy.dot用于ndarray时，出现了一个奇怪的现象。下面是安装在我的numpy上的MacOSX的信息：
{'define_macros': [('NO_ATLAS_INFO', 3), ('HAVE_CBLAS', None)],
'extra_compile_args': ['-msse3',
'-I/System/Library/Frameworks/vecLib.framework/Headers'],
'extra_link_args': ['-
浏览 3提问于2016-03-01得票数 1
4回答
我想从两个不同长度的数组a和b创建一个成对数组的矩阵：
a = np.array([1,2,3])
b = np.array([4,5,6,7])
因此，c矩阵应该如下所示：
[[1,4], [1,5], [1,6], ..., [3,7]]
浏览 0提问于2017-11-08得票数 2
1回答
我有两个密集的矩阵M，N，并且我想计算M*N，但在对元素进行运算时，我不使用规则乘法，我希望*表示取这两个矩阵中的最小值(但+仍然是加法)
例如，如果
M =
(1 2)
(3 4)
然后
M*M =
(1+2 1+2)
(1+3 2+4)
=
(3 3)
(4 6)
因为根据点积的执行方式，(1 2)*(1 3) = min(1,1) + min( 2，3) =1+2，依此类推。
这有可能吗？也许有一个模板类，其中定义了二元运算*，或者别的什么？
谢谢。
编辑:按照建议，我尝试定义我自己的标量类型，但它不起作用。
试试这个：
#include <iostream>
#i
浏览 0提问于2017-06-28得票数 0
1回答
我用过一点神经网络，但用得不多。因此，为了提高我的舒适度，我决定使用其中一个来解决我最喜欢的数学问题:快速矩阵乘法。标准算法将两个n×n矩阵相乘需要O(n^3)。Strassen算法只需要O(n^2.8)。基于Coppersmith和Winograd的工作的算法降到O(n^2.373)，但由于大的常数因子而不切实际。在后两者之间有很大的回旋余地。特别是，如果您可以使用48个或更少的乘法运算将两个4x4矩阵相乘，那么您已经比Strassen做得更好了。
这里是我的设置:我有两个(伪随机生成的) nxn矩阵，A和B。一个神经网络采用A的元素的NMULT线性组合和B的NMULT线性组合，逐点将它们相
浏览 1提问于2017-10-28得票数 1
我有两个矩阵，分别是大小为mxn的T1和T2。我想找出两个矩阵之间的相关系数。
到目前为止，我还没有为它使用任何内置库函数。我正在为它执行以下步骤：
首先，我计算两个矩阵的平均值为：
M1 = T1.mean()
M2 = T2.mean()
然后从相应的矩阵中减去平均值，如下所示：
A = np.subtract(T1, M1)
B = np.subtract(T2, M2)
其中，np是numpy库，A和B是减法后的结果矩阵。
现在，我计算相关系数为：
alpha = np.sum(A*B) / (np.sqrt((np.sum(A))*np.sum(B)))
然而，我得到的值远远大于1，
浏览 3提问于2016-12-01得票数 3
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、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
当取两个向量的点积时，我观察到一个wierd输出。我的代码是
a1 = np.array([[1], [2], [3]])
a2 = np.array([[1, 2, 3]])
print(a1*a2)
print(np.dot(a1, a2))
两者的输出都是相同的，我不明白为什么当被问到点时，它会乘以两个矩阵
产品。对于形状为(x，1)和(1，y)的任意矩阵也是如此。
谢谢
浏览 4提问于2020-10-02得票数 1
回答已采纳
2回答
在我的工作流程中，我遇到了一些阻碍，因为我试图做一个双循环来填充矩阵工作。在我的脚本开始时，我定义：
Tijd_SAT_i <- c("08", "09", "10", "11", "12", "13", "14")
D <- c(0, 0.01, 0.02, 0.03, 0.04, 0.05)
POD <- matrix(1:(length(Tijd_SAT_i)*length(D)), c(length(Tijd_SAT_i), length(D)))
FAR
浏览 1提问于2013-10-02得票数 0
回答已采纳
我自己会说不止一种语言。但我想问一下，如果你有一段视频，你可能会对说印地语的人感兴趣，那么用印地语和英语发布标签会帮助他们更容易地找到视频吗？这就是我的全部要求。
我的视频是印地语。印度人用英语和印地语(拼音)搜索。但他们想要印地语的录像带。
所以写关键字英语和印地语(拼音)好吗？另外，写标题英文和印地语(拼音)也很好？
或仅用英语/印地语书写
浏览 3提问于2017-09-18得票数 0
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回答
我正在为稀疏线性代数计算编写一个库，作为我论文工作的后端，我已经到了一个十字路口。我正在使用现代的Fortran (不要呻吟，它已经继承和多态以及所有的爵士乐已经10年了)。
从软件设计的角度来看，我的主要问题是使迭代求解器能够使用稀疏矩阵，而不知道它们的存储格式。迭代求解者需要知道的唯一功能是如何将矩阵乘以向量。
为此，我使用了一个抽象类sparse_matrix和用于矩阵向量乘法的虚拟方法matvec；然后有几个子类，表示每种存储格式，它们用自己的实现覆盖父类matvec。我相信这就是所谓的“模板”模式，是吗？
我正在考虑重构我的代码以使用组合而不是继承。为了达到这个目的，稀疏矩阵由一个有
浏览 0提问于2013-07-29得票数 3
我有一个矩阵m，它有一些float数。当我把它与恒等矩阵相乘时，当得到的矩阵变成零的时候，一些元素。实际上，矩阵不应该有变化。有谁可以帮我？
在上面的图像中，在矩阵中没有0.013。
浏览 1提问于2021-01-03得票数 0
1回答
因此，这是一个数学事实，如果一个矩阵的行列式等于零，那么这个矩阵一定是奇异的(不可逆)。现在，我遇到的问题是，当我计算矩阵的行列式时，它等于0，然而，当我计算逆时，它是存在的。我认为这与R计算行列式的方式有关，两者并不一致。下面是我正在尝试的代码(我不会打印solve的结果，因为矩阵是100x100)。
> Rinv = solve(R)
>
> det(R)
[1] 0
>
> #Using a Cholesky Factorization
> L = chol(R)
> Q = t(L)
>
> det(L)*det(Q)
[1]
浏览 2提问于2013-08-08得票数 4
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
2回答
我有一个网站，有几页，它将有两种主要语言:西班牙语和英语。然而，在西班牙语中，我正在创建其他语言：
西班牙-阿根廷
西班牙-哥伦比亚
西班牙-秘鲁
西班牙-智利
..。
问题是:我不想重复西班牙语国家的所有内容，而是更改2或3页:主要是因为我想在某些页面中替换国家名称。
在不重复每个西班牙语国家的所有内容的情况下，创建特定国家版本的最佳方法是什么？谢谢!
注意:创建一个“语言中立”页面是个好主意吗？然后在必要时创建英文版本和修改过的西班牙文版本？还不确定“语言中立”是干什么用的。
浏览 0提问于2014-04-17得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
