Question

我有3个矩阵：

T_01 = ['cosd*t1', '-sind*t1', '0', 'd1*cosd*t1'; 'sind*t1', 'cosd*t1', '0', 'd1*sind*t1'; '0', '1', '1', '0'; '0', '0', '0', '1']

T_12 = ['cosd*t2', '-sind*t2', '0', 'd2*cosd*t2'; 'sind*t2', 'cosd*t2', '0', 'd2*sind*t2'; '0', '1', '1', '0'; '0', '0', '0', '1']

T_23 = ['cosd*t3', '-sind*t3', '0', 'd3*cosd*t3'; 'sind*t3', 'cosd*t3', '0', 'd3*sind*t3'; '0', '1', '1', '0'; '0', '0', '0', '1']

我需要进行符号乘法，所以我正在尝试：

mulf(T_01,T_12,T_23)

但是我收到了这个错误：

                     !--error 39 
mulf: Quantidade incorreta de argumentos de entrada: esperava-se 2.

发生了什么事？

Obs：抱歉我的英文。

Answer 1

如果你想要的是获得两个矩阵的符号乘法，你必须实现这样的功能。在这里，我实现了三个功能，它们可以一起执行您想要的任务：

function s = scaProd(a,b)
    //escalar product of two vectors
    //using recursion

    if (a == [] | b == []) then
       s = ""

    elseif (max(size(a)) ~= max(size(b))) | ...
           (min(size(a)) ~= min(size(b))) | ...
           (min(size(a)) ~= 1) then
        error("vectorMulf: Wrong dimensions")

    else
        s = addf( mulf(a(1), b(1)) , scaProd(a(2:$), b(2:$)) )

    end
endfunction

function s = matrixMulf(a,b)
    //matrix multiplication

    acols = size(a,'c');
    brows = size(b,'r');
    if acols ~= brows then
        error("matrixMulf: Wrong dimensions")
    end

    arows = size(a,'r');
    bcols = size(b,'c');
    s = string(zeros(arows,bcols));

    for i = 1 : arows
        for j = 1 : bcols
            s(i,j) = scaProd(a(i,:),b(:,j)');
        end
    end
endfunction

function s = addP(a)
    //encolses each element of a in a pair of parenthesis
    s = string(zeros(a));

    for i = 1 : size(a,'r')
        for j = 1 : size(a,'c')
            s(i,j) = "(" + a(i,j) + ")"
        end
    end
endfunction

以下是它的输出示例。测试代码：

A = [1 2; 3 4];
B = [5 6; 7 8];
C = [9 0; 1 2];
disp(A*B*C)
As = string(A);
Bs = string(B);
Cs = string(C);
disp(matrixMulf(As,addP(matrixMulf(Bs, Cs))))

控制台输出：

   193.   44. 
   437.   100.

!1*(5*9+6*1)+2*(7*9+8*1)  1*(5*0+6*2)+2*(7*0+8*2)  !
!                                                  !
!3*(5*9+6*1)+4*(7*9+8*1)  3*(5*0+6*2)+4*(7*0+8*2)  !

对于你想要的结果，你应该这样做：

使用addP()
执行符号乘法，例如matrixMulf(t1,addP(matrixMulf(t2,t3)))，其中t1，t2，t3是矩阵的随附版本。

最后两个笔记：

在每个乘法步骤中使用addP以获得正确的结果非常重要。您可以通过删除我给出的示例中的(和)来检查这一点：结果不会是正确的。
Scilab 6.0.0中没有mulf和addf功能。所以请记住，如果将Scilab升级到当前的稳定版本，您将无法使用它们。

Scilab符号矩阵乘法

1 个答案: