下面的代码显示了一个矩阵的奇点问题，因为我在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)
e = KtransInv * g * KInv
return e
file = open(‘transformed.txt’, ‘r’)
n = 4
transformed = zeros(n)
for counter, line in enumerate(file):
if counter == n:
break
transformed[counter] = float(line)
CSI = zeros((n, n))
for i in range(n):
for j in range(n):
CSI[i][j] = transformed[abs(i-j)]
A = getA(1)
F = getF(CSI, A)
G = getG(CSI, F, A)
K = zeros((n, n), float)
for j in range(n):
K[0][j] = 0.0001
for i in range(1, n):
for j in range(n):
K[i][j] = ((3.0*70.0*70.0*0.3)/(2.0*.0*.0))*((j*(i-j))/i)*(1.0+(70.0/.0)*j)
E = getE(G, K)
print G
print K
有人有什么建议要修吗？谢谢你
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/232686.html原文链接：https://javaforall.net
