matrixmul.m

来自「各种矩阵相称的算法时间比较」· M 代码 · 共 26 行

function C=matrixmul(A,B)
%此程序用于实现任意两个规模为m*n和n*m的矩阵相乘
%格式为C=matrixmul(A,B)
%A和B分别为两个任意阶矩阵，但要确保阶数符合矩阵相乘的规则
[m,n]=size(A);
[c,d]=size(B);
if (c==n)==1
    for i=1:m
        for j=1:m
            temp=0;
            for k=1:n
                temp=temp+A(i,k)*B(k,j);
            end
            C(i,j)=temp;
        end
    end
    if m>=n
        C; 
    else
        C=C(1:m,1:d);
    end
else
    error('Inner matrix dimensions must agree.')
end