lcmpoly.m

Software for design and tuninig of SISO and MIMO contol systems

function x=lcmpoly(A,y)
%function x=lcmpoly(A,y)
%Find the least common multiple of polynomials in a cell array A
if isempty(A)
   x=1;
   return
end
x=1;
[n]=length(A);
for i=1:n
   x=lcm(x,mt2poly(A{i},y));
end
function x=lcm(a,b)
x=1;
temp=gcdpoly(a,b);
while length(temp) > 1
   x=conv(x,temp);
   a=deconv(a,temp);
   b=deconv(b,temp);
   temp=gcdpoly(a,b);
end
x=conv(x,a);
x=conv(x,b);