建模计算机模拟框图,数学建模报童问题

《数学建模报童问题》由会员分享,可在线阅读,更多相关《数学建模报童问题(5页珍藏版)》请在读根文库上搜索。

1、数 学 实 验 报 告 1 计 算 机 模 拟实 验 目 的 :1.学 习 计 算 机 模 拟 的 基 本 过 程 与 方 法 ;2.会 做 简 单 的 计 算 机 模 拟 。实 验 内 容 :一 、 了 解 什 么 是 模 拟模 拟 就 是 利 用 物 理 的 、 数 学 的 模 型 来 类 比 、 模 仿 现 实 系 统 及 其 演 变 过 程 , 以 寻 求 过 程 规 律 的 一 种 方 法 。 模 拟 的 基 本 思 想 是 建 立 一 个 试 验 模 型 , 这 个 模 型 包 含 所 研 究 系 统 的 主 要 特 点 通 过 对 这 个 实 验 模 型的 运 行 , 获 得 所。

2、 要 研 究 系 统 的 必 要 信 息模 拟 的 方 法 :1、 物 理 模 拟 :对 实 际 系 统 及 其 过 程 用 功 能 相 似 的 实 物 系 统 去 模 仿 。 例 如 , 军 事 演 习 、 船 艇 实 验 、 沙 盘 作 业 等 。物 理 模 拟 通 常 花 费 较 大 、 周 期 较 长 , 且 在 物 理 模 型 上 改 变 系 统 结 构 和 系 数 都 较 困 难 。 而 且 , 许 多 系统 无 法 进 行 物 理 模 拟 , 如 社 会 经 济 系 统 、 生 态 系 统 等 。2、 数 学 模 拟在 一 定 的 假 设 条 件 下 , 运 用 数 学 运 算 。

3、模 拟 系 统 的 运 行 , 称 为 数 学 模 拟 。 现 代 的 数 学 模 拟 都 是 在 计 算 机 上 进 行 的 , 称 为 计 算 机 模 拟 。 计 算 机 模 拟 可 以 反 复 进 行 , 改 变 系 统 的 结 构 和 系 数 都 比 较 容 易 。在 实 际 问 题 中 , 面 对 一 些 带 随 机 因 素 的 复 杂 系 统 , 用 分 析 方 法 建 模 常 常 需 要 作 许 多 简 化 假 设 , 与 面临 的 实 际 问 题 可 能 相 差 甚 远 , 以 致 解 答 根 本 无 法 应 用 。 这 时 , 计 算 机 模 拟 几 乎 成 为 唯 一 的 。

4、选 择 。二 、 报 童 问 题某 报 童 以 每 份 0.03元 的 价 格 买 进 报 纸 ,以 0.05元 的 价 格 出 售 . 根 据 长 期 统 计 ,报 纸 每 天 的 销 售 量 及百 分 率 为 销 售 量 200 210 220 230 240 250 百 分 率 0.10 0.20 0.40 0.15 0.10 0.05已 知 当 天 销 售 不 出 去 的 报 纸 , 将 以 每 份 0.02元 的 价 格 退 还 报 社 .试 用 模 拟 方 法 确 定 报 童 每 天 买 进 报 纸 数 数 学 实 验 报 告 2 量 ,使 报 童 的 平 均 总 收 入 为 最 。

5、大 ?1 系 统 的 假 设 :( 1) 模 拟 时 间 充 分 大 ;( 2) 报 童 购 买 报 纸 量 介 于 销 售 量 最 小 值 与 最 大 值 之 间 ;( 3) 不 考 虑 有 重 大 事 件 发 生 时 卖 报 的 高 峰 期 , 也 不 考 虑 风 雨 天 气 时 卖 报 的 低 谷 期 。2 问 题 分 析报 童 购 进 数 量 应 根 据 需 求 量 确 定 , 但 需 求 量 是 随 机 的 , 所 以 报 童 每 天 如 果 购 进 的 报 纸 太 少 , 不 够 买的 , 会 少 赚 钱 ; 如 果 购 进 太 多 , 卖 不 完 就 要 赔 钱 , 这 样 由 。

6、于 每 天 报 纸 的 需 求 量 是 随 机 的 , 致 使 报 童 每 天的 收 入 也 是 随 机 的 , 因 此 衡 量 报 童 的 收 入 , 不 能 是 报 童 每 天 的 收 入 , 而 应 该 是 他 长 期 ( 几 个 月 、 一 年 )卖 报 的 日 平 均 收 入 。 我 们 可 以 应 用 计 算 机 模 拟 的 方 法 在 模 拟 时 间 充 分 大 的 条 件 下 ( 例 如 10000 天 ) , 模 拟 每 天 的 销 售 量 , 因 而 确 定 每 天 应 买 进 多 少 报 纸 才 能 使 平 均 总 收 入 达 到 最 大 值 。设 该 问 题 的 销 售。

7、 量 为 离 散 型 随 机 变 量 , 有 分 布 列 ii pxP )( )62,1(i 。 所 谓 对 离 散 型随 机 变 量 的 模 拟 , 其 实 质 就 是 通 过 计 算 机 上 的 模 拟 试 验 来 取 得 在 真 实 试 验 下 的 样 本 值 U ,由 分 布 列知 U 只 取 数 值 1x ,2x 中 的 某 一 个 .对 离 散 型 随 机 变 量 的 一 般 模 拟 方 法 如 下 :(1) 在 直 线 上 将 0 ,1区 间 划 分 成 若 干 个 子 区 间 , 其 分 点 记 为 1A ,2A 6A ,分 点 iA 的 坐 标kik k pp 1 )6,2,。

8、1(i如 下 图 所 示 : 0 1A 2A 4A 5A 1 1p 21 pp 41k kp 51k kp x( 2) 在 计 算 机 上 取 伪 随 机 数 rand ;( 3) 判 定 rand 在 x轴 上 所 落 的 位 置 , 若 ),( 111 ik kik k pprand , 则 可 以 认 为 的 一 次 真 实 试 验 样 本值 取 ix 。3 符 号 假 设BUYMIN: 每 天 的 最 小 购 买 量 BUYMAX: 每 天 的 最 大 购 买 量SIMUDAY: 模 拟 时 间sell_amount: 报 童 销 售 量 buy_amount: 报 童 购 买 量pe。

9、rcentage: 销 售 百 分 率 ave_profit: 总 平 均 利 润loop_buy : 当 天 购 买 量 loop_day : 当 天 时 间 数 学 实 验 报 告 3 4 模 拟 框 图 初 始 化 : BUYMIN;BUYMAX;SIMUDAY ;percentagebuy_amount=0;ave_profit=0;loop_buy=BUYMIN; sum_profit=0;loop_day=1产 生 随 机 数 rand , 根 据 rand 选 取 sell_amount loop_day+loop_buySIMUDAY?loop_buyBUYMAX? 输 出 m。

10、 ax(buy_amount), m ax(ave_profit);loop_buy+ YNN Y Y N END 数 学 实 验 报 告 4 5 计 算 机 程 序 : 在 Matlab软 件 包 中 编 程 , 共 需 两 个 文 件 :main.m,Getprofit.m, 主 程 序 为 main.m.% 主 文 件 main.m:BUYMIN=200; % 每 天 的 最 小 购 买 量BUYMAX=250; % 每 天 的 最 大 购 买 量SIMUDAY=1.0e+5; % 模 拟 时 间sell_amount=200:10:250; % 销 售 量percentage=0.1 。

11、0.3 0.7 0.85 0.95 1; % 百 分 率buy_amount=0;ave_profit=0;for loop_buy=BUYMIN:BUYMAXsum_profit=0; for loop_day=1:SIMUDAYindex=find(percentage=rand); % 产 生 随 机 数 , 用 于 决 定 当 天 的 销 售 量sum_profit=sum_profit+GetProfit(loop_buy,sell_amount(index(1);endbuy_amount=buy_amount,loop_buy; % 循 环 嵌 套ave_profit=ave_p。

12、rofit,sum_profit/SIMUDAY; % 循 环 嵌 套endbuy_amount(1)=; % 第 一 个 元 素 置 空ave_profit(1)=;val,id=max(ave_profit) % 显 示 最 大 平 均 收 入 valbuy=buy_amount(id) % 显 示 在 平 均 收 入 最 大 情 况 下 的 每 天 的 购 买 量 buyxlabel=每 天 的 购 买 量 ;ylabel=平 均 利 润 ;plot(buy_amount,ave_profit,*:);% 函 数 GetProfit.m 代 码 :function re=GetProfi。

13、t(a,b)if ab % 供 不 应 求 : 报 童 购 买 量 小 于 销 售 量re=a*(0.05-0.03);else % 供 过 于 求 : 报 童 购 买 量 大 于 销 售 量re=b*(0.05-0.03)+(a-b)*(0.02-0.03);end 运 行 结 果 :val =4.2801 id =21 buy = 220 数 学 实 验 报 告 5 该 结 果 说 明 当 报 童 每 天 买 进 报 纸 数 量 为 220, 报 童 的 平 均 总 收 入 为 最 大 , 且 最 大 为 4.2801。三 、 结 果 分 析使 用 计 算 机 模 拟 法 解 决 卖 报 。

14、童 问 题 ,在 计 算 精 度 以 及 收 敛 速 度 问 题 上 均 能 达 到 一 个 较 好 的 效 果 ,能 有 效 地 解 决 此 问 题 。 计 算 机 模 拟 是 运 用 计 算 机 语 言 程 序 来 模 拟 现 实 世 界 , 为 随 机 性 库 存 问 题 提 供 了 较合 理 的 解 决 方 法 。 通 过 本 次 实 验 , 提 高 了 我 们 对 随 机 性 的 认 知 , 对 现 实 、 模 拟 、 理 论 推 导 等 的 了 解 。 学习 计 算 机 模 拟 的 基 本 过 程 与 方 法 , 会 做 简 单 的 计 算 机 模 拟 , 提 高 了 析 问 题 与 解 决 问 题 的 能 力 , 激 发 我 们对 于 探 寻 科 学 真 理 的 兴 趣 报 告 人 : 郭 舒报 告 时 间 : 2010年 12月 7日 。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值