如何用householder變換求矩陣的qr分解例子

2021-03-22 06:53:03 字數 3227 閱讀 1994

1樓:匿名使用者

||[householder陣]

(1) 設a rn, = ||a||2,通常取 與a1同號,記h=i-2vvt,(v= ),

則ha= - e1. h=i -2vvt稱為householder陣。

(2) 更一般地,對a=(a1,a2,…am,am+1,…,an)t,記 = ,可求出h,使

ha=(a1,a2,…am, ,0,…,0)t。

為此,先在rn-m中求 使 滿足

=(am+1,…,an)t=(- ,0,…,0,0)t,

再作h= ,則ha= (a1,a2,…am,am+1,…,an)t =( a1,a2,…am,- ,0,…,0,0)t

[用householder方法求矩陣的qr分解]

記a=(aij)n*n,由1可知,存在h1=i -2v1v1t,使

h1(a11,a21,…,an1)t=(a11(1),0,…,0)t,

於是 h1a=

又由1知,存在h2= ,使 ,於是

h1a= =

類似地依次進行n-1次,得出

hn-1hn-2…h1a= 。

記r=hn-1hn-2…h1a,q=hnhn-1…h1,得a=q*r

急求matlab複數矩陣qr分解**我知道matlab本身有qr函式 但是我想知道利用householder變換遞迴實現的原理。

2樓:匿名使用者

function [q,r]=qrhs(a)n=size(a,1);

r=a;

q=eye(n);

for i=1:n-1

x=r(i:n,i);

y=[1;zeros(n-i,1)];

ht=householder(x,y);

h=blkdiag(eye(i-1),ht);

q=q*h;

r=h*r;end

求助householde矩陣qr分解

3樓:匿名使用者

|[householder陣]

(1) 設a rn, = ||a||2,通常取 與a1同號,記h=i-2vvt,(v= ),

則ha= - e1. h=i -2vvt稱為householder陣。

(2) 更一般地,對a=(a1,a2,…am,am+1,…,an)t,記 = ,可求出h,使

ha=(a1,a2,…am, ,0,…,0)t。

為此,先在rn-m中求 使 滿足

=(am+1,…,an)t=(- ,0,…,0,0)t,

再作h= ,則ha= (a1,a2,…am,am+1,…,an)t =( a1,a2,…am,- ,0,…,0,0)t

[用householder方法求矩陣的qr分解]

記a=(aij)n*n,由1可知,存在h1=i -2v1v1t,使

h1(a11,a21,…,an1)t=(a11(1),0,…,0)t,

於是 h1a=

又由1知,存在h2= ,使 ,於是

h1a= =

類似地依次進行n-1次,得出

hn-1hn-2…h1a= 。

記r=hn-1hn-2…h1a,q=hnhn-1…h1,得a=q*r

如何用householder變換求hessenberg矩陣

4樓:q我

[householder陣]

(1) 設a rn, = ||a||2,通常取 與a1同號,記h=i-2vvt,(v= ),

則ha= - e1. h=i -2vvt稱為householder陣。

(2) 更一般地,對a=(a1,a2,…am,am+1,…,an)t,記 = ,可求出h,使

ha=(a1,a2,…am, ,0,…,0)t。

為此,先在rn-m中求 使 滿足

=(am+1,…,an)t=(- ,0,…,0,0)t,

再作h= ,則ha= (a1,a2,…am,am+1,…,an)t =( a1,a2,…am,- ,0,…,0,0)t

[用householder方法求矩陣的qr分解]

5樓:深圳華邦瀛

你描述的我根本不明白你在說什麼 你遇到什麼問題?想達到什麼功能?

通過householder變換實現復矩陣的qr變換和實矩陣有什麼區別

6樓:電燈劍客

matlab裡面的一撇運算表示的就是共軛轉置, 所以理論上講householder變換的**基本上是不需要做什麼改動的.

你要小心的是生成householder變換向量那一步, 在實數域當中只有v=x+||x||_2e_1和v=x-||x||_2e_1這兩種選擇, 但是在複數域上不僅僅是1,-1這兩個選擇, 我估計你這步沒寫好, 別的地方應該都不用改**.

qr分解的介紹

7樓:匿名使用者

這裡給出乙個(2×2)矩陣a,在qr分解後用迭代法求解特徵值的過程,僅供參考。

8樓:紫月軍團

qr分解法是目前求一般矩陣全部特徵值的最有效並廣泛應用的方法,一般矩陣先經過正交相似變化成為hessenberg矩陣,然後再應用qr方法求特徵值和特徵向量。它是將矩陣分解成乙個正規正交矩陣q與上三角形矩陣r,所以稱為qr分解法,與此正規正交矩陣的通用符號q有關。

9樓:安徽新華電腦專修學院

function l = rqrtz(a,m)%qr演算法求矩陣全部特徵值

%已知矩陣:a

%迭代步數:m

%求得的矩陣特徵值:l

求乙個用c語言實現的複數矩陣的qr分解,需要完整程式,十分感謝,拜託各位

10樓:明通造化

樓主找到沒,同求……

乙個矩陣可qr分解的充要條件?如何進行qr分解?

11樓:夢想隊員

任何乙個矩陣都可以進行qr分解。有兩種方式:施密特正交化;householder矩陣法

用哪個方法做矩陣qr分解比較好?

12樓:電燈劍客

取決於矩陣的性質

如果沒什麼特殊條件的話householder變換最好, 既穩定工作量又小

如何用matlab將矩陣中的數字變顏色

matlab讀取之後是復乙個rgb三原色組成的三維矩 制陣,現在需要將圖中黑 綠 黃 紅 圖上有兩種紅色需要視作一種 四種顏色的畫素點分別以0,0.01,1和2來表示,希望得到畫素點的二維矩陣,維數等於解析度,求具體實現的 望採納 matlab實現矩陣裡數值的大小對應顏色的深淺 資料如果存在矩陣a裡...

ps如何把證件照變漂亮,如何用PS處理證件照的面部變白

要想拍出漂亮點的 也有幾個小竅門的 1.頭稍微往下傾,眼睛正視攝像頭 2.嘴型要自然,稍微微微笑 3.然後就是要借助photoshop的力量了,是現在普遍用的方法。直接美顏相機,一步搞定 按快捷鍵ctrl m 拉曲線,看自己的感覺調。如何用ps處 照的面部變白 操作步驟如下 來 1 首先開啟自ps軟...

如何用iphone拍出電影,如何用iphone拍出單反效果

因為iphone的感光元件真的很小,因此想要拍出超淺景深 是不可能的,但是利用這個小技巧,就能在拍攝小物體 拍攝食物時增加一些 的可看性喔 1.開啟iphone相機,然後在畫面中長按你想要對焦的物體,大約按三秒後畫面就會出現帶有字眼為 ae af鎖定 的黃色標籤,當這個畫面出現之後,對焦好的距離設定...