MATLAB入门学习-#6-Jacobi、Gauss-Seidel、SOR迭代法编程练习

2024-02-03 08:38

本文主要是介绍MATLAB入门学习-#6-Jacobi、Gauss-Seidel、SOR迭代法编程练习,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

MATLAB入门学习-#6-Jacobi、Gauss-Seidel、SOR迭代法编程练习

  • 1.Jacobi迭代法
  • 2.Gauss-Seidel迭代法
  • 3.SOR迭代法(松弛法)

这三种迭代法是在数值分析课程里学到的,都是求解线性方程组用的,相关的知识可以在数值计算方法的第三章里面找,可以说已经把最核心的怎么算和为什么这么算给交代了,剩下的任务就仅仅是编程实现就完了。
毕竟自己也没怎么编过程序,更是很久以来第一次写function,所以还是比较简陋的,但是也能完成相应的操作。在此就仅仅是记录一下自己学习matlab的历程把…

1.Jacobi迭代法


function [B,g,x]=jacobi(A,b,times,x0)
%[B,g,x]=jacobi(A,b,times,x0)
%
%B为jacobi迭代矩阵
%g为jacobi变形后的b
%x为方程组的解
%A为系数矩阵
%b为矩阵右侧
%times为迭代次数
%x0为初始x
D=diag(diag(A));
c=size(A);
I=eye(c);
B=I-D\A;          %使用左除或者右除要比使用inv计算快
g=D\b;
if vrho(B)>=1    %求B的谱error('jacobi迭代法不收敛!');
end
i=1;
while i<=timesx=B*x0+g;       %来自书上x0=x;disp(['第',num2str(i),'次']);disp(x);i=i+1;
end

2.Gauss-Seidel迭代法

function [M,x]=gauss_seidel(A,b,x0,times)
%[M,x]=gauss_seidel(A,b,x0,times)
%
%M为jacobi迭代矩阵
%x为方程组的解
%A为系数矩阵
%b为矩阵右侧
%times为迭代次数
%x0为初始x0
D=diag(diag(A));
L=D-tril(A);
U=D-triu(A);
M=(D-L)\U;
if vrho(M)>=1     %求B的谱error('jacobi迭代法不收敛!');
end
i=1;
while i<=timesx=((D-L)\U)*x0+(D-L)\b;    %来自书上x0=x;disp(['第',num2str(i),'次']);disp(x);i=i+1;
end

3.SOR迭代法(松弛法)

function [M,x]=SOR(A,b,x0,times,w)
%[M,x]=SOR(A,b,x0,times,w)
%
%M为jacobi迭代矩阵
%x为方程组的解
%A为系数矩阵
%b为矩阵右侧
%times为迭代次数
%x0为初始x0
%w为参数因子
D=diag(diag(A));
L=D-tril(A);
U=D-triu(A);
M=(D-w.*L)\((1-w).*D+w.*U);
if vrho(M)>=1     %求B的谱error('jacobi迭代法不收敛!');
end
i=1;
while i<=timesx=M*x0+(D-w.*L)\(w.*b);    %来自书上x0=x;disp(['第',num2str(i),'次']);disp(x);i=i+1;
end

这篇关于MATLAB入门学习-#6-Jacobi、Gauss-Seidel、SOR迭代法编程练习的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/673578

相关文章

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.

Spring Boot + MyBatis Plus 高效开发实战从入门到进阶优化(推荐)

《SpringBoot+MyBatisPlus高效开发实战从入门到进阶优化(推荐)》本文将详细介绍SpringBoot+MyBatisPlus的完整开发流程,并深入剖析分页查询、批量操作、动... 目录Spring Boot + MyBATis Plus 高效开发实战:从入门到进阶优化1. MyBatis

Java并发编程必备之Synchronized关键字深入解析

《Java并发编程必备之Synchronized关键字深入解析》本文我们深入探索了Java中的Synchronized关键字,包括其互斥性和可重入性的特性,文章详细介绍了Synchronized的三种... 目录一、前言二、Synchronized关键字2.1 Synchronized的特性1. 互斥2.

Python异步编程中asyncio.gather的并发控制详解

《Python异步编程中asyncio.gather的并发控制详解》在Python异步编程生态中,asyncio.gather是并发任务调度的核心工具,本文将通过实际场景和代码示例,展示如何结合信号量... 目录一、asyncio.gather的原始行为解析二、信号量控制法:给并发装上"节流阀"三、进阶控制

Java进阶学习之如何开启远程调式

《Java进阶学习之如何开启远程调式》Java开发中的远程调试是一项至关重要的技能,特别是在处理生产环境的问题或者协作开发时,:本文主要介绍Java进阶学习之如何开启远程调式的相关资料,需要的朋友... 目录概述Java远程调试的开启与底层原理开启Java远程调试底层原理JVM参数总结&nbsMbKKXJx

Python FastAPI入门安装使用

《PythonFastAPI入门安装使用》FastAPI是一个现代、快速的PythonWeb框架,用于构建API,它基于Python3.6+的类型提示特性,使得代码更加简洁且易于绶护,这篇文章主要介... 目录第一节:FastAPI入门一、FastAPI框架介绍什么是ASGI服务(WSGI)二、FastAP

Java深度学习库DJL实现Python的NumPy方式

《Java深度学习库DJL实现Python的NumPy方式》本文介绍了DJL库的背景和基本功能,包括NDArray的创建、数学运算、数据获取和设置等,同时,还展示了如何使用NDArray进行数据预处理... 目录1 NDArray 的背景介绍1.1 架构2 JavaDJL使用2.1 安装DJL2.2 基本操

C#多线程编程中导致死锁的常见陷阱和避免方法

《C#多线程编程中导致死锁的常见陷阱和避免方法》在C#多线程编程中,死锁(Deadlock)是一种常见的、令人头疼的错误,死锁通常发生在多个线程试图获取多个资源的锁时,导致相互等待对方释放资源,最终形... 目录引言1. 什么是死锁?死锁的典型条件:2. 导致死锁的常见原因2.1 锁的顺序问题错误示例:不同

PyCharm接入DeepSeek实现AI编程的操作流程

《PyCharm接入DeepSeek实现AI编程的操作流程》DeepSeek是一家专注于人工智能技术研发的公司,致力于开发高性能、低成本的AI模型,接下来,我们把DeepSeek接入到PyCharm中... 目录引言效果演示创建API key在PyCharm中下载Continue插件配置Continue引言

C#反射编程之GetConstructor()方法解读

《C#反射编程之GetConstructor()方法解读》C#中Type类的GetConstructor()方法用于获取指定类型的构造函数,该方法有多个重载版本,可以根据不同的参数获取不同特性的构造函... 目录C# GetConstructor()方法有4个重载以GetConstructor(Type[]