pq2ss.m

来自「A Matlab toolbox for exact linear time-i」· M 代码 · 共 28 行

% PQ2SS - from LMF to I/S/O representation.
%
% sys = pq2ss(p,q)
%
% P - PxMxL array of coefficients for the polynomial P(z)
% Q - PxPxL array of coefficients for the polynomial Q(z)
%     Note: for SISO system, P and Q can be vectors 
% SYS - ss object of a minimal I/S/O representation
%
% Note: requires the Polynomial Toolbox Version 1.5

function sys = pq2ss(P,Q)

if length(size(P)) == 2 % SISO
  Q = Q(end:-1:1);
  P = P(end:-1:1);
  sys = ss(tf(Q,P,-1));      
else % MIMO
  l1 = size(P,3); l = l1 - 1;
  m  = size(Q,1);
  p  = size(P,1);
  Q  = ppck(reshape(Q,p,l1*m),l);
  P  = ppck(reshape(P,p,l1*p),l);
  % call the Polynomial Toolbox
  [a,b,c,d] = lmf2ss(Q,P);
  sys = ss(a,b,c,d,-1); 
end