大气统计作业探索.docVIP

上机实践 相关回归分析 （1）计算北京夏季降水与气温的相关系数及秩相关系数，并进行显著性检验。 定义计算秩相关系数函数： {function coeff = mySpearman(X,Y); % 本函数用于实现斯皮尔曼等级相关系数的计算操作 % % 输入： % X：输入的数值序列 % Y：输入的数值序列 % % 输出： % coeff：两个输入数值序列X，Y的相关系数 if length(X)~=length(Y) error(两个数值数列的维数不相等); return; end N = length(X); %得到序列的长度 Xrank = zeros(1 , N); %存储X中各元素的排行 Yrank = zeros(1 , N); %存储Y中各元素的排行 %计算Xrank中的各个值 for i = 1 : N cont1 = 1; %记录大于特定元素的元素个数 cont2 = -1; %记录与特定元素相同的元素个数 for j = 1 : N if X(i) X(j) cont1 = cont1 + 1; els