%采樣數據
junzhi=mean(sx);
changdu=length(sx);
t=1:changdu;
h=figure(1);
plot(t,sx)
%去常數項
danweizhen=ones(1,changdu);
chazhi=sx-junzhi*danweizhen';
h=figure(2);
plot(t,chazhi)
% 顯示趨勢項
a=polyfit(t,chazhi',2);    %最小二乘擬合2階擬合%
z=polyval(a,t);
h=figure(3);
plot(t,chazhi,':r',t,z,'b')
%去趨勢項
ququshixiang=chazhi-z';
h=figure(4);
plot(t,ququshixiang)
%平穩性檢驗 ,采用輪次檢驗
for i=0:1:199;        %分組為200個區間,一個區間220個數據
    Xi=ququshixiang(i*220+1:1:(i+1)*220);
    XXi(i+1)=std(Xi);    %對每個分區取標準差
end
 zhongzhi=median(XXi);
h=figure(5);
   plot(1:200,XXi,1:200,zhongzhi*ones(1,200),'r')
   %查交點,對照表A6看是否滿足平穩性假設
%正態性檢驗
h=figure(6) ;
histfit(ququshixiang)   %觀察數據組的頻數直方圖與正態數據的對比
%未進行正態性檢驗
%---------------------
%-----------------------
%選AR(2)模型
%模型參數估計