m-k检验matlab程序出错，求大神指点！

DGT

% Mann-Kendall突变检测 % 数据序列y% 结果序列UFk，UBk2%--------------------------------------------%读取excel中的数据，赋给矩阵y%获取y的样本数%A为时间和气温数据列A=xlswrite('d:\data.xls')x=A(:,1);%时间序列y=A(:,2);%气温数据列N=length(y);n=length(y);% 正序列计算---------------------------------% 定义累计量序列Sk，长度=y，初始值=0Sk=zeros(size(y));% 定义统计量UFk，长度=y，初始值=0UFk=zeros(size(y));% 定义Sk序列元素ss = 0;% i从2开始，因为根据统计量UFk公式，i=1时，Sk(1)、E(1)、Var(1)均为0% 此时UFk无意义，因此公式中，令UFk(1)=0for i=2:n   for j=1:i         if y(i)>y(j)           s=s+1;         else           s=s+0;         end;   end;   Sk(i)=s;   E=i*(i-1)/4; % Sk(i)的均值  Var=i*(i-1)*(2*i+5)/72; % Sk(i)的方差  UFk(i)=(Sk(i)-E)/sqrt(Var);end;% ------------------------------正序列计算end% 逆序列计算---------------------------------% 构造逆序列y2，长度=y，初始值=0y2=zeros(size(y));% 定义逆序累计量序列Sk2，长度=y，初始值=0Sk2=zeros(size(y));% 定义逆序统计量UBk，长度=y，初始值=0UBk=zeros(size(y));% s归0s=0;% 按时间序列逆转样本y% 也可以使用y2=flipud(y);或者y2=flipdim(y,1);for i=1:n    y2(i)=y(n-i+1);end;% i从2开始，因为根据统计量UBk公式，i=1时，Sk2(1)、E(1)、Var(1)均为0% 此时UBk无意义，因此公式中，令UBk(1)=0for i=2:n   for j=1:i         if y2(i)>y2(j)           s=s+1;         else           s=s+0;         end;   end;   Sk2(i)=s;   E=i*(i-1)/4; % Sk2(i)的均值  Var=i*(i-1)*(2*i+5)/72; % Sk2(i)的方差% 由于对逆序序列的累计量Sk2的构建中，依然用的是累加法，即后者大于前者时s加1，% 则s的大小表征了一种上升的趋势的大小，而序列逆序以后，应当表现出与原序列相反% 的趋势表现，因此，用累加法统计Sk2序列，统计量公式(S(i)-E(i))/sqrt(Var(i))% 也不应改变，但统计量UBk应取相反数以表征正确的逆序序列的趋势  UBk(i)=0-(Sk2(i)-E)/sqrt(Var);end;% ------------------------------逆序列计算end% 此时上一步的到UBk表现的是逆序列在逆序时间上的趋势统计量% 与UFk做图寻找突变点时，2条曲线应具有同样的时间轴，因此% 再按时间序列逆转结果统计量UBk，得到时间正序的UBk2，做图用UBk2=zeros(size(y));% 也可以使用UBk2=flipud(UBk);或者UBk2=flipdim(UBk,1);for i=1:n   UBk2(i)=UBk(n-i+1);end;% 做突变检测图时，使用UFk和UBk2% 写入目标xls文件：d:\test2.xls% 目标表单：Sheet1% 目标区域：UFk从A1开始，UBk2从B1开始B=xlswrite('d:\test2.xls',UFk,'Sheet1','A1');B=xlswrite('d:\test2.xls',UBk2,'Sheet1','B1');figure(3)%画图plot(x,UFk,'r-','linewidth',1.5);hold onplot(x,UBk2,'b-.','linewidth',1.5);plot(x,1.96*ones(N,1),':','linewidth',1);axis([min(x),max(x),-5,5]);legend('UF统计量','UB统计量','0.05显著水平');xlabel('t (year)','FontName','TimesNewRoman','FontSize',12);ylabel('统计量','FontName','TimesNewRoman','Fontsize',12);%grid onhold onplot(x,0*ones(N,1),'-.','linewidth',1);plot(x,1.96*ones(N,1),':','linewidth',1);plot(x,-1.96*ones(N,1),':','linewidth',1)

出错提示为Error using ==> xlswriteNot enough input arguments.

DGT

大家还是不要看这个帖子了，我自己看了都受不了，不知道为什么会变成这个样子，好好的一个程序，唉！

kongfeng0824

那是因为你没有排版好，直接摘过来的吧？

DGT

kongfeng0824 发表于 2014-5-12 20:17
那是因为你没有排版好，直接摘过来的吧？

从txt粘过来的，不过问题已经解决了！

小鸽子1991125

DGT 发表于 2014-5-12 21:01
从txt粘过来的，不过问题已经解决了！

弱弱的问一句，您这个程序没有错误，可以直接用吗？我也要用M-K检验，我是新手

小鸽子1991125

可以直接用么？我也要用这个检验，我是新手，不太懂

DGT

小鸽子1991125 发表于 2014-5-13 08:56
可以直接用么？我也要用这个检验，我是新手，不太懂

这个程序有个错，改对了应该可以直接出图，但我不知道错在哪，后来用别的方法做的

amanda_rui

也好想直接用啊

kevinka

国标《农业干旱预警等级》（征求意见稿）附录B有标准

木林森(lynn)

重新排版吧