数值分析题目

求方程组
{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实现)相关推荐

  1. MATLAB用超松弛迭代法求接地金属槽内的电位分布

    工程电磁场仿真实验-用超松弛法求接地金属槽内电位分布 一.实验目的 1.熟悉利用超松弛法描绘出接地金属槽的电位分布情况: 2.学会使用MATLAB进行数值计算,并绘出相应的图形: 二.实验原理 1.有 ...

  2. python不是内部或外部命令的解决方法_详解python常见报错--NO MODULE NAMED _SQLITE3解决方法...

    概述 在安装一些基于python的程序时,经常遇到"ImportError: No module named _sqlite3"问题. 下面记录下解决方法. 解决办法 1.编译sq ...

  3. python random库生成伯努利随机数的方法_详解Python基础random模块随机数的生成

    随机数参与的应用场景大家一定不会陌生,比如密码加盐时会在原密码上关联一串随机数,蒙特卡洛算法会通过随机数采样等等.Python内置的random模块提供了生成随机数的方法,使用这些方法时需要导入ran ...

  4. python中的str方法和repr方法_详解Python中__str__和__repr__方法的区别

    对我当前工程进行全部测试需要花费不少时间.既然有 26 GB 空闲内存,为何不让其发挥余热呢? tmpfs 可以通过把文件系统保存在大内存中来加速测试的执行效率. 但优点也是缺点,tmpfs 只把结果 ...

  5. C++实现牛顿迭代法求一元二次方程

    牛顿迭代法定义如下(来自百度百科): 用牛顿迭代法小试牛刀用来求解一元二次方程的根(工程下载地址[注:不好意思,下载题目写成了二元一次方程,见谅.],release下的应用程序可以直接运行),代码如下 ...

  6. 高斯--塞德尔迭代法求方程组的解(Python实现)

    数值分析题目 求方程组 {5x1+2x2+x3=−12−x1+4x2+2x3=202x1+−3x2+10x3=2\left\{ \begin{array}{c} 5x_1+2x_2 + x_3 = - ...

  7. MATLAB用SOR方法解方程组,SOR迭代法 求MATLAB程序

    满意答案 vwbog 2014.11.02 采纳率:41%    等级:6 已帮助:1507人 function [n,x]=sor22(A,b,X,nm,w,ww) %用超松弛迭代法求解方程组Ax= ...

  8. python迭代法求解方程_第一部分:趣味算法入门;第六题牛顿迭代法求一元三次方程的根...

    100个不同类型的python语言趣味编程题 在求解的过程中培养编程兴趣,拓展编程思维,提高编程能力. 第一部分:趣味算法入门:第六题SRE实战 互联网时代守护先锋,助力企业售后服务体系运筹帷幄!一键 ...

  9. c语言迭代法解线性代数方程,迭代法求线性代数方程组程序.docx

    迭代法求线性代数方程组 计算物理 辽宁科技大学 理学院 骆宾祥 学号:120123802038 迭代法求线性代数方程组 骆宾祥 摘 要:目前在许多实际应用领域,诸如航空.造船以及其它结构工程中,常遇到 ...

最新文章

  1. Elsevier的Greg Landrum访谈 | 成功的开源化学信息软(RDKit)的要素是什么?
  2. css如何让图片不平铺,css怎么设置图片平铺方式?
  3. 【caffe解读】 caffe从数学公式到代码实现5-caffe中的卷积
  4. 还是畅通工程(思想+代码)
  5. Android Studio 2.2 HTTP proxy setting 提示异常
  6. Oracle用户角色权限管理
  7. 如何在 Active Directory 中还原已删除的用户帐户及其组成员身份
  8. 小米android点击,自动点击器小米版
  9. 干货!10分钟,用Python生成图文并茂的PDF报告!
  10. 联想笔记本电脑BIOS中英文对照
  11. 40岁计算机研究生能干什么工作,快四十岁了,还能干点什么?
  12. 我读过的(股票)书籍
  13. 固态硬盘和机械硬盘区别 固态硬盘和机械硬盘的优缺点
  14. 计算机加内存还是固态硬盘,电脑慢加内存还是固态硬盘好
  15. 安卓紧急警报_我们的紧急警报系统依赖于重叠的私有服务混乱局面
  16. HDU 6069 题解
  17. 日本亚马逊海淘经验-2015
  18. 11,MSI文件简介
  19. 前端学习笔记DAY1
  20. SpringBoot中MybatisX插件的简单使用教程(超详细!!)

热门文章

  1. 技术团队的情绪与效率
  2. 如何才能写出一手高质量优美的代码
  3. struct 类型指针技巧
  4. 读《Oracle DBA工作笔记》知识点-获取创建语句
  5. RSS - 简单方便的follow资讯
  6. linux系统常用快捷键,文件类型,后缀名,语言修改,转义
  7. dtree.js树的使用
  8. linux服务器ftp上传文件为空,Linux 服务器 ftp上传文件出现的问题
  9. Nutch爬虫引擎使用分析
  10. 数据结构源码笔记(C语言):英文单词按字典序排序的基数排序