分享一个气候突变点检验的matlab程序

desertsky

猫的他 发表于 2018-8-19 07:41
下载JMA的台风数据，然后取OND季节的，文章的资料和方法里都有
数据的网址
http://www.jma.go.jp/jma/j ...

谢谢，已经下载到，但PANG-CHI HSU论文“An Abrupt Decrease in the Late-Season Typhoon Activity over the
Western North Pacific” 中所用台风标准好像和JMA的标准不一样，所以OND的数据对不起来呢 :(
另，请问纵坐标为PDF的那张图  您用附件程序能生成吗？  我的没有呢

desertsky

猫的他 发表于 2018-8-19 07:41
下载JMA的台风数据，然后取OND季节的，文章的资料和方法里都有
数据的网址
http://www.jma.go.jp/jma/j ...

谢谢，已经下载到，但PANG-CHI HSU论文“An Abrupt Decrease in the Late-Season Typhoon Activity over the
Western North Pacific” 中所用台风标准好像和JMA的标准不一样，所以OND的数据对不起来呢 :(
另，请问纵坐标为PDF的那张图  您用附件程序能生成吗？  我的没有呢

猫的他

desertsky 发表于 2018-8-19 12:31
谢谢，已经下载到，但PANG-CHI HSU论文“An Abrupt Decrease in the Late-Season Typhoon Activity over  ...

不一致，只要能有年代际的突变也行，
那个PDF的图是可以生成的，你在看看吧

猫的他

为了方便测试，我把数据也传上来吧
我大概介绍一下这个的用法（其实很简单，就是一个matlab函数的调用）
第一步，先在matlab中设置好对应的这个外部函数的路径（不会的自行百度）
第二步，导入数据 (在命令行里面输入，count = load('JMA_OND_TY_COUNT.txt');第三步，调用这个函数(在命令行里输入，prob = MultipleChangePointDetection(count);)
当函数执行完毕，就会自己生成贝叶斯断点的后验概率图。

desertsky

猫的他 发表于 2018-8-23 09:40
为了方便测试，我把数据也传上来吧
我大概介绍一下这个的用法（其实很简单，就是一个matlab函数的调用）
...

首先谢谢lz，按照您的流程和数据，还是没有完全出来您所显示的结果，一张图一致，但另一张折线图没有出来，而且纵坐标为PDF的图没有 :(

猫的他

desertsky 发表于 2018-8-23 12:55
首先谢谢lz，按照您的流程和数据，还是没有完全出来您所显示的结果，一张图一致，但另一张折线图没有出来 ...

折线图的确不会自己出，折线图是我为了验证数据的正确性另外用ncl画的；
我matlab出的图也和你的一样是，那个PMF在统计里的全称是probability mass function
如果仔细看Hsu Pang Chi 老师论文里的图注的话，可以发现她画的也是posterior probability mass function，只不过她可能后期改了y轴的标题。
所以你出的图是对的！

紫藤萝

desertsky 发表于 2018-8-23 12:55
首先谢谢lz，按照您的流程和数据，还是没有完全出来您所显示的结果，一张图一致，但另一张折线图没有出来 ...

请问你画出图的横坐标怎么不是时间，只有6个数字？？？

紫藤萝

猫的他 发表于 2018-8-23 16:15
折线图的确不会自己出，折线图是我为了验证数据的正确性另外用ncl画的；
我matlab出的图也和你的一样是 ...

楼主 你好 。请问我们画出图的横坐标是“number of change point”，而不是时间，代表什么意思啊？
为什么和你的横坐标不一样呢？谢谢了

紫藤萝

我的年份是从1960-2013年 54年，横坐标怎么设置为1960-2013啊，即显示出54年每年出现突变的概率