滑动相关

吴东佑

最近，用到了滑动相关分析，学习了论坛里各位前辈发的帖子。最后，参考ncl的滑动相关脚本，改写了matlab的计算函数，以方便使用。
所谓滑动相关，其实就是在给定滑动窗口的前提下，计算两段特定子序列的相关系数，并且这两段子序列会随之发生变化，但是其窗口长度是不变的，因此便得到了滑动相关系数。主要代码如下，具体参数说明及使用可以见附件！
程序如有错误，还望大家可以指出，谢谢！

1. function [cor,date]=run_corr(x,y,time,wSize)
   
2.     N=length(x);
   
3.     if N~=length(y) || N~=length(time)
   
4.         print('run_cor: input array are different sizes: fatal');
   
5.         return
   
6.     end
   
7.    
   
8.     if isa(x,'double') ||  isa(y,'double')
   
9.         cor=ones((N-(wSize-1)),1).*nan;
   
10.         cor=double(cor);
    
11.     else
    
12.         cor=ones((N-(wSize-1)),1).*nan;
    
13.         cor=float(cor);
    
14.     end
    
15.    
    
16.     date=ones((N-(wSize-1)),1).*nan;
    
17.    
    
18.     k=(wSize-1)/2;
    
19.     for j=1:N-wSize+1
    
20.         date(j,1)=time(k)+j;
    
21.     end
    
22.    
    
23.     for i=1:N-wSize+1
    
24.         aa=x(i:i+(wSize-1));
    
25.         bb=y(i:i+(wSize-1));
    
26.         cor(i,1)=corr(aa,bb,'type','Pearson');
    
27.     end
    
28. end

复制代码

kingjames

滑动相关系数  好好学习  谢谢

supermann

楼主你好，请问您参考的是哪里的NCL脚本呢

吴东佑

supermann 发表于 2018-11-27 09:47
楼主你好，请问您参考的是哪里的NCL脚本呢

安装文件下面的脚本

supermann

吴东佑 发表于 2018-11-27 15:02
安装文件下面的脚本

感谢您的回复

常山赵子龙

为什么程序会报错？两个同时间序列？这是什么情况？？？

qsl950915

求ncl脚本哈哈哈

风叶ele

请问NCL中计算滑动相关是哪个函数？ 有这个内置函数吗？

伍孚ms

风叶ele 发表于 2019-8-5 11:24
请问NCL中计算滑动相关是哪个函数？ 有这个内置函数吗？

http://www.ncl.ucar.edu/Document ... buted/run_cor.shtml

微醉的阳光

楼主求教下滑动平均和互相关的区别xcorr？