神剑山庄资源网 Design By www.hcban.com
import numpy as np import time
1.1 Gauss-Seidel迭代算法
def GaussSeidel_tensor_V2(A,b,Delta,m,n,M): start=time.perf_counter() find=0 X=np.ones(n) d=np.ones(n) m1=m-1 m2=2-m for i in range(M): print('X',X) x=np.copy(X) #迭代更新 for j in range(n): a=np.copy(A) for k in range(m-2): a=np.dot(a,x) for k in range(n): d[k]=a[k,k] a[k,k]=m2*a[k,k] 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_2(): 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) GaussSeidel_tensor_V2(A,b,Delta,m,n)
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
神剑山庄资源网 Design By www.hcban.com
神剑山庄资源网
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
神剑山庄资源网 Design By www.hcban.com
暂无Gauss-Seidel迭代算法的Python实现详解的评论...
RTX 5090要首发 性能要翻倍!三星展示GDDR7显存
三星在GTC上展示了专为下一代游戏GPU设计的GDDR7内存。
首次推出的GDDR7内存模块密度为16GB,每个模块容量为2GB。其速度预设为32 Gbps(PAM3),但也可以降至28 Gbps,以提高产量和初始阶段的整体性能和成本效益。
据三星表示,GDDR7内存的能效将提高20%,同时工作电压仅为1.1V,低于标准的1.2V。通过采用更新的封装材料和优化的电路设计,使得在高速运行时的发热量降低,GDDR7的热阻比GDDR6降低了70%。
更新日志
2024年10月05日
2024年10月05日
- 群星《前途海量 电影原声专辑》[FLAC/分轨][227.78MB]
- 张信哲.1992-知道新曲与精丫巨石】【WAV+CUE】
- 王翠玲.1995-ANGEL【新艺宝】【WAV+CUE】
- 景冈山.1996-我的眼里只有你【大地唱片】【WAV+CUE】
- 群星《八戒 电影原声带》[320K/MP3][188.97MB]
- 群星《我的阿勒泰 影视原声带》[320K/MP3][139.47MB]
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[320K/MP3][148.91MB]
- 刘雅丽.2001-丽花皇后·EMI精选王【EMI百代】【FLAC分轨】
- 齐秦.1994-黄金十年1981-1990CHINA.TOUR.LIVE精丫上华】【WAV+CUE】
- 群星.2008-本色·百代音乐人创作专辑【EMI百代】【WAV+CUE】
- 群星.2001-同步过冬AVCD【环球】【WAV+CUE】
- 群星.2020-同步过冬2020冀待晴空【环球】【WAV+CUE】
- 沈雁.1986-四季(2012梦田复刻版)【白云唱片】【WAV+CUE】
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[FLAC/分轨][257.88MB]
- 《国语老歌 怀旧篇 3CD》[WAV/分轨][1.6GB]