[原创]Matlab气候突变检测程序合集（更新了）

why110sc

为楼主点个赞！很厉害

初学的小时

谢谢楼主的慷慨分享

ainibaichi

谢谢楼主的慷慨分享

zhoubaojia

楼主您好，小白提问。我在用yamamoto检测突变的时候，取N=10，P=0.05,画出的图中显示P=0.05的线大约为0.66，并不等于1，但很多文献都表信噪比大于1.0才显著啊。这是为何呢？

sfhstcn2

zhoubaojia 发表于 2019-1-15 12:38
楼主您好，小白提问。我在用yamamoto检测突变的时候，取N=10，P=0.05,画出的图中显示P=0.05的线大约为0.66 ...

书上是说，IH=10，Rsn=1时，通过了0.01显著水平的t检验，但这并不代表0.01显著水平下的t值对应的Rsn=1。我这儿是基于p值反推的Rsn阈值。

taeyeon

怎么把小波分析的实部和小波方差值导出来呢？

sfhstcn2

taeyeon 发表于 2019-1-28 16:31
怎么把小波分析的实部和小波方差值导出来呢？

这个直接调用wavelet那几个函数，都是在那个网站上下载的（帖子上有链接），可以参考上面学习。

twok11

厉害。。。。。。

漠小北

前辈加楼主：

         您好！
          我学习了您在气象家园上传的气候突变检测程序代码，您相关帖子的主题是‘Matlab气候突变检测程序合集（更新了）’。关于突变点的检测有点疑惑，想请教一下您。如果方便，还得麻烦您有空的时候帮我分析一下具体错误的原因在哪里，希望没有打扰到您，谢谢前辈了。
         我主要的问题在于pettitt检验和滑动t检验，先看pettitt检验的结果吧，使用的是您上传的那个代码。原始数据序列的突变点通过肉眼可以大致看到在1998年左右，而且通过pettitt检验之后，确实检测出1998年是突变点，但是问题在于最后统计量p值的计算，我查了很多文献，资料，书籍，都说最后计算出来的p值小于0.05就代表检测出来的突变点具有统计意义，但是我计算出来的p值等于1.9161。不知道到底问题出在了哪里，还得烦请您帮我看看，不胜感激。
         我的学习数据具体情况是1982年到2015年的一组长度为34的数据，具体数据和计算结果在附件excel里面，所使用的程序和matlab计算的结果也截图粘贴在excel表格里面。
         关于滑动t检验的计算结果也粘贴在excel中，我个人感觉滑动t检验计算结果里面，和0值线相交的那个点才是突变点，但是书籍上面说是和上线两条临界线相交的地方才是突变点的位置。但是很明显，原始数据序列不可能存在两个突变点，所以感觉也有点问题，我查了很久的资料，还是想不明白，还得麻烦前辈有空帮我看看这到底该如何解释呢？
         非常感谢前辈阅读这封邮件，祝您在新年‘猪’事顺利，一帆风顺，阖家欢乐！

祝好，
小北同学。

漠小北

附上数据结果，计算结果和所使用的代码截图，希望前辈和园子里面的各位一起帮忙分析看看，不胜感激各位的帮忙。