example7_3.m

系统辨识的M文件

clear all;
N=300;                                                %观测数据长度
M1 = idpoly([1 -2.851 2.717 -0.865],[0 1 1 1],[1 0.7 0.2]); % ARMAX模型
u = idinput(N,'prbs',[0 1]);                          % 输入伪随机信号
e = NORMRND(0,1,N,1);                                 % 随机噪声
y1 = sim(M1,[u e]);                                   % 模型仿真
Q= IDDATA(y1,u);                                      % 获得一批输入输出数据   
AIC=zeros(5,5);  
for i=1:5                                             % ARMA(na,nb)模型中的na
    for j=1:i                                         % ARMA(na,nb)模型中的nb
        Model=armax(Q,'na',i,'nb',j,'nc',i,'nk',1);   % 根据模型类和输入输出数据实现不同阶次系统的辨识
        AIC(i,j)=aic(Model);                          % 辨识结构的AIC准则
    end
    Model                                             % 输出模型
end
AIC                                                   % 输出AIC矩阵