商人过河matlab程序以及解析.pdfVIP

数 学 建 模 作 业 班级：数学 131 姓名：丁延辉 学号： （二）商人过河 Matlab 代码 三个商人三个随从 z=zeros(30,3);%z 为由 (a,b,c) 的列向量组成的 3 行 30 列数组，初 始化为 0 矩阵， a,b,c 代表此刻此岸的商人，仆人数量以及船的运行 状态， c=1 表示即将向彼岸运行 m=zeros(1,20);%m 为一维行向量，初始化为 1 矩阵，用于在后面的 程序中判断第 k 次选择的乘船方案 d=[0,1,1;0,2,1;1,0,1;1,1,1;2,0,1];% 共有 5 种可以选择的乘船方 案，最后面一列全为 1，即用于在后面表示使得 z(k,3) 的取值保持随 着 k 的奇偶性保持着 0-1 变换 . z(1,:)=[3,3,1];% 初始状态为 [3,3,1] k=1; m(k)=1;%第一次默认的乘船方案为决策 1—— d(1) flag=1;% 用于在后面判断是否成功找到方案 answer=0;%用于在后面判断是否找到答案 whilek0%保持 k0 ifm(k)5 flag=0; break; end p=0; z(k+1,:)=z(k,:)+(-1)^k*d(m(k),:);% 每一次的运算规则都是 z(k+1)=z(k)-(-1)^k*d(m(k),:),d(m(k),:) 表示决策方案 a=z(k+1,1);% 将当前情况的矩阵数值复制给 a 商人 ,b 仆人 b=z(k+1,2); c=z(k+1,3); if(a==3(b==0||b==1||b==2||b==3))||(a==1b==1)||(a==2b= =2)||(a==0(b==0||b==1||b==2||b==3))% 判断 (a,b) 是否符合限定 情况 forj=1:k% 判断是否此岸 a,b,c 与之前有重复，如果是，结束此次循 环，重新选择乘船方案 ifa==z(j,1)b==z(j,2)c==z(j,3) ifm(k)~=5% 决策方案只有 5 种，所以 m(k)=5 ， m(k)=m(k)+1;% 因为有重复，所以换下一种决策方案 else while(m(k)==5)(k1) k=k-1;% 回溯，这一步骤已经把所有决策取尽，无可用解法，于是将 后退一步，同时换下一种决策方案 end%while 循环的目的是防止前面几步的决策都是 5，导致 k=k-1,m(k)=m(k)+1 后数组越界，一直找到前面不是 m(k)=5 的步骤 m(k)=m(k)+1; end p=1; break; else p=0; end end ifp==1%程序在跳出内层 for 循环之后，因为要换成决策方案，所以 同时跳出，直接进入下一次 while 循环， continue; end ifa==0b==0%判断是否达到目标情况 answer=1; fprintf(Successfullyfound!\n 每一次的此岸人员分布：商人仆人 \n) fori=1:100 fprintf( 第%2d次%d%d\n,i,z(i,1),z(i,2)) ifz(i,1)==0z(i,2)==0 break; end end%如果不是，进入下一步骤，计算 z(k+2) ifm(k)~=5 m(k)=m(k)+1;%这是正常的进入下一次， 所以仍从 d1 乘船决策 1 开始 els