example5_1.m

来自「一个黑盒对输入输出矩阵进行辨识模型并输出系数」· M 代码 · 共 23 行

clear all; hold off;   N=400; 
% A=[ 1  -1.5  0.7 ];   B=[0 1 0.5];    C=[1  0.5];  D=[1  0.5];   
n=4;
a=1;
n1=4;
n2=1;
mprbs(n,a,n1,n2);
U =idinput(N,'prbs',[0 1])*1;                               % 输入伪随机信号
subplot(1,2,1);stem(U);title('系统输入序列');
E =NORMRND(0,0.9,N,1);                                             % 随机噪声
 ek=zeros(N, 1);   y=zeros(N, 1);                                   % 变量初始化
for i=3:N
    ek(i)=[E(i)]*[1]';                                % 模拟观测噪声
    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
subplot(1,2,2);stem(y);title('系统输出Y');
z=iddata(y, U);                                                     % 输入输出数据组
Thita=rels(z, 2, 2,0);                                             % 调用rels 函数
x=1:N;figure
plot(x,Thita(1,x),'r'); hold on;  plot(x,Thita(2,x),'m');           % 作图显示参数递推值
plot(x,Thita(3,x),'b');  plot(x,Thita(4,x),'y');
  hold off;
  grid on