《计算方法实验》由会员分享，可在线阅读，更多相关《计算方法实验（33页珍藏版）》请在金锄头文库上搜索。
1、计算方法实验指导姓名学号院系 专业哈尔滨工业大学计算方法实验指导根据实际问题建立的数学模型，一般不能求出所谓的解析解，必须针对数学模型的特 点确定适当的计算方法，编制出计算机能够执行的计算程序，输入计算机，进行调试，完 成运算，如果计算结果存在问题或不知是否正确，还需要重新确定新的计算方法，再编制 出计算程序，输入计算机，重新调试，完成运算，直至获得正确的计算结果，这就是数值 计算的全部过程。学生在学习“计算方法”和“高级语言”等课程时普遍存在的问题是：只会套用教科 书中的标准程序进行数值计算， 很少有人能够独立地将学过的数值算法编制成计算机程序， 至于灵活应用已经掌握的算法求解综合性较大的课
2、题，则更是困难的事情。编写计算方法实验指导的目的是：突出数值计算程序结构化的思想。提高学生的 编程能力，加深对“计算方法”课程内容的理解和掌握，为”计算方法“课程的教学服务， 进一步奠定从事数值计算工作的基础。具体地1. 根据“计算方法”课程内容的特点，给出五个典型算法的分析流程，学生可以利用 所掌握的“高级语言”顺利地编制出计算机程序，上机实习，完成实验环节的教学要求。2. 所有的计算实习题目都经过任课教师逐一检验，准确无误。3. 充分利用循环的思想、迭代的思想，给出算法结构描述和程序语言的对应关系，有 利于学生编制相应的程序。4. 结合实习题目，提出实验要求，要求学生按规范格式写出相应的实
3、验报告，实验报 告成绩记入期末总成绩。需要提醒学生：不能简单地套用现成的标准程序完成实验题目， 应当把重点放在对算法的理解、程序的优化设计、上机调试和计算结果分析上，否则就失 去实验课的目的啦。5. 五个具体的实验题目是：实验题目 1 拉格朗日 (Lagrange) 插值实验题目2龙贝格(Romberg)积分法实验题目 3 四阶龙格库塔 (RungeKutta) 方法实验题目 4 牛顿 (Newton) 迭代法实验题目5高斯(Gauss)列主元消去法要求必须完成其中三个(如果全部完成更好) 。实验题目 1 拉格朗日 (Lagrange) 插值方法概要：j ;则满足条件a,b时，有误差给定平面上
4、n 1个不同的数据点(xk,f(xk) , k 0,1,L ,n ,人Xj , iPn(Xk)f (Xk) , k 0,1 丄,n的n次拉格朗日插值多项式是存在唯一的。若xk a, b, k 0,1,L , n，且函数f (x)充分光滑，则当x估计式f(n 1()f (x) Fn(x)(x xo)(x xJL (x Xn)，a,b(n 1)!拉格朗日插值算法实验实验目的：利用拉格朗日插值多项式Fn(x)求f (x)的近似值输 入：n 1个数据点(x-f&k)， k 0,1,L , n ;插值点x输 出：f (x)在插值点x的近似值巳(x)程序流程：1 置 y 0.0 ； k 02 当 k n
5、时，做 2.1 2.42.1 置 I 1.0 ；2.2 对 j 0,1,L ,k 1,k 1,L ,n，置 l l(x 百)/仪 Xj)2.3 置 y y l f (xj2.4 置 k k 13输出x, y4停机问题1拉格朗日插值多项式的次数n越大越好吗？考虑下面两个拉格朗日插值问题：1(1) 设f(x) 二，x 5,5，考虑等距节点的拉格朗日插值多项式巳(x),即将区1 x间5,5进行n等分，记h 100， Xk5.0 k h， k 0,1,L , n，构造Pn(x)，利用拉格朗n日插值多项式Pn(x)作为f (x)的近似值。分别取n 5，n 10，n 20，同时计算Pn(x)在x 0.75
6、，x 1.75，x 2.75，x 3.75，x 4.75处的函数值。(2) 设f (x) ex，x 1,1，考虑等距节点的拉格朗日插值多项式Pn(x)，即将区间2 01,1进行n等分，记h 一，Xk1.0 k h， k 0,1,L ,n，构造Pn(x)，利用拉格朗日插n值多项式Pn(x)作为f(x)的近似值。分别取 n 5, n 10，n 20 ,同时计算Pn(x)在x 0.95，x 0.05，x 0.05，x 0.95处的函数值。问题2插值区间越小越好吗？考虑下面两个拉格朗日插值问题：1(1)设f (x) r，x 1,1，考虑等距节点的拉格朗日插值多项式R(x)，即将区1 x2 0间1,1进
7、行n等分，记h 2，Xk1.0 k h，k 0,1丄，n,构造Pn(x)，利用拉格朗日n插值多项式Pn(x)作为f (x)的近似值。分别取n 5, n 10，n 20 ,同时计算Pn(x)在0.95处的函数值。x 0.95, x 0.05, x 0.05, x(2)设f(x) ex, x 5,5，考虑等距节点的拉格朗日插值多项式Pn(x)，即将区间一 2 05,5进行n等分，记h ，Xk1.0 k h，k 0,1,L ,n，构造Pn(x)，利用拉格朗日n插值多项式Pn(x)作为f (x)的近似值。分别取n 5，n 10，n 20，同时计算Pn(x)在x 4.75，x 0.25，x 0.25，x
8、 4.75处的函数值。问题3在区间1,1考虑拉格朗日插值问题，为了使得插值误差较小，应如何选取插值节点？考虑下面两个拉格朗日插值问题：(1)设f (x)，x 1,1，考虑非等距节点的拉格朗日插值多项式Pn(x)，记1 xXk cos(2k 1)，k 0,1,L ,n，构造R(x)，利用拉格朗日插值多项式 R(x)作为f(x)的近2(n 1)似值。分别取 n 5 ,n 10, n 20，同时计算 Pn(x)在 x0.95 ,x 0.05 ,x 0.05 ,x 0.95处的函数值。(2)设f(x) ex , x 1,1，考虑非等距节点的拉格朗日插值多项式R(x)，记Xk cos(2k 1) , k
9、 0,1,L ,n，构造R(x)，利用拉格朗日插值多项式 巳(x)作为f(x)的近 2(n 1)似值。分别取 n 5 ,n 10, n 20，同时计算 Pn(x)在 x0.95 ,x 0.05 ,x 0.05 ,x 0.95处的函数值。问题4考虑拉格朗日插值问题，内插比外推更可靠吗?考虑下面两个拉格朗日插值问题：(1)设f(x)J，关于以X。1 ,治4 , X2 9为节点的拉格朗日插值多项式P2(x),利用拉格朗日插值多项式P2(x)作为f(x)的近似值。同时计算F2(X)在x 5 , x 50 ,x 115， x 185处的函数值。(2) 设f(x)浪，关于以X。36，人49，X2 64为节
10、点的拉格朗日插值多项式 F2(x)，利用拉格朗日插值多项式P2(x)作为f (x)的近似值。同时计算P2(x)在x 5，x 50， x 115， x 185处的函数值。(3) 设f(x)上，关于以xo 100，为121，X2 144为节点的拉格朗日插值多项式 F2(x)，利用拉格朗日插值多项式P2(x)作为f (x)的近似值。同时计算P2(x)在x 5，x 50， x 115， x 185处的函数值。(4) 设f(x) ,x，关于以xo 169，捲196，X2 225为节点的拉格朗日插值多项 式P2(x)，利用拉格朗日插值多项式 P2(x)作为f(x)的近似值。同时计算 P2(x)在x 5，
11、x 50， x 115， x 185处的函数值。思考题：1. 对实验1存在的问题，应如何解决？2. 对实验2存在的问题的回答，试加以说明3. 对实验3存在的问题的回答，试加以说明4. 如何理解插值问题中的内插和外推？写出实验报告实验题目2龙贝格(Romberg)积分法方法概要：利用复化梯形求积公式、复化辛普生求积公式、复化柯特斯求积公式的误差估计式计算积分f(x)dx。记h -_-， xk a k h， k 0,1,L ,n，其计算公式：an一般地，利用龙贝格算法计算积分，要输出所谓的T数表龙贝格(Romberg)积分法实验b实验目的：利用龙贝格(Romberg)积分法计算积分f(x)dxa输
12、入：a,b, N,输出：龙贝格T数表程序流程：1 置 h 口，m 1n2输出Ti3对i2,3丄,N，做 3.1 3,53.1置ii21 1置t211 T1h f (a (k22 k 12)h)输出T23.2置S21-(4T2 T1)3输出S23.3对m1,置 C216S2S)输出C2，转 3.63.4对m2，置 R2163(64C2G)输出R2，转 3.63.5对m3，置 tolR2R如果tol ，则停机，否则转3.6h3.6RiR2,CiC2 , SS2,TiT2, h , m24停机问题1:利用龙贝格(Romberg)积分法计算积分1 oX2exdx ,103(2)1 ex sin xdx
13、 ,10(3)011十，101 1(4)0dx，1010问题2：被积函数无界，如何处理?提示：f(0) limsinx 1x 0 X10 6提示：引进变换t/1 X(3)1 COSX ,0.XdX，106提示:利用等式?0SXXdX1 1 dx0-.x1 COSX 10dX，第一个积分值等于2,第二个积分,X利用f(0)lim COS10 ；X 0.X也可以考虑利用分部积分1 COSXcosxdC- x)X，1(4)110 6提示：利用第一类 Gauss-Chebyshev求积公式问题3：积分区间无限，如何处理?2(1) e x dx,，10 610 2提示：利用weXdx作近似1冷X，106
14、提示：利用变换t .T（3）x236e cos xdx ，10提示：Gauss-Hermite 求积公式（4）e xsin2 xdx，10 60 7提示：Gauss-Lagurre 求积公式思考题:1. 输入的参数N有什么意义？2. 在实验1中二分次数和精度的关系如何？3. 在实验2中给出的提示具有普遍性吗？存在其它的方法吗？试加以说明4. 在实验3中给出的提示具有普遍性吗？存在其它的方法吗？试加以说明写出实验报告实验题目3四阶龙格一库塔（Runge Kutta）方法方法概要：给定常微分方程初值问题记Xn a n h， n 0,1 L ,N，利用四阶龙格一库塔方法可逐次求出微分方程初值问题的数值解yn，n 1,2,L ,N四阶龙格一库塔(Runge Kutta)方法实验实验目的：利用四阶龙格一库塔(Runge Ku
