数值分析上机题3

上传人:ni****g 文档编号:447911142 上传时间:2024-01-21 格式:DOC 页数:14 大小:286.50KB
返回 下载 相关 举报
数值分析上机题3_第1页
第1页 / 共14页
数值分析上机题3_第2页
第2页 / 共14页
数值分析上机题3_第3页
第3页 / 共14页
数值分析上机题3_第4页
第4页 / 共14页
数值分析上机题3_第5页
第5页 / 共14页
点击查看更多>>
资源描述

《数值分析上机题3》由会员分享,可在线阅读,更多相关《数值分析上机题3(14页珍藏版)》请在金锄头文库上搜索。

1、数值分析上机题目3实验一1. 根据Matlab语言特点,描述 Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法。2. 编写Jacobi迭代法、Gauss-Seidel迭代法和 SOR迭代法的 M文件。3给定AR20 20为五对角矩阵312412321412314121214312141231214123(1)选取不同的初始向量 x(0)及右端面项向量b,给定迭代误差要求,分别用编写的 Jacobi迭代法和Gauss-Seidel迭代法程序求解,观察得到的序列是否收敛?若收敛,通过迭代次数分析计算结果并得出你的结论。用编写的SOR迭代法程序,对于(1)所选取的初始向量 x(0)及

2、右端面项向量b进行求解,松驰系数3取13 10A(-7 )&k10A(-7 )&k10A(-7 )&k10A(-5) & k10A(-2) & k100000k=k+1;g=y1-y0;y=x1-x0;A1=A0+(g-A0*y)/(y*y)*y;x0=x1;x1=x0-A1y1;A0=A1;y0=f(x0);y1=f(x1);endx=x1k=k运行结果,0.3973-0.3592-0.1240 k =53. 尝试不同初值。x0=0.1,0.1,0.1采用Newt on法0.2459-0.6348-0.084868Broyden 秩 1 法0.0213-0.4117-0.2104 k =7该结果说明,采用不同的初值,对结果没有影响。4、计算两种方法的用时采用Newt on法t 1=0.7623Broyden 秩 1 法t2 =0.2444采用 Broyden 秩 1 法运算比较快。

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 办公文档 > 活动策划

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号