第二節(jié)QR分解知識發(fā)現(xiàn)

1、第一節(jié)第一節(jié) qrqr分解分解qrqr分解也稱為正交三角分解分解也稱為正交三角分解 矩陣矩陣qrqr分解是一種特殊的三角分解，在解決分解是一種特殊的三角分解，在解決矩陣特征值的計算、最小二乘法等問題中起到重矩陣特征值的計算、最小二乘法等問題中起到重要作用。要作用。主要內(nèi)容：主要內(nèi)容：1 1矩陣的矩陣的qrqr分解分解- schmidt- schmidt正交化方法正交化方法2 2矩陣的矩陣的qrqr分解分解- householder- householder變換、變換、 givensgivens變換變換1學習園地qrqr分解定理分解定理任意一個滿秩實任意一個滿秩實( (復）矩陣復）矩陣a a，都

2、可唯一地分解，都可唯一地分解a = qr a = qr ，其中其中q q為為正交（酉）矩陣，正交（酉）矩陣，r是具有正是具有正對角元的上三角矩陣。對角元的上三角矩陣。由于由于x x 1 1, ,x x 2 2, , , ,x x n n 線性無關，將它們用線性無關，將它們用schmidtschmidt正交正交證明證明設設a a是一個實滿秩矩陣是一個實滿秩矩陣, a, a的的n n個列向量為個列向量為 x x 1 1, ,x x 2 2, , , ,x x n n 定義定義:設設.nnca如果存在如果存在n階酉矩陣階酉矩陣q和和n階上三角矩陣階上三角矩陣r r，使得，使得qra 則稱之為則稱之為

3、a a的的qrqr分解或酉三角分解分解或酉三角分解當當 時，則稱為時，則稱為a的正三角分解的正三角分解nnra化方法得標準正交向量化方法得標準正交向量e e 1 1, ,e e 2 2, , , ,e e n n2學習園地nnnnnnebebebxebebxebx221122211221111其中其中nibii, 2 , 1,0從而有從而有nnnnnnbbbbbbeeexxx2221121121213學習園地nnnnnbbbbbbreeeq2221121121,令令iqqt則則則則如果如果再證唯一性再證唯一性,11rqqra由此得由此得dqrrqq1111式中式中d=rd=r1 1r r-1-

4、1仍為具有正對角元的上三角矩陣。由于仍為具有正對角元的上三角矩陣。由于 dddqdqqqittt11即即d d為正交矩陣，因此為正交矩陣，因此d d為單位矩陣（正規(guī)上三角為對角陣）為單位矩陣（正規(guī)上三角為對角陣）故故rdrrqdqq111,4學習園地說明：說明：1若不要求若不要求r具有正對角元，則具有正對角元，則a的不同的不同qr分解僅在正交分解僅在正交矩陣的列和上三角矩陣矩陣的列和上三角矩陣r的對應行相差模為的對應行相差模為1的因子。的因子。該定理的證明過程給出了利用該定理的證明過程給出了利用schmidtschmidt正交化方法求可逆矩陣正交化方法求可逆矩陣qrqr分解的方法。分解的方法。

5、例例 求矩陣求矩陣a a的的qrqr分解分解110201221a解解，則，則記記122,102,011321xxx2 2若若a a為滿秩復矩陣，則存在酉矩陣為滿秩復矩陣，則存在酉矩陣q q與復非奇異上三角矩與復非奇異上三角矩陣陣r r，使，使a = qr a = qr 5學習園地tyyyxyyyxtyyyxyyxyyxyyxyxyxy2 , 1 , 121 , 1, 131231132),(),(1),(),(33121),(),(2211222311131112將將 正交化正交化321,xxxtyytyytyyeee2 , 1 , 11 , 1, 10 , 1 , 1663332221332

6、211單位化單位化6學習園地336233132121122322eeexeexex整理得整理得,03633663322663322q令令363300302222rqra則則7學習園地例例1 1：利用：利用schmidtschmidt正交化方法求矩陣的正交化方法求矩陣的qrqr分解分解212240130a設設,2 , 2, 1,1 , 4 , 3,2 , 0 , 0321tttxxx則則 321,xxx線性無關，首先將它們正交化得：線性無關，首先將它們正交化得：,2 , 0 , 011txy1),(),(221112yxyyyyx2),(),(1),(),(3322231113yyxyyyyxy

7、yyxtyyx0 ,56,5851213tyx0 , 4 , 31212再單位化再單位化：,1 , 0 , 02111tye,0 ,54,535122tye8學習園地,0 ,53,542133tye于是：于是：1112eyx21212521eeyyx32132132251eeeyyyx從而從而 qra00153540545302150212,1 , 0 , 02111tye,0 ,54,535122tye9學習園地householderhouseholder變換變換o+otihr2)(3)(h則則記記即：該變換將向量即：該變換將向量 變成了以變成了以 為法向量為法向量的平面的對稱向量的平面的對

8、稱向量 。householderhouseholder變換又稱為反射變換或鏡像變換，有明變換又稱為反射變換或鏡像變換，有明顯的幾何意義。在顯的幾何意義。在 中，給定一個向量中，給定一個向量 ，令，令 表示表示 關于平面關于平面 （以（以 為法向量）為法向量）的反射變換所得像，的反射變換所得像，如圖所示，如圖所示，3r10學習園地定義定義 設設 是一個單位向量，令是一個單位向量，令nchih2)(則稱則稱h h是一個是一個householderhouseholder矩陣或矩陣或householderhouseholder變換。變換。性質(zhì)性質(zhì)5.1.1 5.1.1 設設h h是一個是一個house

9、holderhouseholder矩陣，則矩陣，則（1 1）h h是是hermitehermite矩陣，矩陣， ；（2 2）h h是酉矩陣，是酉矩陣， ；（3 3）h h是對合矩陣，是對合矩陣， ；（4 4）h h是自逆矩陣是自逆矩陣（5 5）diagdiag( (i i, ,h h ) ) 也是一個也是一個householderhouseholder矩陣矩陣; ;（6 6）det hdet h = -1 = -1。hhhihhhih2hh111學習園地其中其中 為實數(shù)。為實數(shù)。定理定理 設設 是一個單位向量，則對于任意的是一個單位向量，則對于任意的ncu ncxauhx uaxxah,2nc

10、當當 時，取單位向量時，取單位向量 使使0 auxnc0 xhauxxxxixhhh)(22)(存在存在householderhouseholder矩陣矩陣h h，使得，使得證明證明 當當x=0 x=0時，任取單位向量時，任取單位向量則則則則002)(hixh12學習園地所以所以 當當 時，取時，取aux ,2auxauxxauxauxauxixixhht22)(22)(uuaxuauaxxxauxauxhhhhh2)()(由于由于auauxauxauxxauxxxhhh)()()()(2)()()()()(2auxauxauxxauxxhhxauxxuaxxxxuauaxxxhhhhhhh)

11、(2)(22213學習園地推論推論1 1 對于任意的對于任意的 ，存在，存在householderhouseholder矩陣矩陣h h，使，使ncx1aehx其中其中 為實數(shù)。為實數(shù)。12,eaxxah) 1,(,2)(uuruuuihtnt1aehx2xa 推論推論2 2 對于任意的對于任意的 ，存在，存在householderhouseholder矩陣矩陣h hnrx上述結論表明，可以利用上述結論表明，可以利用householderhouseholder變換將任意向量變換將任意向量化為與第一自然基向量化為與第一自然基向量 平行的向量（共線）平行的向量（共線）。 nrx1e，其中，其中使得使

12、得得得14學習園地例例2 2 用用householderhouseholder變換將向量變換將向量化為與化為與 平行的向量。平行的向量。tiix2,232xte0, 0, 11iexh21iaeaxxah2,12ia325301211iiaexaex13iehx 因此因此解解 由于由于為了使為了使為實數(shù)，取為實數(shù)，取令令112102145105101512iiiiihh則則也可取也可取 或或3aia3說明說明15學習園地1 1 將矩陣將矩陣a a按列分塊按列分塊 , ,取取na,2121121111111,aeaeahih111200*,11121111bahhhahn利用利用househol

13、derhouseholder矩陣求矩陣的矩陣求矩陣的qrqr分解的步驟：分解的步驟：則則16學習園地2 2 將矩陣將矩陣 按列分塊，按列分塊，)1()1(1nncbnb,32122221221222,bebebuhuuih222222001hht2211200*0*)(caaahh)2()2(2nncc取取則則其中其中17學習園地121nhhhq則則 a=qra=qr依次進行下去，得到第依次進行下去，得到第n-1n-1個個n n階的階的householdhousehold矩陣矩陣h hn n-1-1，使得，使得raaaahhhnn*2112133因因 為自逆矩陣，令為自逆矩陣，令 ih18學習

14、園地例例2：已知矩陣：已知矩陣,112240130a利用利用householderhouseholder變變換求換求a a的的qrqr分解分解因為因為,2 , 0 , 01t記記, 2211a令令21111111eaeat1 , 0 , 121則則hih1112,001010100從而從而1302402121ah記記,3 , 4t則則, 5222b令令22222221ebeb,3 , 1101thih2222,43345119學習園地記記,43034000100122hht則則rahh20015021212取取0053404305121hhq則則qra20學習園地givensgivens變換變

15、換x 2yx o我們知道，平面坐標系我們知道，平面坐標系 中的旋轉角為中的旋轉角為 變換可變換可表示為表示為2rt t是正交矩陣，稱為平面旋轉矩陣。是正交矩陣，稱為平面旋轉矩陣。將其推廣到一般的將其推廣到一般的n n維酉空間中，維酉空間中，可以得到初等旋轉變換，也稱為可以得到初等旋轉變換，也稱為givensgivens變換。變換。cossinsincos,2121txxtyy21學習園地定義定義 設設記記n n階矩陣階矩陣ncsc,122 sc)()()()(111111lklkcssctkl由由 所確定的線性變換稱為所確定的線性變換稱為givensgivens變換或初等旋轉變換。變換或初等旋

16、轉變換。klt稱稱 為為givensgivens矩陣或初等旋轉矩陣；矩陣或初等旋轉矩陣；klt容易驗證，容易驗證，givensgivens矩陣是矩陣是酉矩陣酉矩陣，且，且 。 1detklt22學習園地定理定理 對于任意向量對于任意向量 ，存在，存在givensgivens變換變換 ，使得，使得 的第的第l l個分量為個分量為0 0，第，第k k個分量為非負實數(shù)，其個分量為非負實數(shù)，其余分量不變。余分量不變。ncxkltxtkltnkltnyyyxtxxxx,2121),( ,lkjxycxsxyxsxcyjjlkllkk證明證明 記記由由givensgivens矩陣的定義可得矩陣的定義可得2

17、3學習園地當當 時，取時，取c c=1,=1,s s=0=0，則，則t tkl kl = = i i, ,此時此時022lkxx),(, 0lkjxyyyjjlk當當 時，取時，取022lkxx2222,lkllkkxxxsxxxc),(002222222222lkjxyxxxxxxxxyxxxxxxxxxxyjjlklklklkllklklllkkkk, ,結論成立。結論成立。則則24學習園地與第一自然基向量與第一自然基向量推論推論 給定一個向量給定一個向量 ，則存在一組，則存在一組givensgivens矩陣矩陣 ， 使得使得ncxnttt11312,1212131exxtttnncx1e

18、tnxxxx,2112ttnxxxxxt, 0 ,3222112稱為用稱為用givensgivens變換化向量變換化向量證明證明 設設由上述定理存在由上述定理存在givensgivens矩陣矩陣使得使得共線。共線。25學習園地依此繼續(xù)下去，可以得出依此繼續(xù)下去，可以得出tnxxxxxxtt, 0, 0,433222112131222221121310, 0,exxxxxttttnn對于對于 又存在又存在givensgivens矩陣矩陣 ，使得，使得xt1213t26學習園地例例3 3 用用givensgivens變換化向量變換化向量 與第一自然基向量與第一自然基向量共線共線 tiix2,25,2222121xxixix5,5211isic1000525055212iiiit20512xt解解 由于由于取取則構造則構造givensgivens矩陣矩陣27學習園地3, 2, 523