先将矩阵打包成一个类,重载加法和减法运算,就可以很简单地实现了,目前这个程序中还有一些冗余的地方,比如将行数和列数分别讨论,其实可以归为一个矩阵size。类的定义感觉需要改进,不过懒得做下去了
# -*- coding: utf-8 -*-
def matrixproduct(a, b):
def matrixproductmask(mat_a, mat_b):
if mat_a.row == 1:
c11 = [[mat_a.mat_list[mat_a.row_list[0]][mat_a.col_list[0]] *
mat_b.mat_list[mat_b.row_list[0]][mat_b.col_list[0]]]]
return Martrix(c11)
else:
mat_a11 = mat_a.divide('11')
mat_a12 = mat_a.divide('12')
mat_a21 = mat_a.divide('21')
mat_a22 = mat_a.divide('22')
mat_b11 = mat_b.divide('11')
mat_b12 = mat_b.divide('12')
mat_b21 = mat_b.divide('21')
mat_b22 = mat_b.divide('22')
s1 = mat_b12 - mat_b22
s2 =