雅可比与高斯赛德尔方法求解线性代数方程（Julia）

用Julia语言写的迭代解法，求解如下方程

PS：该方法用雅可比方法收敛而高斯赛德尔方法不收敛

using LinearAlgebra
using Statistics# 构建方程
A = [1.0 2 -21 1 12 2 1
]
b = [1.035
]# 雅可比迭代矩阵B和g
D = Diagonal(A)
E = -(UpperTriangular(A)-D)
F = -(LowerTriangular(A)-D)
B = D\(E + F)
g = D\b# 求解函数
function solve(B,g,mode)size = length(g)res = zeros(size)res_last = ones(size)while(abs(mean(res - res_last))>1.0E-5)res_last = resres = B * res + gif mode == 1println(res)endendreturn res
end# 求解
solve(B,g,1)# 构建高斯赛德尔方法矩阵并求解
B = (D-E)\F
g = (D-E)\b
solve(B,g,1)

雅可比与高斯赛德尔方法求解线性代数方程（Julia）相关推荐

用matlab求解线性代数方程
用matlab求解线性代数方程 % 功能:线性代数方程求解:AX=B % 解的判定矩阵:C=[A B] % 求解方法分三种情况讨论: % 唯一解:A为非奇异方阵,x=inv(A)*B % 无穷解:ra ...
【算法设计zxd】第3章迭代法03 线性代数方程
线性代数方程:  设线性代数方程组具有如下特征:  算法框架 (1)设置线性代数方程组的初值X=(x1,-,xn-1 ,xn): (2)构造迭代方程xi=gi (X) (i=1,-,n-1,n ...
求解线性最小二乘问题的奇异值分解及广义逆法的C++实现
求解线性最小二乘问题的广义逆法的C++实现 1,功能 2,方法说明 3,函数语句与形参说明第一步,求对系数矩阵进行奇异值分解(muav函数) #include "stdlib.h" ...
2021-01-13 Matlab求解微分代数方程 (DAE)
Matlab求解微分代数方程 (DAE) 什么是微分代数方程? 微分代数方程是一类微分方程,其中一个或多个因变量导数未出现在方程中.方程中出现的未包含其导数的变量称为代数变量,代数变量的存在意味着不能 ...
奇异值分解（SVD）方法求解最小二乘问题
奇异值分解(SVD)方法求解最小二乘问题 1 奇异值分解(SVD)原理 1.1 回顾特征值和特征向量 1.2 SVD的定义 1.3 求出SVD分解后的U,Σ,V矩阵 1.4 SVD的一些性质 2 线性 ...
c语言编程非线性方程求解,c语言计算机编程三种方法求解非线性方程
c语言计算机编程三种方法求解非线性方程本科专业学年论文题目:非线性方程求解比较姓名: 何娟专业: 计算机科学技术系班级: 08 级本科(2)班指导老师: 刘晓 ...
c语言编程非线性方程求解,c语言计算机编程三种方法求解非线性方程.doc
c语言计算机编程三种方法求解非线性方程.doc 本科专业学年论文题目非线性方程求解比较姓名何娟专业计算机科学技术系班级 08 级本科(2)班指导老师刘晓娜 ...
求解欧拉方程的c语言,用有限体积方法求解欧拉方程
<用有限体积方法求解欧拉方程>由会员分享,可在线阅读,更多相关<用有限体积方法求解欧拉方程(12页珍藏版)>请在人人文库网上搜索. 1.有限体积法求解二维可压缩Euler方程计 ...
哈希查找解决地址冲突的两种最常见方法（线性探测再散列，链地址法）C++实现
哈希查找解决地址冲突的两种最常见方法(线性探测再散列,链地址法)C++实现参考文章: (1)哈希查找解决地址冲突的两种最常见方法(线性探测再散列,链地址法)C++实现 (2)https://www. ...

雅可比与高斯赛德尔方法求解线性代数方程（Julia）

雅可比与高斯赛德尔方法求解线性代数方程（Julia）相关推荐

最新文章

热门文章