compensator.m

several numerical method code

numg=1000;
deng=poly([0 -5 -20]);
G=tf(numg,deng);
numc=[1 3.51272];
denc=[1 37.7144];
Gc=tf(numc,denc);
Ge=G*10.7366*Gc;
T=feedback(Ge,1);
step(T)