ck1k5.m

来自「状态反馈matlab仿真程序」· M 代码 · 共 26 行

N=1;

a_ITAE(1,1)=[1.41];
a_ITAE(2,1:2)=[2.15 1.75];
a_ITAE(3,1:3)=[2.7 3.4 2.1];
a_ITAE(4,1:4)=[3.4 5.5 5.0 2.8];
a_ITAE(5,1:5)=[3.95 7.45 8.6 6.6 3.25];     
a_ITAE(6,1:6)=[4.58 10.64 13.54 15.05 10.42 4.475];
a_ITAE(7,1:7)=[5.15 13.30 22.20 25.75 21.60 12.80 5.20];

aop=a_ITAE(N,1:N);
l0=1;
for (i=1:1:N)
   l(i)=tao2^i/fac(i);
end

%k(N+1)=(aop(N)*T-l(N-1)*T/l(N)-1)/K,
%for (i=N:-1:3)
%   k(i)=(aop(i-1)*l(N)*T-K*k(N+1)*l(i-1)-T*l(i-2)-l(i-1))/K;
%end
%k(2)=(aop(1)*l(N)*T-K*k(N+1)*l(1)-T*l0-l(i-1))/K;
k(N+1)=(aop(N)*T-l0*T/l(N)-1)/K,
k(1)=(l(N)*T-K*k(N+1)*l0-l0)/K,