计算方法-Jacobi（雅可比）迭代法
太累了，不想写太多说明了，看ppt看代码吧，推荐华东理工大学的慕课
链接：https://www.bilibili.com/video/av83437134
计算方法-Jacobi（雅可比）迭代法
针对上面的课堂作业，我写了下面的代码，运行成功
"""
@Date:2020/2/24
@Author:indigoDeveloper
@E-mail：changwenhan0530@gmail.com
"""
import numpy as npdef Jacob(a,b):x1 = x2 = x3 = 1 #设置(1,1,1)为初始点count = 0X = [[x1],[x2],[x3]]Xarray = np.array(X) #转为矩阵barray = np.array(b)array = np.array(a) #读入矩阵L = U = H = [[0,0,0],[0,0,0],[0,0,0]]#创建同等大小的矩阵L UD = [[4,0,0],[0,5,0],[0,0,7]]Darray = np.array(D)Darray_N = np.linalg.inv(Darray) #取逆矩阵for i in range(len(a)):for j in range(len(a[0])):#为L添加内容if i>j:L[i][j] = a[i][j]#为U添加内容elif j>i:U[i][j] = a[i][j]#相加生成H阵for i in range(len(a)):for j in range(len(a[0])):H[i][j] = (L[i][j] + H[i][j])/2Harray = np.array(H)#开始迭代while count<20:a_1 = np.dot(Darray_N,Harray)a_2 = np.dot(a_1,Xarray)a_3 = -1*a_2a_4 = np.dot(Darray_N,barray)Xarray = a_3 + a_4count += 1print('第{}次'.format(count))print(Xarray)
def main():A=[[4,2,1],[2,5,-2],[1,-2,7]]B=[[23],[29],[0]]Jacob(A,B)if __name__ == "__main__":main()
给你们看看效果
就到这儿吧，每次都是深夜写代码才写的下去，或者该改改习惯了
2020/2/24
计算方法-Jacobi（雅可比）迭代法相关推荐
- 用matlab求雅可比迭代法,基于matlab的jacobi(雅可比)迭代法求解线性方程组
说明推导见此博客:https://blog..net/zengxyuyu/article/details/53054880 源码见下面: main.m clear clc A = [8 -3 2;4 ...
- 雅可比（Jacobi）迭代法解线性方程组的Matlab实现
雅可比(Jacobi)迭代法解线性方程组的Matlab实现 代码 运行 代码 迭代法解线性方程组的基本思想是构造一串收敛到解的序列,即建立一种从已有近似解计算新的近似解的规则,有不同的计算规则得到不同 ...
- 计算方法：列主元消去法，LU分解法， 雅可比迭代法，高斯塞德尔迭代法 解线性方程（C++）
Matrix.h包括矩阵类Matrix的定义,Matrix.cpp包括该类成员函数的实现,LinearEqu.h包括线性方程类LinearEqu的定义,继承自Matrix类,其中solve()方法为列 ...
- C++使用雅可比迭代法(jacobi)和高斯赛德尔迭代法(G-S)解线性方程组
理论部分: 简言之 jacobi迭代法就是把线性方程组化为关于第i行xi的等式,然后利用等式循环计算xi Gauss-Seidel迭代法就是在jacobi迭代的基础上,在计算第i+1行的xi+1时,带 ...
- 雅可比迭代法的收敛性及谱半径（自己总结心得）
说到雅可比迭代法收敛性,首先引入矩阵谱半径的概念 **谱半径** 当谱半径ρ小于1,即收敛,且半径越少,收敛速度越快. 雅可比迭代法 #include<iostream> #include<cstdio> #include<cmath> using name ...
最新文章
- 熊出没之伐木机器人_熊出没第一首富之争：李老板真的是第一吗？难道不是他第一？...
- iOS 工程师被要求写小程序，过分吗？| 每日趣闻
- 【童心制物】一篇很硬的标新立异级别的体验测评——聊新版造物编程盒
- 京东「卖家日志」系统的构建 | 流式计算日志系统应用实践
- 年前最后一期！知识图谱理论与实战（线下集训）
- C语言试题五十九之请编写一个函数fun，它的功能时：求fibonacci数列中大于t的最小的一个数，结果由函数返回。
- 如果你需要在 XHTML 中声明 DOCTYPE,必须使用到jsp:text动作元素
- C++ 对比 C语言
- JavaScript学习（十六）—实现购物车加减数量，计算总金额
- string的operate+=
- 双击java安装包没反应_wps安装包双击没反应
- Django ---uploads files
- Oracle drop table
- 当面试官问到以后的职业规划，这样回答绝对称赞你是高手！（二）
- 非常详细的STM32 CAN通信的贴子，从总线细节到编程实现
- php内部网,php – WordPress内部网址重写是如何工作的？
- 爬取qq音乐的评论并生成词云——以《听妈妈的话》为例
- 来看看怎么通过a标签打开一个对话框
- 达梦数据库迁移数据报错 字符串截断
- 联想win8改win7