- 积分
- 80
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-7-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
M=load('D:\text3.txt');
x=M(:,1);y=M(:,2);t=zeros(55,1);y1=zeros(55,1);
%MMT方法检验%
%先设定子序列的长度(取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:50
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:55),t,'r-');
y1(:)=2; %n=30-2时,显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1961,2015,-5,5]);
xlabel('年');ylabel('t');title('MMT(n=5)检验');
所用程序时根据网上别人的改的,在此有1961年到2015年降水数据,画出的图 2000年以后为0,Matlab小白想请问该怎么修改程序呢,还是哪里出错了
谢谢
|
|