- n=length(x);
- for i=1:fix(n/2)
- a(i)=floor(n/i);
- y=reshape(x(1:a(i)*i),i,a(i));
- mx=mean(y);
- cm=y-ones(i,1)*mx;
- cum=cumsum(cm);
- s=std(y);
- rs(i)=(max(cum)-min(cum))/s;
- end
- loglog(a,rs,'*')
|
楼主: 盗梦人
|
1146
0
R/S运行的小程序 |
加好友,备注jltj京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明


