2022年随机模型 .pdf
1随机模型在现实世界中 , 不确定现象是普遍存在的. 例如 , 漂浮在液面上的微小粒子不断地进行着杂乱无章运动, 粒子在任一时刻的位置是不确定的 ; 又如公共汽车站等车的人数在任一时刻也是不确定的, 因为随时都可能有乘客的到来和离去. 这类不确定现象, 表面看来无法把握 , 其实 , 在其不确定的背后, 往往隐藏着某种确定的概率规律, 因此 , 以概率和数理统计为基础的随机模型就成为解决此类问题最有效的工具之一. 依随机规律是否随时间的变化而变化, 随机时模型可分为静态和动态两类 , 前者只涉及到随机变量(向量 )的概率分布及其数字特征, 后者则要处理随机过程和随机微分方程, 本讲章主要讨论前者. 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页,共 38 页 - - - - - - - - - 21 电梯问题有r个人在一楼进入电梯, 楼上有 n 层, 设每个乘客在任何一层出电梯的概率相同, 试建立一个概率模型, 求直到电梯中的乘客下完时 , 电梯需停次数的数学期望. 分析: 对于此问题 , 容易想到该问题为离散随机变量求数学期望的问题 , 既然求电梯需停次数的数学期望, 那么电梯需停次数便为要定义的随机变量, 用 X 表示, 显然 X 的取值范围为,min,2 ,1nr, 然后需计算iXP=, 那么所求的答案即为,min1iXpinri=?=, 用上述思路求解电梯问题理论上完全正确 , 然而我们却很难得到问题的一个简洁的表达式结果 , 原因是古典概率iXP=的计算相当复杂 , 而且要依据r, n的不同情况具体来求. 下面我们看一些具体的例子, 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页,共 38 页 - - - - - - - - - 3例 1: 取2=r, 3=n时, 此时 X 的取值范围为2, 113211133CP X?=23222233CP X?=从而所求的解为35322311=?+?例 2: 取3=r, 2=n时, 此时 X 的取值范围为2, 1123111 24CP X?=43222233=-=XP从而所求的解为名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页,共 38 页 - - - - - - - - - 447432411=?+?例 3: 取3=r, 4=n时, 此时 X 的取值范围为3 ,2, 1143111416CP X?=3433!63416CP X?=313443413!92416CCP X-? -?=从而所求的解为1637166316921611=?+?+?上述三例是在r, n给出具体数据的情况下的计算, 可以看出 , 随着r, n数据的增大 , 计算变的愈加复杂 , 且没有明显的规律可言 . 而当r, n未给出具体数据时 , 用上述思路求名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页,共 38 页 - - - - - - - - - 5解问题 , 想要得到具体的表达式就更为困难. 下面我们换个角度考虑该问题, 我们将电梯在第i层是否停下来这一事件作为随机变量iY , 1=iY表示停下来, 0=iY表示电梯未停 , 其中i取值为,2, 1n, 这样问题便转化为求iY 的期望之和 , 由题意容易得知电梯在任何楼层上是否停留这一概率完全相同即?-=rrinnnnY)1(0)1(11从而iY 的期望即为rinnYE)1(1)(-=那么原问题的解即为名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页,共 38 页 - - - - - - - - - 6)1(1)(1rniinnnYE-?=显然要比刚才的方法来得简单, 而且得到了统一的表达式结果 , 避免了用第一种方法在计算古典概率时对r, n 大小的具体讨论 . 我们用前面的三个例子验证上面的结论: 例 1 中, 取2=r, 3=n时, 得到的结果为35, 而35)313(132=-?例 2 中, 取3=r, 2=n时, 得到的结果为47, 而47)212(123=-?名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页,共 38 页 - - - - - - - - - 7例 3 中, 取3=r, 4=n时, 得到的结果为1637, 而1637)414(1 43=-?通过电梯问题求解的讨论, 可以看出在解决带随机性现象的问题中 , 方法的选取是非常重要的, 只有采取合适的方法才会事半功倍 . 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页,共 38 页 - - - - - - - - - 82 钓鱼问题为了估计湖中鱼的数量, 先从湖中钓出r条鱼做上记号后又放回湖中 , 然后再从湖中钓出S条鱼 , 结果发现 S条鱼中有 x条鱼标有记号 . 问应该如何估计湖中鱼的数量N? 分析与求解该问题就是要从第二次钓出的标有记号的鱼所占的比例估计出湖中鱼的数量 . 首先我们假设放回湖中的鱼在湖中的分布是均匀的 . 则第二次钓出的标有记号的鱼数X 是一个随机变量 , X 服从超几何分布xs xrNrsNCCP XxC-?=(*) 其 中 x 为 整 数 , 且,min)(,0maxsrxrNs-. 用),(NxL表示(*) 式的右端 , 则取使),(NxL达到极大值的N作为N的估计量 . 直接对N求导考察极值比较困难 , 我们用比名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页,共 38 页 - - - - - - - - - 9值法来研究),(NxL的变化( ,)( ,)( ,1)()()L x NNrNsA x NL x NNNrsx-=?-22()()Nrs NrsNrs NNx-+=-+(*) 从(*) 式看出 , 当且仅当xrsN NxLNxL, 而当且仅当xrsN 时 , )1,(),(-时, 求得*11ln(19118 )ln22plll= -+ (3.8) 将120.0007,0.004pp=代入(3.8)式得到*0.49p =(百万元 ). 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页,共 38 页 - - - - - - - - - 174 报童的策略背景简介报童问题即单周期库存问题(Single-Period Problem), 是供应链管理中最重要的模型之一, 其历史可以追溯到 1888年 著名 经济 学家 Edgeworth应 用 它 解决 银 行的 现金 流(cash-flow)问题, 1955 年, Whitin首次建立了受价格影响的报童问题模型 . 目前 , 报童问题在生产、服务、管理和金融等领域成功地取得了广泛的应用. 现在报童问题衍生出许多扩展问题 . 与众多扩展模型相比 , 经典报童问题模型是最简单最基本的问题 , 它可以描述为 : 报童每天早晨以单位批发价b 从报社买进报纸 , 然后以单位零售价a出售 , 晚上将没有卖掉的报纸当作废品以价格c(cba)处理掉 . 同时假设 : (1) 报童拥有购买足够多报纸的资金; (2) 报纸过剩只能以低于零售价的价格v处理; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页,共 38 页 - - - - - - - - - 18(3) 报纸供应不足 , 会遭受缺货惩罚 ; (4) 不计其他费用 (如交通费、摊位费等 ). 报童应该如何确定订购量而获得最高的利润呢? 显然 , 报童应该根据市场需求量来确定订购量, 而市场需求量是随机的. 假设报童通过经验已经掌握了市场需求量的随机规律, 我们就可以建立随机优化模型来求解报童问题了. 报童每天清晨从报社购进报纸零售, 晚上将没有买掉的报纸退回 . 每份报纸的购进价为元b, 零售价为 a元, 退回价为 c 元,. 报童售出一份报纸赚ba -元, 退回一份报纸赔cb -元. 报童每天如果购进的报纸太少, 不够卖时会少赚钱 , 如果购进的报纸太多, 卖不完时会赔钱 . 试为报童筹划每天应如何确定购进的报纸数使得收益最大? 模型一问题分析报童应该根据需求量确定购进量, 而需求量是随机的, 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页,共 38 页 - - - - - - - - - 19所以这是一个风险决策问题. 假定报童已经通过自己每天的卖报经验或其它渠道掌握了需求量的分布规律, 需求量r为一连续型随机变量, 密度函数为)(rf,假设每天的购进量为n, 由于需求量r是随机的 , 可以小于 n、 等于 n或大于 n, 这就导致报童每天的收入也是随机的, 所以作为优化模型的目标函数 , 不能是报童每天的收入函数, 而应该是他长期卖报的日平均收入 . 从概率论大数定律的观点看这相当于报童每天收入的期望值 , 以下简称它为平均利润. 模型建立显然 , 若 nr, 则以 a 价售出r份报纸 , 以 c 价售出nr-份报纸 , 若nr, 则全部 n以 a价售出 , 故平均利润为0( )()()()( )nF nab rbc nrf r dr=-()( )nab nf r dr+-0()()() ( )nab nacnr f r dr=-名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页,共 38 页 - - - - - - - - - 200( )()()( )ndF nabacf r drdn=-由一元函数极值存在的必要条件可知( )0dF ndn=, 得0( )nabf r drac-=-而22( )()( )0d F nac f ndn= -. 从而知满足上式的n 可以使平均利润达到最大 . 易知上式等价于0( )( )nnf r drabbcf r dr+-=-上式左边是报童订购n份报纸时 , 不能将它卖完的概率与能将它卖完的概率之比, 右边则表示卖出一份报纸的盈利与退回一份报纸亏损之比, 该式表明最优购进量是使这两个名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页,共 38 页 - - - - - - - - - 21比相等的购进量 . 实例 : 如3 .0,6 .0, 1=cba, 需求量r服从正态分布)10,100(2N, 则当报童的订报量 n满足0( )43( )nnf r drabbcf r dr+-=-或等价于04( )7nabf r drac-=-不难计算得到102n =时长期平均收益最大 . 若3 .0,6 .0,2 .1=cba此 时02( )3nabf r drac-=-, 104n =时长期平均收益最大 . 若2.0,6.0,8.0=cba此 时01( )3nabf r drac-=-, 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页,共 38 页 - - - - - - - - - 2296n =时长期平均收益最大 . 再如, ab-=0.3, bc-=0.1, 需求量r为服从 2000,4000上均匀分布的连续型随机变量, 其密度函数为1 200020004000( )0 xf r ?=?其它由0( )nabf r drac-=-知, 2000( )nf r drabac-=-=0.30.30.1+=0.75, 于是 , 有 n-2000=1500, 因此 , 报童的最优策略是订购3500份. 模型二(报童问题 )报童每天要到邮局去订报, 出售一份报纸可获得利润()Aab=-元, 但如卖不出退回邮局, 每份报纸要名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页,共 38 页 - - - - - - - - - 23损失()Bbc=-元. 根据以往经验 , 得知每天需求量为k 份的概率为kp . 问报童每天应订购多少份报纸, 才能使它获利的期望值最大 . 如果卖报童问题中的顾客每天需求量X 是一个离散型随 机 变 量 , 设 报 童 每 天 订 购 的 份 数 为 n份 , 于 是 有()kP Xkp=, 记出售一份报纸可获得利润()Aab=-(元) 退回邮局一份报纸要损失()Bbc=-(元) 则报童每天的利润()f X可用下列公式来表示 : ()()AnXnf XAXnX BXn?=?-, 可解得售出概率p应满足下述不等式BpAB+(4.2) 其中100(1)()1()1nniiniipP XnP XiP Xip= += +=-=-. 即001nniiiiBAppABAB=-?+, (4.3) 其中BAB+称为临界值 . 那么, 最好的 n 应满足名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 25 页,共 38 页 - - - - - - - - - 26(1)( ), (1)( )c nc n c nc n-+即100nniiiiAppAB-=+例如, 已知某种报纸每天需求量Nj的概率分布如下 : 表 1 需求量 Nj0 1 2 3 4 5 6 7 8 P(Nj) 0.05 0.10 0.100.250.200.150.050.05 0.05 每出售一份报纸 , 可获利4 角; 若当天卖不掉 , 每份报纸将损失 3 角. 试问每日应进多少份报纸 ? 因为我们引进了需求量随机变量X, 所以我们将表1 作一点修改并就在表中进行计算, 如表 2: 表 2 n 0 1 2 3 4 5 6 7 8 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 26 页,共 38 页 - - - - - - - - - 27P(X =n) 0.05 0.10 0.100.250.200.150.050.05 0.05 P(X n) 0.05 0.15 0.250.500.700.850.900.95 0.10 P(X n+1) 0.95 0.85 0.750.500.300.150.100.05 0.00 现在本问题中, A=4, B=3, 所以BAB+0.43, 在 n=3 时, P(X 3+1)= 0.50 满足式 (3), 所以最佳订购量 n*=4 份. 该模型也可以理解为单周期随机库存问题, 即假定在一个周期末库存的货物对下一个周期没有任何价值, 即模型适用于仅有一次机会存贮以供需求的产品如时装、新鲜食品、月饼等 . 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 27 页,共 38 页 - - - - - - - - - 285 最佳订票问题一问题提出在激烈的市场竞争中, 航空公司为争取更多的客源而开展的一个优质服务项目是预订票业务. 公司承诺 , 预先订购机票的乘客如果未能按时前来登机, 可以乘坐下一班飞机或退票 , 无需附加任何费用. 当然也可以订票时只订座, 登机时才付款 , 这两种办法对于下面的讨论是等价的. 设某种型号的飞机容量为n, 若公司限制预定n张机票 , 那么 , 由于总会有一些订了机票的乘客不按时来登机, 致使飞机因不满员飞行而利润降低, 甚至亏本 , 如果不限制订票数量呢 , 那么当持票按时前来登机的乘客超过飞机容量时, 必然会引起那些不能登机飞走的乘客(以下称被挤掉者)的抱怨 . 公司不管以什么方式予以补救, 都会导致受到一定的经济损失 , 如客源减少 , 或挤到以后班机的乘客, 公司要无名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 28 页,共 38 页 - - - - - - - - - 29偿供应食突宿或者付给一定的赔偿金等. 这样, 综合考虑公司的经济利益 , 必然存在一个恰当的订票数量和限额. 假设飞机容量为300, 乘客准时到达机场而未乘上飞机的赔偿费是机票价格的%10, 飞行费用与飞机容量、机票价格成正比(由统计资料知 , 比例系数为 0.6, 乘客不按时前来登机的概率为 0.03), 请你: (1)建立一个数学模型 , 给出衡量公司经济利益和社会声誉的指标 , 对上述预定票业务确定最佳的预定票数量. (2)考虑不同客源的不同需要, 如商人喜欢上述这种无约束的预定票业务 , 他们宁愿接受较高的票价; 而按时上下班的雇员或游客 , 愿意以若不能按时前来登机, 则机票失效为代价 , 换取较低额的票价. 公司为降低风险 , 可以把后者作为基本客源 . 根据这种实际情况 , 制定更好的预订票策略 . 二.模型的假设及符号说明1、模型的假设名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 29 页,共 38 页 - - - - - - - - - 30 假设预订票的乘客是否按时前来登机是随机的. 假设已预订票的乘客不能前来登机的乘客数是一个随机变量 . 假设飞机的飞行费用与乘客的多少无关. 2、符号说明n: 飞机的座位数 , 即飞机的容量 ; g: 机票的价格 ; f : 飞行的费用 ; b: 乘客准时到达机场而未乘上飞机的赔偿费; m: 售出的机票数 ; k: 已预订票的乘客不能前来登机乘客数, 即迟到的乘客数 , 它是一个随机变量 ; kp : 已预订票的m个乘客中有k个乘客不能按时前来登机的概率 ; p : 每位乘客迟到的概率 ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 30 页,共 38 页 - - - - - - - - - 31( )mPj: 已预订票前来登机的乘客中至少挤掉j 人的概率 , 即社会声誉指标 ; S: 公司的利润 ; ES: 公司的平均利润 . 三.问题的分析及数学模型1、问题的分析通过上面引进的符号易知, 赔偿费gb1.0=, 习行费用ngf6. 0=, 每位乘客迟到概率03.0=p, 已预订票的m个乘客中 , 恰有k个乘客不能按时前来登机, 即迟到的乘客数k服从二项分布()pmB, 此时, (1)(0,1,2,)kkm kkmpC ppkm-=-=当nkm-时, 说明km -个乘客全部登机 , 此时利润fgkmS-=)(当nkm-时, 说明有n个乘客登机 , 有nkm-个乘名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 31 页,共 38 页 - - - - - - - - - 32客没有登上飞机 , 即被挤掉了 , 此时利润bnkmfngS)(-=根据以上的分析 , 利润S可表示为 : ?-=)()()()(nmknkmbnkmfngnmknkmfgkmS迟 到 的 乘 客 数1,2, 1 , 0-=nmk时 , 说 明 有nkm-个乘客被挤掉了 ; 迟到的乘客数mnmnmk, 1,+-=时, 说明已来的km -个乘客全部登机了 . 于是平均利润-=-=-+-=mnmkknmkkpfgkmpbnkmfngES)()(10因为名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 32 页,共 38 页 - - - - - - - - - 33-=-=-=-=-=+-=-=-101010010)()()()()1)()(nmkknmkknmkkmkknmkkmnmkkkpgkgEpfmgfmgkpkpgpfmgpfgkm所以101100() )()()( )m nkkm nm nkkkkESngfmkn b pmgfmgfpgE kgkp-=-=-+-+1010( )()()( )()()m nkkm nkkmE kgfngfmkn bmgfgk pmE k gfbgm kn p-=-=-+-+=-+-由于kmkkmkppCppnBk-=)1(),(, 可知, 随机变量k名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 33 页,共 38 页 - - - - - - - - - 34的数学期望mpkE=)(, 此时, -=-+-=10)1()()()1 (nmkkmkkmppCnkmgbfmgpES2、数学模型通过以上对问题的分析, 可以在一定的社会声誉指标( )mPj范围内 , 寻求合适的m, 根据Ngf6.0=的关系 , 使得目标函数fES达到最大 , 即10101max(1)(1)()(1)10.610.971.1(300)(1)1180m nkkm kmkm nkkm kmkESbp mmkn C ppfNgmmkC pp-=-=-+-=-下面考虑社会声誉指标 . 由于jknm+=, 所以jnmk-=, 即当被挤掉的乘客数为 j 时, 等价的说法是恰有jnm-个迟到的乘客 . 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 34 页,共 38 页 - - - - - - - - - 35公司希望被挤掉的乘客人数不要太多, 被挤掉的概率不要太大 , 可用至少挤掉j 人的概率作为声誉指标, 相应地k的取值范围为jnmk-=,2, 1 , 0, 社会声誉指标0()(1)m njkkm kjmkP mC pp-=-四、模型求解为 了 对 模 型 进 行 求 解 , 可 以 分 别 给 定m, 比 如,350,306,305=m, 计算fES/, 同时 , 给定 j , 比如 取5=j, 计算社会声誉指标( )mPj, 从中选取使fES/最大 , 且社会声誉指标( )mPj小于等于某个 (比如取05.0=)最佳订票数m. 下面给出 MATLAB计算程序 . % 飞机最佳订票策略ch43 % 文件名 : ch43.m 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 35 页,共 38 页 - - - - - - - - - 36%m 表示售出的票数; Es 表示平均利润 ; p 表示声誉指标 ; for m=305: 325 sm=0; p=0; for k=0;m-305 pp=(prod(m-k+1:m)/prod(1:k)*0.03k*0.97(m-k); p=p+pp; sm= sm+(m-k-300)*pp/prod(1:k); end Es= (1/180)*0.97*m-1.1*sm-1; m Es p end 执行后可输出以下结果: m ES P 305 0.6436 9.2338e-005 306 0.6490 9.3723e-004 307 0.6543 0.0048 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 36 页,共 38 页 - - - - - - - - - 37308 0.6596 0.0167 309 0.6649 0.0442 310 0.6703 0.0952 311 0.6756 0.1742 312 0.6810 0.2796 313 0.6864 0.4028 314 0.6917 0.5314 315 0.6971 0.6525 316 0.7024 0.7566 317 0.7078 0.8388 318 0.7132 0.8890 319 0.7185 0.9399 320 0.7239 0.9661 321 0.7293 0.9818 322 0.7347 0.9907 323 0.7400 0.9954 324 0.7454 0.9979 325 0.7508 0.9990 从 计 算 结 果 易 见 , 当m=309 时 , 社 会 声 誉 指 标名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 37 页,共 38 页 - - - - - - - - - 38()53090.04420.05,p=所以为了使尽/ESf 量大, 且要满足社会声誉指标( )50.05,pm 则最佳订票数可取m=309. 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 38 页,共 38 页 - - - - - - - - -