jordan.m

System will automatically delete the directory of debug and release, so please do not put files on t

function [A,B,C,D]=jordan(num,den)
if length(num)>length(den)
    deg=length(num)-1;
else
    deg=length(den)-1;
end

[r,p,k]=residue(num,den);
r=flipud(r);
p=flipud(p);

A=zeros(deg,deg);
c=0;
for i=1:deg
    A(i,i)=p(i);
    if i>1&&p(i)==p(i-1)
        A(i-1,i)=1;
        c=c+1;
    end
end

B=ones(deg,1);
B(1:c)=0;

C=r';

D=k;