写了一个M-K趋势分析的M程序，大家帮我看看

lml7980

1. <blockquote>%% M-K趋势分析

复制代码

有点疑问：得到的Z值说明了增加或减少的趋势，我试着用纯增加趋势的数据来验证，结果
只能通过90%信度检验，更高的却无法通过了，为什么会这样呢？？

lml7980

%% M-K趋势分析
clear,clc
% x=rand(30,1);
x=data;          %写入data,要求为一列数据

n=length(x);
for i=1:(length(x)-1);
  if (x(i+1)-x(i))>0;
      b(i)=1;
          else if (x(i+1)-x(i))<0;
           b(i)=-1;
             else
          b(i)=0;
         end
      end
  end
S=sum(b);

   if S>0;
   Z=(S-1)./sqrt((n*(n-1)*(2*n+5)./18))
       else if S==0;
           Z=0
             else
                  Z=(S+1)./(n*(n-1)*(2*n+5)./18)^(1/2)
         end
    end
  if abs(Z)>=1.28;
      disp('通过信度90%显著性检验')
  else if abs(Z)>=1.64;
     disp('通过信度95%显著性检验')
   else if abs(Z)>=2.32;
     disp('通过信度99%显著性检验')  
   else  
     disp('无法通过信度显著性检验')
         end
      end
  end
  
  


                  

denghaijun

M-k检验，1.64对应0.1, 1.96对应0.05, 2.576对应0.01；
“if abs(Z)>=1.28;
      disp('通过信度90%显著性检验')
  else if abs(Z)>=1.64;
     disp('通过信度95%显著性检验')
   else if abs(Z)>=2.32;
     disp('通过信度99%显著性检验')  
   else  
     disp('无法通过信度显著性检验')
         end
      end
  end”
改成：
if  abs(Z)>=1.64 && abs(Z)<1.96 ;
      disp('通过信度90%显著性检验')
  else if abs(Z)>=1.96 && abs(Z)<2.576 ;
     disp('通过信度95%显著性检验')
   else if abs(Z)>=2.576 ;
     disp('通过信度99%显著性检验')  
   else  
     disp('无法通过信度显著性检验')
         end
      end
  end

lml7980

denghaijun 发表于 2015-7-29 22:41
M-k检验，1.64对应0.1, 1.96对应0.05, 2.576对应0.01；
“if abs(Z)>=1.28;
      disp('通过信度90%显著 ...

非常感谢！！

忍冬。

谢谢楼主分享！等我试验后再反馈~

大花园

表示看不懂

qinjin@

怎么没有最后出图的代码

qinjin@

qinjin@ 发表于 2018-4-12 22:20
怎么没有最后出图的代码

同求出图代码

sciren

楼主  为啥我输入1-30  不管怎么输入 都是不通过检验啊