下一节:调和Ritz值 上一级:基本概念 上一节:正交投影法

斜投影方法

在斜投影方法中,我们给定两个子空间 \mathcal{L}\mathcal{K}, 并寻求一个近似向量 \tilde{u} \in \mathcal{K} 和一个标量 \tilde \lambda \in \mathcal{C} 满足Petrov-Galerkin条件,
v^{\ast} (A - \tilde \lambda I) \tilde{u} = 0 \quad \forall \ v \in \mathcal{L} \ .
子空间 \mathcal{K} 被称为右子空间, \mathcal{L} 被称为左子空间。 可以设计一种类似于Rayleigh-Ritz过程的方法, 并且可以通过矩阵形式方便地描述, 即将近似特征向量 \tilde{u} 用某些基表示, 并为 \mathcal{L} 的基制定Petrov-Galerkin条件。 这次我们需要两个基,一个用 V 表示的子空间 \mathcal{K} 的基, 另一个用 W 表示的子空间 \mathcal{L} 的基。我们假设这两个基是双正交的,即 w^{\ast}_i v_j = \delta_{ij}
W^{\ast} V\ =\ I ,
其中 I 是单位矩阵。 然后,像之前一样写成 \tilde{u} = V y , 上述Petrov-Galerkin条件产生与正交投影相同的近似问题, 只是矩阵 B_m 现在定义为
B_m \ =\ W ^{\ast} A V .
为了使双正交对 V, W 存在, \mathcal{L}\mathcal{K} 必须满足以下附加假设。

对于 \mathcal{K}\mathcal{L} 的任意两个基 VW

\det ( W ^{\ast} V ) \ne 0 .

显然,这个条件不依赖于所选的特定基, 它等价于要求 \mathcal{K} 中的任何向量都不与 \mathcal{L} 正交。

斜投影方法得到的近似结果可能比正交投影方法的条件更差。 因此,人们可能会怀疑是否有必要使用斜投影方法。 然而,基于斜投影的方法可以提供一些优势。 特别是,它们可能能够同时计算出良好的左特征向量和右特征向量的近似。 正如后面将看到的,还有一些基于斜投影技术的方法, 需要的存储量远少于类似的正交投影方法。

通过将这种技术与一些(更昂贵的)正交投影方法的步骤相结合, 可以减少使用非正交 VW 的缺点。



下一节:调和Ritz值 上一级:基本概念 上一节:正交投影法
Susan Blackford 2000-11-20