开发工具:
文件大小: 2kb
下载次数: 0
上传时间: 2019-08-13
详细说明: MK检验程序出错有请大侠-MK.m % Time Series Trend Detection Tests % % [ z, sl, lcl, ucl ] = trend % % where z = Mann-Kendall Statistic % sl = Sens Slope Estimate % lcl = Lower Confidence Limit of sl % ucl = Upper Confidence Limit of sl % y = Time Series of Data % dt = Time Int erval of Data % % Bob Newell, February 1996 % %-------------------------------------------------- % load unnamed1.mat; y=unnamed1; function [ z, sl, lcl, ucl ] = trendMK % n = length; %-------------------------------------------------- % Mann-Kendall Test for N > 40 % disp; if n < 41, disp; end; % calculate statistic s = 0; for k = 1:n-1, for j = k 1:n, s = s sign - y ); end; end; % variance v = * ) / 18; % test statistic if s == 0, z = 0; elseif s > 0, z = / sqrt; else z = / sqrt; end; % should calculate Normal value here nor = 1.96; % results disp ] ); disp ) ] ); disp ] ); if abs < nor, disp; z = 0; elseif z > 0, disp; else disp; end; %---------------------------------------------------- % disp; % calculate slopes ndash = n * / 2; s = zeros; i = 1; for k = 1:n-1, for j = k 1:n, s = - y ) / / dt; i = i 1; end; end; % the estimate sl = median; disp ] ); % variance v = * ) / 18; m1 = fix ) / 2 ); m2 = fix ) / 2 ); s = sort; lcl = s; ucl = s; disp( [ Lower Confidence Limit = ... num2str ] ); disp( [ Upper Confidence Limit = ... num2str ] ); %---------------------------------------------------- % the end
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台,提供交流渠道,下载内容来自于网络,除下载问题外,其它问题请自行百度。
- 本站已设置防盗链,请勿用迅雷、QQ旋风等多线程下载软件下载资源,下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载,请稍后再次尝试;或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎,请到消费记录里找到下载记录反馈给我们,经确认后退回积分.
- 如下载前有疑问,可以通过点击"提供者"的名字,查看对方的联系方式,联系对方咨询.
相关搜索: