超松弛迭代法求方法组的解(Python实现)
数值分析题目
求方程组
{5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2\left\{ \begin{array}{c} 5x_1+2x_2 + x_3 = -12 \\ -x_1 + 4x_2+2x_3 = 20 \\ 2x_1+-3x_2+10x_3=2 \end{array} \right. ⎩⎨⎧5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2
超松弛迭代其实基于高斯-塞德尔迭代的
可以对比着之前的高斯迭代法对比一下效果
高斯–塞德尔迭代法求方程组的解(Python实现)
代码
import numpy as npA = np.array([[5.0, 2, 1], [-1, 4, 2], [2, -3, 10]])
B = np.array([-12.0, 20, 3])
x0 = np.array([1.0, 1, 1])
x = np.array([0.0, 0, 0])
w = 1.5
times, MT = 0, 1000while times < MT:tempx = x0.copy()for i in range(3):temp = 0for j in range(3):if i != j:temp += x0[j] * A[i][j]x[i] = (B[i] - temp) / A[i][i]x0[i] = x[i]x = w * x + (1-w) * tempxcalTemp = max(abs(x - tempx))times += 1if calTemp < 1e-4:break
print(times)
print(x)
- 结果算出来的X为[-4, 3 , 2]
超松弛迭代法求方法组的解(Python实现)相关推荐
- MATLAB用超松弛迭代法求接地金属槽内的电位分布
工程电磁场仿真实验-用超松弛法求接地金属槽内电位分布 一.实验目的 1.熟悉利用超松弛法描绘出接地金属槽的电位分布情况: 2.学会使用MATLAB进行数值计算,并绘出相应的图形: 二.实验原理 1.有 ...
- python不是内部或外部命令的解决方法_详解python常见报错--NO MODULE NAMED _SQLITE3解决方法...
概述 在安装一些基于python的程序时,经常遇到"ImportError: No module named _sqlite3"问题. 下面记录下解决方法. 解决办法 1.编译sq ...
- python random库生成伯努利随机数的方法_详解Python基础random模块随机数的生成
随机数参与的应用场景大家一定不会陌生,比如密码加盐时会在原密码上关联一串随机数,蒙特卡洛算法会通过随机数采样等等.Python内置的random模块提供了生成随机数的方法,使用这些方法时需要导入ran ...
- python中的str方法和repr方法_详解Python中__str__和__repr__方法的区别
对我当前工程进行全部测试需要花费不少时间.既然有 26 GB 空闲内存,为何不让其发挥余热呢? tmpfs 可以通过把文件系统保存在大内存中来加速测试的执行效率. 但优点也是缺点,tmpfs 只把结果 ...
- C++实现牛顿迭代法求一元二次方程
牛顿迭代法定义如下(来自百度百科): 用牛顿迭代法小试牛刀用来求解一元二次方程的根(工程下载地址[注:不好意思,下载题目写成了二元一次方程,见谅.],release下的应用程序可以直接运行),代码如下 ...
- 高斯--塞德尔迭代法求方程组的解(Python实现)
数值分析题目 求方程组 {5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2\left\{ \begin{array}{c} 5x_1+2x_2 + x_3 = - ...
- MATLAB用SOR方法解方程组,SOR迭代法 求MATLAB程序
满意答案 vwbog 2014.11.02 采纳率:41% 等级:6 已帮助:1507人 function [n,x]=sor22(A,b,X,nm,w,ww) %用超松弛迭代法求解方程组Ax= ...
- python迭代法求解方程_第一部分:趣味算法入门;第六题牛顿迭代法求一元三次方程的根...
100个不同类型的python语言趣味编程题 在求解的过程中培养编程兴趣,拓展编程思维,提高编程能力. 第一部分:趣味算法入门:第六题SRE实战 互联网时代守护先锋,助力企业售后服务体系运筹帷幄!一键 ...
- c语言迭代法解线性代数方程,迭代法求线性代数方程组程序.docx
迭代法求线性代数方程组 计算物理 辽宁科技大学 理学院 骆宾祥 学号:120123802038 迭代法求线性代数方程组 骆宾祥 摘 要:目前在许多实际应用领域,诸如航空.造船以及其它结构工程中,常遇到 ...
最新文章
- Elsevier的Greg Landrum访谈 | 成功的开源化学信息软(RDKit)的要素是什么?
- css如何让图片不平铺,css怎么设置图片平铺方式?
- 【caffe解读】 caffe从数学公式到代码实现5-caffe中的卷积
- 还是畅通工程(思想+代码)
- Android Studio 2.2 HTTP proxy setting 提示异常
- Oracle用户角色权限管理
- 如何在 Active Directory 中还原已删除的用户帐户及其组成员身份
- 小米android点击,自动点击器小米版
- 干货!10分钟,用Python生成图文并茂的PDF报告!
- 联想笔记本电脑BIOS中英文对照
- 40岁计算机研究生能干什么工作,快四十岁了,还能干点什么?
- 我读过的(股票)书籍
- 固态硬盘和机械硬盘区别 固态硬盘和机械硬盘的优缺点
- 计算机加内存还是固态硬盘,电脑慢加内存还是固态硬盘好
- 安卓紧急警报_我们的紧急警报系统依赖于重叠的私有服务混乱局面
- HDU 6069 题解
- 日本亚马逊海淘经验-2015
- 11,MSI文件简介
- 前端学习笔记DAY1
- SpringBoot中MybatisX插件的简单使用教程(超详细!!)