eige1.m

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

num=[0 1 3];
den=[1 3 2];
[A B C D]=tf2ss(num, den);
syms x
[m,n]=size(A);
eq=det(x*eye(m,n)-A);
s=solve(eq,x);
s=s';
P=ones(m,n);
for i=2:m
    for j=1:n
           P(i,j)=(s(j))^(i-1);
    end
end

A_new=inv(P)*A*P;
B_new=inv(P)*B;
C_new=C*P;
D_new=D;

[num1,den1]=ss2tf(A,B,C,D);

fprintf('transfer function for new state space is equal to\n')
printsys(num1,den1)
fprintf('\n original transfer function was\n')
printsys(num,den)
step(A,B,C,D)