数值计算基础

1、精品文档精品文档数值计算基础实验指导书2010年目录实验一直接法解线性方程组的1实验二插值方法10实验三数值积分4实验四常微分方程的数值解6实验五迭代法解线性方程组与非线性方程8实验一直接法解线性方程组、实验目的掌握列选主元消去法与追赶法解线性方程组。二、实验内容分别写出Gauss列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。1、用Gauss列选主元消去法求解方程组2.55.32.3 -5.1 Xi9.61.5X2'I1.7-4.31X3-3.713.85.5J2、

2、用追赶法求解方程组-2-2三、实验仪器设备与材料主流微型计算机四、实验原理1、Gauss列选主元消去法（A|B）,其中A是上三角矩阵。即:对于AX=Ba11a12.anb1”1ai2一«anb1”a21a22,a2nb201-a2nb2-aT-m-一-an1an2,annbnj<00-annbn)（A|B）进行变换为1）、消元过程:将k从1至Un-1a、列选主元选取第k列中绝对值最大元素max ak至囱ik作为主元。b、换行akj=a。，j=k1,nbk=bic、归一化akj /akk = akj, j =k 1, ,n bk/akk = bkd、消元aj -aikakj -

3、aj ,i = k 1,n; j = k 1, , nb -aikbk - bi,i = k 1, ,n丁 . .2）、回代过程：由（A|B）解出Xn,Xn,，X1。bn / ann =' Xnnbk. akj Xj =,Xk,k =n -1, ,2,1j =k 12、追赶法线性方程组为:ab2CiXia2b3C2X2a3C3bn JandbnCn J做LU分解为:：12分解公式:：1(i =2,3, ,n)二bi 一 i - jX3Xnjan八Xn(i =2,3, ,n)Ciai(i =1,2, ,n-1)f1f2f3fn< fnLy=fUx=y精品文档回代公式：'_f

4、ly1=«1(i =2,3,n)(i = n -1,n - 2, ,1)fi-iyiJy二二ixn=yn=X=VPiXi书五、实验步骤1、理解并掌握全选主元消去法与高斯-塞德尔迭代法公式；2、画出全选主元消去法与高斯-塞德尔迭代法的流程图3、使用C语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用武汉科技大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序写在实验报告册内；3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注意事项注意如何定义数据结构以保存矩阵和解以降低算法的复杂性。八、思

5、考题若使用全主元消去法，在编程中应如何记录保存对于未知数的调换。精品文档实验二插值方法一、实验目的掌握拉格郎日插值法与牛顿插值法构造插值多项式。二、实验内容分别写出拉格郎日插值法与牛顿插值法的算法，编写程序上机调试出结果，要求所编程序适用于任何一组插值节点，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。已知下列函数表Xi0.561600.562800.564010.56521yi0.827410.826590.825770.82495求x=0.5635时的函数值。三、实验仪器设备与材料主流微型计算机四、实验原理已知n个插值节点的函数值，则可由拉格郎日插值公式与牛顿插值

6、公式构造出插值多项式,从而由该插值多项式求出所要求点的函数值。拉格郎日插值公式与牛顿插值公式如下：1、Lagrange插值公式nLn(x)=lo(x)y°li(x)yi.ln(x)yn=、y*(x)kz01k(x)=(x -x0)(x -x1)(x -x-)(x- xk i) (x -xn)(xk -xo)(xk -xi)(xk -xk)(xk -xk 1) (xk -xn)n=nj=sj=kx fjxk - xj2、Newton插值公式Nn(x)-f(xo)fxo,xi(x-xo)fxo,xi,x2(x-xo)(x-xi)fx0,x1,xn(xx°)(xXi)(xxn)五

7、、实验步骤1、理解并掌握拉格郎日插值法与牛顿插值法的公式;2、画出拉格郎日插值法与牛顿插值法算法的流程图;3、使用C语言编写出相应的程序并调试验证通过。六、实验报告要求1、统一使用武汉科技大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。精品文档2、源程序需打印后粘贴在实验报告册内；3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。七、实验注意事项Newton插值法在编程时应注意定义何种数据结构以保存差商。八、思考题比较Lagrange插值法与Newton插值法的异同。精品文档实验三数值积分、实验目的掌握梯形复合求积法与Rom

8、berg法计算定积分。二、实验内容分别写出变步长梯形复合求积法法与Romberg法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。1sinx,0dx,；m0.00001。三、实验仪器设备与材料主流微型计算机四、实验原理通过变步长梯形法与龙贝格法，我们只要知道已知n个求积节点的函数值，则可由相应的精品文档公式求出该函数的积分值，从而不需要求该函数的原函数。变步长梯形法与龙贝格法公式如下:1、变步长梯形法nJ-hTn八-f(Xi)-f(Xi1)y2h,nIf2、f(Xi)f(b)2y1hnT2n=

9、-Tn-f(Xii/2)22if用 T2n n名来控制精度2、龙贝格法T2n =Tn2h n二 ' f (Xi 1/2)2 ifSnCn=4T2n315- 3Tn-SRn-64c一 63c2n15I63Cn用R2n-Rn|-名来控制精度五、实验步骤1、理解并掌握变步长梯形法与龙贝格法的公式；2、画出变步长梯形法与龙贝格法的流程图3、使用C语言编写出相应的程序并调试验证通过六、实验报告要求1、统一使用武汉科技大学实验报告本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。2、源程序需打印后粘贴在实验报告册内；3、运行结果以屏幕截图形式保存并打

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

11、器设备与材料主流微型计算机四、实验原理常微分方程的数值解主要采用“步进式”，即求解过程顺着节点排列次序一步一步向前推进,在单步法中改进欧拉法和四阶龙格-库塔法公式如下：1、改进欧拉法Vn1=ynhf(Xn，yn)hyn1=ynf(Xn,Yn)f(Xn1,丫口1)22、四阶龙格-库塔法.hyn+=yn+(kl+2k2+2k3+k4)6kl=f(Xn,yn)心2 = f (XnK1h - 2十ny,h - 2十卜3=f(Xn+2，yn+”2)k4=f(Xnh,ynhk3)五、实验步骤1、理解并掌握改进欧拉法与四阶龙格-库塔法的公式;2、画出改进欧拉法与四阶龙格-库塔法的流程图3、使用C语言编写出相

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

13、高斯-塞德尔迭代法与牛顿迭代法的算法，编写程序上机调试出结果，要求所编程序适用于任何一个方程的求根，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。1、高斯-塞德尔迭代法求解线性方程组一 7219153-2-211J32-2区-2X25X313-/41一4171-1R 12、用牛顿迭代法求方程X3-X-1=0的近似根，名E0.00001,牛顿法的初始值为1。三、实验仪器设备与材料主流微型计算机四、实验原理二分法通过将含根区间逐步二分，从而将根的区间缩小到容许误差范围。牛顿通过迭代的方法逐步趋进于精确解，该两种方法的公式如下：1、高斯-塞德尔迭代法1)判断线性方程组是否主对角占优naiaijaii,i1,2,nj1j-i2)直接分离xi,即nXi=(di-%bijXj)/aii,i=1,2,nj1建立高斯-塞德尔迭代格式为：i1n(k1)(k1)v(k)、Xi=(di-ajXj-aijXj)/aH,i=1,2,nj1j丑13)取初值迭代求解至所要求的精度为止。2、牛顿法f(Xk)f(Xk)五、实验步骤1、理解并掌握二分法与牛顿法的公式；2、画出二分法与牛顿法的流程图3、使用C语言