【趋势分析方法三】MATLAB代码实现TFPW-MK检验

news2024/9/28 21:26:06

目前水文时间序列趋势分析的方法很多,主要分为参数检验非参数检验两大类:

  • 参数检验中常用的有线性回归法、滑动平均法、累积距平法等
  • 非参数检验则主要包括Mann-Kendal(MK)法和 Spearman 秩次相关法等

虽然从理论上讲,参数检验法较非参数检验可获得更有效的检验结果,但由于水文数据序列存在非同一分布、缺失值或异常值、季节性变化、自相关性等诸多问题,使得参数检验法的使用在水文趋势检验中受到了诸多限制。而非参数检验法凭借其不受样本值分布类型影响等特点,目前已被广泛应用于水文时间序列的趋势检验领域。其中最为常用的检验方法,即Mann-Kendal(MK)法。
利用Mann-Kendall(MK)方法进行水文序列趋势检验时, 去趋势预置白(Trend-free pre-whitening, TFPW) 作为处理水文序列自相关性影响的重要方法之一,其处理过程的合理性及方法的适用性在变化环境下备受关注。

1 原理

在TFPW-MK方法中,为了消除序列自相关性的影响,需要对序列预先进行去趋势处理,即实现趋势成分和非趋势成分的有效分离,以保证预置白处理不致趋势成分丢失。
在这里插入图片描述

1.1 TFPW-MK趋势检验法

TFPW-MK趋势检验法是针对待检序列的自相关性问题,提出的一种前置移除型的MK 趋势检验方法,前置程序包括去趋势预置白两个处理过程,可有效降低序列中自相关性对检验结果的影响,避免检验结果失真。

1.2 去趋势方法

在传统的TFPW-MK 趋势检验法中,采用了Theil-Sen 法,即线性去趋势法—坡度法(SM 法)进行处理。

1.2.1 经验模态分解法(Empirical Mode Decomposition, EMD)

经验模态分解法简称EMD 法,是一种针对非线性数据的处理方法。其基本思想是将原始序列不同尺度的波动信号逐步分解,则得到一系列平稳性较好的固有模态函数(Intrinsic Mode Function,IMF)和具有单一性的残余项(Residual)两部分。在水文领域中,一般近似认为IMFs 表征了水文序列中的随机成分和周期成分的变化,而常将残余项视为趋势成分。

1.2.2 一阶差分法(First Difference, FD)

一阶差分法是统计学中最基础的去趋势方法,简称FD法。一阶差分即序列中连续相邻的两项之差。

1.2.3 对数线性法(Log-linear detrending, LLD)

1.2.4 直线滑动平均法(Linear moving average, LMA)

1.2.5 直线滑动回归平均法(Linear regressive moving average, LRMA)

参考

1.论文-J2002-The influence of autocorrelation on the ability to detect trend in hydrological series
2.论文-J2016-去趋势预置白方法对径流序列趋势检验的影响
3.Mann-Kendall(MK)去趋势预置白(TFPW)水文气象序列的趋势和突变分析
4.R工具箱-modifiedmk
modifiedmk
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/366508.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

2023年,什么行业更有发展前景?

关于有前景有发展的行业推荐,小课今天还是推荐咱们IT互联网行业。 很多人会说现在懂电脑的那么多,这个行业都饱和了,很多学电脑的找不到工作都改行了。但事实是现在每个行各业都需要互联网,需要懂电脑的技术人才,尤其是在云计算、大数据到来…

了解一下TCP/IP协议族

在《简单说说OSI网络七层模型》中讲到,目前实际使用的网络模型是 TCP/IP 模型,它对 OSI 模型进行了简化,只包含了四层,从上到下分别是应用层、传输层、网络层和链路层(网络接口层),每一层都包含…

《安富莱嵌入式周报》第304期:开源硬件耳机设计,AI单片机STM32N6已确定为M55内核,另外还有新品STM32H5, H50X, H7R, H7S发布

往期周报汇总地址:嵌入式周报 - uCOS & uCGUI & emWin & embOS & TouchGFX & ThreadX - 硬汉嵌入式论坛 - Powered by Discuz! 更新一期视频教程: 第6期ThreadX视频教程:图文并茂吃透RTOS运行机制,任务管理&…

微信小程序的全局弹窗以及全局实例

全局组件 微信小程序组件关系中,父组件使用子组件需要在父组件index.json中引入子组件,然后在父组件页面中使用,这种组件的对应状态是一对一的,一个组件对应一个页面。如果有一个全局弹窗(登录)&#xff0…

拒绝select *,Navicat快捷查询表中所有字段以,隔开两种方式【附带阿里开发手册】

一、前言 今天公司规定了不要使用SELECT *进行查询,让我想到阿里的《Java 开发手册》 中的ORM 映射规范也是这样的,于是翻出来看看,刚刚好重温一下! 规范看了一定要实践,严格规范自己的代码风格,做一个优…

TEMU联盟计划用意何在?做推广达人真的能收入满满吗?

据东哥近期了解,Temu在北美市场表现十分火爆,甚至冲上了AppStore下载榜第一名。Temu在美国市场上采用了类似PDD的病毒式传播策略,以实惠的产品和折扣吸引消费者并动员普通人大量传播链接和App下载,所以有了TEMU联盟计划&#xff0…

Sprng依赖注入(三):构造方法注入是如何工作的?

前言这是Spring依赖注入系列的第三篇,前两篇主要分析了Spring bean依赖属性注入的两种方式,是字段注入和setter方法注入,单独比较这两种方式,会发现其过程和工作原理非常类似,那么构造方法注入会不会也和前两种比较类似…

搭建个人博客保姆级教程(二)

文章目录一. Springboot项目打包成jar包二. mySql部署三. UI 项目部署一. Springboot项目打包成jar包 使用IDEA进行打包,当然也有其他部署方式,如使用maven进行打包,可自行查询资料。 1.打开项目,右击项目选中Open Module Settin…

【C++的OpenCV】第五课-OpenCV图像常用操作(二):OpenCV的基本绘图、平滑滤波(模糊)处理

让我们继续一、OpenCV基本绘图1.1 OpenCV关于绘图的操作1.1.1 cv::Point()1.1.2 cv::Scalar()1.1.3 cv::line()画线1.1.4 cv::rectangle()画矩形1.1.5 cv::circle()画圆二、图像的平滑滤波处理2.1 概念2.2 OpenCV关于图像模糊的操作2.2.1 常用滤波器的分类2.2.2 各种滤波方法具…

【4】linux命令每日分享——cd切换路径

大家好,这里是sdust-vrlab,Linux是一种免费使用和自由传播的 类UNIX操作系统,Linux的基本思想有两点:一切都是文件;每个文件都有确定的用途;linux涉及到IT行业的方方面面,在我们日常的学习中&am…

云计算学习课程——越来越重要的云安全

2023,越来越多的企业和组织正在或即将把核心系统和数据迁移上云端,其中以公有云和服务居多,那么就意味着在数据迁移的过程中会出现安全问题的几率更大。企业也越来越注重云安全体系,对我们云计算运维工程师来说,也是一…

Linux 日志讲解

目录 日志系统的分类 常见日志文件 日志等级 journalctl查看指定日志信息 日志系统的分类 系统日志 记录系统的运行情况和内核信息 用户日志 记录用户的访问信息 程序日志 记录程序运行的各种时间(每个服务程序都有自己独立的日志文件&#x…

UE4 手把手教你做插件(4)做一个可以拖入场景的小方块插件

0,前言 学习一下如何编写一个有actor的插件,能够将它拖入场景里面。 1,创建插件 (1)创建一个空白插件 添加之后vs里面就会多出一个目录: (2)给插件添加一个actor 创建完actor之后&am…

渗透测试之交换式网络嗅探实验

渗透测试之交换式网络嗅探实验实验目的一、实验原理1.1 网络嗅探器Sniffer的工作原理1.2 网络嗅探器的分类1.3 网络嗅探器Sniffer的作用二、实验环境2.1 操作机器2.2 实验工具Sniffer2.3 安装工具Sniffer三、实验步骤1. 熟悉Sniffer工具的启动2. 进行监听3. 熟悉Sniffer工具的介…

Python数据结构与算法篇(四)-- 链表的实现

实现线性表的另一种常用方式就是基于链接结构,用链接关系显式表示元素之间的顺序关联。基于链接技术实现的线性表称为链接表或者链表。 采用链接方式实现线性表的基本想法如下: 把表中的元素分别存储在一批独立的存储块(称为表的结点)里。保…

刷题记录:牛客NC51112Stars in Your Window 扫描线

传送门:牛客 题目描述: Here comes the problem: Assume the sky is a flat plane. All the stars lie on it with a location (x, y). for each star, there is a grade ranging from 1 to 100, representing its brightness, where 100 is the brightest and 1 is the we…

数据库大量数据导出导入的操作

基于上一篇文章,我想到如果在数据库之间拥有大量数据的表数据的导入导出,该如何快速完成表的导入导出工作呢? 思路一:使用db软件工具导出数据,然后向新数据库导入数据。我用的是dbserver22.3.4,不出意外的是&#xff0…

【读书笔记】《深入浅出数据分析》第二章 检验你的理论

文章目录一,相关分析方法1,相关系数二,相关性不等于因果关系三,证明因果关系,“控制变量法”?本章主要说明了两个问题: 1,相关性不等于因果关系 2,如何判断两种数据之间是相关性&am…

深圳/东莞/惠州师资比较强的CPDA数据分析认证

深圳/东莞/惠州师资比较强的CPDA数据分析认证培训机构 CPDA数据分析师认证是中国大数据领域有一定权威度的中高端人才认证,它不仅是中国较早大数据专业技术人才认证、更是中国大数据时代先行者,具有广泛的社会认知度和权威性。 无论是地方政府引进人才、…

Android 高性能列表:RecyclerView + DiffUtil

文章目录背景介绍一般刷新 notifyDataSetChanged()局部刷新实现调用代码准备工作创建 MyDiffUtilCallback 类继承 DiffUtil.Callback 抽象类MyAdpter 类代码实现步骤总结通过 log 证实 diffutil 的局部刷新diffutil 优化后台线程参考主线程参考diff 更新优化后写法相关参考背景…