请问Matlab如何生成一个随机的半正定矩阵?

N=5;
a=eye(N)*sprandsym(N,3);
while(prod((1:N)'.*(eig(a)>0))==0)
a=eye(N)*sprandsym(N,3);
end
a
eig(a)
a一定可以满足要求,至于运算时间么,看你运气了