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

再举一个简单的M文件实例：sunooy 发表于 2006-6-22 0:07:00 ? 日本数学家角谷提出一个著明的猜想: 任何一个大于1的自然数经过如下的有限步运算,最后结果为1 ? (1)如果是偶数,将其除以2 ? (2)如过是奇数,将其乘以3,再加上1 例如 n＝1010/2=5 5*3+1=1616/2=8 8/2＝44/2＝22/2＝1 根据上面提供的猜想：建立M文件，以便方便的验证对不同自然数，角谷猜想是否正确。程序可以参考这里：user2/sunooy/archives/2006/1275717.shtml? 第一步：建立M文件 % M－file, jiaogu.m%验证对不同的自然数n，角谷猜想是否成立n=input(Type n:);? ?%从键盘输入n，并赋值给ni=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;endf???? %输出每次运算后的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 Type n:10f =10???? 5??? 16???? 8???? 4???? 2???? 1I =?? 6n = 1 －－－－－－－－－－－－－－ n＝55时，经过20步可将55化为1??f =? 55???????? 166????????? 83???????? 250???????? 125???????? 376???????? 188????????? 94????????? 47???????? 142 .....? 4?????????? 2?????????? 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???? 1I =25n =1 还可给n取更多的值或更大的值，总之，建立M文件够，验证工作变得十分方便了！ --------------------------------------------- 特别说明：如果把while n~=1改为while n~=1n0程序会更完美。写成while n~=1表面上看是没有错的。但因为角谷猜想毕竟是个猜想，我们无法验证所有的数，如果存在某个特别大的数使猜想不成立，它会从n逐渐的变小，直到负的无穷大，而不能被化为1。这样程序就会进入死循环！所以应该写成 while n~=1n0 ，这就保证如果有某个数不满足角谷猜想，一旦发现它不是正数了，就停止循环，从而避免了死循环。 ??????? 如何建立和调用M文件(M-File)呢?(原创)sunooy 发表于 2006-6-21 21:17:00 1.建立M文件的必要性: 对于常用的或很长的程序,如果每次都从头输入,则非常费时间,也容易出错.将这些程序做成M文件就方便多了--只要输入自定义的命令,就可以方便的调用,编程的效率大大提高! 2.如何建立M文件 (1)在记事本(.txt)或maltab命令窗口中编辑好M文件 例:编写一个输入半径,计算圆的面积和球的体积 的M文件. 前两行是文件头,一般说明文件的 类型(M-file),文件名(tutex2),和文件的功能(输入半径,计算圆的面积和球的体积),后面是文件体,主要是解决特定问题的程序代码. % M-file, tutex2.m? % 输入半径,计算圆的面积和球的体积 r = input(Type radius:); ％从键盘输入Type radius的值，并将它赋值给r. area=pi*r^2; volume=(4/