要理解Gram-Schmidt过程是如何被翻译成这段代码作为实现的

试图从以下解释中理解Gram-Schmidt过程:

http://mlwiki.org/index.php/Gram-Schmidt_Process

计算步骤对我来说是有意义的。然而,同一文章中包含的Python实现似乎并不一致。

def normalize(v):    return v / np.sqrt(v.dot(v))n = len(A)A[:, 0] = normalize(A[:, 0])for i in range(1, n):    Ai = A[:, i]    for j in range(0, i):        Aj = A[:, j]        t = Ai.dot(Aj)        Ai = Ai - t * Aj    A[:, i] = normalize(Ai)

从上面的代码中,我们看到它对V1和b进行了点积运算,但是(V1,V1)部分并未作为分母使用(参见下面的方程)。我想知道下面的方程是如何被翻译成代码中的for循环的?

enter image description here


回答:

这就是代码所做的

enter image description here

基本上,它对前一个向量(A中的列)进行归一化,并将当前向量投影到它上面,然后从当前向量中减去这个投影。

为了计算的整洁,每个向量都会进行归一化处理。

上面的V2方程没有对前一个向量进行归一化,因此存在差异。

Related Posts

L1-L2正则化的不同系数

我想对网络的权重同时应用L1和L2正则化。然而,我找不…

使用scikit-learn的无监督方法将列表分类成不同组别,有没有办法?

我有一系列实例,每个实例都有一份列表,代表它所遵循的不…

f1_score metric in lightgbm

我想使用自定义指标f1_score来训练一个lgb模型…

通过相关系数矩阵进行特征选择

我在测试不同的算法时,如逻辑回归、高斯朴素贝叶斯、随机…

可以将机器学习库用于流式输入和输出吗?

已关闭。此问题需要更加聚焦。目前不接受回答。 想要改进…

在TensorFlow中,queue.dequeue_up_to()方法的用途是什么?

我对这个方法感到非常困惑,特别是当我发现这个令人费解的…

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注