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②高斯赛德尔迭代法:建立函数文件。
随机阅读
- 从上沙到彩田南路坐什么车 下沙到彩田路海鹰大厦
- cet4各题型的题量及分值? cet4每小题分值
- 银川大兴镇卫生院地址 被狗咬了,在哪里可以打疫苗
- 求大神写一篇去凤凰古城的攻略 ? 天门剧院《水乡情》
- 小学阶段的最后一次运动会为题 一年一度的运动会即将召开,这也是我们初中生活中最后一次运动会,对此,你有什么感想,请写一篇800字的文章.(文体不限,诗歌、戏剧除外)
- 走马观花进口博览会 去威海旅游哪个季节比较合适?
- 昆明汇仁堂是药店吗 谁能提供国企医药公司排名~谢谢~~
- 气动阀无法关闭怎么处理 博雷气动头更换密封圈
- 他卡西醇乳膏作用 他卡西醇软膏作用
- 朋友的记忆七年级叙事作文 关于沉着冷静的作文叙事
- 教师有那些心理健康问题 教师有哪些心理健康的标准?
- 福建樱之缘生态旅游开发有限公司怎么样? 绍武市大竹镇大竹村村务
- iris火焰检测控制器 什么是显卡?作用是什么?
- 谁能发我铜雀台高清资源(720P或者1080P)麻烦发邮箱? 铜雀台正片资源
- 万能试验机的夹头能夹非圆柱型拉伸试样吗?比方说板状拉伸试样,还是得换专用夹头?机器原配的夹头可以吗? 标准板状拉伸试验
- 为什么怀孕期间容易得胆囊炎? 女性激素影响胆囊的是
- 幼儿园安吉玩具教案 幼儿园物质设备配备的分析与思考
- 内心很纠结,该怎么办? 纠 四风 不能止步 作风建设永远在路上
- 哈尔滨新区 航空气动院 请问你认识哈尔滨气动院的学生是哈工程的么?我本科和研究生都是哈工程的,想去气动院,我是机械电子的,
- 强调理性认识依赖于感性认识,这是因为 强调理性认识 否认感性认识是