mathematica实现非负矩阵分解

Mathematica中可以使用NonnegativeMatrixFactorization函数实现非负矩阵分解。该函数需要两个参数：要分解的矩阵和分解后的矩阵的秩。以下是使用NonnegativeMatrixFactorization函数进行非负矩阵分解的示例：

假设我们有一个3x4的非负矩阵A：

A={{1,2,3,4},{5,6,7,8},{9,10,11,12}}

我们可以使用以下代码进行非负矩阵分解：

k=2;
{W,H}=NonnegativeMatrixFactorization[A,k];

其中，k=2表示我们希望将矩阵A分解为两个非负矩阵W和H的乘积，{W,H}表示分解后得到的两个矩阵。我们可以使用MatrixForm函数将它们打印出来：

MatrixForm[W]

MatrixForm[H]

注意，非负矩阵分解的结果不一定唯一，因此，不同的初始值可能会得到不同的分解结果。