一个LU分解存在但不唯一的例子

最新推荐文章于 2023-11-28 08:57:02 发布

InactionWorld

最新推荐文章于 2023-11-28 08:57:02 发布

阅读量1.2k

点赞数 1

文章标签： python

本文链接：https://blog.csdn.net/weixin_43793510/article/details/121110492

版权

如题，对于矩阵A，只要保证L32=U33-2即可。

这是测试代码:

import numpy as np
import random

def isEqual(x):
    A = np.array([[1, 1, 1], [2, 2, 1], [3, 3, 1]])
    L = np.array([[1, 0, 0], [2, 1, 0], [3, x, 1]])
    U = np.array([[1, 1, 1], [0, 0, -1], [0, 0, x-2]])
    return A == np.matmul(L, U)

print(isEqual(random.randint(-10000, 10000)))
'''
[[ True  True  True]
 [ True  True  True]
 [ True  True  True]]
'''