亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

Jacobi迭代算法的Python實現詳解

系統 1592 0
            
import numpy as np
import time
          

1.1 Jacobi迭代算法

            
def Jacobi_tensor_V2(A,b,Delta,m,n,M):
start=time.perf_counter()#開始計時
find=0#用于標記是否在規定步數內收斂
X=np.ones(n)#迭代起始點
x=np.ones(n)#用于存儲迭代的中間結果
d=np.ones(n)#用于存儲Ax**(m-2)的對角線部分
m1=m-1
m2=2-m
for i in range(M):
print('X',X)
a=np.copy(A)
#得Ax**(m-2)
for j in range(m-2):
a=np.dot(a,X)
#得d 和 (2-m)Dx**(m-2)+(L'+U')x**(m-2)
for j in range(n):
d[j]=a[j,j]
a[j,j]=m2*a[j,j]
#迭代更新
for j in range(n):
x[j]=(b[j]-np.dot(a[j],X))/(m1*d[j])
#判斷是否滿足精度要求
if np.max(np.fabs(X-x))
            
          

1.2 張量A的生成函數和向量b的生成函數:

            
def Creat_A(m,n):#生成張量A
size=np.full(m, n)
X=np.ones(n)
while 1:
#隨機生成給定形狀的張量A
A=np.random.randint(-49,50,size=size)
#判斷Dx**(m-2)是否非奇異,如果是,則滿足要求,跳出循環
D=np.copy(A)
for i1 in range(n):
for i2 in range(n):
if i1!=i2:
D[i1,i2]=0
for i in range(m-2):
D=np.dot(D,X)
det=np.linalg.det(D)
if det!=0:
break
#將A的對角面張量擴大十倍,使對角面占優
for i1 in range(n):
for i2 in range(n):
if i1==i2:
A[i1,i2]=A[i1,i2]*10
print('A:')
print(A)
return A
#由A和給定的X根據Ax**(m-1)=b生成向量b
def Creat_b(A,X,m):
a=np.copy(A)
for i in range(m-1):
a=np.dot(a,X)
print('b:')
print(a)
return a
          

1.3 對稱張量S的生成函數:

            
def Creat_S(m,n):#生成對稱張量B
size=np.full(m, n)
S=np.zeros(size)
print('S',S)
for i in range(4):
#生成n為向量a
a=np.random.random(n)*np.random.randint(-5,6)
b=np.copy(a)
#對a進行m-1次外積,得到秩1對稱張量b
for j in range(m-1):
b=outer(b,a)
#將不同的b疊加得到低秩對稱張量S
S=S+b
print('S:')
print(S)
return S
def outer(a,b):
c=[]
for i in b:
c.append(i*a)
return np.array(c)
return a
          

1.4 實驗一

            
def test_1():
Delta=0.01#精度
m=3#A的階數
n=3#A的維數
M=200#最大迭代步數
X_real=np.array( [2,3,4])
A=Creat_A(m,n) 
b=Creat_b(A,X_real,m)
Jacobi_tensor_V2(A,b,Delta,m,n)
          

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 久久久精品国产四虎影视 | 亚洲国产精品a一区 | 四虎国产精品4hu永久 | 亚洲精品日本高清中文字幕 | 国产精品亚洲综合一区在线观看 | 久久综合玖玖爱中文字幕 | 岛国不卡 | 欧美视频一区在线 | 91精品国产91久久久久久麻豆 | 成人三级视频在线观看 | 四虎影城| 爱搞逼综合 | 亚洲波多野结衣日韩在线 | 亚洲高清二区 | 亚洲国产精品免费 | 精品久久久久久影院免费 | 在线亚洲成人 | 老色鬼a∨在线视频在线观看 | 天天尻逼 | 欧洲美女a视频一级毛片 | 四虎影院国产精品 | 99亚洲精品视频 | 亚洲视频综合网 | 国内精品久久影视 | 色国产精品一区在线观看 | 亚洲一区二区三区国产精品 | 2021国产精品系列一区二区 | 97视频在线观看视频在线精品 | 国产成人综合亚洲 | 亚洲精品www久久久久久 | 亚洲欧美一区二区三区二厂 | 日韩一区二区免费看 | 波多野结衣手机视频一区 | 国产成人免费观看在线视频 | 久久综合狠狠综合久久 | 欧美jizzhd精品欧美另类 | 亚洲最新在线 | 欧美精品亚洲精品日韩专 | 国产高清精品一级毛片 | 亚洲高清国产一线久久 | 亚洲国产成人久久综合一区77 |