lqr.m

control souce code with LQR

A=[0 1;0 -1];
B=[0;1];
F=[10 7.426];
syms s;
G=F*inv(s*eye(2)-A)*B;
[d n]=numden(G);
num=sym2poly(d);
den=sym2poly(n);
H = tf(num,den);
nyquist(H);
%bode(H);
figure(2)
margin(H);
grid on;