Python中的矩阵乘法

时间:2013-02-21 13:50:17

标签: python matrix multiplication

你好(请原谅我的英文),我对使用矩阵乘法的python有一个很大的疑问,我创建了一个列表列表并乘以一个缩放矩阵,这就是我所做的并且我不能执行乘法运算索引的问题,我用纸和笔检查它是否有效,我做了一些不好的事情来容纳索引,或者我从一开始就错误地容纳矩阵了?

def main():
if len(sys.argv) > 1:
    v = int(sys.argv[1])
else:
    print "error python exe:"
    print "\tpython <programa.py> <num_vertices>"
A = []
for i in range(v):
    A.append([0]*2)
for i in range(v):
    for j in range(2):
        A[i][j] = input("v: ")  


print A 
    Escala(A)

def Escala(A):
    print "Escala"
    sx = input("Sx: ")
    sy = input("Sy: ")
    S = [(sx,0),(0,sy)]
    print S
    M = mult(S,A)
    print M



def mult(m1,m2):
M = zero(len(m1),len(m2[0])) 
    for i in range(len(m2)):   
        for j in range(len(m2[0])): 
            for k in range(len(m1)): 
                M[i][j] += m1[k][j]*m2[k][j]
    print M                 
    return M

def zero(m,n):
# Create zero matrix
new_matrix = [[0 for row in range(n)] for col in range(m)]
return new_matrix

1 个答案:

答案 0 :(得分:5)

这对我来说似乎不对:

M[i][j] += m1[k][j]*m2[k][j]

不应该是:

M[i][j] += m1[i][k]*m2[k][j]