statespace.m

clc,clear all
disp('CALCULATION OF STATE TRANSITION MATRIX, STATE VECTOR AND OUTPUT')
a=input('Matrix [A]=');
b=input('Matrix [B]=');
c=input('Matrix [C]=');
d=input('Matrix [D]=');
xo=input('Initial States Vector=');
disp('Eigenvalues')
autov=eig(a)
disp('State Transition Matrix')
syms t tao s
mte=expm(a*t);
pretty(mte)
ft=input('Input function or Input Vector f(t)=');
%fts=laplace(ft);
mtetao=subs(mte,t,t-tao);
fttao=subs(ft,t,tao);
xt=mte*xo+int(mtetao*b*fttao,tao,0,t);
disp('State Matrix')
xt=simplify(xt)
disp('Output')
yt=c*xt+d*ft;
yt=simplify(yt)