用Matlab进行MK趋势分析与突变检验.pdf

用Matlab进行MK趋势分析与突变检验.pdf

ID:56755657

大小:65.85 KB

页数:4页

时间:2020-07-07

用Matlab进行MK趋势分析与突变检验.pdf_第1页
用Matlab进行MK趋势分析与突变检验.pdf_第2页
用Matlab进行MK趋势分析与突变检验.pdf_第3页
用Matlab进行MK趋势分析与突变检验.pdf_第4页
资源描述:

《用Matlab进行MK趋势分析与突变检验.pdf》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、%M-K趋势检定clear,closeall,clc%clear:清变数closeall:清图面clc:清画面%defineandassignthefullfilepathusing"fileopen"dialog[filenamefilepath]=uigetfile('data1.xls');full_filepath=[filepathfilename];[X,TXTX,RAWX]=xlsread(full_filepath,1);%数据全部读入,数据缺失不影响结果x=X(:,1);%x时间y=X(:,2);%y数据%计算Sn=size(y,1);%数据个

2、数S=0;fori=1:n-1S=S+sum(sign(y(i+1:n)-y(i)));%S计算式end%计算VarSVarS=n*(n-1)*(2*n+5)/18;%计算ZifS>0Z=(S-1)/sqrt(VarS);elseZ=(S+1)/sqrt(VarS);end%计算Zabsalpha1=0.05;%信度95%的显著水平alpha2=0.01;%信度99%的显著水平PZ1=norminv(1-alpha1/2,0,1);PZ2=norminv(1-alpha2/2,0,1);H=0;%虚无假设Zabs=abs(Z);ifZabs>=PZ1H=1;elseH=0;endP

3、_value=2*(1-normcdf(abs(Z),0,1));%若P_value比alpha1小,则否定虚无假设%计算倾斜度ndash=n*(n-1)/2;%对称矩阵上半部slope1=zeros(ndash,1);%起始归零m=0;fork=1:n-1,forj=k+1:n,m=m+1;slope1(m)=(y(j)-y(k))/(x(j)-x(k));%分母非(j-k)end;end;slope=median(slope1);%中位数%历线绘图yd=max(y)-min(y);figureplot(x,y,'b-o','linewidth',1.5);axis([min(x

4、),max(x),min(y)-0.2*yd,max(y)+0.2*yd]);%全距外扩20%xlabel('时间','FontName','TimesNewRoman','FontSize',12);ylabel('数据','FontName','TimesNewRoman','Fontsize',12);title('数据历线图')%添加标题gridonoutput='数据历线图';saveas(gcf,output,'jpg')%M-K突变检定Sk=zeros(size(y));%起始归零UFk=zeros(size(y));%起始归零s1=0;fori=2:nforj=1:

5、iify(i)>y(j)s1=s1+1;elses1=s1+0;end;end;Sk(i)=s1;E=i*(i-1)/4;%均值Var=i*(i-1)*(2*i+5)/72;%方差UFk(i)=(Sk(i)-E)/sqrt(Var);end;%起始归零y2=zeros(size(y));Sk2=zeros(size(y));UBk=zeros(size(y));s2=0;fori=1:ny2(i)=y(n-i+1);%逆序end;fori=2:nforj=1:iify2(i)>y2(j)s2=s2+1;elses2=s2+0;end;end;Sk2(i)=s2;E=i*(i-1)/

6、4;%Sk2(i)的均值Var=i*(i-1)*(2*i+5)/72;%Sk2(i)的方差UBk(i)=-(Sk2(i)-E)/sqrt(Var);end;UBk2=zeros(size(y));fori=1:nUBk2(i)=UBk(n-i+1);%逆序end;%线性回归x1=x-x(1)+1;%x1可为非连续时间序列,有缺失数据无所谓x1,非x1=[1:n]'r=corrcoef(x1,y)%相关系数R2=r(1,2)^2C=polyfit(x1,y,1)%C(1):一次项系数C(2):常数项系数%画UFk,UBkM-K统计量曲线图dFB=max(max(UFk)-m

7、in(UFk),max(UBk2)-min(UBk2));dFB1=min(min(UFk),min(UBk2))-0.2*dFB;%全距外扩20%dFB2=max(max(UFk),max(UBk2))+0.2*dFB;ifdFB1>-PZ2dFB1=-5;endifdFB2

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。