Matlab_AR模型阶数确定(6)

附录MATLAB 代码

（1）

% AR model order clc;clear;

addpath(genpath(pwd),1);

data = load('Case2Damage0_1_0.01_0.001_40_10_150_123_1.mat'); dofy = [2,4,6,8,10,12,14,16];

x = zscore(data.acc(:,dofy)); order = [10:10:130]; for i = 1:length(order) for j =1:size(x,2)

m = ar(x(:,j),order(i),'yw'); am(i,j) = aic(m); end end

plot(order,am,'LineWidth',2,'Marker','*') xlabel('AR order'),ylabel('AIC')

legend('sensor 2','sensor 4','sensor 6','sensor 8',... 'sensor 10','sensor 12','sensor 14','sensor 16')

（2）

clc;clear;

addpath(genpath(pwd),1);

data = load('Case2Damage0_1_0.01_0.001_40_10_150_123_1.mat'); dofy = [2;4;6;8;10;12;14;16]; order = 80;

ts = zscore(data.acc(:,dofy)); k = 1;

name = {'sensor 2','sensor 4','sensor 6','sensor 8',... 'sensor 10','sensor 12','sensor 14','sensor 16'}; Ounit = repmat({'m/s^2'},length(dofy),1);