MATLAB实现矩阵快速幂.docVIP

MATLAB实现矩阵快速幂 武汉大学 MATLAB及其应用课程 课程论文 学生姓名： 学 号： 电话： 班 级： 指导教师： 2016年4月4日 目录 一、实验题目 ................................................................................. 3 二、实验目的 ................................................................................. 3 三、实验设备和条件 .................................................................... 3 四、实验原理 ................................................................................. 3 五、实验内容与过程 .................................................................... 6 六、程序代码截图 ........................................................................ 8 七、实验结果 ............................................................................... 10 八、实验收获与总结 .................................................................. 12 一、实验题目 利用矩阵快速幂算法与符号运算求解大型斐波那契数。 二、实验目的 lt;1gt;实现在较短的时间内完成大规模运算。 lt;2gt;充分利用MATLAB方便的矩阵运算以及超高精度的符号运算的特色，进行对大规模数值的求解运算。 三、实验设备和条件 lt;1gt;设备：MATLAB R2015a软件 lt;2gt;条件： ①MATLAB拥有着极其方便的矩阵运算，每一个数组都可以表示为相应维度的矩阵，可以通过两个数组的直接相乘来获得“矩阵1×矩阵2”运算的结果矩阵。 ②MATLAB拥有着精度超高的符号运算，而且使用起来非常方便。 四、实验原理 根据题目，“利用矩阵快速幂算法与符号运算求解大型斐波那契数”，需要将该实验的题目拆成“矩阵快速幂算法”、“大型”和“斐波那契数”几个关键点进行原理的阐述。 为了方便描述，稍排了个序进行原理讲解： lt;1gt;“斐波那契数” 斐波那契数列（Fibonacci sequence），又称黄金分割数列、因数学家列昂纳多·斐波那契（Leonardoda Fibonacci ）以兔子繁殖为例子而引入，故又称为“兔子数列”，指的是这样一个数列：0、1、1、2、3、5、8、13、21、34、??在数学上，斐波纳契数列以如下被以递推的方??定义：F（0）=0，F（1）=1，F（n）=F(n-1)+F(n-2)（n≥2，n∈N*）在现代物理、准晶体结构、化学等领域，斐波纳契数列都有直接的应用。 而斐波那契数，就是指在已经定义了递推初始值的前提下，持续递推斐波那契数列所得的其中一个数。 lt;2gt;“大型” 对于诸如C，C++，JAVA等语言，进行高精度数字的求解总是非常具有挑战性的，即使只是进行大规模整数运算，C++也需要铺写大量的代码定义大数类和重载操作符。如果再加入矩阵运算，使用起来极其不方便。甚至在大数类定义不够严谨的情况下，轻则造成运行效率极低，数组空间溢出等问题，重则造成数据溢出还找半天找不到bug的“惨案”。 而MATLAB的符号类精度之高非常恐怖，且操作异常方便，可以像处理普通整型一样地 处理符号类。因此，使用MATLAB可以让“大型”运算变得具有非常强的可行性。 所谓“大型”，当然不会是求第10个、第20个、第100个斐波那契数这等简单循环便 20可以轻易得到结果的数据（当然，其实第100个斐波那契数已经达到了3.54×10的整型存 放不下的大数）。 所谓“大型”，应该是指第一万个、第十万、百万，乃至第千万、第一亿个斐波那契数。根据下图所示，斐波那契数列有着跟指数