numpy中的QuiRk-使用反射器在QR分解中将所有归一化元素分配为0

2024-04-20 • 问答

数组中所有归一化的元素都设置为0。

我正在尝试为数组创建反射器。当我用Python中的第一个元素(x[1:] /= x[0],x[0] = 1)对数组中的所有元素进行归一化时，我在第一个元素中得到的都是零和1。

import numpy as np
#project x via householder
u = np.array([9,4,1,3,1])
x = np.array([3,1])

for i in range(1,len(u)):
  u[i] /= u[0]
u[0] = 1

结果应该像这样(u = [1,4/9,1/9,3/9,1/9])

lizihao3hao 回答：numpy中的QuiRk-使用反射器在QR分解中将所有归一化元素分配为0

当我们定义一个numpy数组时，我们可以通过在np.array中指定dtype来定义它的类型。

import numpy as np
#project x via householder
u = np.array([9,4,1,3,1],dtype = np.float64)

for i in range(1,len(u)):
  u[i] /= u[0]
u[0] = 1
print(u)

给我们

[1.         0.44444444 0.11111111 0.33333333 0.11111111]

linear-algebra numerical-methods numpy qr-decomposition

本文链接：https://www.f2er.com/3169601.html

numpy中的QuiRk-使用反射器在QR分解中将所有归一化元素分配为0

lizihao3hao 回答：numpy中的QuiRk-使用反射器在QR分解中将所有归一化元素分配为0

大家都在问