数值计算基础实验指导书1

1、数值计算基础实验指导书2012年目录实验一直接法解线性方程组的1实验二插值方法10实验三数值积分4实验四常微分方程的数值解6实验五迭代法解线性方程组与非线性方程81、用一-21I 00-00-210000-210实验一直接法解线性方程组、实验目的掌握列选主元消去法与追赶法解线性方程组。二、实验内容分别写出Guass列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。Guass列选主元消去法求解方程组三、实验仪器设备与材料主流微型计算机四、实验原理1、Guass列选主元消去法对于A

2、X=B一1）、消元过程：将（A|B）进行变换为（A|B）,其中是上三角矩阵。即:a11a12anb1”1a12,-anb1”a21a22一.a2nb201-a2nb2-aaTa-aaaa1an1an2.annbnJ<00-annbn/k从1至Un-1a、列选主元选取第k列中绝对值最大元素maxHij作为主元。b、换行akj二a。，j=k1,nbk=bic、归一化akj /akk - akj, j =k 1, ,n b"akk = bkd、消元aj -aikakj - aj ,i = k 1,n; j = k 1, ,nb -aikbk - bi,i = k 1, ,n. .2）

3、、回代过程：由（A|B）解出xn%,，人。bn / ann =' xnnbk. akj xj =,Xk,k =n -1, ,2,1j =k 12、追赶法线性方程组为:ab2、1a2b3c2x2a3C3bn Janbnan做LU分解为:：12分解公式:Yi:1=b1,二 i(i =2,3, ,n)- i-4x3f1f2f3fnfn(i =23 ,n)Cii(i =1,2, ,n-1)Ly=fUx=y回代公式:(i =23,n)(i = n -1,n -2, ,1)、，fi-iyi1yi二二ixn=yn=港=V一Rxi卡五、实验步骤1、理解并掌握列选主元消去法与追赶法；2、画出列选主元消去

4、法与追赶法的流程图3、使用C语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用武汉科技大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序需打印后粘贴在实验报告册内；3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注意事项注意如何定义数据结构以保存矩阵和解以降低算法的复杂性。八、思考题若使用全主元消去法，在编程中应如何记录保存对于未知数的调换。实验二插值方法一、实验目的掌握拉格郎日插值法与牛顿插值法构造插值多项式。二、实验内容分别写出拉格郎日插值法与牛顿插值法的算法，编写程序上机调试出结果，要求所

5、编程序适用于任何一组插值节点，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。已知下列函数表Xi0.561600.562800.564010.56521yi0.827410.826590.825770.82495求x=0.5635时的函数值。三、实验仪器设备与材料主流微型计算机四、实验原理已知n个插值节点的函数值，则可由拉格郎日插值公式与牛顿插值公式构造出插值多项式,从而由该插值多项式求出所要求点的函数值。拉格郎日插值公式与牛顿插值公式如下：1、Lagrange插值公式nLn(x)=lo(x)y°li(x)y1.ln(x)yny*(x)kz01k(x)=(x

6、-x0)(x -x1) (x-xy)(x-xk i) (x-xn)(xk -xO)(xk -xi) (xk -xk)(xk -xk i) (xk -xn)n=nj=sj：*x fjxk -xj2、Newton插值公式Nn(x)=f(x°)fx°,xi(x-x0)fx°,xi，x2(x-x0)(x-x1)fx0,x1,xn(x-x0)(x-Xi)(x-xn4)五、实验步骤1、理解并掌握拉格郎日插值法与牛顿插值法的公式;2、画出拉格郎日插值法与牛顿插值法算法的流程图;3、使用C语言编写出相应的程序并调试验证通过。六、实验报告要求1、统一使用武汉科技大学实验报告本书写，

7、实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序需打印后粘贴在实验报告册内；3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注意事项Newton插值法在编程时应注意定义何种数据结构以保存差商。八、思考题比较Lagrange插值法与Newton插值法的异同。实验三数值积分一、实验目的掌握复化梯形法与龙贝格法计算定积分。二、实验内容分别写出变步长梯形法与Romberge法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。,、1si

8、nx求dx,名<0.00001。0x三、实验仪器设备与材料主流微型计算机四、实验原理通过变步长梯形法与龙贝格法，我们只要知道已知n个求积节点的函数值，则可由相应的公式求出该函数的积分值，从而不需要求该函数的原函数。变步长梯形法与龙贝格法公式如下:1、变步长梯形法nJ-hTn八-f(Xi)f(Xii)i=02h门=-f(a)2、f(x)f(b)2i11 hnT2n=-Tn一二.f(xii/2)2 2»用T2n-Tn<名来控制精度2、龙贝格法梯形法则二阶公式四阶公式六阶公式八阶公式T0(0)T°Ti(0)T0(2)Ti(1)T2(0)T0(3)Ti(2)T2(1)T

9、3(0)T°(4)1(3)T2(2)T3(1)T4(0)24681nO(h2)O(h4)O(h6)O(h8)O(h).4mTm(k)1)Tm(kI)7mL4-1用Tm(0)-TmA(1)W8来控制精度五、实验步骤1、理解并掌握变步长梯形法与龙贝格法的公式；2、画出变步长梯形法与龙贝格法的流程图3、使用C语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用武汉科技大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序需打印后粘贴在实验报告册内；3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注

10、意事项1sinx在sdx积分中，被积函数在x=0点函数值为1,对该点在程序设计中应注意对其的0x定义。八、思考题使用复化梯形法与复化Simpson法来计算该问题有何缺点？实验四常微分方程的数值解一、实验目的掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题。二、实验内容分别写出改进欧拉法与四阶龙格-库塔求解的算法，编写程序上机调试出结果，要求所编程序适用于任何一阶常微分方程的数值解问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。，一2求1y-xy步长h=0.25。y(0)=2(0<x<5)三、实验仪器设备与材料主流微型计算机四、实验原理常微分方

11、程的数值解主要采用“步进式”，即求解过程顺着节点排列次序一步一步向前推进,在单步法中改进欧拉法和四阶龙格-库塔法公式如下：1、改进欧拉法Vn1=ynhf(%,yn)hVn1=Vnf(Xn,Vn)f(Xn1，Vn1)22、四阶龙格-库塔法.hVn+=yn+(kl+2k2+2k3+k4)6kl=f(Xn,yn)一hh4k2=f(Xn+2,yn+-k1)k3=f(Xn+2，yn+hk2)k4-f(Xnh,ynhk3)五、实验步骤1、理解并掌握改进欧拉法与四阶龙格-库塔法的公式；2、画出改进欧拉法与四阶龙格-库塔法的流程图3、使用C语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用武汉科技

12、大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序需打印后粘贴在实验报告册内;3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注意事项=2-xyy(0) = 2的精确解为y = 2/(1 + x2),通过调整步长，观察结果的精度(0 < x < 5)的变化八、思考题如何对四阶龙格-库塔法进行改进，以保证结果的精度。实验五迭代法解线性方程组与非线性方程一、实验目的掌握高斯-塞德尔迭代法求解线性方程组与牛顿迭代法求方程根。二、实验内容分别写出高斯-塞德尔迭代法与牛顿迭代法的算法，编写程序上机调试出结

13、果，要求所编程序适用于任何一个方程的求根，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。1、高斯-塞德尔迭代法求解线性方程组7219153-2 -2 11J 32一2x1-2 x25 X313 X4一4171-1-0 132、用牛顿迭代法求方程x-x-1=0的近似根，名0.00001,牛顿法的初始值为1。三、实验仪器设备与材料主流微型计算机四、实验原理二分法通过将含根区间逐步二分，从而将根的区间缩小到容许误差范围。牛顿通过迭代的方法逐步趋进于精确解，该两种方法的公式如下：1、高斯-塞德尔迭代法1)判断线性方程组是否主对角占优nZau-aii,i1,2,nj1j-i2)直接分离xi,即nxi=(di-%bijxj)/aii,i=1,2,nj1建立高斯-塞德尔迭代格式为：i1n(k1)(k1)、(k)、x二(di-、a/j-'aijxj)/a",i=1,2,nj1j丑13)取初值迭代求解至所要求的精度为止。2、牛顿法xk1=xkf(xk)f(xk)五、实验步骤1、理解并掌握二分法与牛顿法的公式；2、画出二分法与牛顿法的流程图3、使用C语言编