第1个回答 2015-01-01
例一
y = sin([1:300]') + 0.5*randn(300,1);
y = iddata(y);
mb = ar(y,4,'burg')
mfb = ar(y,4)
bode(mb,mfb)
例二
% 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')本回答被网友采纳