代码：
A=[1,2,3,4;
3,4,3,6;
6,6,8,9];%输入任意待算矩阵
s=size(A);
cm=zeros(s(2));
for i=1:s(1,2)-1;
for j=i+1:s(1,2);
c=corrcoef(A(:,i),A(:,j));
cm(i,j)=c(1,2);
end
cm运行结果：
