赞
踩
下面的代码显示了矩阵的奇点问题,因为我在Pycharm中工作了
raise LinAlgError("Singular matrix")
numpy.linalg.linalg.LinAlgError: Singular matrix
我想问题是K但我无法理解究竟是怎么回事:
from numpy import zeros
from numpy.linalg import linalg
import math
def getA(kappa):
matrix = zeros((n, n), float)
for i in range(n):
for j in range(n):
matrix[i][j] = 2*math.cos((2*math.pi/n)*(abs(j-i))*kappa)
return matrix
def getF(csi, a):
csiInv = linalg.inv(csi)
valueF = csiInv * a * csiInv * a
traceF = valueF.trace()
return 0.5 * traceF
def getG(csi, f, a):
csiInv = linalg.inv(csi)
valueG = (csiInv * a * csiInv) / (2 * f)
return valueG
def getE(g, k):
KInv = linalg.inv(k)
Ktrans = linalg.transpose(k)
KtransInv = linalg.inv(Ktrans)
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。