数学实验Fiboncci数列.docxVIP

Fibonacci数列实验内容1.（1）首先建立一个名叫shulie.m的m文件，依据条件建立an的数组。然后显示出an的图形。functionshulie(n)an=[1]; fori=2:n an=[an,an(i-1)+1./an(i-1)]; endplot(x,an) 1．（2）所得图形如上，初步观察，应该是一个幂函数，幂为1/2.而且通过观察图形，an的极限不存在。1．（3）因为推测an是n的1/2次幂，所以运用matlab的plot函数以n为横坐标，an^2为纵坐标画图，结果如下：functionshulie(n)an=[1]; fori=2:n an=[an,an(i-1)+1./an(i-1)]; endxn=1:n; an=an.^2 ;plot(an); 由上图可知，图像接近于直线，猜想的到认证。为了得到直线的斜率和截距，可以运用matlab中的polyfit函数进行求解，语句如下polyfit(xn,an,1) %将an^2与n进行线性拟合结果如下：ans = 2.0148 0.7931所以线性拟合的结果为：an^2=2.0148x+0.7931为了的到线性拟合与实际情况的近似情况，将线性拟合的直线与实际的线画于同一坐标系下，functionshulie(n)fn1=[];