华工数学实验报告-斐波那契数列.docVIP

《数学实验》报告

学院：电子信息学院

专业班级：信息工程电联班

学号：

姓名：

实验名称：实验二斐波那契数列

实验日期：2016/04/05

实验目的

认识Fibonacci数列，体验发现其通项公式的过程。了解matlab软件中，进行数据显示与数据拟合的方式提高对数据进行分析与处理的能力。

实验任务

1.讨论调和级数的变化规律

（1）画出部分和数列变化的折线图，观察变化规律；

（2）引入数列，作图观察其变化，猜测是否有极限；

（3）引入数列，作图观察其变化，寻找恰当的函数拟合；

（4）讨论调和级数的部分和数列的变化规律。

2.人口问题是我国最大社会问题之一，估计人口数量和发展趋势是我们制定一系列相关政策的基础。从人口统计年鉴，可查我国从1990年至2010年人口数据资料如下，试根据表中数据，分析人口增长的规律，并以此预测2011年和2012年的人口数量，然后与实际人口数量做对比评价模型的优劣，并对我国人口政策提出建议。

表1不同年份我国的人口数量（万）

年份

1990

1991

1992

1993

1994

1995

1996

1997

数量

114333

115823

117171

118517

119850

121121

122389

123626

年份

1998

1999

2000

2001

2002

2003

2004

2005

数量

124761

125786

126743

127627

128453

129227

129988

130756

年份

2006

2007

2008

2009

2010

数量

131448

132129

132802

133450

134091

实验过程

3.1实验原理

3.1.1任务一

通过用for循环语句来进行操作，用plot进行画图，通过看图猜测函数的类型，判断是对数函数，取指数后，利用polyfit进行拟合，判断猜测成立。

3.1.2任务二

用polyfit进行拟合，R5=dot(y-polyval(p5,t),y-polyval(p5,t))计算拟合残差，

再用polyval预测2011和2012年的人口。

3.2算法与编程

3.2.1任务一

functionfib1(n)

s=[1,1.5];

fori=3:n

s=[s,s(i-1)+1/i];

end

plot(s);

functionfib2(n)

s=[1];

fori=2:2*n

s=[s,s(i-1)+1/i];

end

h=[0.5];

forj=2:n

h=[h,s(2*j)-s(j)];

end

plot(h)

functionfib3(n)

s=[1];

fori=2:2*n

s=[s,s(i-1)+1/i];

end

g=[1.5];

forj=2:n

g=[g,s(2*j)];

end

plot(g)

y=exp(g);

x=[1:n];

p=polyfit(x,y,3)

plot(x,y,ro,x,polyval(p,x),b)

legend(数据点,拟合曲线);

3.2.2任务二

functionfib4()

t=[0:20];

y=[114333,115823,117171,118517,119850,121121,122389,123626];

y=[y,124761,125786,126743,127627,128453,129227,129988,130756,131448,132129,132802,133450,134091];

p2=polyfit(t,y,2)

p4=polyfit(t,y,4)

p5=polyfit(t,y,5)

R5=dot(y-polyval(p5,t),y-polyval(p5,t))%计算拟合残差

R4=dot(y-polyval(p4,t),y-polyval(p4,t))

R2=dot(y-polyval(p2,t),y-polyval(p2,t))

plot(t,y,r+,t,polyval(p2,t),t,polyval(p4,t),t,polyval(p5,t))

legend(我国人口数量,2阶拟合,4阶拟合,5阶拟合)

a=polyval(p5,21)

b=polyval(p5,22)

3.3计算结果或图形

3.3.1任务一

fib3(50)

p=

3.561