ts = iddata(zscore(data.acc(:,dofy)),[],data.dt,... 'OutputName',name,'OutputUnit',Ounit,... 'Name','ASCE-benchmark 120DOF');
m = ar(ts(:,k),order,'yw');
compare(ts(:,k),m,'r-.',1,1000:1300);
(3)
resid(m,ts(:,k),'corr',25)
[e,r] = resid(m,ts(:,k),'corr',40);
resid(r)
plot(e(1000:1300))
参考文献
Matlab help
[1] DA SILVA S, DIAS J NIOR M, LOPES JUNIOR V. Damage detection in a benchmark structure using AR-ARX models and statistical pattern recognition [J]. J Braz Soc Mech Sci Eng, 2007, 29(2): 174-84.