av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁(yè)技術(shù)文章
文章詳情頁(yè)

python共軛梯度法特征值迭代次數(shù)討論

瀏覽:22日期:2022-06-23 11:15:43
共軛梯度法,特征值聚堆情況下迭代次數(shù)討論

輸入各種特征值聚堆與分散時(shí)的矩陣,并應(yīng)用共軛梯度法,觀察迭代次數(shù)與聚堆情況的關(guān)系。

python共軛梯度法特征值迭代次數(shù)討論

因?yàn)閷?duì)角矩陣的對(duì)角線元素為其特征值,則用對(duì)角矩陣討論較為方便代碼

import numpy as npdef cg(x0, A, b): r0 = np.dot(A, x0) - b p0 = -r0 rk = r0 pk = p0 xk = x0 t = 0 #記錄迭代次數(shù) while np.linalg.norm(rk) >= 1e-6: rr = np.dot(rk.T, rk) ak = rr / np.dot(np.dot(pk.T, A), pk) xk = xk + ak * pk rk = rk + ak * np.dot(A, pk) bk = np.dot(rk.T, rk) / rr pk = -rk + bk * pk t += 1 return xk, t#輸入列表,生成以列表為對(duì)角元素的對(duì)角矩陣def Diagonal_matrix(D): n = len(D) diag = np.zeros((n,n)) for i in range(n): diag[i][i] = D[i] return diag#矩陣對(duì)角線元素D_1 = [1, 1, 1, 1, 1, 6, 7, 8, 9, 10]D_2 = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]D_3 = [0.8, 0.9, 1, 1.1, 1.2, 6, 7, 8, 9, 10]D_4 = [1 - 2*1e-7, 1 - 1e-7, 1, 1 + 1e-7, 1 + 2*1e-7, 6, 7, 8, 9, 10]D_5 = [1, 1, 1, 2, 2, 2, 3, 3, 3, 10]#初始值x0 = np.zeros((10,1))b = np.ones((10,1)) #生成對(duì)角矩陣diag1 = Diagonal_matrix(D_1)diag2 = Diagonal_matrix(D_2)diag3 = Diagonal_matrix(D_3)diag4 = Diagonal_matrix(D_4)diag5 = Diagonal_matrix(D_5)#共軛梯度法迭代x_1, n_1 = cg(x0, diag1, b)x_2, n_2 = cg(x0, diag2, b)x_3, n_3 = cg(x0, diag3, b)x_4, n_4 = cg(x0, diag4, b)x_5, n_5 = cg(x0, diag5, b)n = [n_1, n_2, n_3, n_4, n_5]#輸出for i in range(5): print(’矩陣’,i + 1 ,’的迭代次數(shù)為: ’, n[i])

python共軛梯度法特征值迭代次數(shù)討論

矩陣1,前5個(gè)元素聚堆且都為相同元素

python共軛梯度法特征值迭代次數(shù)討論

矩陣2,特征值分散

python共軛梯度法特征值迭代次數(shù)討論

矩陣3,前5個(gè)特征值聚堆,但是最大差為0.4 ,而cg法精度為1e-6

python共軛梯度法特征值迭代次數(shù)討論

矩陣4,前5個(gè)特征值聚堆,且相差最大小于1e-6

python共軛梯度法特征值迭代次數(shù)討論

矩陣5,三聚堆輸出:

python共軛梯度法特征值迭代次數(shù)討論

分析:

聚堆特征值可看作一個(gè)特征值 特征值差小于迭代精度時(shí)被看作聚堆 例如矩陣5,前三個(gè)對(duì)角元素看作一個(gè),4-6元素看作一個(gè),7-9看作一個(gè) 一共4個(gè)元素,則需要迭代4次

以上就是python共軛梯度法特征值迭代次數(shù)討論的詳細(xì)內(nèi)容,更多關(guān)于python共軛梯度法迭代的資料請(qǐng)關(guān)注好吧啦網(wǎng)其它相關(guān)文章!

標(biāo)簽: Python 編程
相關(guān)文章:
主站蜘蛛池模板: 欧美一级欧美三级在线观看 | 午夜性视频 | 久精品久久| 国产激情视频网址 | 成人欧美一区二区三区在线播放 | 国产成人99久久亚洲综合精品 | 国产四区 | 99精品国产在热久久 | 久久视频免费看 | 欧美一区二区免费 | 亚洲天堂免费在线 | 久久成人精品视频 | 一级毛片在线播放 | 亚洲中午字幕 | 美女激情av| 国产在线播放一区二区三区 | 中文字幕乱码一区二区三区 | 日韩国产精品一区二区三区 | 天天干干| 不卡一区二区三区四区 | 在线观看国产网站 | 毛片在线看片 | 国产无套一区二区三区久久 | 久久精品亚洲精品国产欧美 | 欧美一级黑人aaaaaaa做受 | 国产视频一视频二 | 91精品国产色综合久久不卡98 | 欧美日韩综合一区 | 成人免费一级 | 国产精品日韩在线观看 | 久久国产精品一区 | 国产色在线 | h在线免费观看 | 日本天堂视频在线观看 | 国产在线精品一区二区三区 | 国产在线www| 欧美精品久久 | 日韩在线小视频 | 亚洲国产一区二区三区四区 | 国产精品久久久久无码av | 四虎影视 |