1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
function s1=%r_m_r(s1,s2)
//s1=s1*s2
//!
[s1,s2]=sysconv(s1,s2),
[n1,d1]=s1(["num","den"]);
[n2,d2]=s2(["num","den"]);
sz1=size(n1)
sz2=size(n2)
indef=%f
if or(sz1==-1) then
n1=n1+0;d1=d1+0;sz1=[1 1]
if prod(sz2)==1 then indef=%t,else error(14),end
end
if or(sz2==-1) then
n2=n2+0;d2=d2+0;sz2=[1 1]
if prod(sz1)==1 then indef=%t,else error(14),end
end
if prod(sz1)==1 then
[num,den]=simp(n1.*n2(:),d1.*d2(:)),
num=matrix(num,sz2)
den=matrix(den,sz2)
elseif prod(sz2)==1 then
[num,den]=simp(n1(:).*n2,d1(:).*d2),
num=matrix(num,sz1)
den=matrix(den,sz1)
else
if size(sz1,"*")>2|size(sz2,"*")>2 then error(10),end
if sz1(2)<>sz2(1) then error(10),end,
l1=sz1(1);m1=sz1(2);m2=sz2(2);
for i=1:l1, pp(i)=lcm(d1(i,:)),end,
for j=1:m2,
y=lcm(d2(:,j)),
for i=1:l1,
yij=y*pp(i),
x=0;
for k=1:m1,
x=x+n1(i,k)*n2(k,j)*pppdiv(y,d2(k,j))*pppdiv(pp(i),d1(i,k)),
end
num(i,j)=x,den(i,j)=yij,
end
end
[num,den]=simp(num,den),
end
if indef then
num=num*eye()
den=den*eye()
end
s1=rlist(num,den,s1.dt)
endfunction
|
