mk检验的uf和ub表示什么意思
sisi韩
2024-05-10 16:31:06
最佳回答
uf,ub]=mannkendall(y);% % y即为待检测的时间序列figure,hold on;plot(uf,'-b');plot(ub,'-r');function [uf,ub]=mannkendall(timeseries)n=length(timeseries);uf=smk(timeseries,n);for i=1:nyy(i)=timeseries(n+1-i);endu_res=smk(yy,n);for i=1:nub(i)=-u_res(n+1-i);endfunction u_res=smk(y,n)m_res=zeros(n,1);md_res=zeros(n,1);u_res=zeros(n,1);m_res(1)=0;for i=2:nm_res(i)=0;md_res(i)=0;for j=1:i-1if y(i)<y(j)m_res(i)=m_res(i)+0;elsem_res(i)=m_res(i)+1;endmd_res(i)=md_res(i-1)+m_res(i);endendu_res(1)=0;for i=2:ne=i*(i-1)/4;var=i*(i-1)*(2*i+5)/72;u_res(i)=(md_res(i)-e)/sqrt(var);end 20210311