ndarray列修改异常：Gram-Schmidt正交化过程中赋值失效

核心原因 #

你的数组a大概率是整数类型（如int32/int64）。当你将计算得到的浮点数结果赋值给整数数组的列时，浮点数会被自动截断为整数，像0.22这类小于1的数值会直接变成0，最终导致赋值后列全为0。

验证方式 #

执行以下代码查看数组类型：

print(a.dtype)

如果输出是int32、int64等整数类型，即可确认此问题。

解决方案 #

1. 提前指定浮点数类型：初始化数组时直接设置 dtype 为浮点数，比如：

   a = np.array(your_initial_data, dtype=np.float64)
   

2. 转换数组类型：在处理前将整数数组转为浮点数数组：

   a = a.astype(np.float64)
   

3. 简化冗余代码：你代码中连续的reshape(-1,1).reshape(-1)完全多余，同时要注意循环范围应取列数而非行数，可直接简化为：

   for j in range(a.shape[1]):  # 用a.shape[1]获取列数，替代len(a)（len(a)返回行数）
       norm = np.linalg.norm(a[:, j])
       a[:, j] = a[:, j] / norm
   

内容的提问来源于stack exchange，提问作者noor napso