我想生成正随机半定矩阵。我正在寻找一种算法,或更理想的是使用C,matlab,java或任何语言的简单算法实现。
示例代码(Python):
from scipy import random, linalg matrixSize = 10 A = random.rand(matrixSize,matrixSize) B = numpy.dot(A,A.transpose()) print 'random positive semi-define matrix for today is', B