用Matlab进行MK检验
6页1、用 Matlab 进行 Mann-Kendall 分析一、M-K趋势分析法与M-K突变检验的原理介绍1、Mann-Kendall趋势分析法非参数检验,又称任意分布检验, 它不对变量的分布做严格假定,检验不针对特定的参数,而是模糊地对变量分布的中心位置或分布状态做检验,由于其不对总体分布做严格假定,因而适用性强12。因此,本文采用非参数的 Mann-Kendall检验法对昌马河流域近 50年的气 候水文要素时间序列显著性检验,定量反映变化趋更计算公式如下:S 1 .var(S), Zc0,(1.1)S 1 var( S),Xi)(1.2)在公式中, n 1 n Ssgn(Xki 1 k i 11,0sgn()0,01,0(1.3)varS n(n1)(2n5)t(tt1)(2t5)18(1.4)式中:Xk, Xi为连续的气候、水文等数据序列,n为数据集合的总长度,t为每个单位的宽精品资料度,2表示所有单位的总和。衡量趋势大小的指标为:Xi Median xj(1.5)式中:1jix(m)s=s+1;elseif x(n)=x(m)s=s+0;elses=s-1;endendendvarS
2、 n(n 1)(2n 5) t(t 1)(2t 5) 18%ci算 vars 和 e%t/vars=len*(len-1)*(2*len+5)/18;ZcS 1 , var( S)0, S 1 . var( S)if s0zc=(s-1)/sqrt(vars);elsezc=(s+1)/sqrt(vars);end%十算zaza=var(x);sign=0;zc1=abs(zc);if zc1 = zasign=1;elsesign=0;end% 十算倾斜度ndash = len * ( len - 1 ) / 2; slope1= zeros( ndash, 1 ); m=1;for k = 1:len-1,for j = k+1:len,Xi XjMedian , j islope1(m) = ( x(j) - x(k) ) / ( j - k ) ; % m = m + 1;end;end;slope= median( slopel );2、M-K突变检验clc;y=data(1:50,6);%输入数据n=length(y);Sk=zeros(size(y); UFk=zeros
3、(size(y);s1=0;kSkri (k 2,3,n)i 1for i=2:n%for j=1:iif y(i)y(j) s1=s1+1;elseri1,XiXjs1=s1+0;%end;0,XXj(j 1,2,L,i)n (n 1)E(Sk)end;Sk(i尸s1;%E=i*(i-1)/4;Var=i*(i-1)*(2*i+5)/72;UFk(i)=(Sk(i)-E)/sqrt(Var);UF k回E(Sk)Var (sQ(k 1,2,n)end;y2=zeros(size(y);Sk2=zeros(size(y);UBk=zeros(size(y);s2=0;for i=1:ny2(i)=y(n-i+1);end;for i=2:nfor j=1:iif y2(i)y2(j)s2=s2+1;elses2=s2+0;end;end;Sk2(i)=s2;的方差 Va%)也”E=i*(i-1)/4; % Sk2(i) 的均值Var=i*(i-1)*(2*i+5)/72; % Sk2(i)UBk(i)=0-(Sk2(i)-E)/sqrt(Var); end;UBk2=zeros(size(y);for i=1:nUBk2(i)=UBk(n-i+1);end;year=1961:2010;year=year;xlswrite(D:test.xls,year,Sheet1,A1);xlswrite(D:test.xls,UFk,Sheet1,B1);xlswrite(D:test.xls,UBk2,Sheet1,C1);Welcome ToDownload !欢迎您的下载,资料仅供参考!
《用Matlab进行MK检验》由会员枫**分享,可在线阅读,更多相关《用Matlab进行MK检验》请在金锄头文库上搜索。
小学生关于我的中国梦演讲稿500字
互感器伏安特性测试仪
诸城电网运行维护分界点及运行维护管理责任划分规定
六年级数学下册 第一单元《圆柱和圆锥》单元综合测试题1(无答案)北师大版
X实业股份有限公司课心得
金蝶协同办公软件V62用户手册-知识管理篇
信阳操作系统服务项目申请报告
音乐教案《我们祖家是歌乡》
模拟卷厦门市湖里科技创业服务中心入驻协议书
衡水关于成立CPU公司可行性报告(模板)
第四代移动通信技术的特点
个人的工作安排三篇
2023年工程工作计划范文(七篇).doc
论如何提高我国物流企业的竞争能力
小学春季开学典礼学生讲话稿(2022年开学典礼学生讲话稿6篇)
Word绝招和Excel常用函数大全
学生安全教育记录表格参考
公司前台年终工作总结范文(4篇).doc
英语教学工作总结范本(2篇).doc
土地增值税鉴证报告
2023-08-01 27页
2022-11-24 4页
2023-01-04 3页
2023-01-17 11页
2022-11-29 4页
2022-12-27 16页
2022-10-15 3页
2022-08-30 9页
2022-10-10 9页
2023-01-14 54页