数学实验4 - 副本.docVIP

数学与统计学院 数 学 实 验 报 告 姓 名：苟于翠 专 业：信息与计算科学 学 号：201370020107 数列与级数 一、实验目的 1、通过使用编程复习并巩固以前学过的数列与级数的知识； 2、通过编程演示Fabonacci数列、调和级数以及3n+1问题的函数图像及函数关系式； 3、通过图示的方法发现数列与级数的规律及其极限行为，并体会数列与级数在理论与实际应用中的差距； 二、实验的内容、步骤和结果分析 （1） Fibonacci数列 练习1 实验内容：分别取N=20，50，100，200，500，观察Fibonacci数列的折线图。Fibonacci数列是否单调增？它是否趋于无穷？它增加的速度是快还是慢？你能否证实你的观察？ 实验步骤： 方法一：画Fibonacci数列的折线图 语句 FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined\ rightarrow True] ] FibShow[20] 运行结果： FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined\ rightarrow True] ] FibShow[50] 运行结果： FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined\ rightarrow True] ] FibShow[100] 运行结果： FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined\ rightarrow True] ] FibShow[200] 运行结果 FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined\ rightarrow True] ] FibShow[500] 运行结果 结果分析：从实验得出的五个图像可以看出，Fibonacci数列得变化速度非常快，数列单调递增而且趋于无穷大。n的取值越大，图像越陡峭，即递增越快。 练习2：用直线去拟合（）， 实验内容：分别取N=2000，5000，10000，用直线去拟合数据（），，由此求数列的近似表示。注意观察的线性项的系数，它与黄金分割数有何联系？ 实验步骤： n=2000; FibFit[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Log[Fibonacci[i]]}]]; Fit[t,{1,x},x] ] 当n=2000时,结果为Out[1]=-0.80431+0.481212 x n=5000; FibFit[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Log[Fibonacci[i]]}]]; Fit[t,{1,x},x] ] 当n=5000时,结果为Out[2]= -0.804555+0.481212 x n=10000; FibFit[n_Integer]:= Module[ {t={},i}, For[i=1,i?n,i++,AppendTo[t,{i,Log[Fibonacci[i]]}]]; Fit[t,{1,x},x] ] 当n=10000时,结果为Out[3]= -0.804637+0.481212 x 结果分析：从实验结果可以看出，当给点n的值越大，线性拟合的结果越趋于稳定，而且的线性项的系数与黄金分割数的和近似等于1。 （2）调和级数