zdata=[a,b,c,d];
a=[0.3;3;4;.3;.4;.3];
x=zdata;
k=length(zdata(:,1));
n=length(zdata(1,:));
copula_cdf=a(4)*(sum((x.^(-a(1)))')'-n+1).^(-1/a(1))+a(5)*exp(-sum(((-log(x)).^a(2))')'.^(1/a(2)))-a(6)/a(3)*log(1+prod((exp(-a(3)*x)-1)')'/(exp(-a(3))-1)^(n-1));
for i=1:k
    empirical_cdf(i)=mean((x(:,1)<=x(i,1)).*(x(:,2)<=x(i,2)).*(x(:,3)<=x(i,3)).*(x(:,4)<=x(i,4)));
end
empirical_cdf=empirical_cdf';
[h,p,ksstat,cv]=kstest(empirical_cdf,copula_cdf)
程序是这样的，最后结果说列向量必须是个矩阵，那位大神帮忙看看怎么改啊