作者iyenn (晓风)
看板Grad-ProbAsk
标题Re: [理工] [线代]-QR分解
时间Sun Feb 7 15:51:26 2010
※ 引述《lovefo (lovefo)》之铭言：
: ※ 引述《atled (乔巴)》之铭言：
: : 想请问一下
: : 一个没有满秩的矩阵
: : 例如
: : [ 1 0 0 0 ]
: : [ 0 1 0 0 ]
: : [ 0 0 3 6 ]
: : [ 0 0 2 4 ]
: : 该怎麽做QR分解呢
: : GS正交化要取三个还是四个阿?
: V1=[1 0 0 0]^T
: V2=[0 1 0 0]^T
: V3=[0 0 3 2]^T
: V4=[0 0 6 4]^T - 26/13 [0 0 3 2]^T
: =[0 0 0 0]
: A=[U1 U2 U3 U4]
: V1=U1 =>U1 = V1
: V2=U2 =>U2 = V2
: V3=U3 =>U3 = V3
: V4=U4 - 2U3 =>U4 = V4 + 2U3
: 1 0 0 0 1 0 0 0
: 0 1 0 0 0 1 0 0
: A=[ 0 0 3/√13 0] [0 0 √13 2√13]
: 0 0 2/√13 0 0 0 0 0
: Q R
: Q:4*3 R:3*4
: 不知道对不对
换我有问题了= =a
我在书上查了一下
theorem.
if A is a mxn matrix with linearly independent columns,
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
then A can be factored as A=QR,where Q is an mxn matrix
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
whose columns form an orthonormal basis for the column
space of A,
Proof,
let u1,u2,..,un be the LI columns of A,
by G-S process,
v1=u1
<ui,v1> <ui,v2>
vi=ui- -------v1 - -------v2 -...
<v1,v1> <v2,v2>
1
finally, wi=------vi
||vi||
ui can be rewritten as a linear combination of the w-vectors
u1=r11w1+r21w2+...+rn1wn
u2=r12w1+r22w2+...+rn2wn
...
un=r1nw1+r2nw2+...+rnnwn
span{v1,v2,...,vi}=span{w1,w2,..,wi}
since wj is ortheogonal to span{w1,w2,..,wi}
for j > i,it is orthogonal to ui
->hence rji=0 for j>i
let Qbe the matrix whose columns are w1,w2,...,wj.
let rj=[r1j]
[r2j]
[...]
[rnj]
A=[u1 u2 ... un]
=[Qr1 Qr2 ... Qrn]=QR
R=[r11 r12 ... r1n]
[0 r22 ... r2n]
[0 0 ... ...]
[0 0 ... rnn]
不是行独立不能拆吧??
而且硬拆R也不是上三角矩阵吧...
In the mathematical discipline of linear algebra,
a triangular matrix is a special kind of square matrix
where the entries either below or above
the main diagonal are zero.
By wiki.
感觉只是类似QR拆解的东西.但不是QR?
请高手救救我的线代Orz
--
为者常成,行者常至.
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 123.193.214.165
※ 编辑: iyenn 来自: 123.193.214.165 (02/07 15:52)
※ 编辑: iyenn 来自: 123.193.214.165 (02/07 16:06)
1F：推 CRAZYAWIND:我是顷向0的那行删掉= =97台大生医96 97中山电机己的 02/07 16:09
2F：→ CRAZYAWIND:都是0那行删掉= = 总不会这几张考卷都出错吧 02/07 16:09
4F：→ iyenn:第一题, 02/07 16:09
5F：推 CRAZYAWIND:烦脑了= = 到底该怎做呢............ 02/07 16:11
6F：→ iyenn:我认为不能拆 Orz. 但是题目真的出计算,我还是会拆... 02/07 16:11
7F：→ iyenn:毕竟计算题写不能拆当答案...感觉不好XDD 02/07 16:12
8F：→ iyenn:我觉的这指是类QR拆解 02/07 16:12
9F：→ iyenn:有没有高手说明一下>"< 02/07 16:12
10F：→ CRAZYAWIND:我也是觉得如果写不能拆........很危险 02/07 16:12
11F：→ iyenn:但是像96中正资工那题 我就会写B非LI不能拆 QR 02/07 16:13
12F：→ iyenn:算了,等高手说明吧...继续大扫除 02/07 16:13
14F：→ CRAZYAWIND:f 台大的一题QR 20分耶= = 如果写不能拆20分很重 02/07 16:16
15F：推 lovefo:我回去翻小黄的讲义 我上面有把他删掉 可能当时有讲吧XD 02/07 16:20
16F：→ iyenn:所以是= =? 02/07 16:21
17F：推 lovefo:因为题库班讲义是别人给我的 他是上小黄的 他有划掉 02/07 16:27
18F：推 luckysky1:orthonormal basis 行长度要为1 所以要划掉 02/07 17:44
19F：推 cherry770519:没错 小黄上课有强调 把0那行划掉 变成一个梯型矩阵 02/08 13:37