mk突变点检测_MK突变检测程序
MannKendall
突变检测算法
MATLAB
源码
题目:
MannKendall
突变检测算法
MATLAB
源码
%%
以下是单边
MannKendall
突变检测算法
% GreenSim
团队原创作品,转载请注明
% Email:greensim@163.com
% GreenSim
团队主页:
http://blog.sina.com.cn/greensim
% [color=red]
欢迎访问
GreenSim
——算法仿真团队
→[url=http://blog.sina.com.cn/greensim]http://blog.sina.com.cn/green
sim[/url][/color]
clear
load data
X=MJL;
N=length(X);
U=zeros(N-1,1);
for t=2:N
x=X(1:t);
S=0;
n=length(x);
for k=1:(n-1)
for j=(k+1):n
S=S+sign(x(j)-x(k));
end
end
VarS=n*(n-1)*(2*n+5)/18;
if S>0
Z=(S+1)/sqrt(VarS);
elseif S==0
Z=0;
else
Z=(S-1)/sqrt(VarS);
end
U(t-1)=Z;
end
figure(1)
plot(1:(N-1),U,'linewidth',1.5);
hold on
plot(1:(N-1),1.96*ones(N-1,1),':','linewidth',1);
legend('
统计量
','0.05
显著水平
');
hold on
plot(1:(N-1),0*ones(N-1,1),'-.','linewidth',1);
plot(1:(N-1),1.96*ones(N-1,1),':','linewidth',1);
mk突变点检测_MK突变检测程序相关推荐
- mk突变点检测_mk突变检测
寄信人: fig ([晋]『康杰中学』) 标 题: 没主题 发信站: 南京大学小百合站 (Sun May 23 19:18:25 2010) 来 源: 172.16.74.204 Fortran: p ...
- mk突变点检测_Mann-Kendall突变检测(mk突变检测)
本帖最后由 vb1987 于 2013-6-12 23:27 编辑 %最近写论文需要用到MK检验法,网上收集到大量的matlab代码,但是没有一个代码能够 %完全正确运行或者分析信息不全,结合多位网友 ...
- mk突变点检测_科学网—从网上找的M-K突变检验的程序 - 张乐乐的博文
%从matlab论坛上找的MK突变检验的程序,这个程序运行的结果跟我自己编写程序运行出来的结果一样,但是跟魏凤英老师书上的例子出图结果不一样 A=xlsread('test-mk.xlsx'); x= ...
- mk突变点检测_MK检验突变分析 matlab
% Mann-Kendall突变检测 % 数据序列y % 结果序列UFk,UBk2 %读取excel中的数据,赋给矩阵y %获取y的样本数 %A为时间和降水数据列 x=降水(:,1);%时间序列 y= ...
- mk突变点检测_气候突变检测
1.气候突变是什么意思 科学家们所说的气候变化,通常指的是"气候渐变".也就是我们经常说的,"地球在逐渐变暖",或者是"气候在慢慢变暖".但 ...
- python:遥感时间序列处理——Pettitt突变点检测(突变时间/年份)
本文记录了使用python语言和Pettitt突变检测算法检测遥感数据30年的NDVI时间序列中的突变时间(年份)的代码. Pettitt可以检测时间序列中最显著的突变点,其检测结果只有一个,相对MK ...
- python:遥感时间序列处理——Mann-Kendall(MK)突变检测(突变次数)
作者:CSDN @ _养乐多_ 本文记录了使用python语言读取数据和逐像素求MK突变次数的方法和代码. 并以2001年到2022年的NDVI时间序列数据为例,展示了研究区内NDVI发生突变的次数. ...
- pettitt检验r语言_[原创]Matlab气候突变检测程序合集(更新了)
登录后查看更多精彩内容~ 您需要 登录 才可以下载或查看,没有帐号?立即注册 x 本帖最后由 sfhstcn2 于 2017-10-19 07:29 编辑 mutation_test这是一包代码,是基 ...
- 【突变检验方法二】MATLAB实现贝叶斯突变检测
MATLAB实现贝叶斯突变检测 1 贝叶斯突变检测 2 原理 3 MATLAB相关代码 3.1 调用函数 3.2 案例 参考 另:其它语言实现贝叶斯突变检测 1 贝叶斯突变检测 贝叶斯突变检测属于概率 ...
最新文章
- 奇葩错误:不同变量名称相同
- 使用Maven插件构建SpringBoot项目,生成Docker镜像push到DockerHub上
- kickstart中ks.cfg指定目标机ip的小备忘
- ​​​​​​​2016最新CocoaPods安装与使用
- 【云计算虚拟化】docker在ubuntu14.04下的安装
- js 常用的正则表达式
- Java学习小程序(4)数列求和
- 3.js模式-策略模式
- Delphi使用server组件控制excel的一点经验
- golang切片转为字符串,字符串转切片
- ORA-03113: end-of-file on communication channel
- r 64位连接mysql_R与Mysql数据库连接的两种方法
- win7(32位)U盘安装、卸载ubuntu(64位)双系统
- 时域、频域与时域分析、频域分析简单理解
- 云尚发卡PHP,PHP云尚发卡系统V1.5.7源码,专门为个人或小型企业提供在线售卡
- ajax 网页加速,20 种提升网页速度的技巧
- Rhino学习教程——1.2
- 摇一摇 周边 与 ibeacon
- 杂七杂八(9): IDEA初始化配置 插件收集
- Javescript第二周学习