实验三程序的编辑及调试..docxVIP

实验三 程序的编辑及调试姓名： 许梦阳 班级： 学号： 2010301743 实验目的掌握MATLAB程序编辑、运行及调试方法。实验程序1、启动MATLAB后，点击File|New|M-File，启动MATLAB的程序编辑及调试器（Editor/Debugger），编辑以下程序，点击File|Save保存程序，注意文件名最好用英文字符。点击Debug|Run运行程序，在命令窗口查看运行结果，程序如有错误则改正。注：数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如：21310516842163105168421运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。%classic 3n+1 problem from number theory.while 1 n=input(Enter n,negative quits:); if n=0 break end a=n; while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end aend d1Enter n,negative quits:1a = 1Enter n,negative quits:2a = 2 1Enter n,negative quits:3a = 3 10 5 16 8 4 2 1Enter n,negative quits:5a = 5 16 8 4 2 1Enter n,negative quits:7a = 7 22 11 34 17 52 26 13 40 20 10 5 16 8 4 2 12、编程求满足的最小m值。s=0;for i=1:1:100; s=s+2^i; if s10000 x=i; break endendx3、输出所有的“水仙花数”，该类数是一个3位数，其各位数字立方和等于该数本身。方法一x=[];for i=100:1:999;a=fix(i/100);b=fix(i/10)-a*10;c=i-a*100-b*10;if i==a^3+b^3+c^3 x=[x i];endendx方法二x=[];for i=100:1:999;a=rem(i,10);%个位数b=(rem(i,100)-a)/10;%十位数c=(i-a-b*10)/100;%百位数if i==a^3+b^3+c^3;x=[x i];endendx sxh1x = 153 370 371 407实验结果及分析1、启动MATLAB后，点击File|New|M-File，启动MATLAB的程序编辑及调试器（Editor/Debugger），编辑以下程序，点击File|Save保存程序，注意文件名最好用英文字符。点击Debug|Run运行程序，在命令窗口查看运行结果，程序如有错误则改正。注：数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。实验验证得到的结果是正确的，但是仔细分析过程，发现while语句中的条件是n1，所以把程序进一步完善如下：while 1 n=input(Enter n,negative quits:) if n=0 break end a=n; if n==1 a=4; n=a; end while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end aend2、编程求满足的最小m值。 d2x = 13该程序需要使用循环语句，由于要求的结果相当于是循环的次数，但是次数未知。因此在程序中取一个较大的值使得结果能在该区间内，再根据循环即可得出结果。3、输出所有的“水仙花数”，该类数是一个3位数，其各位数字立方和等于该数本身。 sxhx = 153 370 371 407 sxh1x = 153 370 371 407求水仙花数的基本思路是把三位数的个位、十位和百位三位分离开，利用if语句对水仙花数的条件加以判断，如果是的话就存入数组，不是的话用break语句跳出该次循环，继续下一次循环，依次对每个三位数进行判断。方法一采用了fix函数，利用取整的方法分理出三位数，而方法二则采取的取余的方法对三位数进行分离。但是整体的判断思路相同。思考题用对分法求解方程在[0，1]内的解，并验证，在程序中统计出对分次数。