用QR方法解決特徵值問題寫出matlab程序原始碼

時間 2022-10-07 02:30:04

1樓:匿名使用者

前面的函式是求特徵值的函式,後面的是呼叫的例子和驗證function l = rqrtz(a,m)%qr演算法求矩陣全部特徵值

%已知矩陣:a

%迭代步數:m

%求得的矩陣特徵值:l

a = hess(a);

for i=1:m

n = size(a);

n = n(1,1);

u = a(n,n);

[q,r]=qr(a-u*eye(n,n));

a = r*q+u*eye(n,n);

l = diag(a);

enda=[0 5 0 0 0 0;1 0 4 0 0 0;0 1 0 3 0 0;0 0 1 0 2 0;0 0 0 1 0 1;0 0 0 0 1 0]

a =0 5 0 0 0 01 0 4 0 0 00 1 0 3 0 00 0 1 0 2 00 0 0 1 0 10 0 0 0 1 0>> rqrtz(a,50)

ans =

-3.2030

3.2030

-1.8837

1.8837

-0.6167

0.6167

>> eig(a)

ans =

-3.3243

3.3243

-1.8892

-0.6167

1.8892

0.6167

2樓:匿名使用者

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

%已知矩陣:a

%迭代步數:m

%求得的矩陣特徵值:l

a = hess(a);

for i=1:m

n = size(a);

n = n(1,1);

u = a(n,n);

[q,r]=qr(a-u*eye(n,n));

a = r*q+u*eye(n,n);

l = diag(a);

end------------------------------------

a=[0 5 0 0 0 0;1 0 4 0 0 0;0 1 0 3 0 0;0 0 1 0 2 0;0 0 0 1 0 1;0 0 0 0 1 0]

a =0 5 0 0 0 01 0 4 0 0 00 1 0 3 0 00 0 1 0 2 00 0 0 1 0 10 0 0 0 1 0>> rqrtz(a,50)

ans =

-3.2030

3.2030

-1.8837

1.8837

-0.6167

0.6167

>> eig(a)

ans =

-3.3243

3.3243

-1.8892

-0.6167

1.8892

0.6167

確定地基承載力特徵值的方法包括哪些

左岸居東 地基承載力特徵值確定方法可歸納為三種 1.按土的抗剪強度指標以理 式計算。2.按地基荷載實驗或觸探實驗確定。3.按有關規範提供的承載力或經驗公式確定。地基承載力特徵值的存在,其目的很明確,就是為了掌握地基的承載規律,發揮地基的承載能力,合理確定地基承載力,讓地基不致於因荷載作用而發生剪下破...

確定地基承載力有哪些方法,地基承載力特徵值是什麼,確定方法有哪些

封疆大吏 1 原位試驗法 是一種通過現場直接試驗確定承載力的方法。包括 靜 載荷試驗 靜力觸探試驗 標準貫入試驗 旁壓試驗等,其中以載荷試驗法為最可靠的基本的原位測試法。2 理 式法 是根據土的抗剪強度指標計算的理 式確定承載力的方法。3 規範 法 是根據室內試驗指標 現場測試指標或野外鑑別指標,通...

初中最值問題解決方法,初中數學求最值問題的方法

肖斌綿陽 學習中沒有高手,學無分先後,達者即可為師。交流而已 初中涉及的數學求最值問題,複雜點就是二次函式在區間 t1,t2 內求最大值或最小值 最值與極值的區別就是,極大值可能是最大值,可能不是最大值,與誰比較?端點函式值 極小值可能是最小值,也可能不是最小值,與誰比較?端點函式值 所以,知識點要...