test_rels.m

系统辨识的最小二乘递推算法、辅助变量法、增广最小二乘法及偏差补偿法的matlab程序设计实例。

load input data_in;
load output data_out;

data_sum = size(data_in);
U0 = data_in(2,1:data_sum(2));
Z0 = data_out(2,1:data_sum(2));

m(1)=input('order of A(z)');
m(2)=input('order of B(z)');
m(3) =input('order of D(z)');                            
n = max(m);

Len=size(Z0);  
L = Len(2) - n;
P = 10e6*eye(m(1)+m(2)+m(3));                       
sita = 0.001*ones(m(1)+m(2) +m(3),1);                  
h = zeros(m(1)+m(2) +m(3),1);                         
z = 0;
MatixI = eye(m(1)+m(2) +m(3));                       
V = zeros(1, Len(2));

for i=1:L
    for j=1:m(1)
        h(j,1) = -Z0(n+i-j);
    end;
    for j=1:m(2)
        h(j+m(1),1) = U0(n+i-j);
    end;
    for j=1:m(3)
        h(j+m(1)+m(2),1) = V(n+i-j);                         
    end;

z = Z0(n+i);
    K = P*h/(h'*P*h+1);                                
    sita = sita + K*(z - h'*sita);
    P = (MatixI - K*h') * P;              
V(n+i) = z - h'*sita;                               
end

sita