example7.m

一个黑盒对输入输出矩阵进行辨识模型并输出系数

clear all;
N=300;       %观测长度
A=[1 -1.5 0.705];
B=[0 0.065 0.048 0.008];
C=[1 -0.527 0.0695];%ABC为模型的参数
M1=idpoly(A,B,C);         %ARMAX模型
u=idinput(N,'prbs',[0 1]);      %输入伪随机信号
e=NORMRND(0,1,N,1);%随机噪声
y1=sim(M1,[u e]);%模型仿真
t=1:300;
subplot(2,2,1),plot(t,u);title('M');
subplot(2,2,2),plot(t,e);title('随机噪声');
subplot(2,2,3),plot(t,y1);title('输出信号');
Q=IDDATA(y1,u);%获得一批输入输出数据
AIC=zeros(5,5);
for i=1:5    %ARMA(na,nc)模型中的na
    for j=1:i%ARMA(na,nc)模型中的na
        Model=armax(Q,'na',i,'nb',j,'nc',i,'nk',1);%不同阶次系统的辨识
        AIC(i,j)=aic(Model);%辨识结果的AIC准则
    end
    Model     %输出模型
end
AIC            %输出AIC矩阵