《轮次检验法》由会员分享，可在线阅读，更多相关《轮次检验法（4页珍藏版）》请在装配图网上搜索。
1、轮次检验法原理轮次检验法是平稳性的非参数检验法，该方法只涉及一组实测数 据，而不需要假设数据的分布规律，因此本方法具有良好的实用性。 轮次（游程）检验法的工作步骤：1）将待检验数据序列xi分成m个等分，并计算各等分的方差q22）计算O = bmax。_2min23）当叶 02时将该等分记为“ + ”，当叶 02时将该等分记为“-”4）将 m 个等分按“+”“-”排成观察序列，并将符号相同的连续序列定义为一个轮次（游程）。统计轮次（游程）数 r 作为检验统计量。5） N1 = N （“+”的总数）， N2 = N- （“-”的总数）1 2 -6）当叫和N2 15时，认为是小样本量，查轮次（游程）
2、检验分布丄厶表，可得在显著性水平a =0.05的上下限r u，rl。7）若r , r15时，可认为是大样本量，这时可用正态来近似，丄厶可用正态分布表来定出检验的接受域和否定域。计算统计量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)
3、; %输出每一组的方差值 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; %统计“-”的总数end
4、fprintf(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&N215|N215 %大样本量，用正太分布近似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
轮次检验法
