两个序列的互相关计算(基于R语言实现)

基于R语言实现两个时间序列c1和c2的互相关计算，假设c1和c2都是一维数据，它们分别保存在c1.txt和c2.txt文件中，在R界面终端需要用到read.table()函数将数据导入，并用ccf()函数计算互相关，具体代码如下：

>R
>c1 <- read.table('c1.txt')
>c2 <- read.table('c2.txt')
>y13 <- ccf(c1,c2) #结果见下图y13 <- ccf(c1,c2,lag.max=x),lag_max=x表示最大延迟x

#将计算结果保存到本地，需要把ccf()函数返回的数据结构适当转化再保存
>ccf_values <- y13$acf
>lag_values <- y13$lag
>ccf_df <- data.frame(Lag=lag_values,CCF=ccf_values)
>write.table(ccf_df,file='ccf.txt',sep='\t',row.names=FALSE)