shujuyuchuli.m

用于随机数据的数据预处理

%采样数据
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)模型
%模型参数估计