数值分析与实验

1/12数值分析与实验论文关键词列主元高斯消去法雅可比法高斯赛德尔迭代法幂法论文摘要本文通过实例对线性方程组数值解法和矩阵的特征值及特向量的计算进行了探讨。在对线性方程组数值解法的讨论下用到了列主元高斯消去法、雅可比法和高斯赛德尔迭代法。正是高斯消去法在消元时存在一些必须的条件，才启发我们通过列主元高斯消去法来对线性方程组数值解法作进一步的研究，达到了很好的的效果。同时用雅可比法和高斯赛德尔迭代法对相类似的问题的探讨来比较它们的优劣，使我们在分析问题时能更好的把握方法。在求矩阵按模最大的特征值及对应特征向量时，本文用到了幂法，可以使现实中很多复杂的计算简单。第一章线性方程组数值解法实验目的熟悉求解线性方程组的有关理论和方法；会编制列主元消去法，雅可比及高斯赛德尔迭代法的程序；通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。实验内容列主元高斯消去法求解线形方程组；2/12雅可比法和高斯赛德尔迭代法解方程组；题目列主元高斯消去法求解线形方程组方程组为列主元高斯消去法算法将方程用增广矩阵表示1）消元过程对K1，2，,N11选主元，找使得如果则矩阵A奇异，程序结束；否则执行3如果则交换第K行与第行对应元素位置，JK,N1消元，对IK1,N计算对JK1,N1计算）回代过程1若则矩阵A奇异，程序结束；否则执行2；对IN1,2,1计算INCLUDEINCLUDEVOIDCOLPIVOTFLOATC,INTN,FLOATXINTI,J,T,KFLOATPFORI0IKI3/12FORJI1JIFFABSCJN1IFABSCKN1IKJIFKIFORJIJPCIN1JCIN1JCKN1JCKN1JPFORJI1JPCJN1I/CIN1IFORTITFORIN1I0IFORJN1JI1JCIN1NXJCIN1JXICIN1N/CIN1IVOIDMAIN4/12VOIDCOLPIVOTFLOAT,INT,FLOATINTIFLOATX4FLOATC451,1,2,1,8,2,2,3,3,20,1,1,1,0,2,1,1,4,3,4,COLPIVOTC0,4,XFORI0I输出结果结果分析从输出结果可以得到，从结果和过程可以知道这种方法一般能保证舍入误差不扩散，这个方法基本上是稳定的。LOCALHOST1题目雅可比法解方程组方程组为雅可比迭代法算法设方程组AXB的系数矩阵的对角线元素，M为迭代次数容许的最大值为容许误差。1取初始向量令K0对I1,25/12,N计算如果则输出结果；否则执行4如果则不收敛，终止程序；否则，转2程序INCLUDEINCLUDEDEFINEEPS1E6DEFINEMAX100VOIDJACOBIFLOATA,INTN,FLOATXINTI,J,K0FLOATEPSILON,SFLOATYNEWFLOATNFORI0IWHILE1EPSILON0KFORI0IS0FORJ0JIFJICONTINUE6/12SAIN1JXJYIAIN1NS/AIN1IEPSILONFABSYIXIFORI0IIFEPSILONPRINTF“DIEDAICISHUWEIDN“,KRETURNIFKMAXPRINTF“DIEDAIFASAN“RETURNDELETEYVOIDMAINSINTIFLOATA4510,1,2,0,11,0,8,1,3,11,2,1,10,0,6,1,3,1,11,25FLOATX4JACOBIA0,4,XFORI0I输出结果7/12结果分析迭代次数增加时，精度越高。从输出结果可以看出此方程组的迭代次数为17，迭代结果越来越接近精确解了，于是,1题目高斯赛德尔迭代法解方程组方程组为高斯赛德尔迭代法算法设方程组AXB的系数矩阵的对角线元素，M为迭代次数容许的最大值为容许误差。1取初始向量令K0对I1,2,N计算如果则输出结束；否则执行4如果则不收敛，终止程序；否则，转2程序INCLUDEINCLUDEDEFINEN00VOIDMAININTIFLOATX48/12FLOATC4510,1,2,0,11,0,8,1,3,11,2,1,10,0,6,1,3,1,11,25VOIDGAUSSSEIDELFLOAT,INT,FLOATGAUSSSEIDELC0,4,XFORI0IVOIDGAUSSSEIDELFLOATA,INTN,FLOATXINTI,J,K1FLOATD,DX,EPSFORI0IWHILE1EPS0FORI0ID0FORJ0JIFJICONTINUEDAIN1JXJDXAIN1ND/AIN1IEPSFABSDXXIXIDXIFEPSPRINTF“迭代次数是DN“,K9/12RETURNIFKNPRINTF“迭代发散NN“RETURNK输出结果结果分析从输出结果可以看出此方程组的迭代次数为7，此时能得到精确结果是,107374176从结果和原有知识可以知道其系数矩阵是严格对角占优的。所以此迭代解法有很好的收敛性。1方法比较雅可比法和高斯赛德尔迭代法解方程组两种方法的比较。由于此题的系数矩阵是严格对角占优的，所以雅克比迭代法和高斯赛德尔迭代法都是收敛的，这两种迭代法没迭代一步均是作一次矩阵和向量的乘法，但前者需要10/122组工作单元分别存放和，而后者只需要1组工作单元。对于同一个线性方程组，这两种方法可能同时收敛，也可能同时发散，也可能其一收敛，而另一发散。但当两者皆收敛时，一般来说高斯赛德尔迭代法比雅克比迭代法收敛快。实际中更多的是使用逐次超松弛迭代法。第二章矩阵的特征值及特征向量的计算实验目的在数学和物理中，很多问题都需要计算矩阵的特征值及特征向量，它们是线性代数中的一个重要课题，而在实际问题中，这样的计算是很复杂的，有的要求矩阵按模最大特征值及相应的特征向量，有些则要求全部特征值及特征向量，根据不同的要求计算方法大体上可分为2种类型。本实验用的是幂法求矩阵按模最大的特征值及对应特征向量，要求领会求矩阵特征值及特征向量的幂法的方法，并要求会编制幂法的计算程序，来计算有关问题。实验内容利用幂法求矩阵按模最大的特征值及对应特征向量。1幂法求矩阵按模最大的特征值及对应特征向量用幂法求矩阵按模最大的特征值及其相应的特征向量，使，11/12幂法算法幂法是求矩阵主特征值的一种迭代方法。设有N个线性无关的特征向量，而相应的特征值满足，则对任意非零初始向量按下述公式构造向量序列其中表示中最大的分量，并且有，。用幂法计算实对称矩阵的特征值时，可用RAYLEIGH商作加速。设的RAYLEIGH商为则，当时，将比更快趋于。1程序1输出结果结果分析主特征值为；相应的特征向量为。幂法是求矩阵主特征值的一种有效方法，特别当矩阵为大型稀疏时，更显得如此。但由于特征值的分布无法事先预测，因此不能控制收敛速度，往往需要利用某些加速技巧。所以计算时我们要根据需要选择计算方法来计算矩阵的特征值及特征向量。12/12参考文献1袁