协方差函数matlab代码

z=[1,5;3,4;4,6;5,3];

%%
%matlab函数结果
% b=cov(z)
% 
%     2.9167   -0.8333
%    -0.8333    1.6667
%%
z_mean=mean(z);
result = zeros(2,2);
for i = 1:2
    for j=1:2
        result(i,j)=(z(:,i)-z_mean(i))'*(z(:,j)-z_mean(j));
        %result(i,j)=sum((z(:,i)-z_mean(i))'*(z(:,j)-z_mean(j)))/(size(z,1)-1);%用这个公式得到与matlab相同的结果
    end
end
%自编函数结果 
%     8.7500   -2.5000
%    -2.5000    5.0000