原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
matlab黄金比例
黄金分割
Fibonacci数列Fn通过如下递推格式定义：
Fn=Fn-1+Fn-2，其中F0=F1=1
要求：
计算前51项Fibonacci数，并存入一个向量
利用上述向量计算比值Fn/Fn-1
验证该比值收敛到黄金比例 (1+sqrt(5))/2
二、黄金比例
（1）建立f.m文件在黄金比例文件夹，其功能是计算第n个Fibonacci数
f.m源代码如下（截图如图4）
%计算第n项Fibonacci数
function y=f(n)
f=[1 1];
for i=3:100
f(i)=f(i-1)+f(i-2);
end
y=f(n)
end
图4：f.m源代码截图
计算前51项的值，将其存入向量a，代码如下，截图如图5
n=1:51;
a=f(n);
图5：计算过称截图 （2）计算比值Fn/Fn-1，源代码如下，截图如图6
n=2:51;
y=f(n)/f(n-1);
图6：计算比值过程
（3验证该比值收敛到黄金比例
建立lim.m文件在黄金比例文件夹，其功能是观察f(n)/f(n-1)的曲线和黄金比例的曲线区别。
Lim.m源代码如下（截图如图7）
%计算前50项Fibonacci数，并图表表示
bl=(1+sqrt(5))/2;
a=0;
for n=3:50
a=f(n)/f(n-1);
plot(n,bl,bx);
hold on;
plot(n,a,r+);
end
图7：lim.m源代码截图
在黄金比例文件夹路径下的Command Window输入lim，显示曲线，红色+线为f（n）/f（n-1）的变化曲线，蓝色+线为黄金比例 (1+sqrt(5))/2的直线,如图8.
图8：曲线图
由图像可以得到结论f（n）与f（n-1）比值收敛到黄金比例 (1+sqrt(5))/2。
文档评论（0）