数值计算基础实验指导 部分实验源代码 复习指导 三套试题及其答案.docVIP

数值计算基础实验指导书2010 年目录实验一 直接法解线性方程组的 3 实验二 插值方法 12 实验三 数值积分 6 实验四 常微分方程的数值解 8 实验五 迭代法解线性方程组与非线性方程 10实验一 直接法解线性方程组一、实验目的掌握全选主元消去法与高斯-塞德尔法解线性方程组。二、实验内容分别写出 Guass 列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编 程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以 下列数据验证程序的正确性。1、用 Guass 列选主元消去法求解方程组?2.52.3? 5.1? ? x1 ??3.7??5.39.61.5 ? ? x ? ? ?3.8????8.11.7? 4.3? ??? ??2 ?x3 ???? ???5.5??2、用追赶法求解方程组?? 2 0 0 00 ? ? x1 ??? 10?? 1 ? 2 0 00 ? ? x ?? 0 ?? ? ? 2 ? ? ?? 0 1? 2 00 ? ? x3 ? ? ? 0 ?? ? ? ? ? ?? 0 01 ? 20 ? ? x4 ?? 0 ?0 0 01 ? 2 x5 0三、实验仪器设备与材料主流微型计算机四、实验原理1、Guass 列选主元消去法对于 AX =B~A B 是上三角矩阵。即：1）、消元过程：将（A|B）进行变换为 ( ~ | ~) ，其中 A? a a ? ab ?? ? 1 a ? a b ?? 11 12 1n 1 ? ? 12 1n 1 ?? a21 a22 ? a2 nb2 ?? 0 1? a2 nb2 ?? ? ?? ? ? ? ? ? ?? ? ????? an1an 2?b?? ann n ?b???? 0 0??b?? ann n ?b?k 从 1 到 n-1 a、 列选主元选取第 k 列中绝对值最大元素 max aik k ?i?n作为主元。b、 换行akjbkc、 归一化? aij , j ? k ? 1,?, n? bid、 消元akj / akkbk / akk? akj , j ? k ? 1,?, n? bkaij ? aik akj ? aij , i ? k ? 1,?, n; j ? k ? 1,?, n bi ? aik bk ? bi , i ? k ? 1,?, n2）、回代过程：由 ( ~ | ~) 解出 x , x,?, x 。A Bbn / ann ? xn nn n?1 1bk ?? akj x j ? xk , k ? n ? 1,?,2,1j ?k ?12、追赶法线性方程组为：? a c?? x ?? f ?? 1 1?? 1 ?? 1 ?? b2??a2 c2?b3 a3 c3??? x2 ?3?? x ?3?? ?? f 2 ?3? f ?3? ??? ? ? ???? ?? ? ? ? ???? bn ?1??an ?1cn?1???????? x???????n?1?? ??? f???? f??n?1?? bnan ?? xn ?? f n ?做 LU 分解为：??? ? 1 ? ?? 1 ? ? 1 ?? ? 2 ? 2?? ? ??? ? 1 ? 2 ?? ? ? ? ?L ? ? 3 3?, R ? ? ??? ? ? ? ????? ???? ? ? ?? ? ?? ? ??1 n?1 ??? ? n? n ? ? 1 ?分解公式：??? ? a?? i i(i ? 2,3,?, n)??1 ? b1 ,? i ? bi ? ? i ? i ?1?(i ? 2,3,?, n)i?? ? cii?? ? i(i ? 1,2,?, n ? 1)则Ax ? f? LUx ? f?Ly ? f? ??Ux ? y回代公式：y?? f1y?? 1??y?y?? i?1? f i ? ? i yi ?1? i(i ? 2,3,?, n)?xn ? yn??xi ? yi ? ? i xi ?1(i ? n ? 1, n ? 2,?,1)五、实验步骤1、理解并掌握全选主元消去法与高斯-塞德尔迭代法公式；2、画出全选主元消去法与高斯-塞德尔迭代法的流程图3、使用 C 语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用《武汉科技大学实验报告》本书写，实验报告