ck1k2.m

状态反馈matlab仿真程序

omiga0=sqrt(0.1767);
aop1=1.9222/omiga0;
printsys([1],[1 aop1*omiga0 omiga0^2])
l0=1;
l1=tao3;

k(2)=(aop1*omiga0*T-l0*T/l1-1)/K,
k(1)=(omiga0^2*l1*T-K*k(2)*l0-l0)/K,