实验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红星
微信“图像处理与模式识别研究所”关注我呦