- 积分
- 180
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-7-24
- 最后登录
- 1970-1-1

|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
M=load('wyj.txt');
x=M(:,1);y=M(:,2);t=zeros(30,1);y1=zeros(30,1);
%MMT方法检验%
%先设定子序列的长度(取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:25
averx1=mean(y(i-4:i));
averx2=mean(y(i+1:i+5));
s1=var(y(i-4:i));
s2=var(y(i+1:i+5));
s=sqrt((n1*s1+n2*s2)/(n1+n2-2));
t(i-4)=(averx2-averx1)/(s*m);
end
plot(x(1:30),t,'r-');
y1(:)=2; %n=30-2时,显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1971,2000,-6.5,3]);
xlabel('年');ylabel('t');title('MMT(n=30)检验');
我这是1971-2000年30年的气象资料,请各位大神帮我看看第二行括号里的冒号那怎么改,还有检验水平怎么看,非常感谢!!!
|
评分
-
查看全部评分
|