测试矩阵是对称还是斜对称的？

我正在做我的Python实验室练习。我做了代码，但我不知道如果矩阵是对称或扭曲对称的，我将如何进行测试。问题是： 1.如果矩阵是对称的或者是偏斜对称的，则进行测试。 2.我将如何使用断言来确保矩阵是方阵？ 3.有人可以重新检查我的代码吗？谢谢！测试矩阵是对称还是斜对称的？

def symmetric(X): 
H = X.copy() 
D = len(H) 
F = len(H[0]) 
listx = [] 
listy = [] 
for d in range(0,D): 
    listx.append(0) 
for f in range(0,F): 
    listy.append(listx) 
K = np.array(listy) 
for i in range(0,D): 
    for j in range(0,F): 
     K[j][i] = H[i][j] 
return K 

def skewsymmetric(Y): 
H = Y.copy() 
D = len(H) 
F = len(H[0]) 
listx = [] 
listy = [] 
R = -1 
for d in range(0,D): 
    listx.append(0) 
for f in range(0,F): 
    listy.append(listx) 
K = np.array(listy) 
for i in range(0,D): 
    for j in range(0,F): 
     K[j][i] = H[i][j]*R 
return K 

if __name__ == "__main__": 
X = np.array([[20,120,200],[120,10,150],[200,150,30]]) 
Y = np.array([[0,1,-3],[-1,0,-2],[3,2,0]]) 
print "X =", X 
print "Y =", Y 
print "This is Symmetric matrix:", symmetric(X) 
print "This is Skew-symmetric matrix:", skewsymmetric(Y)

来源

2014-09-01 Brylle Cuetara Garcia

有一个简单的方法来检查，如果矩阵是对称的使用numpy的：

(arr.transpose() == arr).all()

如果你需要检查它是否是斜对称，然后在比较更改第二个“改编”的符号：

(arr.transpose() == -arr).all()

来源

2014-09-01 09:37:28 masteusz

我不能使用转置我的代码。因为我的实验练习需要使用“for循环”和“断言”。所以我需要坚持我的代码，只需要改变一些东西。也可以将这两种方法合而为一。 – 2014-09-01 09:44:08

我收到“AttributeError：all not found”。我是否缺少导入或其他内容？ – brandones 2016-09-03 19:26:39

我无法重现您的错误。也许你的进口有错误？你可以声明数组没有问题吗？ – masteusz 2016-09-10 17:20:46

测试矩阵是对称还是斜对称的？

回答

相关问题