els.m

关于系统辨识的一些仿真程序

num=[1 6.5 14 11.5 3];
den=[1 20 120 232 320];
u=mseries(60);
[A,B]=getval(num,den);
[e,dif,c,z,zm,zs]=elsm(A,B,u);
i=1:60;
figure(1);
plot(i,dif,'red')
figure(2);
for a=1:14
    subplot(4,4,a);
    plot(i,c(a,:),'black')
end