施密特正交化过程编程c语言,利用C程序编写格拉姆-施密特正交化的过程.docx

利用C程序编写格拉姆-施密特正交化的过程

格拉姆-施密特正交化

在 HYPERLINK "/wiki/%E7%BA%BF%E6%80%A7%E4%BB%A3%E6%95%B0" \o "线性代数" 线性代数中，如果HYPERLINK "/wiki/????§ˉ??oé?′" \o "内积空间" 内积空间上的一组向量能够组成一个HYPERLINK "/wiki/?-???oé?′" \o "子空间" 子空间，那么这一组向量就称为这个子空间的一个基。Gram－Schmidt正交化提供了一种方法，能够通过这一子空间上的一个基得出子空间的一个HYPERLINK "/wiki/?-￡?o¤??o" \o "正交基" 正交基，并可进一步求出对应的HYPERLINK "/wiki/???????-￡?o¤??o" \o "标准正交基" 标准正交基。

这种正交化方法以HYPERLINK "/wiki/J%C3%B8rgen_Pedersen_Gram" \o "en:J?rgen Pedersen Gram"J?rgen Pedersen Gram和 HYPERLINK "/wiki/Erhard_Schmidt" \o "en:Erhard Schmidt" Erhard Schmidt命名，然而比他们更早的HYPERLINK "/wiki/%E6%8B%89%E6%99%AE%E6%8B%89%E6%96%AF" \o "拉普拉斯"拉普拉斯(Laplace)和 HYPERLINK "/wiki/%E6%9F%AF%E8%A5%BF" \o "柯西" 柯西(Cauchy)已经发现了这一方法。在HYPERLINK "/w/index.php?title=%E6%9D%8E%E7%BE%A4%E5%88%86%E8%A7%A3&action=edit&redlink=1" \o "李群分解(页面不存在)" 李群分解中，这种方法被推广为HYPERLINK "/wiki/?2??3????è§￡" \o "岩泽分解" 岩泽分解(HYPERLINK "/wiki/Iwasawa_decomposition" \o "en:Iwasawa decomposition" Iwasawa decomposition)。

在 HYPERLINK "/wiki/%E6%95%B0%E5%80%BC%E8%AE%A1%E7%AE%97" \o "数值计算" 数值计算中，Gram－Schmidt正交化是HYPERLINK "/wiki/??°????¨3?????§" \o "数值稳定性" 数值不稳定的，计算中累积的舍入误差会使最终结果的正交性变得很差。因此在实际应用中通常使用HYPERLINK "/wiki/è±a??ˉé???°???·?????￠" \o "豪斯霍尔德变换" 豪斯霍尔德变换或HYPERLINK "/wiki/Givens???è??" \o "Givens旋转" Givens旋转进行正交化。

记法

： HYPERLINK "/wiki/%E7%BB%B4%E6%95%B0" \o "维数" 维数为n?的内积空间

：中的元素，可以是向量、HYPERLINK "/wiki/?????°" \o "函数" 函数，等等

：与的 HYPERLINK "/wiki/%E5%86%85%E7%A7%AF" \o "内积" 内积

：、……张成的 HYPERLINK "/wiki/%E5%AD%90%E7%A9%BA%E9%97%B4" \o "子空间" 子空间

：在上的 HYPERLINK "/wiki/%E6%8A%95%E5%BD%B1" \o "投影" 投影

基本思想

Gram-Schmidt正交化的基本想法，是利用 HYPERLINK "/w/index.php?title=%E6%8A%95%E5%BD%B1%E5%8E%9F%E7%90%86&action=edit&redlink=1" \o "投影原理(页面不存在)" 投影原理在已有正交基的基础上构造一个新的正交基。

设。是上的维子空间，其标准正交基为，且不在上。由投影原理知，与其在上的投影之差

是正交于子空间的，亦即正交于的正交基。因此只要将单位化，即

那么就是在上扩展的子空间的标准正交基。

根据上述分析，对于向量组张成的空间?()，只要从其中一个向量(不妨设为)所张成的一维子空间开始(注意到就是的正交基)，重复上述扩展构造正交基的过程，就能够得到?的一组正交基。这就是Gram-Schmidt正交化。

算法

首先需要确定已有基底向量的顺序，不妨设为。Gram-Schmidt正交化的过程如下：

这样就得到上的一组正交基，以及相应的标准正交基。

例

考察如下 HYPERLINK "/wik