数学建模强化宝典(13)M-K检验法

前言

       M-K检验法,全称为Mann-Kendall检验法,是一种非参数的假设检验方法,广泛应用于时间序列数据的趋势性变化检验,特别是气候序列中的趋势分析和突变点检测。以下是对M-K检验法的详细介绍:

一、定义与背景

       M-K检验法由Frank Mann和Donald R. Kendall提出并发展,因此得名Mann-Kendall检验法。该方法不依赖于数据的具体分布形式,能够处理含有异常值和非正态分布的数据,非常适合于气候、水文、环境等领域的时间序列数据分析。

二、基本原理

       M-K检验法的基本原理是通过比较每个数据点与其之前数据点的大小,来检测时间序列数据中的单调趋势(上升、下降或没有趋势)。具体步骤如下:

  1. 构造秩序列:对于时间序列中的每一个数据点,计算其之前所有数据中比它大或小的数据点个数,从而构造出一个秩序列。
  2. 计算统计量:基于秩序列,计算统计量UFk(正向序列)和UBk(反向序列),这两个统计量用于评估时间序列的趋势性和突变点。
  3. 显著性检验:将UFk和UBk与设定的显著性水平(如0.05)对应的临界值进行比较,判断趋势是否显著,以及是否存在突变点。

三、应用领域

M-K检验法广泛应用于以下领域:

  1. 气候变化研究:用于检测气候序列中的趋势变化和突变点,如气温、降水等气候要素的变化趋势。
  2. 水文水资源研究:分析河流流量、地下水位等水文要素的时间序列数据,评估其变化趋势和突变情况。
  3. 环境科学研究:在环境监测和评估中,利用M-K检验法分析污染物浓度、环境质量等指标的变化趋势。

四、优点与局限性

优点:
  1. 非参数检验:不依赖于数据的具体分布形式,适用于各种类型的时间序列数据。
  2. 抗干扰能力强:对异常值不敏感,能够较好地处理非正态分布的数据。
  3. 计算简便:计算过程相对简单,易于实现自动化处理。
局限性:
  1. 无法确定趋势形式:M-K检验法只能判断时间序列是否存在单调趋势,但无法确定趋势的具体形式(如线性、非线性等)。
  2. 对周期性变化不敏感:对于存在周期性变化的时间序列数据,M-K检验法的检测效果可能不佳。

五、实现 

       下面是一个简单的MATLAB函数,用于执行Mann-Kendall趋势检验,并计算统计量UFk(正向序列)以及可选的UBk(反向序列,用于确定突变点): 

function [UF, UB, pValue] = mannKendall(x)  % 输入:  % x - 时间序列数据向量  %  % 输出:  % UF - 正向序列的Mann-Kendall统计量  % UB - 反向序列的Mann-Kendall统计量(可选)  % pValue - 对应的p值(此函数不直接计算,通常需要查表或使用其他方法)  n = length(x);  s = 0;  varS = 0;  UF = zeros(1, n);  % 初始化  for k = 2:n  for j = 1:k-1  if x(k) > x(j)  s = s + 1;  elseif x(k) < x(j)  s = s - 1;  end  end  UF(k) = (s - k * (k - 1) / 2) / sqrt(k * (k - 1) * (2 * k + 5) / 18);  varS = varS + s^2;  s = 0;  end  % UB的计算需要反向序列  UB = zeros(1, n);  if nargout > 2  % 计算方差  varS = (n * (n - 1) * (2 * n + 5) / 18) - varS;  if varS == 0  varS = 1; % 避免除以零  end  % 反向序列  y = fliplr(x);  for k = 2:n  UB(k) = -UF(n+1-k); % 利用正向序列的结果  % 注意:这里直接利用UF的结果进行转换,实际上UB的计算需要完整重复上述过程,  % 但对于标准Mann-Kendall测试,通常只关注UF,UB用于确定突变点  end  % 标准化UB  UB = UB / sqrt(varS / 18);  end  % 注意:此函数不直接计算p值,通常需要查表或使用统计软件  pValue = NaN; % 占位符  
end

注意

  1. 上述代码主要计算了UF统计量,并可选地计算了UB统计量(虽然UB的计算是基于UF的简化,实际中可能需要独立计算)。
  2. p值的计算通常不是直接通过M-K统计量本身进行的,而是需要参考相应的统计表或使用专门的统计软件/函数。
  3. 对于趋势的显著性判断,通常需要将UF或UB的值与一定的显著性水平(如α=0.05对应的临界值)进行比较。
  4. 在实际应用中,你可能还需要对时间序列进行预处理,如去趋势、季节性调整等,以提高M-K检验的准确性。

六、结论

       M-K检验法作为一种非参数的假设检验方法,在气候、水文、环境等领域的时间序列数据分析中发挥着重要作用。其优点在于不依赖于数据的具体分布形式且抗干扰能力强,但也存在无法确定趋势形式和对周期性变化不敏感的局限性。在实际应用中,应根据具体研究目的和数据特点选择合适的分析方法。

 结语 

成功是优点的发挥

失败是缺点的累积

!!!

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

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

相关文章

嵌入式初学-C语言-数据结构--四

栈 1. 基本概念 栈是一种逻辑结构&#xff0c;是特殊的线性表。特殊在&#xff1a; 只能在固定的一端操作 只要满足上述条件&#xff0c;那么这种特殊的线性表就会呈现一种“后进先出”的逻辑&#xff0c;这种逻辑就被称为栈。栈 在生活中到处可见&#xff0c;比如堆叠的盘子…

Kubernetes部署(haproxy+keepalived)高可用环境和办公网络打通

HAProxy Keepalived 部署高可用性入口&#xff1a; 部署两台或多台节点运行 HAProxy 作为负载均衡器。使用 Keepalived 实现 VIP&#xff08;虚拟 IP&#xff09;&#xff0c;为 HAProxy 提供高可用性。Keepalived 会监控 HAProxy 的状态&#xff0c;如果主节点失效&#xff0…

Spring AOP,通知使用,spring事务管理,spring_web搭建

spring AOP AOP概述 AOP面向切面编程是对面向对象编程的延续&#xff08;AOP &#xff08;Aspect Orient Programming&#xff09;,直译过来就是 面向切面编程,AOP 是一种编程思想&#xff0c;是面向对象编程&#xff08;OOP&#xff09;的一种补充。&#xff09; 面向切面编…

【C++二分查找】1482. 制作 m 束花所需的最少天数

本文涉及的基础知识点 C二分查找 LeetCode1482. 制作 m 束花所需的最少天数 给你一个整数数组 bloomDay&#xff0c;以及两个整数 m 和 k 。 现需要制作 m 束花。制作花束时&#xff0c;需要使用花园中 相邻的 k 朵花 。 花园中有 n 朵花&#xff0c;第 i 朵花会在 bloomDay…

南京网站设计手机用的网站

近年来&#xff0c;随着移动互联网的快速发展&#xff0c;越来越多的用户通过手机浏览网页&#xff0c;这使得网站设计逐渐向移动端倾斜。在南京&#xff0c;网站设计特别注重适配手机端&#xff0c;这不仅是用户体验的提升&#xff0c;也是市场竞争的需要。一个响应式的网站能…

vue3 v-bind=“$attrs“ 的一些理解,透传 Attributes相关说明及事例说明

1、可能小伙伴们经常会在自己的项目中看到v-bind"$attrs"&#xff0c;这个一般是在自定义组件中看到。 比如&#xff1a; <template><BasicModalv-bind"$attrs"register"registerModal":title"getTitle"ok"handleSubm…

Python进程间网络远程通讯方式:socket、pipe、RPC详解!

背景 最近在进行开发工作的时候&#xff0c;遇到了一个场景&#xff1a; pc程序需要和安卓设备进行通讯和接口调用。 此时就需要进行远程调用方法。然而大学时代有关于远程过程调用的知识都还给了老师……所以在此进行一个复习&#xff0c;并进行实战演练&#xff01; 网络…

UE驻网失败问题(三)

这个问题是lab问题&#xff0c;现象如下&#xff1a; 期望UE注册在SA网络下&#xff0c;咋一看没有5G MIB/SIB1打印&#xff0c;好像是没搜到5G小区&#xff0c;而实际上并不是这样。 在查看搜网过程时会发现如下log打印&#xff1a; [I nr5g_rrc_acq.c 3544] RRC ACQ: Band 41…

一个全面、先进、公平且模块化的开源RAG框架

两个关键问题限制了 RAG 的发展&#xff1a; 新型 RAG 算法之间缺乏全面和公平的比较。 像 LlamaIndex 和 LangChain 这样的开源工具使用了高级抽象&#xff0c;这导致了透明度的缺失&#xff0c;并限制了开****发新算法和评估指标的能力。 RAGLAB&#xff1a;是一个模块化的…

QML学习三:qml设计器报错 Line: 0: The Design Mode requires a valid Qt kit

开发环境:Qt 6.5.3 LTS 1、Qt 6.5.3 LTS 2、Pyside6 3、Python 3.11.4 4、win11 默认不打开设计器的时候可以看到我们默认是有Python的环境,而且点击运行是可以运行的。但是当打开qml设计器时提示下面这个错误,提示需要一个可用的套件。 解决办法: 点击+号创建新的套件…

OpenAI官宣ChatGPT用户破2亿,新模型比GPT-4强100倍!

在刚刚召开的“KDDI SUMMIT 2024”大会上&#xff0c;OpenAI日本首席执行官長崎忠雄表示&#xff0c;截止至今年8月底&#xff0c;ChatGPT的月活用户数量突破2亿大关&#xff0c;成为史上最快达到这一成就的软件产品。 同时&#xff0c;ChatGPT企业级用户也获得大幅度增长达到…

Windows自动化程序开发指南

自动化程序的概念 “自动化程序”指的是通过电脑编程来代替人类手工操作的一类程序或软件。这类程序具有智能性高、应用范围广的优点&#xff0c;但是自动化程序的开发难度大、所用技术杂。 本文对自动化程序开发的各个方面进行讲解。 常见的处理对象 自动化程序要处理的对…

STM32关于keil使用过程中遇到的问题

1.设备管理器STlink驱动确认安装完成&#xff0c;但是keil里一直识别不到&#xff0c;换下载器也没用 &#xff08;1&#xff09;问题描述 我的问题是这样产生的&#xff1a;之前用标准库开发STM32的时候&#xff0c;STLink能够正常使用&#xff0c;然后使用HAL库开发的时候出…

声音克隆工具CosyVoice

阿里的免费声音克隆工具CosyVoice CosyVoice 是阿里通义实验室在七月初开源的一款专注于自然语音合成的语音大模型&#xff0c;它具备多语言、多音色和细腻的情感控制能力。这个系统支持中文、英文、日文、粤语和韩语五种语言的语音生成&#xff0c;并且在语音合成的效果上远超…

WebShell流量特征检测_哥斯拉篇

90后用菜刀&#xff0c;95后用蚁剑&#xff0c;00后用冰蝎和哥斯拉&#xff0c;以phpshell连接为例&#xff0c;本文主要是对后三款经典的webshell管理工具进行流量分析和检测。 什么是一句话木马&#xff1f; 1、定义 顾名思义就是执行恶意指令的木马&#xff0c;通过技术手…

STM32G474之DAC

STM32G474分别使用CORDIC硬件和“math.h”的正弦值&#xff0c;从DAC1和DAC2输出。 1、DAC特点 PA4的附加功能为DAC1_OUT1&#xff0c;无需映射&#xff0c;直接将它配置为模拟功能&#xff0c;就可以使用了。 PA6的附加功能为DAC2_OUT1&#xff0c;无需映射&#xff0c;直接将…

物料主数据的分类及编码规则概述

物料主数据分布在设计、工艺、生产、采购、销售、存储、物流过程中的多个业务部门中&#xff0c;类型繁杂且不易于管理&#xff0c;必须对于物料主数据进行分类管理。 本文将从物料主数据的分类原则、编码样例、编码规范以及编码设计等方面进行详细介绍&#xff0c;帮助企业更…

STM32高级定时器生成互补PWM的原理与代码实现

文章目录 前言一 CubeMx配置1.1 TIM1 Mode and Configuration1.2 Paramter Settings 二 程序代码三 仿真分析总结 前言 互补 PWM&#xff08;Complementary PWM&#xff09;是指一对逻辑状态互为反相的 PWM&#xff08;脉冲宽度调制&#xff09;信号。这种信号配置常见于电机控…

数字人实战第三天——wav2lip部署教程

wav2lip数字人 一、wav2lip简介 Wav2Lip 是一种通过将音频与视频中的嘴唇动作同步的技术&#xff0c;旨在生成与音频内容高度匹配的口型动画。其主要应用是让视频中的人物嘴唇动作与配音或其他音频输入精确同步&#xff0c;这在电影配音、虚拟主持人、在线教学、影视后期处理等…

秋招最新大模型算法面试,熬夜都要肝完它

&#x1f4a5;大家在面试大模型LLM这个板块的时候&#xff0c;不知道面试完会不会复盘、总结&#xff0c;做笔记的习惯&#xff0c;这份大模型算法岗面试八股笔记也帮助不少人拿到过offer ✨对于面试大模型算法工程师会有一定的帮助&#xff0c;都附有完整答案&#xff0c;熬夜…