需求

计算不同站点，两个变量月尺度序列上的相关系数

代码

function [R2,P2] = cor_p(data1,data2)
% By Yang,2019/8/15
% 计算两个变量间的相关性和P值
% 输入：data1,data2两个变量矩阵（行列数一致），每一行（m）代表一个样本变量
% 输出: R2一列（m）的相关系数
% 输出：P2一列（m）的P值
R2=[]; P2=[];
m = size(data1,1);  %有多少个变量
for i = 1:m
    [R,P] = corrcoef(data1(i,:),data2(i,:)); %计算一个变量的相关性和P值
    R2 = [R2;R(1,2)];                        %累加
    P2 = [P2;P(1,2)];
end
end