example5_1.m

来自「系统辨识的M文件」· M 代码 · 共 17 行

clear all; hold off;   N=1000; 
% A=[ 1  -1.5  0.7 ];   B=[0 1 0.5];    C=[1  0.5];  D=[1  0.5];   
U = idinput(N,'prbs',[0 1])*2;                                      % 输入伪随机信号
E = NORMRND(0,0.9,N,1);                                             % 随机噪声
 ek=zeros(N, 1);   y=zeros(N, 1);                                   % 变量初始化
for i=3:N
    ek(i)=[E(i) ek(i-1)]*[1  -0.8]';                                % 模拟观测噪声
    y(i)=[-y(i-1) -y(i-2) U(i-1) U(i-2)]*[-1.5 0.7 1.0 0.5]'+ek(i); % 模拟观测值  
end
z=iddata(y, U);                                                     % 输入输出数据组
Thita=rels(z, 2, 2, 2);                                             % 调用rels 函数
x=1:N;
plot(x,Thita(1,x),'r'); hold on;  plot(x,Thita(2,x),'r');           % 作图显示参数递推值
plot(x,Thita(3,x),'b');  plot(x,Thita(4,x),'b');
plot(x,Thita(5,x),'g');  hold off;