eige.m

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

numerator=[0 1 3];
denominator=[1 3 2];
[A B C D]=tf2ss(numerator, denominator);
[m,n]=size(A);
eig=roots(poly(A));
P=ones(m,n);
for i=2:m
    for j=1:n
           P(i,j)=(eig(j))^(i-1);
    end
end

A1=inv(P)*A*P;
B1=inv(P)*B;
C1=C*P;
D1=D;

[numerator1,denominator1]=ss2tf(A,B,C,D);
fprintf('transfer function for new state space is equal to\n')
printsys(numerator1,denominator1)
fprintf('\n original transfer function was\n')
printsys(numerator,denominator)
%step response is:
step(A,B,C,D)