博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
矩阵论——正交矩阵&Gram-Schimidt正交化
阅读量:3960 次
发布时间:2019-05-24

本文共 1293 字,大约阅读时间需要 4 分钟。

标准正交基:

q i T q j = { 0 if     i   / = j 1 if     i = j q_i^Tq_j=\begin{cases}0&\text{if }\;i \mathrlap{\,/}{=}j\\1&\text{if }\;i=j\end{cases} qiTqj={
01if i/=jif i=j
优点:1. 计算方便
   2. 不会上溢,不会下溢
   3.所有方程使用标准正交基后,都变得非常简单。

对于标准正交列向量组成的矩阵 Q = [ q 1 , q 2 , . . . , q n ] Q=\begin{bmatrix}q_1,q_2,...,q_n\end{bmatrix} Q=[q1,q2,...,qn](当 Q Q Q是方阵时,称为正交矩阵),有

  1. Q T Q = I Q^TQ=I QTQ=I
  2. 其列空间的投影矩阵 P = Q Q T { = I if     Q    i s    s q u a r e P=QQ^T\begin{cases}=I&\text{if }\;Q\;is\;square\end{cases} P=QQT{
    =Iif Qissquare

注:即使 q i q_i qi不是 n 维向量。

A T A x ^ = A T b A^TA\hat x=A^Tb ATAx^=ATb

Now A is Q
x ^ = Q T b \hat x=Q^Tb x^=QTb
⟹ x ^ i = q i T b \Longrightarrow \hat x_i=q^T_ib x^i=qiTb
即已知标准正交基Q, b b b 在 第 i i i个 基方向上的投影的数量积 x ^ i \hat x_i x^i等于 q i T b q^T_ib qiTb.

Gram-Schmidt

目的

使线性无关向量组标准正交化.

算法

设有线性无关向量组 a , b , c a,b,c a,b,c

A = a A=a A=a

在这里插入图片描述

B = b − A T b A T A a B=b-\dfrac{A^Tb}{A^TA}a B=bATAATba

A T b A T A \dfrac{A^Tb}{A^TA} ATAATb b b b a a a 上的。

C = c − A T c A T A A − B T c B T B B C=c-\dfrac{A^Tc}{A^TA}A-\dfrac{B^Tc}{B^TB}B C=cATAATcABTBBTcB

q 1 = A ∣ ∣ A ∣ ∣ , q 2 = B ∣ ∣ B ∣ ∣ , q 3 = C ∣ ∣ C ∣ ∣ q_1=\dfrac{A}{||A||},q_2=\dfrac{B}{||B||},q_3=\dfrac{C}{||C||} q1=AA,q2=BB,q3=CC即为标准正交向量组。

特点
  1. A → Q A\rightarrow Q AQ,其列空间不变.
  2. A=QR,R是上三角矩阵.
缺点

计算时,可能会出现 根号

转载地址:http://vymzi.baihongyu.com/

你可能感兴趣的文章
C 指针解读
查看>>
有关乱码的处理---中国程序员永远无法避免的话题
查看>>
JSP的运行内幕
查看>>
python超简单的web服务器
查看>>
代理模式、静态代理、动态代理、aop
查看>>
Struts1.x Spring2.x Hibernate3.x DWR2.x整合工具文档v1.00
查看>>
大型Web2.0站点构建技术初探
查看>>
机器学习算法汇总:人工神经网络、深度学习及其它
查看>>
解决Spring中AOP不能切入Struts的DispatchAction方法的问题
查看>>
出国以后才知道英语应该怎么学
查看>>
计算机专业权威期刊投稿经验总结
查看>>
如何在三个月内学会一门外语?
查看>>
看看你对Linux到底了解多少?
查看>>
网上看到的:ARM入门最好的文章(转)
查看>>
中国最美情诗100句
查看>>
javascript注册window的onload事件问题研究
查看>>
客户端技术分页控件javascript+css,可用于任何服务器端技术
查看>>
学习Swing 的网站[转]
查看>>
Google App engine 的第一个应用 midispot
查看>>
提问的智慧
查看>>