python实现高斯(Gauss)迭代法

最新推荐文章于 2024-09-10 08:45:16 发布

此生无分起相思

最新推荐文章于 2024-09-10 08:45:16 发布

阅读量7.4k

点赞数

分类专栏： python 数学

本文链接：https://blog.csdn.net/cswfqxs_/article/details/84068660

版权

python 同时被 2 个专栏收录

2 篇文章 1 订阅

订阅专栏

数学

2 篇文章 0 订阅

订阅专栏

#Gauss迭代法 输入系数矩阵mx、值矩阵mr、迭代次数n(以list模拟矩阵 行优先)
def Gauss(mx,mr,n=100):
    if len(mx) == len(mr):  #若mx和mr长度相等则开始迭代 否则方程无解
        x = [] #迭代初值 初始化为单行全0矩阵
        for i in range(len(mr)):
            x.append([0])
        count = 0 #迭代次数计数
        while count < n:
            for i in range(len(x)):
                nxi = mr[i][0]
                for j in range(len(mx[i])):
                    if j!=i:
                        nxi = nxi+(-mx[i][j])*x[j][0]
                nxi = nxi/mx[i][i]
                x[i][0] = nxi
            count = count + 1
        return x 
    else:
        return False

#调用 Gauss(mx,mr,n=100) 示例
    

mx = [[8,-3,2],[4,11,-1],[6,3,12]]
mr = [[20],[33],[36]]
print(Gauss(mx,mr,20))