rls.m

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

num=[1 6.5 14 11.5 3];
den=[1 20 120 232 320];
digits(64);
u=mseries(200);
[A,B]=getval(num,den);
x1=length(A);
x2=length(B);
[e,c,z,dif]=rlsm(u,A,B);
i=1:200;
figure(1);
plot(i,dif,'black')
figure(2);
for j=1:9
    subplot(3,3,j);
    plot(i,c(j,:),'red');
end