jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析

41528d3028836879cd698677c3999917.gifjacobi迭代法实验MATLAB程序数值分析

例1. 求线性方程组 得近似解。精确解为x*=[3,2,1]’。 解:对方程进行移项就得 记为Ax=b,或写为x=B0 x+f,其中 取初始值,代入原方程组可得再将把它代入可得.反复利用这个计算过程,得到一向量序列和一般的计算公式(迭代公式) 简写为迭代到第10次有 从此例看出,由迭代法产生的向量序列x(k)逐步逼近方程组的精确解x*. 6.1常用迭代法 定义1 (ⅰ)对于给定的方程组x=Bx+f ,用公式逐步代入 求近似解的方法称为迭代法(或称为一阶定常迭代法,这里B 与k无关). (ⅱ)如果 存在(记为x*),称此迭代法收敛,显然x*就是方 程组的解,否则称此迭代法发散. 迭代法的流程图为: ①为初始向量, ②是判断条件, 即时停止运行 ③k 是循环次数。 ④中带入初始值,然后赋给 ① Jacobi迭代法 对一般方程组 从第一个方程解出x1,第二个方程解出x2,…,记成 用矩阵写法即x=Bx+g,B的对角元皆零,可拆成 B=L+UL是B下三角部分,U是B上三角部分Jacobi迭代法如下述. 任取初始近似x(0),对k=1,2,…计算 直至║x(k+1)-x(k)║≤ε,预定的精度.用矩阵记号,即任取初始近似x(0),对k=1,2,…计算 x(k+1)=Bx(k)+g, 直至║x(k+1)-x(k)║≤ε(通常对迭代法限定最大迭代次数也是必要的).Jacobi迭代法的流程图为: 在以上的流程图中,先读入数据,即先输入系数矩阵A,常数向量b, 初始值,停止条件和最大循环次数。图中是,在我们迭代公式中的。k 是循环次数,N 是最大循环次数。 例2. 利用Jacobi方法求方程组 的近似解。 解 把原方程改为 任取初始近似x(0),对k=1,2,…计算 直至║x(k+1)-x(k)║≤ε,预定的精度. 此即Jacobi迭代法.计算结果如下表. k x1(k) x2(k) x3(k) 0 0 0 0 1 0.720000000000 0.830000000000 0.840000000000 2 0.971000000000 1.070000000000 1.150000000000 3 1.057000000000 1.157100000000 1.248200000000 4 1.085350000000 1.185340000000 1.282820000000 5 1.095098000000 1.195099000000 1.294138000000 6 1.098337500000 1.198337400000 1.298039400000 7 1.099441620000 1.199441630000 1.299334980000 8 1.099811159000 1.199811158000 1.299776650000 9 1.099936445800 1.199936445900 1.299924463400 10 1.099978537270 1.199978537260 1.299974578340 11 1.099992769394 1.199992769395 1.299991414906 二 实验部分 本章实验内容: 实验题目:Jacobi迭代法,Gauss-Saidel迭代法,SOR迭代法。 实验内容:利用MATLAB ,编制求Ax=b的各迭代计算方法的程序。 实验目的:了解迭代法的运用性,进行各迭代法数值结果的比较,并找出一个计算量小的,使迭代法加速收敛的迭代方法。 编程要求:①利用迭代法,初始向量为x(0) ②同时利用Jacobi法和Gauss-Seidel法来进行对比。 ③利用SOR迭代法来进行对比。 计算算法:①Jacobi迭代法的算法为: ②Gauss-Saidel迭代法的算法为: ③SOR迭代法的算法为: 实验例题⑴: 条件:取 实验例题⑵: 条件:取选择适当的松弛因子。 程序①: function [X,Y]=JacobiGS(A,b,p,p1,del,max) % A为线性方程组的系数矩阵,b为自由项,p和p1为两种迭代法的初始解,del为限制数,max为循环的限制次数。 n=length(b); for k1=1:max for j=1:n Y(j)=(b(j)-A(j,[1:j-1])*p1(1:j-1)-A(j,[j+1:n])*p1(j+1:n))/A(j,j); if j==1 X(1)=(b(1)-A(1,[2:n])*p(2:n))/A(1,1); elseif j==n X(n)=(b(n)-A(n,[1:n-1])*(X(1:n-1)) )/A(n,n); else X(j)=(b(j)-A(j,[1:j-1])*(X(1:j-1)) -A(j,[j+1:n])*p(j+1:n))/A(j,j); end err=abs(norm(X -p)); reerr=err/(norm(X)+eps); p=X ; if (err

下载提示(请认真阅读)1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。

2.下载的文档,不会出现我们的网址水印。

下载文档到电脑,查找使用更方便

10 积分

还剩页未读,继续阅读 关 键 词:jacobi迭代法 MATLAB程序 jacobi迭代法MATLAB程序 matlab编程 Jacobi 迭代法matlab程序 数值分析 matlab程序 数值分析 雅可比迭代法

50bcf9764bbd4714f9fe0eb9e92f395a.gif  蚂蚁文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/566921.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

docker部署python web应用_安装docker并部署web项目

一.docker简介1、docker定义:docker是一个用来装应用的容器,就像杯子可以装水,笔筒可以装笔,书包可以放书一样。你可以把“hello world!”放到docker中,也可以把网站放到docker中,你可以把任何你…

java 水印乱码,linux java程序加水印及中文乱码方案

在linux,centos环境下的,生成的带文字的水印图片在显示为方框乱码。img.setFont(new Font("宋体", Font.BOLD, 20));写中文进入图片之前必须设置字体,而且这个字体必须支持中文,否则就会出现乱码或者方框、问号等等。Bu…

resultmap拿不到数据_阿里巴巴国际站每日电商运营工作数据表格

一、日常运营工作表1.数据日报统计每天的流量数据:2.数据周报统计每周的流量数据:3.数据月报统计每月流量数据:前半部分后半部分4.P4P月数据总览统计P4P的流量数据:大图包含上面所有的方案5.P4P日数据总览6.爆款数据统计每月统计爆…

php oracle新增数据类型,oracle数据类型

数据类型字符串类型字符串类型用来存储字符串数据char用于描述定长的字符型数据0~2000nchar用来存储Unicode字符集的定长字符型数据0~1000varchar2用于描述可变长的字符型数据0~4000nvarchar2用来存储unicode字符集的可变长字符型数据0~1000long用来存储边长的字符串0~2GBvarch…

jenkins ssh 远程部署_Jenkins部署jar到远程服务器

首选得确定Jenkins的安装完整,还有插件的安装,除默认插件,此处必备插件:Publish Over SSH,SSH plugin,因为我源码在gitlab所以关于gitlab的插件我也安装了,看各自的情况选择安装。然后到系统管理…

jieba分词_Jieba.el – 在Emacs中使用jieba中文分词

jieba.el在Emacs中使用jieba中文分词众所周知, Emacs并没有内置中文分词系统, 以至于 forward-word 和 backward-word 以及 kill-word 等以单词为单位的操作只能粗暴的标点符号进行确定所谓”词汇”, (其实是中文里的句子).jieba.el 利用nodejieba对buffer中的中文句子进行分割…

python的request请求401_Python模拟HTTPS请求返回HTTP 401 unauthorized错误

Python模拟HTTPS请求返回HTTP 401 unauthorized错误开始是使用的 httplib模块,代码如下:header {"Content-type": "application/json", "Accept": "*/*" }params { ‘source‘:‘en‘, ‘target‘:‘es‘, ‘…

java php 女生数量,萌妹子告诉你php和java如何选

原标题:萌妹子告诉你php和java如何选phpPHP 独特的语法混合了C、Java、Perl 以及 PHP 自创新的语法。它可以比 CGI 或者 Perl 更快速的执行动态网页PHP具有非常强大的功能,所有的CGI的功能PHP都能实现,而且支持几乎所有流行的数据库以及操作系…

new 结构体指针_Go:我应该用指针替代结构体的副本吗?

logo对于许多 golang 开发者来说,考虑到性能,最佳实践是系统地使用指针而非结构体副本。我们将回顾两个用例,来理解使用指针而非结构体副本的影响。1. 数据分配密集型让我们举一个简单的例子,说明何时要为使用值而共享结构体&…

Oracle 同音字查询,汉字的演变过程100字,汉字的演变图片

中国的汉字博大精深,历史悠久。早期甲骨文主要是线条图,统一而优美。后来,它逐渐演变成一幅由笔画组成的图画,这是当代汉字的基础商周时期,甲骨文逐渐成为金文。西周时期,青铜器铸造,铭文用于记…

python音频聚类_python实现鸢尾花三种聚类算法(K-means,AGNES,DBScan)

python实现鸢尾花三种聚类算法(K-means,AGNES,DBScan)发布时间:2020-08-31 21:23:24

matlab 实验数据 传递函数,《传递函数MATLAB实验》.ppt

传递函数MATLAB实验 练习: 各种典型环节的阶跃响应曲线 1、比例环节(K) 从图形库浏览器中拖曳Step(阶跃输入)、Gain(增益模块)、Scope模块到仿真操作画面,连接成仿真框图。 结论:比例环节 K 改变,则放大倍数改变。 2、积分环节( 1…

localstorage存储大小_Cookie 已凉,Web 存储该这么做!

本文经授权转自公众号CSDN(ID:CSDNnews)作者 | 浪里行舟责编 | 郭芮随着移动网络的发展与演化,我们手机上现在除了有原生 App,还能跑“WebApp”——它即开即用,用完即走。一个优秀的 WebApp 甚至可以拥有和原生 App 媲美的功能和体…

三折线弹塑性滞回模型matlab,动力弹塑性滞回模型-迈达斯汇总.doc

9-1 概要非线性抗震分析方法可分为非线性静力分析方法和非线性动力分析方法。其中非线性静力分析方法(静力弹塑性分析)因其理论概念易于理解、计算效率高、整理结果较为容易等原因为设计人员所广泛使用。但是由于静力弹塑性分析存在反映结构动力特性方面的缺陷、使用的能力谱是…

python实现knn算法鸢尾花_Python学习之knn实现鸢尾花分类

# K近邻算法# 导入相关库文件import numpy as npimport matplotlib.pyplot as plt#import pandas as pdfrom sklearn import neighbors, datasets# 导入数据集,数据集sklearn自带,X与y一一对应dataset datasets.load_iris()# 获取鸢尾花前两列花萼长度和…

php获取变量数据类型,php如何确定变量的数据类型

在php中,数据类型有:Boolean 布尔类型、Integer 整型、Float 浮点型、String 字符串、Array 数组、Object 对象、Resource 资源类型、NULL;知道一个数据的类型,能够更加有效地进行代码逻辑处理。1、使用 var_dump() 函数,可以获取…

深入理解java虚拟机 - jvm高级特性与最佳实践(第三版)_JVM虚拟机面试指南:年薪30W以上高薪岗位需求的JVM,你必须要懂!...

JVM的重要性很多人对于为什么要学JVM这个问题,他们的答案都是:因为面试。无论什么级别的Java从业者,JVM都是进阶时必须迈过的坎。不管是工作还是面试中,JVM都是必考题。如果不懂JVM的话,薪酬会非常吃亏。其实学习JVM并…

php ajax xmlhttpreq 上传文件 get,使用Ajax XmlHttpRequest上传文件

使用Ajax XmlHttpRequest上传文件嗨,我正在尝试使用此代码发送带有xmlhttprequest的文件。var url "http://localhost:80/....";$(document).ready(function(){document.getElementById(upload).addEventListener(change, function(e) {var file this.fi…

linux 针对目录空间配额,linux磁盘配额quota

Linux是一个多用户多任务的操作系统,在使用中可能会有几个人对服务器有操作,几个用户共同使用一个共享磁盘的情况,因为我们的硬盘是有限的,我们需要对用户的空间进行限制。这里使用磁盘配额,可以很方便的对用户的空间进…

label居中_表格固定列宽时如何居中?

列宽固定居中的设置的时候,我们通常使用 p{宽度} 来指定固定的列宽,这时单元格会自动换行,换行之后是左对齐的,如何获得居中对齐呢?\begin{tabular}{|p{54pt}l|p{71pt}c|p{71pt}c|}\hline Method& Train set&T…