(转)
ALPHA=.01; % Pr[type 1 error] for threshold line to test acf for N(0,1/sqrt(n))
pvec=[0,1,2]; % orders of AR
%pvec=[4,5,6];
d=0; % number of differences
qvec=[0,1,2]; % orders of MA
%qvec=[2,3];
RESIDPLOTS=0; % > 0 to plot residual acfs
%dataname='data2(:,1)'; % assumed loaded
%descriptor='US Production'; % strings to be used for plotting
%xvalues='quarters';
%dindex=[1:120]; % select which of the input data to use
brillx=textread('brillx.txt','%f'); % this puts 912 elements in a column
%brilly=textread('brilly.txt','%f');
dataname='brillx'; %%%%%%%%%%%%% dataname statement %%%%%%%%%%%%%%%%%%%
descriptor='WOBBLE';
xvalues='time';
dindex=[