| 查看: 630 | 回复: 0 | ||
[求助]
【请教】延时相关分析问题
|
|
各位虫友，大家好！ 最近想做两个变量的延时相关分析，请人编了个matlab程序代码，但是感觉结果有些问题，比如，在计算延迟0-12个月的径流与降水量二者的相关系数时，延迟0个月（也就是两个变量同期，不延迟），按理论，同期的径流和降水量的相关系数应该最高，可0个月的延迟相关系数只有0.38，我怀疑是程序的算法有些问题，请高手帮我看看，这个程序到底问题出在哪里，不胜感激！ 附程序代码和数据如下 close all; clear all; clc; fid1=fopen('G:\科研工具\Spearman互相关系数\data.txt','r'); data=fscanf(fid1,'%d %f %d %f',[4 inf]); data=data'; fclose(fid1); N=length(data(:,2)); tongji=zeros(13,1); x=zeros(12,1); y=zeros(12,1); indx=0; for i=1:12:N X=data(i:i+11,2); Y=data(i:i+11,4); x=x+X; y=y+Y; indx=indx+1; end x=x/indx; y=y/indx; x=x-mean(x); y=y-mean(y); [c,LAGS]=xcorr(x,y,'coeff'); tongji=tongji+c(11:end); m=0:12; plot(m,c(11:end),'-bo'); hold on; save tongji.mat tongji; fuzhi1=ones(13,1)*0.01; fuzhi2=ones(13,1)*0.05; plot(m,fuzhi1,'--r'); hold on; plot(m,fuzhi2,'-k'); hold on; set(gca,'xtick',0:1:12); xlabel('滞后时间/月'); ylabel('相关系数'); axis([0,12,-1,1]); title('spearman相关系数'); legend('r0.01临界值','r0.05临界值'); A=[m',tongji,fuzhi1,fuzhi2]; fid2=fopen('G:\科研工具\Spearman互相关系数\data_deal.txt','w'); fprintf(fid2,'%d %f %f %f \r\n',A'); fclose(fid2); 原数据网盘下载链接：http://pan.baidu.com/s/1pKFtj7d 其中，第一列和第三列是年月时间，第二列是径流数据，第四列是降水量数据。 |
回复此楼» 猜你喜欢
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
基金申报
已经有3人回复
-
国自然申请面上模板最新2026版出了吗？
已经有9人回复
-
溴的反应液脱色
已经有6人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
推荐一本书
已经有10人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
