数值计算基础.docx

数值运算基础 试验指导书 2021 年 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 目录 试验一 试验一 直接法解线性方程组的 ................................ 1 试验二 插值方法 ........................................... 10 试验三 数值积分 ............................................ 4 试验四 常微分方程的数值解 .................................. 6 试验五 迭代法解线性方程组与非线性方程 ...................... 8 |品. |可. |编. |学. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 试验一 直接法解线性方程组 一、试验目的 |精. |品. |可. |编. 把握列选主元消去法与追逐法解线性方程组； 二、试验内容 分别写出 Gauss 列选主元消去法与追逐法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题；试验中以以下 数据验证程序的正确性； 1、用 Gauss列选主元消去法求解方程组 |辑. |学. |习. |资. |料. * | * | * 2.5 5.3 8.1 2.3 9.6 1.7 5.1 x1 1.5 x2 4.3 x3 3.7 3.8 5.5 | * | |欢. |迎. |下. |载. 2、用追逐法求解方程组 20 2 0 0 0 0 x1 10 1 2 0 0 0 x 0 0 1 2 0 0 x3 0 0 0 1 2 0 x4 0 0 0 0 1 2 x5 0 三、试验仪器设备与材料 主流微型运算机 四、试验原理 ~ ~中是上a11a12a1nb1 ~ ~ 中是上 a11 a12 a1n b1 1 a12 a1n b1 a21 a 22 a2 n b2 0 1 a 2n b2 1）、消元过程：将 （A|B ）进行变换为  ( A | B) ，其  三角矩阵；即： an 1 an 2 ann bn 0 0 ann bn k 从 1 到 n-1 a、 列选主元 选取第 k 列中肯定值最大元素  max  aik  作为主元； b、 换行  akj bk  aij , j bi  k 1, k i n , n akj akj / akk bk / akk akj , j bk k 1, ,n d、 消元 aij bi aik akj aik bk aij , i bi , i k 1, 1, ,n; j ,n k 1, , n k 2）、回代过程：由 ( A | B) 解出 xn , xn 1, ~ ~ , x1； bn / ann n bk xn akj x j xk , k n 1, ,2,1 j k 1 2、追逐法 线性方程组为： a1 b2 c1 a2 b3 c2 a3 c3 x1 x2 x3 f 1 f 2 f 3 bn 1 an bn 1 cn 1 an xn 1 xn f n 1 f n 做 LU 分解为： 1 1 1 2 2 1 2 L 3 3 , R 1 n n n 1 1 分解公式： i 1 ai b1, ci i (i 2,3, , n) i bi i i 1 (i 2,3, , n) i (i 1,2, , n 1) 就 Ax f LUx f Ly Ux f y |精. |可. |可. |编. |辑. |辑. |学. |习. |资. *|料. * | * * | * *| * | |欢. |欢. |迎. |迎. |下. |载. |载. 回代公式： y1 y 1 f1 1 y f i i yi 1 i (i 2,3, , n) i xn xi yn yi i xi 1 (i n 1, n 2, 五、试验步骤 |精. |品. |可. |编. |辑. |学. |习. |资. |料. * | * | * | * | |欢. |迎. |下. |载. 1、懂得并把握全选主元消去法与高斯 -塞德尔迭代法公式； 2、画出全选主元消去法与高斯 -塞德尔迭代法的流程图 3、使用 C 语言编写出相应的程序并调试验证通过 六、试验报告要求 1、统一使用《武汉科技高校试验报告》本书写，试验报告的内容要求有：试验目的、试验内容、程序流程图、源程序、运行结果及试验小结六个部分； 2、源程序写在试验报告册内； 3、运行结果以屏幕截图形式储存并打印后粘贴在试验报告册内； 七、试验留意