商人过河问题-MATLAB实现

function draw(pm,qianan,n) clf; axis equal axis(0.5+[0,14,0,11]) set(gca,'xtick',0:1:14,'ytick',0:1:11,'xcolor','w','ycolor','w') set(gca,'color','y') hold on river=[7,1;7,2;7,3;7,4;7,5;7,6;7,7;7,8;7,9;7,10;7,11;8,1;8,2;8,3;8,4;8,5;8,6;8,7;8,8;8,9;8,10;8,11]; scatter(gca,river(:,1),river(:,2),900,'bs','filled') cs=[pm,pm]; qa=qianan(n,:); %红色为商人，绿色为仆人 for i=1:2 if qa(i)==0 continue elseif mod(qa(i),2)~=0 side=(qa(i)-1)/2; people=[2+i,6]; for j=1:side people=[people;2+i,6+j;2+i,6-j]; end else side=qa(i)/2; people=[2+i,6]; for j=1:side people=[people;2+i,6+j;2+i,6-j]; end people=people(1:qa(i),:); end if i==1 scatter(gca,people(:,1),people(:,2),450,'rs','filled') else scatter(gca,people(:,1),people(:,2),450,'gs','filled') end end da=cs-qa; for i=1:2 if da(i)==0 continue elseif mod(da(i),2)~=0 side=(da(i)-1)/2; people=[10+i,6]; for j=1:side people=[people;10+i,6+j;10+i,6-j]; end else side=da(i)/2; people=[10+i,6]; for j=1:side people=[people;10+i,6+j;10+i,6-j]; end people=people(1:da(i),:); end if i==1 scatter(gca,people(:,1),people(:,2),450,'rs','filled') else scatter(gca,people(:,1),people(:,2),450,'gs','filled') end end title('商仆过小河') end