Matlab作业1..docVIP

喻晓磊 200731510103 电子科学与技术 习题一：编写M文件，计算以下问题： 口口口口×口＝口口口口 以上9个口代表1～9这9个数字，不得遗漏或重复 要求：给出解算思路和M文件代码 注意算法的效率 解题理念： 1，设以上方框中的数字分别是1~9之间的数字a、b、c、d、e、f、g、h、i。 2，我们可以推出，个位数e肯定不是1或者5 ，否则两个四位数的个位必然相等，即d=i；再者，其不可能等于9，由不重复的1~9组成的四位数，最大为9876，最小为1234，故e的最大可能取值为9876/1234=8.003……，即最大取8 。 3，我们看到，e的最小取值为2，那么a肯定不能大于或等于5 ，否则乘法运算后的结果为5位数，所以a 只能取1~4。 4，在此基础上，我们采取“穷举法”，逐个试验，从a 开始，每一位与前面出现位的数字不相等，然后检测其是否满足算式。 程序代码如下： global a global b global c global d global e global f global g global h global i for a=1:4 for b=1:9 if b~=a; for c=1:9 if (c~=b)(c~=a); for d=1:9 if (d~=b)(d~=a)(d~=c); for e=2:8 e~=5; if (e~=a)(e~=b)(e~=c)(e~=d); for f=1:9 if (f~=a)(f~=b)(f~=c)(f~=d)(f~=e); for g=1:9 if (g~=a)(g~=b)(g~=c)(g~=d)(g~=e)(g~=f); for h=1:9 if (h~=a)(h~=b)(h~=c)(h~=d)(h~=e)(h~=f)(h~=g); for i=1:9 if (i~=a)(i~=b)(i~=c)(i~=d)(i~=e)(i~=f)(i~=g)(i~=h); x=1000*a+100*b+10*c+d; y=e*x; z=1000*f+100*g+10*h+i; if y==z; r=x s=e t=y end; end; * * * * * * end; end; 得到的结果如下： 我们看到，结果有两组： 1738 * 4= 6952 1963 * 4= 7852 这两组结果都显示到了MATLAB软件的命令窗（COMMAND WINDOW）中，我们是使用了一个小技巧：将x、e、y的值赋给r、s、t时并没有在句末加分号，所以计算的过程也显示出来了，但是看工作区（work space）中，r、s、t的值却只对应后面一组，这是因为，当循环进行到算出第二组结果时，便替换掉了第一组的值。 习题2：丁字路口红绿灯问题建模 丁字路口，每条车道同时只允许一辆车通过，如图所