4

我找到了一个乘以矩阵的代码。

% SWI-Prolog has transpose/2 in its clpfd library
:- use_module(library(clpfd)).

% N is the dot product of lists V1 and V2.
dot(V1, V2, N) :- maplist(product,V1,V2,P), sumlist(P,N).
product(N1,N2,N3) :- N3 is N1*N2.

% Matrix multiplication with matrices represented
% as lists of lists. M3 is the product of M1 and M2
mmult(M1, M2, M3) :- transpose(M2,MT), maplist(mm_helper(MT), M1, M3).
mm_helper(M2, I1, M3) :- maplist(dot(I1), M2, M3).

如果我输入:mult([[1,2],[3,4]],[[5,6],[7,8]],X).那么我得到 X = [[19, 22], [43, 50]].

但是我怎么能得到一个X = [[1*5+2*7, 1*6+2*8], [3*5+4*7, 3*6+4*8]] .

PS我是序言的新手。谢谢!

4

1 回答 1

7

这很容易:不要使用 is/2 来评估算术表达式,只需不计算它们并使用复合项而不是它们的数值。我为 product/3 做:而不是

product(N1,N2,N3) :- N3 is N1*N2.

我写:

product(N1, N2, N1*N2).

你只需要写一个对应版本的sumlist/2。

于 2012-03-26T11:46:05.903 回答