import numpy as np
with open('matrix.txt', 'r') as f:
x = []
for line in f:
x.append(map(int, line.split()))
f.close()
a = array(x)
l, v = eig(a)
exponent = array(exp(l))
L = identity(len(l))
for i in xrange(len(l)):
L[i][i] = exponent[0][i]
print L
我的代码打开一个包含矩阵的文本文件:
1 2
3 4
并将其x
作为整数放入列表中。然后将列表
x
转换为数组a
。的特征值
a
放在 中l
,特征向量放在 中v
。然后我想把 exp(a) 放在另一个数组
exponent
中。L
然后我创建一个任意长度的单位矩阵l
。我的 for 循环应该在单位矩阵的对角线上取值
exponent
并替换 1,但我收到错误消息invalid index to scalar variable
.
我的代码有什么问题?