sor迭代算法的matlab程序 function[x k]=EqtsSOR(A,b,x0,omiga,eps)%超松弛(SOR,Successive Over-Relaxation)迭代法求解线性方程组Ax=b%[x k]=EqtsSOR(A,b,x0,eps)%x:解向量,列向量%k:迭代次数%A。什么是SOR? SOR是英文Standard Operation Records的缩写,它的中文含意是记录、凭证类文件。SOR包括了:⑴物料管理记录;⑵生产操作记录;⑶质量检验与管理记录;⑷设备运行与管理记录。关于数值分析超松弛(SOR)法的公式问题 因为 A=D-L-U,所以 L,U 里面的元素相当于 A 里元素加负号,即相当于-A 里的元素。用Matlab写的雅各比i和高斯塞德尔以及SOR迭代法 1.用雅克比迭代法和高斯-赛德尔迭代法求解下列方程组,取迭代初值[0;0;0]。(1)编程求解,并与用数学软件求解的结果对比。(2)考察迭代法的收敛性,若均收敛,对比两种方法的收敛速度。解:源程序:①雅克比迭代法:建立函数文件jacobi.mfunction[n,x]=jacobi(A,b,X,nm,w)用雅克比迭代法求解方程组Ax=b输入:A为方程组的系数矩阵,b为方程组右端的列向量,X为迭代初值构成的列向量,nm为最大迭代次数,w为误差精度输出:x为求得的方程组的解构成的列向量,n为迭代次数n=1;m=length(A);D=diag(diag(A));令A=D-L-U,计算矩阵DL=tril(-A)+D;令A=D-L-U,计算矩阵LU=triu(-A)+D;令A=D-L-U,计算矩阵UM=inv(D)*(L+U);计算迭代矩阵g=inv(D)*b;计算迭代格式中的常数项下面是迭代过程while nx=M*X+g;用迭代格式进行迭代if norm(x-X,2)disp('迭代次数为');ndisp('方程组的解为');xreturn;上面:达到精度要求就结束程序,输出迭代次数和方程组的解endX=x;n=n+1;end下面:如果达到最大迭代次数仍不收敛,输出警告语句及迭代的最终结果(并不是方程组的解)disp('在最大迭代次数内不收敛。');disp('最大迭代次数后的结果为');x②高斯赛德尔迭代法:建立函数文件。
随机阅读
- 金华万通大厦会展部门 有哪些有趣的大数据?
- 越来越越来越造句 用越来越什么越来越什么造句
- 建筑施工安全检查中关于基坑支护的五个保证项目是什么? 基槽底有水井
- 迷你世界怎么刷皮肤 迷你世界我想要皮肤
- 从长春桥东到朝阳区大柳树怎么坐公交车,最快需要多久 长春桥站到朝阳门站
- 骑士精神盾 骑士精神是什么意思?
- 山西冶金技师学院我从太原火车站坐几路车去 太原冶金技师学院在吗
- 理想的筛检试验对人体无害 添加抗结剂的食盐对人体有害吗?
- 如何去无锡永丰路212号? 无锡53路永丰路
- 国内有哪些适合一个人去旅行的地方? 浙江平阳县竹林园机械小吃车
- 显存与内存的区别在哪里? tsv芯片封装有哪几种
- 尿妊辰试验金标法原理 请问尿妊娠试验(金标法)是查什么的呢
- 哈尔滨民南街 哈尔滨桥南街二十号在哪
- 如何看待电视剧《杀破狼》官宣檀健次饰演顾昀、陈哲远饰演长庚? 十万买断母子情
- 通州区永乐店镇什么时候拆迁 规划 通州永乐店拆大务村
- 娶个小十岁的妻子 娶一个小自己十几岁的妻子,是种什么体验?一个39岁男人说了实话
- 昌平20路公交换乘879路公交车从哪里下车 南农家园到304医院南门
- 家庭趣事作文:家庭“战争” 作文家庭战争600字
- 作文教学目标和三维目标一样吗 教学三维目标是什么?三个维度的详细解释,确定三维目标的原则
- 河南睢县在什么地方? 睢县在那个地方