基于MATLAB平台的三例数学实验.docVIP

基于MATLAB平台的三例数学实验 杨亚辉 （海南软件职业技术学院 琼海 571400） 摘要：论述了在MATLAB平台上完成关于角谷猜想、隐函数图像绘制和首位数问题的数学实验，说明了实验的具体步骤和结论. 关键词： 数学实验；角谷猜想；隐函数图像的绘制；首位数问题；MATLAB Three mathematics experiments based on MATLAB Yang Ya-Hui ( Hainan Software Vocational---Technical College ,Qionghai,Hainan,571400) Abstract: The paper, based on MATLAB, expounds how to complete the following mathematics experiments, such as Sizuo Kakutani assumption , Implicit Function graph-drawing, The first figure question, narrates the concrete procedure and conclusion of the experiment, and indicates the immense power of MATLAB and important function of mathematics experiment in the mathematics education. Key words: Mathematics experiment, Sizuo Kakutani assumption, Implicit Function graph-drawing, The first figure question, MATLAB 1.角谷猜想的验证 1.1角谷猜想 日本数学家角谷提出一个著明的猜想，其内容为： 任何一个大于1的自然数经过如下的有限步运算,最后结果为1. (1)如果是偶数,将其除以2.(2)如过是奇数,将其乘以3,再加上1. 1.2角谷猜想的验证 1.2.1实验说明用手工计算的方式可以验证比较小的自然数，但对于比较大的自然数（例 大于1亿的自然数）用手工计算起来就很吃力，也很费时间。要想快速的、批量的（即某个给定区间上所有的自然数）验证对于任意给定大小的自然数角谷猜想是否成立，就必须利用计算机强大的计算能力。下面是一个在MATLB平台上完成的验证角谷猜想的实验. 1.2.2实验步骤 第一步：根据角谷猜想的内容,在MATLAB平台上建立验证角谷猜想的M文件. % M－file, jiaogu.m，将文件保存为 jiaogu n=input(Type n:); i=1; f=[n]; while n1 if mod(n,2)==0 n=n/2; else n=n*3+1; end i=i+1; f(i)=n; end f I=i-1 n %若角谷猜想成立，则n＝1，否则n为某一负数。 第二步：调用jiaogu.m文件验证角谷猜想 在matlab窗口输入jiaogu命令.即可调用第一步建立的程序, 回车后显示Type n: 立即给n赋不同的值，马上得到结果。比如： 令n=8975431,经过287步可将8975431化为1.令n=95478631245,经过259可将95478631245化为1. 1.2.3实验结论 通过若干次的验证,发现这些数都满足角谷猜想.给n取更多的值或更大的值，可以进一步增强我们对角谷猜想的信任.对jiaogu.m文件稍作修改，就可以让它逐一验证某个给定区间上每一个自然数是否满足角谷猜想，实现角谷猜想验证的自动化和批量化,大大的提高了实验的效率。限于篇幅，不再赘述。 2.隐函数图像的绘制隐函数即满足方程f(x,y)=0的x,y之间的关系式.隐函数的显化是困难的,甚至 是不可能的,所以往往无法先定义一个x向量再求出相应的y向量,这给绘制隐含数f(x,y)的图像带来了困难.许多学生由于没有见过隐函数的图像,所以感到隐函数很难理解、很抽象.MATLAB提供一个专门绘制隐函数曲线的函数,该函数调用格式为 ezplot(f(x,y)) 下面通过一个实例来演示用该函数绘制隐函数图像的方法 问题:试绘制出隐函数的图像 求解: 利用如下的MATLAB命令,可得到图（a））） 图 2 参考文献: [1] 薛定宇，陈阳泉，高等应用数学问题的MATLAB求解[M]，北京，清华大学出版社，2004，33-34 [2] 杨亚辉，基于MAT