python mk趋势检验_求问！MK趋势检验和突变检验！

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

分享一段用matlab进行MK趋势的代码给你，希望你能用上

%α=0.05

% Time Series Trend Detection Tests

% [ z, sl, lcl, ucl ] = trend( y, dt )

% where z = Mann-Kendall Statistic

% sl = Sen's Slope Estimate

% lcl = Lower Confidence Limit of sl

% ucl = Upper Confidence Limit of sl

% y = Time Series of Data

% dt = Time Interval of Data

%--------------------------------------------------

function [ z, sl, lcl, ucl ] = trend( y, dt )

n = length( y );

%--------------------------------------------------

% Mann-Kendall Test for N > 40

disp( 'Mann-Kendall Test:' );

if n < 41,

disp( 'WARNING - sould be more than 40 points' );

end;

% calculate statistic

s = 0;

for k = 1:n-1,

for j = k+1:n,

s = s + sign( y(j) - y(k) );

end;

% variance ( assuming no tied groups )

v = ( n * ( n - 1 ) * ( 2 * n + 5 ) ) / 18;

% test statistic

if s == 0,

z = 0;

elseif s > 0,

z = ( s - 1 ) / sqrt( v );

else,

z = ( s + 1 ) / sqrt( v );

end;

% should calculate Normal value here

nor = 1.96;

% results

disp( [ ' n = ' num2str( n ) ] );

disp( [ ' Mean Value = ' num2str( mean( y ) ) ] );

disp( [ ' Z statistic = ' num2str( z ) ] );

if abs( z ) < nor,

disp( ' No significant trend' );

z = 0;

elseif z > 0,

disp( ' Upward trend detected' );

else,

disp( ' Downward trend detected' );

end;

%----------------------------------------------------

python mk趋势检验_求问！MK趋势检验和突变检验！相关推荐

python mk趋势检验_时序数据常用趋势检测方法
背景在最近的项目中,需要自动检测某段时间内的某个指标是上升了还是下降了,因此需要研究下常用的时序数据趋势检测方法. 方法一斜率法原理斜率法的原理就是使用最小二乘等方法对时序数据进行拟合,然后根 ...
python mk趋势检验_【C语言】MK趋势检验C语言代码
[C] 纯文本查看复制代码#include #include void main() { FILE *fp; char filename[40] ; int i,j,k;//循环用 int h,l; ...
如何用python画一个心形图_求问怎样用python/python turtle画“心”
展开全部 python turtle画52614个同心圆方法 import turtle #draw first circle turtle.penup() turtle.goto(0,-200) t ...
笨办法学习python应该看第几版_求问：完全小白学习Python看《笨方法学Python3》还是看《Python编程从入门到实践》?...
完全小白,建议看 <Python编程从入门到实践> 或 <像计算机科学家一样思考Python 第2版> 本书以培养读者以计算机科学家一样的思维方式来理解Python语言编程.贯 ...
python批处理文件编写_求大神用python写个批处理程序？
文本如下,只截了一段. 我想要把注释的内容放到正文里,序号去掉,注释用[]括起来. 例:秦兴师临周[姚本续:周显王,后语.]而求九鼎,周君患之,以告颜率[姚本续:率,名也.当如字.或云,力出切,后语注 ...
普通人学python好吗_有人问，普通人学python真的有意义吗？看看大家都是怎么说的吧...
普通人学python有意义吗? 现在随着python越来越火,尤其是它成为了人工智能的第一编程语言,还被纳入了中小学的教育中.并且python的应用范围很广泛,可以解决很多专业或非专业的问题. 但py ...
简单python代码实例_求简洁优美的python代码例子、片段、参考资料
展开全部建2113议你去看一本书:<计算机程序5261的构造与解释>.里面4102用的语言是Scheme,一种Lisp的方言.通1653过这本书学习程序的抽象.封装,以及重要的函数式编程 ...
python log函数_求你别再花大价钱学 Python 之爬虫实战
引子 Python 基本概念 Python 优势和劣势优势 Python 的劣势 Python 安装设置 Python 基本语法程序例子 Python 基本语法 Python 爬虫实现爬虫相关 ...
量化交易python入门书籍推荐_求量化投资入门书籍或课程推介。？
入门,还是要看一个人的知识背景. 1.有一本:西蒙斯打开量化投资的黑箱,还有西蒙斯的文艺复兴公司的介绍(喜马拉雅),另外B站有西门斯的演讲. 通过对西蒙斯的了解,你要对自己的知识结构.背景.资源和量化 ...

python mk趋势检验_求问！MK趋势检验和突变检验！

python mk趋势检验_求问！MK趋势检验和突变检验！相关推荐

最新文章

热门文章