爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: Aires

[源程序] 超前滞后相关+滑动平均

  [复制链接]

新浪微博达人勋

发表于 2016-8-12 09:55:11 | 显示全部楼层
{:eb502:}谢谢分享
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2016-8-15 10:09:08 | 显示全部楼层
楼主,计算相关系数的时候感觉有问题,你先对整个X与Y进行标准化,后面计算的时候直接取他们标准化后一部分的数据进行计算,这样感觉不对,应该是取X,Y相应超前滞后的部分,然后标准化再计算,之前你的方法,算出来的相关系数会出现小于-1的情况,可能问题就是出在这里。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-8-15 18:19:49 | 显示全部楼层
小小书 发表于 2016-8-15 10:09
楼主,计算相关系数的时候感觉有问题,你先对整个X与Y进行标准化,后面计算的时候直接取他们标准化后一部分 ...

哦,这样啊,我现在用matlab自带的xcorr函数,你可以查查怎么用
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-8-16 09:42:27 | 显示全部楼层
Aires 发表于 2016-8-15 18:19
哦,这样啊,我现在用matlab自带的xcorr函数,你可以查查怎么用

我把程序改了下,也用xcorr做过,在样本小的情况下有偏差,样本大的情况下很接近,前辈你是用xcorr(x,y,t,'coeff')这样的语句做的?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-8-16 09:42:30 | 显示全部楼层
Aires 发表于 2016-8-15 18:19
哦,这样啊,我现在用matlab自带的xcorr函数,你可以查查怎么用

我把程序改了下,也用xcorr做过,在样本小的情况下有偏差,样本大的情况下很接近,前辈你是用xcorr(x,y,t,'coeff')这样的语句做的?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-8-16 10:21:16 | 显示全部楼层
小小书 发表于 2016-8-16 09:42
我把程序改了下,也用xcorr做过,在样本小的情况下有偏差,样本大的情况下很接近,前辈你是用xcorr(x,y,t ...

嗯,是的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-8-16 13:42:31 | 显示全部楼层
matlab这个语句做的和根据公式算的有偏差,也不知道为什么。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-8-16 14:28:00 | 显示全部楼层
小小书 发表于 2016-8-16 13:42
matlab这个语句做的和根据公式算的有偏差,也不知道为什么。

它超前滞后的算法不大一样吧,你测试的结果偏差很大么?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-8-16 14:50:55 | 显示全部楼层
Aires 发表于 2016-8-16 14:28
它超前滞后的算法不大一样吧,你测试的结果偏差很大么?

我用的是月数据,当用30年月数据,有微小偏差,当做50年,全部重合,但是当小于20年,偏差就比较明显了,小于10年偏差更大了,通常偏差都是出现在画出图形的波峰处,波形都是一样的,xcorr算的最大值通常是小于公式算的,样本小的时候可达0.1的偏差。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-8-16 16:57:42 | 显示全部楼层
小小书 发表于 2016-8-16 14:50
我用的是月数据,当用30年月数据,有微小偏差,当做50年,全部重合,但是当小于20年,偏差就比较明显了, ...

哦,那应该说明这两种方法算的应该都是可信的。有一个可能性,就是算超前滞后相关时取的数据长度是否一致,比如数据一共5个数,一种是同期相关时是用整个数据(5个数),然后逐渐错开做超前滞后;另一种是比如要求最大滞后2个数,所以采用每次超前滞后的数据长度是一样的,即1~3与1~3同期,1~3与2~4错开一个,1~3与2~5错开两个。不知道你用的是哪种算法算的?我觉得是可能样本量比较大时,即便选取的数据段不一致,自由度还是比较大的,相关系数比较稳定,你test一下。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表