一个简单的算法产生正半定矩阵 [英] A simple algorithm for generating positive-semidefinite matrices

问题描述

我要生成正随机半正定矩阵。我在寻找一种算法或以上preferably一个简单的实现算法的C，MATLAB，java或任何语言的。

解决方案

1. 在生成随机矩阵
2. 将它自己的换位乘以
3. 您已经获得了半正定矩阵。

例如code（蟒蛇）：

 从SciPy的进口随机的，linalg
matrixSize = 10
A = random.rand（matrixSize，matrixSize）
B = numpy.dot（A，A.transpose（））
打印随机正半定义基质今天是'，B
 

I want to generate positive random semi-definite matrices. I am looking for an algorithm or more preferably an simple implementation of the algorithm in C, matlab, java or any language.

解决方案

1. generate random matrix
2. multiply it by it's own transposition
3. you have obtained a positive semi-definite matrix.

Example code (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

这篇关于一个简单的算法产生正半定矩阵的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！