rgls.m

来自「关于系统辨识的一些仿真程序」· M 代码 · 共 22 行

u=mseries(100);
a=length(u);
num=[1 6.5 14 11.5 3];
den=[1 20 120 232 320];
[A,B]=getval(num,den);
num1=[1];
den1=den;
[C,D]=getval(num1,den1);
[nois,e,cf,ce]=rglsm(A,B,C,u);
i=1:a;
figure(1);
plot(i,nois,'red')
figure(2);
for a=1:9
    subplot(3,3,a);
    plot(i,cf(a,:),'black')
end
figure(3);
for b=1:4
    subplot(2,2,b);
    plot(i,ce(b,:),'green')
end