【Matlab代码】斐波拉契数列

实验1：显示数列的函数类型

function fib1(n)%显示数列前n项
fn=[1,1];%将前两项放到数组fn中
for i=3:n
    fn=[fn,fn(i-2)+fn(i-1)];%将第i项添加到数组fn中
end
plot(fn)%前n项显示成曲线

实验2：验证上步得到的结论。
显示取对数后的数列前n项：

function fib2(n)%显示取对数后数列前n项
fn=[1,1];%将前两项放到数组fn中
for i=3:n
    fn=[fn,fn(i-2)+fn(i-1)];%将第i项添加到数组fn中
end
fn=log(fn);%将原来的数据取对数
plot(fn)%前n项显示成曲线

实验3：获取数列的近似函数关系
用一阶多项式做拟合：

function y=fib3(n)%取对数后做线性拟合
fn=[1,1];%将前两项放到数组fn中
for i=3:n
    fn=[fn,fn(i-2)+fn(i-1)];%将第i项添加到数组fn中
end
xn=1:n;%横坐标
fn=log(fn)%将原来的数取对数
y=polyfit(xn,fn,1)%一阶多项式拟合

注意1：

修改为：fn=[fn,fn(i-2)+f(i-1)]
注意2：

修改为：y=polyfit(xn,fn,1)
实验4：观察拟合数据与原始数据的吻合。

function fib4(n)%显示拟合数据与原始数据
fn1=[];%装拟合数据的数组
for i=1:n
    fn1=[fn1,0.4476*1.618^i];%将第i项添加到数组fn1中
end
fn2=[1,1];%装原始数据的数组
for i=3:n
    fn2=[fn2,fn2(i-2)+fn2(i-1)];%将第i项添加到数组fn2中
end
x=1:n;
plot(x,fn1,x,fn2,'r*')%显示，fn1-兰线，fn-2红星

微信“图像处理与模式识别研究所”关注我呦