轮次检验法

轮次检验法原理轮次检验法是平稳性的非参数检验法，该方法只涉及一组实测数 据，而不需要假设数据的分布规律，因此本方法具有良好的实用性。 轮次（游程）检验法的工作步骤：1）将待检验数据序列xi分成m个等分，并计算各等分的方差q22）计算O = bmax。_2min23）当叶> 02时将该等分记为“ + ”，当叶< 02时将该等分记为“-”4）将 m 个等分按“+”“-”排成观察序列，并将符号相同的连续序列定义为一个轮次（游程）。统计轮次（游程）数 r 作为检验统计量。5） N1 = N （“+”的总数）， N2 = N- （“-”的总数）1 2 -6）当叫和N2 «15时，认为是小样本量，查轮次（游程）检验分布丄厶表，可得在显著性水平a =0.05的上下限r u，rl。7）若r , <r<r，则为平稳性数据lu8）当当叫或N2 >15时，可认为是大样本量，这时可用正态来近似，丄厶可用正态分布表来定出检验的接受域和否定域。计算统计量z =二09) N = N1N2口 = 2N1N2 + 1No =»NiN2(2NiN2-N)寸N2(N-1)10）在a =0.05，若|Z|S1.96,则为平稳性数据。流程图轮次检验的 matlab 程序代码clearx=rand(1,100); %取随机数序列for i=1:25a=x(4*i-3):4*i); %平均分组c(i)=var(a);%求每一组的方差fprintf('c(%d)=%fn',i,c(i); %输出每一组的方差值 endc=c(1) c(2) c(3) c(4) c(5) c(6) c(7) c(8) c(9) c(10) c(11) c(12) c(13) c(14) c(15) c(16) c(17) c(18) c(19) c(20) c(21) c(22) c(23) c(24) c(25);cev=(max(c)+min(c)/2 ; %求方差的均值N1=0;N2=0;for i=1:25if c(i)>cevb='+'q(i)=1;%用来计量观察序列正负号的变化N1=N1+1; %统计“+”的总数elseb='-'q(i)=-1;N2=N2+1; %统计“-”的总数endfprintf(b,q(i);r=1;endfor i=2:1:25if q(i)*q(i-1)=1r=r;elser=r+1;%统计游程数rendendfprintf('r=%dn',r); %输出游程数N1,N2if N1<=15&N2<=15fprintf('请查轮次检验分布表n');%通过查轮次检验用分布表判断数据平稳性elseif N1>15|N2>15 %大样本量，用正太分布近似N=N1+N2;u=2*N1*N2/N+1;cigma=sqrt(2*N1*N2*(2*N1*N2-N)/(NA2*(N-1);z=(r-u)/cigma;fprintf('z=%fn',z);z=abs(z); %取 z 的绝对值if z<=1.96fprintf（'此数据为平稳数据'）;elsefprintf（'此数据为不平稳数据'）;endend参考文献：杨位钦时间序列分析与动态数据建模北京：北京理工 大学，1987吉林大学 汽车工程学院 车辆工程 袁中亮研-03yuanzhongliang09