数值分析课程设计报告.docVIP

数值分析课程设计报告 班级 学号 所在单位 指导老师 201 年 实验一 1．1 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只给猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？ 试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题（15621）。 ①程序清单 n=input(input n：); for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p==fix(p),break,end end disp([x,p]) 运行这段程序后，屏幕出现要求从键盘输入 x 数据的信息inpu n，输入1200后，MALAB计算出合适的 x 和 p0 的值为 1023 15621 p0为整数，必须取 (x +1) 为 4 5( =1024)的倍数，一种简单的处理可取 x = 1023。 1．2 设， （1）从尽可能精确的近似值出发，利用递推公式： 计算机从到的近似值； （2）从较粗糙的估计值出发，用递推公式： 计算从到的近似值； 分析所得结果的可靠性以及出现这种现象的原因。 ①程序清单 （1） syms x n; int(1/(x+5),0,1); ans log(2)+log(3)-log(5) eval(ans) ans = 0.1823 s=0.1823 for n=1:20 s=-5*s+1/n; disp(s) end (2) syms x n; int(x^30/(x+5),0,1) //所得结果为 syms x n; int(x^30/(x+5),0,1) ans = - 931322574615478515625*log(5/6) - 79095966183067699902965545527073/465817912560 eval (ans) 0 //即I（30）≈0 s=0; for n=30:-1:2 s=-s/5+1/5*n; disp(s) end ②运行结果 (1)s = 0.1823 0.0885 0.0575 0.0458 0.0208 0.0958 -0.3125 1.7054 -8.4018 42.1200 -210.5002 1.0526e+003 -5.2629e+003 2.6314e+004 -1.3157e+005 6.5786e+005 -3.2893e+006 1.6447e+007 -8.2233e+007 4.1116e+008 -2.0558e+009 (2)s = 0 6 4.6000 4.6800 4.4640 4.3072 4.1386 3.9723 3.8055 3.6389 3.4722 3.3056 3.1389 2.9722 2.8056 2.6389 2.4722 2.3056 2.1389 1.9722 1.8056 1.6389 1.4722 1.3056 1.1389 0.9722 0.8056 0.6389 0.4722 0.3056  ③算法性能分析 第一小题当中首先算出较为精确的I0，之后经过代入的方法计算出 I（20），结果是比较准确的。但是第二小题当中先是算出必是很精确的I（30），本来所得