再举一个简单的M文件实例.docxVIP

8 8 再举一个简单的M文件实例： sunooy 发表于 2006-6-22 0:07:00 日本数学家角谷提出一个著明的猜想： 任何一个大于1的自然数经过如下的有限步运算，最后结果为1 如果是偶数，将其除以2 如过是奇数，将其乘以3,再加上1 例如n = 10 10/2=5 5*3+1=16 16/2=8 8/2 = 4 4/2 = 2 2/2 = 1 根据上面提供的猜想：建立M文件，以便方便的验证对不同自然数，角谷猜想 是否正确。程序可以参考这里： user2/su nooy/archives/2006/1275717.shtml 第一步：建立M文件 % M — file, jiaogu.m %验证对不同的自然数n,角谷猜想是否成立 n=input(Type n:); %从键盘输入n,并赋值给n i=1; %标记运算的步数 f=[n]; %存放每次运算后的n的值 while n~=1 %当n不等于1时就循环，直到n= 1时为止。 if mod(n,2)==0 n=n/2; else n=n*3+1; end i=i+1; f(i)=n; end f %输出每次运算后的n的值组成的集合 I=i-1 %输出运算的步数,f集合中会有i个元素，但实际上是只经过i-1步的运 算 n %输出最后一步所得n的值，显然，一般n = 1，除非猜想不成立时，n会 取到其它的值！ 第二步：保存M文件: 进入File菜单，打开New菜单，选择M — file，在保存，注意为了直观最好将 文件名保存为jiaogu 第三步：调用M文件 在matlab窗口输入调用命令jiaogu回车后显示Type n:立即给n赋不同的 值，马上得到结果。比如： n= 10时，经过6步可将10化为1 16 8 4 2 1 TOC \o 1-5 \h \z Type n:10 f =10 5 16 8 4 2 1 I = 6 n= n= 55时，经过20步可将55化为1 f = 55 166 83 250 4 47 142 ..…4 2 125 376 188 9 1 （为节省空间，中间的数略去了） n= 100时，经过25步可将100化为1 f =100 50 25 76 38 19 58 29 88 44 22 11 34 17 52 26 13 40 20 10 5 16 8 4 2 1 I =25 n =1 还可给n取更多的值或更大的值，总之，建立M文件够，验证工作变得十分方 便了！ 特别说明： 如果把while n~=1改为while n~=1 n0程序会更完美。写成while n~=1表 面上看是没有错的。但因为角谷猜想毕竟是个猜想，我们无法验证所有的数， 如果存在某个特别大的数使猜想不成立，它会从 n逐渐的变小，直到负的无穷 大，而不能被化为1。这样程序就会进入死循环！所以应该写成 while n~=1 n0，这就保证如果有某个数不满足角谷猜想，一旦发现它不是正数了，就停 止循环，从而避免了死循环。 如何建立和调用 M文件(M-File)呢?(原创) sunooy 发表于 2006-6-21 21:17:00 建立M文件的必要性： 对于常用的或很长的程序，如果每次都从头输入，则非常费时间，也容易出错?将 这些程序做成M文件就方便多了 --只要输入自定义的命令，就可以方便的调用, 编程的效率大大提高！ 如何建立M文件 (1)在记事本(.txt)或maltab命令窗口中编辑好 M文件 例:编写一个输入半径，计算圆的面积和球的体积 的M文件. 前两行是文件头，一般说明文件的 类型(M-file),文件名(tutex2),和文件的功能 (输入半径,计算圆的面积和球的体积),后面是文件体，主要是解决特定问题的程 序代码? % M-file, tutex2.m %输入半径，计算圆的面积和球的体积 r = input(Type radius:); %从键盘输入 Type radius的值，并将它赋值给r. area=pi*rA2; volume=(4/3)*pi*rA3; fprintf(The radius is %12.5f\n,r) fprintf(The area of a circle is %12.5f\n,area) fprintf(The volume of a sphere is %12.5f\n,volume) ⑵将上面的代码保存到工作空间里.保存方法是：File/New/F-ile/保存 如何调用M文件？ 在matlab工作窗口直接键入M文件名，即可调用. 以上面的程序为例， tutex2 %键入文件名，回车 Type radius: %提示输入半径,输入半径，比如13.5,回车 显示如下 Type radius:13.5 The radiu