数学建模——商人过河(Matlab)

最新推荐文章于 2025-01-02 14:31:43 发布
置顶 最新推荐文章于 2025-01-02 14:31:43 发布
阅读量1.2w 收藏 131
点赞数 32
分类专栏： 数学建模 文章标签： 数学建模 多步决策问题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zy440458/article/details/114271787
版权
本文通过数学模型解决了商人与随从安全过河的问题，确保在任何时刻任一岸的商人数都不少于随从数，避免发生危险。采用状态转移的方法，通过允许状态集合和允许决策集合来寻找所有可能的安全过河路径。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

问题描述

三名商人各带–个随从乘船渡河,一只小船只能容纳二人,由他们自己划行.随从们密约,在河的任一岸，一旦随从的人数比商人多,就杀人越货.但是如何乘船渡河的大权掌握在商人们手中.商人们怎样才能安全渡河呢?
对于这类智力游戏,经过一番逻辑思索是可以找出解决办法的.这里用数学模型求解，一是为了给出建模的示例,二是因为这类模型可以解决相当广泛的一类问题,比逻辑思索的结果容易推广.

问题分析

问题已经理想化，所以可以抽象成一个多步决策问题，即在往返的过程中都要保证两岸及船上的商人数大于等于随从数。用允许安全状态变量表示岸上的人员状态，允许决策变量表示船上的人员状况，可以找到状态随决策变化的规律，问题转换为在允许状态下的决策达到安全渡河。

建立模型

此处自行参考《数学模型》姜启源第四版的——商人们怎样安全过河的案例

以下给出代码实现的流程图
在这里插入图片描述

代码

main.m

%% 商人安全过河问题
close, clear, clc
%% 声明全局变量
% setOfAllowedStates 允许状态集合
% setOfDecisionVariables 允许决策集合
% collectionOfStatusMarkers 状态标记集合
% numberOfBusinessmen 商人的总人数
% numberOfFollowers 随从的总人数
% maximumLoadCapacity 船的最大承载人数
% totalSetOfDecisionVariables 允许决策个数
% totalSetOfAllowedStates 允许状态个数
% 0ollectionOfStatusMarkersBegin 初始编号
% collectionOfStatusMarkersProcess 中间编号
% collectionOfStatusMarkersEnd 终止编号
global setOfAllowedStates setOfDecisionVariables collectionOfStatusMarkers numberOfBusinessmen collectionOfStatusMarkersProcess cnt k;
global collectionOfStatusMarkersBegin collectionOfStatusMarkersEnd numberOfFollowers maximumLoadCapacity totalSetOfDecisionVariables totalSetOfAllowedStates;
setOfAllowedStates = [];
setOfDecisionVariables = [];
totalSetOfDecisionVariables = 0;
totalSetOfAllowedStates = 0; 

%% 输入
fprintf('输入大于3 3 2时候将不再作图\n')
numberOfBusinessmen = input('请输入商人的人数：');
numberOfFollowers = input('请输入随从的人数：');
maximumLoadCapacity = input('请输入船的最大承载人数：');

%% 作图属性
axis([0 numberOfBusinessmen 0 numberOfBusinessmen]);
set(gca, 'XTick', 0 : 1 : numberOfBusinessmen), set(gca, 'YTick', 0 : 1 : numberOfFollowers);
xlabel('此岸商人数'), ylabel('此岸随从数'), title('允许状态集合')
grid on

%% 允许状态集合
for businessmen = 0 : numberOfBusinessmen
    for followers = 0 : numberOfFollowers
        if(businessmen == 0 || businessmen == numberOfBusinessmen || (businessmen >= followers && numberOfBusinessmen - businessmen >= numberOfFollowers - followers))
            totalSetOfAllowedStates = totalSetOfAllowedStates + 1;
            setOfAllowedStates(totalSetOfAllowedStates, :) = [businessmen, followers];
            hold on
            plot(businessmen, followers, 'b*', 'MarkerSize', 10);
        end
    end
end

setOfAllowedStates = [setOfAllowedStates, ones(totalSetOfAllowedStates, 1); setOfAllowedStates, zeros(totalSetOfAllowedStates, 1)];
totalSetOfAllowedStates = totalSetOfAllowedStates * 2;
fprintf('共有%d种允许状态\n', totalSetOfAllowedStates)
printfAllowedStates = input('是否打印允许状态集合y/n？','s');
if printfAllowedStates == 'y'
    fprintf('允许状态集合：')
    setOfAllowedStates
end
%% 允许决策集合
for businessmen = 0 : numberOfBusinessmen
    for followers = 0 : numberOfFollowers
        if(businessmen + followers >= 1 && businessmen + followers <= maximumLoadCapacity)
            totalSetOfDecisionVariables = totalSetOfDecisionVariables + 1;
            setOfDecisionVariables(totalSetOfDecisionVariables, :) = [businessmen, followers];
        end
    end
end
fprintf('共%d个允许决策\n',totalSetOfDecisionVariables);
printfDecisionVariables = input('是否打印允许决策集合y/n？','s');
if printfDecisionVariables == 'y'
    fprintf('允许决策集合：');
    setOfDecisionVariables
end
%% 状态标记集合
% 初始状态利用哈希表
collectionOfStatusMarkers = zeros(totalSetOfAllowedStates, 1);
collectionOfStatusMarkersBegin = find(ismember(setOfAllowedStates, [numberOfBusinessmen, numberOfBusinessmen, 1], 'rows') == 1); % 初始状态编号
collectionOfStatusMarkersEnd = find(ismember(setOfAllowedStates, [0, 0, 0], 'rows') == 1); % 终止状态编号
cnt = 0;
k = 1;
collectionOfStatusMarkers(collectionOfStatusMarkersBegin) = 1;
collectionOfStatusMarkersProcess(1) = collectionOfStatusMarkersBegin;

%% 模拟过河（递归实现）
crossTheRiver(collectionOfStatusMarkersBegin);
if cnt == 0
    fprintf('This question currently has no answer.\n')
else
    fprintf('共有%d个方法\n', cnt)
end

crossTheRiver.m

function crossTheRiver(present)
    global setOfAllowedStates setOfDecisionVariables collectionOfStatusMarkers collectionOfStatusMarkersProcess k  numberOfBusinessmen;
    global collectionOfStatusMarkersEnd totalSetOfDecisionVariables cnt;
    if present == collectionOfStatusMarkersEnd
        printCurrentResults();
        if numberOfBusinessmen <= 3
            plotCurrentResults(cnt + 1);
        end
    else
        for i = 1 : totalSetOfDecisionVariables
            possible(1, [1 2]) = setOfAllowedStates(present, [1 2]) + ((-1)^(setOfAllowedStates(present, 3))) * setOfDecisionVariables(i, :);
            possible(1, 3) = 1 - setOfAllowedStates(present, 3);
            tag = ismember(setOfAllowedStates, possible, 'rows');
            if sum(tag) == 1
                next = find(tag == 1);
                if collectionOfStatusMarkers(next) == 0
                    collectionOfStatusMarkers(next) = 1;
                    k = k + 1;
                    collectionOfStatusMarkersProcess(k) = next;
                    crossTheRiver(next);
                    collectionOfStatusMarkers(next) = 0;
                    k = k - 1;
                end
            end
        end
    end
end

printCurrentResults.m

function printCurrentResults()
    global setOfAllowedStates numberOfBusinessmen collectionOfStatusMarkersProcess cnt k numberOfFollowers
    cnt = cnt + 1;
    fprintf('方法 %d：\n',cnt);
    fprintf('\t\t\t\t此\t岸\t\t\t\t\n');
    fprintf('\t商人数\t\t\t随从数\t\t\t船的位置\n');
    for i=1:k
        fprintf('%d:\t %d\t\t\t\t %d',i,setOfAllowedStates(collectionOfStatusMarkersProcess(i),1),setOfAllowedStates(collectionOfStatusMarkersProcess(i),2));
        if setOfAllowedStates(collectionOfStatusMarkersProcess(i),3)==1 % 船在此岸
            fprintf('\t\t\t\t\t此岸\n');
        else
            fprintf('\t\t\t\t\t彼岸\n');
        end
    end
    fprintf('\n');
end

plotCurrentResults.m

function plotCurrentResults(cntPicture)
    global setOfAllowedStates collectionOfStatusMarkersProcess k numberOfBusinessmen numberOfFollowers;
    figure(cntPicture)
    axis([0 numberOfBusinessmen 0 numberOfBusinessmen]);
    set(gca, 'XTick', 0 : 1 : numberOfBusinessmen), set(gca, 'YTick', 0 : 1 : numberOfFollowers);
    xlabel('此岸商人数'), ylabel('此岸随从数'), title(['方法', num2str(cntPicture)])
    grid on
    for i = 2 : k
        hold on
        plot([setOfAllowedStates(collectionOfStatusMarkersProcess(i - 1),1), setOfAllowedStates(collectionOfStatusMarkersProcess(i),1)],...
            [setOfAllowedStates(collectionOfStatusMarkersProcess(i - 1),2), setOfAllowedStates(collectionOfStatusMarkersProcess(i),2)], 'LineWidth',3)
        pause(0.2)
    end
end

动态展示

在这里插入图片描述

如有错误以及可以改进的地方欢迎在下方评论区留言！

商人渡河问题(MATLAB版)
Blogs of zcy1221
03-18 1万+
本文将详细介绍商人渡河问题的求解思路和实现方法(MATLAB).
【数学建模】商人随从过河的一般模型
EnkiduGilgamesh的博客
07-23 1351
三名商人各带一个随从乘船渡河，一只小船只能容纳二人，由他们自已划行，随从们密约，在河的任一岸，一旦随从的人数比商人多就杀人越货。但是如果乘船渡河的大权掌握在商人们手中，商人们怎样才能安全渡河呢？又若四名商人各带一个随从过河，商人如何才能安全渡河？一般地，有n名商人带n个随从过河，船每次能渡k人过河，试讨论商人能安全过河时，n与k应满足什么条件？
数学建模商人过河问题
07-13
解决了商人过河问题，并且可推广到多种情形，直接用MATLAB运行即可，
数学建模之初等模型——商人过河问题（Matlab实现）
热门推荐
captaincoke的博客
03-01 2万+
无关前言： 进入大学一年半了，发现许多不久前学过的知识转头就会忘并且自己的学习能力也没那么尽如人意，在知道了程序员的神器小黄鸭后，想着也许把我所学到的知识讲清楚分享出来更有利于自己的深度理解，同时也能帮到他人避开我所走的弯路。来到这个平台的目的就是分享自己的所学所感并与他人交流，大家共同进步。如果大家发现了我的问题（其实我挺菜的），请在评论区畅所欲言！ 问题叙述 三名商人和他们的三个仆从打算乘船过...
数学建模-商人过河问题
weixin_45895718的博客
03-21 1万+
商人过河问题 1.问题重述 三名商人各带一个随从乘船渡河，一只小船只能容纳二人，由他们自己划行。随从们密约，在河的任一岸，一旦随从的人数比商人多，就杀人越货。但是如何乘船渡河的大权掌握在商人们手中，商人们怎样才能安全渡河呢? 2.问题分析 该问题是一个多步决策的问题，商人的下一步决定方案依赖于此时的状态，同时商人可以选择的方案有多种可能，因此，要想顺利将商人送达到对岸，需要对每一步进行决策和假设，通过一步步推理，分析出最后的结果 3.模型假设 将商人和随从的数目抽象成具体数值 将河的此岸和彼岸抽象成两
商人过河问题的Matlab程序
03-09
商人过河问题的Matlab程序，供学习数学建模或者对趣味性数学感兴趣的人参考。
数学建模-数学模型——商人过河问题.zip
01-20
最后，文档"数学建模-数学模型——商人过河问题.doc"很可能会提供更详细的解题步骤和具体算法实现，包括可能的数据分析、模型验证和结果解释。通过对这个问题的深入探讨，不仅可以提升对数学建模的理解，还能锻炼...
数学建模——商人与随从过河问题
weixin_43326910的博客
03-30 3180
# -*- coding: utf-8 -*- """ ①判断当前状态是否合理 - 得到答案输出 ②合理进行 运输 适用于全部 """ ##初始条件， U = input("请输入小船容量：\n")# 小船容量 x0 = input("请输入商人数：\n")# 商人数 x0 = int(x0) U = int(U) y0 = x0 #随从 sk0 = [x0, y0] #初始s...
商人过河问题——MATLAB实现
最新发布
前排观众的博客
01-02 577
【代码】商人过河问题——MATLAB实现。
数学建模之商人过河（编程法）
qq_44322532的博客
04-18 8321
商人安全过河数学建模培训第一天题目描述题目分析商人过河C++源码 数学建模培训第一天 题目描述 第一问：商人和随从共三对，小船容量2人？ 第二问：四名商人各带一名随从，小船容量3人？ 第三问：假设小船容量三人，最多可以有几对商人及随从安全过河？ 题目分析 多步决策问题 方法1：图解法（不做介绍） 方法2：编程法（易于推广） 第一步：表示状态矩阵，三维（商人数量，随从数量，小船所在侧） S={(...
商人过河问题MATLAB程序
07-24
商人过河MATLAB程序，通篇只有一个while循环，可实现n个商人，n个仆人情况下的最优解（步骤最少）
商人过河问题的MATLAB程序
08-30
数学建模中经典问题商人过河问题的MATLAB源代码，绝对是自己想出来的算法，比枚举算法的时间复杂度要小得多，算法中运用了队列作为存储结构，自己感觉非常的巧妙，前天用C语言实现后，一直想着用MATLAB实现，现在将同样算法的MATLAB源代码提供给大家，希望广大朋友给我提出意见和建议，Email：liang.wang.hubei@gmail.com，谢谢
商人过河的MATLAB程序
09-27
此源代码可以解决经典的商人过河问题，对大家应该会有一些帮助。
商人过河问题的MATLAB实现
08-29
商人过河问题的MATLAB实现，MATLAB源代码。
商人过河问题的MATLAB和C语言程序
08-30
数学建模中经典问题商人过河问题的C语言和MATLAB源代码，绝对是自己想出来的算法，比枚举算法的时间复杂度要小得多，算法中运用了队列作为存储结构，自己感觉非常的巧妙！ 希望广大朋友给我提出意见和建议Email：liang.wang.hubei@gmail.com，谢谢
商人过河MATLAB
weixin_44737931的博客
06-22 3015
%商人过河模型MATLAB实现 clc clear n_s=3;%商人数 n_p=3; %随从数 n_r=2; %船最大承载人数 S_start=[n_s,n_p]; %初始状态 S_end=[0,0]; %最终状态 %生成安全集 S=[]; %安全集 for ii=0:n_s; for jj=0:n_p; if ii>=jj&&n_s-ii>=n_p-jj||ii0|...
商人过河问题，MATLAB代码
qq_43561757的博客
05-24 1218
https://blog.csdn.net/weixin_30562977/article/details/116066354 %%%%%%%%%%%%%%%% 开始 %%%%%%%%%%%%%%%%%%%%%% function jueche=guohe clear all clc %%%%%%%%%%程序开始需要知道商人和仆人数；%%%%%%%%%%%%% shangren=input(输入商人数目: ); puren=input(输入仆人数目: ); rongliang=input(输入船的最大容量:
matlab三个商人三个随从,数学建模 商人过河
weixin_35034536的博客
03-17 7631
数学建模课程作业论文题目：对商人过河问题的研究指导教师：黄光辉小组成员：黄志宇(20156260)车辆工程04班牛凯春(20151927)电气工程05班文逸楚(20150382)工商管理02班一、问题重述3名商人带3名随从乘一条小船过河，小船每次只能承载至多两人。随从们密约，在河的任一岸，一旦随从的人数比商人多，就杀人越货。乘船渡河的方案由商人决定，商人们如何才能安全渡河呢？二、问题分析本题针对商...
商人过河程序matlab编程实现
04-06
适合matlab初学者，提供了商人过河程序的实现，适合所有情况
商人过河问题数学建模matlab
09-13
商人过河问题可以使用数学建模来解决，而MATLAB是一个非常适合进行数学建模和计算的工具。下面是一个使用MATLAB进行商人过河问题数学建模的示例： 首先，我们可以使用二进制向量来表示商人、野人和船的位置状态，其中0表示左岸，1表示右岸。假设商人和野人的数量分别为M和N。 下面是MATLAB代码示例： % 商人和野人的数量 M = 3; N = 3; % 初始状态（左岸） initial_state = [ones(1, M+N), 0]; % 目标状态（右岸） target_state = [zeros(1, M+N), 1]; % 状态转移规则函数 state_transition = @(state, action) state + action; % 判断状态是否合法 is_valid_state = @(state) ... all(state(1:M) >= state(M+1:end) | state(1:M) == 0) && ... all(state(M+1:end) >= state(1:M) | state(M+1:end) == 0); % 判断是否达到目标状态 is_goal_state = @(state) all(state == target_state); % 使用递归函数解决问题 solution = recursive_solve(initial_state, []); % 递归函数 function solution = recursive_solve(state, path) % 如果已经达到目标状态，则返回路径 if is_goal_state(state) solution = path; return; end % 遍历所有可能的动作 actions = generate_actions(state); for i = 1:size(actions, 1) action = actions(i, :); % 计算新状态 new_state = state_transition(state, action); % 如果新状态合法，则继续递归求解 if is_valid_state(new_state) solution = recursive_solve(new_state, [path; action]); % 如果找到解，则返回 if ~isempty(solution) return; end end end % 如果找不到解，则返回空 solution = []; end % 生成所有可能的动作 function actions = generate_actions(state) M = sum(state(1:end-1) == 1); N = sum(state(1:end-1) == 0); % 动作格式：[商人移动数 野人移动数 船移动方向] % 商人独自移动 actions = [-1 0 -1; -2 0 -1; 1 0 1; 2 0 1]; % 野人独自移动 actions = [actions; 0 -1 -1; 0 -2 -1; 0 1 1; 0 2 1]; % 商人和野人一起移动 for i = 1:M for j = 1:N if i + j <= 2 actions = [actions; -i -j -1; i j 1]; end end end end 这个示例代码使用了递归求解的方法来找到商人过河问题的解。代码中的state_transition函数定义了状态转移规则，is_valid_state函数判断状态是否合法，is_goal_state函数判断是否达到目标状态。generate_actions函数生成所有可能的动作。 注意：这只是商人过河问题的一个简单数学建模示例，实际问题可能需要更多的约束和复杂的规则。你可以根据具体需求进行修改和扩展。
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值