数值分析上机试题作业.doc

PAGE PAGE 0 数值分析上机实习报告 姓名： ****** 学号： ******* 专业： 电话： 序 言 所用程序语言：本次数值分析上机实习采用Visual Basic.Net和Matlab作为程序设计语言，利用Visual Basic.Net可视化的编程实现方法，采用对话框形式进行设计计算程序界面，并将结果用表格或文档的格式给出。 程序概述： 第一题中在Matlab中用雅格比法与高斯－赛德尔迭代法解方程组，并研究二者的收敛性，上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛影响度。 第二题使用Visual Basic.NET程序设计语言完成了“松弛因子对SOR法收敛速度的影响” ，通过在可视化界面下输入不同的n和w值，点击按钮直接可看到迭代次数及计算结果，观察了不同的松弛因子w对收敛速度的影响。 目 录 TOC \o 1-3 \h \z \u 一．雅格比法与高斯—塞德尔迭代法解方程组 3 1. 计算结果 3 1）雅格比法计算结果 3 2）高斯—塞德尔迭代法计算结果 4 2. 结果分析 5 3. 程序清单 5 1）雅格比法 5 2）高斯—塞德尔迭代法 6 二．松弛因子对SOR法收敛速度的影响 8 1. 迭代次数计算结果 8 1）n=10时迭代次数 8 1）n=20时迭代次数 9 2. 计算X()结果 10 1）n=10时X()计算结果 10 2）n=20时X()计算结果 11 3. 对比分析 12 4. 程序清单： 12 三．实习总结 16 实验课题（一）用雅格比法与高斯－赛德尔迭代法解方程组 一、用雅格比法与高斯－赛德尔迭代法解下列方程组Ax=b，研究其收敛性，上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛有无影响。 (1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T, b2=[100,-200,345]T, (2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1] T, b2=[5,0,-10]T, (3)A行分别为A1=[1,3],A2=[-7,1]；b=[4,6]T。 1、计算结果 1）雅格比法计算结果： 当A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]时： 当b1=[-3,2,4]T 时 结果为： -0.7282 0.8076 0.2548 迭代次数： 9 (2)当b2=[100,-200,345]T 时 结果为： 36.3619 -2.0737 114.0417 迭代次数： 15 当A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]时 (1)当b1=[3,2,1] T时 经程序计算（变换N值），雅克比法发散！ (2)当b2=[5,0,-10]T时 经试算（变换N值），雅克比法发散！ (3) 结果为： 经试算（变换N值），雅克比法发散！ 2）高斯－赛德尔迭代法计算结果 当A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]时： (1) 当b1=[-3,2,4]T 时 结果为： -0.7254 0.8086 0.2538 迭代次数： 6 (2) 当b2=[100,-200,345]T 时 结果为： 36.3637 -2.0751 114.0415 迭代次数： 11 当A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]时 (1) 当b1=[3,2,1] T时 结果为： 5.7854 0.7704 -4.2446 迭代次数： 17 (2) 当b2=[5,0,-10]T时 结果为： 32.6912 7.7129 -42.3233 迭代次数： 23 (3) 当b=[4,6] T时 经试算（变换N值），高斯法发散！ 2、结果分析 通过对雅克比法和高斯法的上机编程实习，分析对比实验结果可得：在方程组Ax=b中，右端