- 积分
- 359
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-2
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
>> M=load('D:\matlab\气温.txt');
x=M(:,1);y=M(:,2);t=zeros(30,1);y1=zeros(38,1);
%MMT方法检验%
%先设定子序列的长度(取n1=5;n2=5)%
n1=5;n2=5;
m=sqrt((1.0/n1)+(1.0/n2));
%开始计算%
for i=5:34
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(5:34),t,'r-');
y1(:)=32; %n=14时,显著性检验标准,99.5%
hold on;
plot(x,y1,'k:',x,-y1,'k:');
axis([1961,2010,-5,5]);
xlabel('年');ylabel('t');title('MMT(n=5)检验');
??? Error using ==> plot
Vectors must be the same lengths.
出现这个错误,是为什么呢?如何解决?
|
|