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))<Delta: find=1 break X=np.copy(x) end=time.perf_counter()#结束计时 print('时间:',end-start) print('迭代',i) return X,find,i,end-start
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)
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
华山资源网 Design By www.eoogi.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
华山资源网 Design By www.eoogi.com
暂无评论...
更新日志
2024年10月02日
2024年10月02日
- 纪钧瀚《钢琴阅读时光 雨中书店聆听轻音乐》[FLAC/分轨][399.62MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[320K/MP3][87.4MB]
- 证声音乐图书馆《走向自然 疗心爵士乐》[FLAC/分轨][184.94MB]
- 陈慧娴.2018-Priscilla-Ism演唱会3CD(2024环球红馆40复刻系列)【环球】【WAV+CUE】
- 郑秀文.1999-我应该得到(国)【华纳】【WAV+CUE】
- 陈家慧.2011-钢琴酒吧2CD【龙吟唱片】【WAV+CUE】
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[320K/MP3][45.01MB]
- 证声音乐图书馆《雨季 蓝调吉他 Rainy Blues》[FLAC/分轨][109.13MB]
- 赞多《序章》[320K/MP3][45.54MB]
- 许巍.2004-每一刻都是崭新的【步升大风】【WAV+CUE】
- 群星.2024-四方馆影视原声带【韶愔音乐】【FLAC分轨】
- 陈雷.1997-安锁咧【金圆唱片】【WAV+CUE】
- 关淑怡.2013-MY.FAVORITE.SK.3CD【环球】【WAV+CUE】
- Sweety.2006-花言乔语【丰华】【WAV+CUE】
- 李恕权.2003-回·20年全精选2CD【SONY】【WAV+CUE】