练习2.37
这道题花了我太长的时间了,一开始我就把题目中的m看成了w。然后题中给出的dot-product的两个参数我还以为一个是向量另一个是矩阵。怎么算都算不出来,直到看到“返回求和...”里的w只有一个i而没有j。好了,那么开始按照题目的要求来做题了。
既然发现了自己的错误,那么就知道了dot-product是干嘛的了,它可以用来求一个矩阵中的一列和一个向量的积。因此定义出matrix-*-vector就不难了。
(define (matrix-*-vector m v)
(map (lambda (col)
(dot-product col v))
m))
下面就来用题目中的数据来测试一下。
(matrix-*-vector ‘( (1 2 3 4) (4 5 6 6)(6 7 8 9) ) )
;Value: (30 56 80)
接下来我们来看看transpose,之所以把matrix-*-matrix放在最后是因为其的定义中需要transpose。
(define (transpose mat)
(accumulate-n cons ‘() mat))
看似这么一行代码,但其功能可强大了。
(transpose ‘((1 2 3 4) (4 5 6 6) (6 7 89)))
;Value: ((1 4 6) (2 5 7) (3 6 8) (4 69))
如果这里accumulate的变化有什么不明白的,可以回过头看看练习2.33的解答中的截图。而accumulate-n只不过是个accumulate的外壳而已,在变换中accumulate-n都将会慢慢变成accumulate。
在线性代数里我们学过,两个矩阵的乘积mn的第一列第一行的值等于m的第一列和n的第一行的点积,mn的第一列第二行的值等于m的第一列和n的第二行的点积……
(define (matrix-*-matrix m n)
(let ((cols (transpose n)))
(map (lambda (col-of-m)
(matrix-*-vector colscol-of-m))
m)))