数学建模案例分析2 随机存储模型--概率统计方法建模
数学建模概率模型

2
记为X ~ N(, 2 )
背景:如果决定试验结果X的是大量随机因素的总和,假设
各个因素之间近似独立,并且每个因素的单独作用相对均匀 地小,那么X的分布近似正态分布。
如:同龄人的身高、体重、考试分数、某地区年降水量等。
3、数学期望的概念和计算 描述了随机变量的概率取值中心—均值
数学期望
Y gX
E( X ) xk pk k 1
E( X ) xf ( x)dx E(Y ) EgX g( xk ) pk k 1
E(Y ) Eg( X )
g( x) f ( x)dx
4、MATLAB中相关的的概率命令
常见的几种分布的命令字符为: 正态分布:norm 指数分布:exp 泊松分布:poiss 二项分布:bino
G(n)
n
0
[(
a
b)r
(b
c)(n
r
)]
p(r
)drຫໍສະໝຸດ n(ab)np(r
)dr
dG (a b)np(n)
n
(b c) p(r)dr
dn
0
(a b)np(n) n (a b) p(r)dr
n
(b c)0 p(r)dr (a b)n p(r)dr
例3 有10台机床,每台发生故障的概率为0.08,而10台机床工作 独立,每台故障只需一个维修工人排除.问至少要配备几个维修 工人,才能保证有故障而不能及时排除的概率不大于5%。
解:随机变量X示发生故障的机床的台数,则 X ~ B(10,0.08)
即P{X n} 0.95
数学建模—概率模型 ppt课件

数学建模—概率模型
v3统计图(examp05-03) v箱线图(判断对称性) v频率直方图(最常用) v经验分布函数图 v正态概率图(+越集中在参考线附近,越近似正态分布)
v4分布检验 vChi2gof,jbtest,kstest,kstest2,lillietest等 vChi2gof卡方拟合优度检验,检验样本是否符合指定分布。它把观测数据分 组,每组包含5个以上的观测值,根据分组结果计算卡方统计量,当样本够 多时,该统计量近似服从卡方分布。 vjbtest,利用峰度和偏度检验。
3 单因素一元方差分析步骤
( example07_01.m 判断不同院系成绩均值是否相等)
数据预处理
正态性检验 lillietest (p>0.05接受)
方差齐性检验 vartestn (p>0.05接受)
方差分析
anoval (p=0 有显著差别)
多重比较:两两比较,找出存在显著差异的学院,multcompare
构造观测值矩阵,每一列对应因素A的一个水平,每一行对应因素B的一个
水平
方差分析
anova2 得到方差分析表
方差分析表把数据差异分为三部分(或四部分): 列均值之间的差异引起的变差 列均值之间的差异引起的变差 行列交互作用引起的变差 (随机误差) 后续可以进行多重比较,multcompare,找出哪种组合是最优的
Computer Science | Software Engineering & Information System
数学建模—概率模型
目的:用一个函数近似表示变量之间的不确定关系。 1 一元线性回归分析 做出散点图,估计趋势;计算相关系数矩阵; regress函数,可以得到回归系数和置信区间,做残差分析,剔除异常点,重 新做回归分析 Regstats 多重线性或广义回归分析,它带有交互式图形用户界面,可以处 理带有常数项、线性项、交叉项、平方项等模型 robustfit函数:稳健回归(加权最小二乘法)
概率统计建模方法

第1章概率方法建模简介第2章数据统计描述和分析第3章方差分析第4章回归分析第5章马氏链模型第6章时间序列模型第7章主成分分析及应用第8章判别分析简介及应用主讲:山东大学数学学院陈建良2第1章概率方法建模简介随机性模型,是指研究的对象包含有随机因素的规律,以概率统计为基本数学工具,其结果通常也是在概率意义下表现出来。
随机因素的影响可以用概率、平均值(即数学期望)等的作用来体现。
自然界中的现象总的来说可以概括为两大现象:确定性现象和随机现象在确定性现象中可以忽略随机因素的影响,在随机现象中必须考虑随机因素的影响。
确定性离散模型,主要使用差分方程方法、层次分析方法以及比较简单的图的方法和逻辑方法等方法建立模型;确定性连续模型,主要使用微积分、微分方程及其稳定性、变分法等方法建立模型;§2 概率方法建模实例分析实例一、报童的策略问题1.问题描述报童每天清晨从报站批发报纸零售,晚上将未卖完的报纸退回。
设每份报纸的批发价为b,零售价为a,退回价为c,且设a>b>c,因此报童每售出一份报纸赚(a-b),退回一份赔(b-c)。
若批少了不够买就会少赚,若批多了买不完就赔钱,报童如何确定每天批发报纸的数量,才能获得最大收入?92. 分析显然应根据需求量来确定批发量。
一种报纸的需求量是一随机变量。
假定报童通过自己的实践经验或其它方式掌握了需求量的随机规律,即在他的销售范围内每天报纸的需求量为X = x 份的概率为P(x),则通过P(x) 和a, b, c 就可建立关于批发量的优化模型。
3.数学模型设每天批发量为n,因需求量x 是随机的,因此x可以小于、等于或大于n,从而报童每天的收入也是随机的,作为优化模型的目标函数,应考虑他长期(半年、一年等)卖报的日平均收入。
据概率论中的大数定律,这相当于报童每天收入的期望值(以下简称平均收入)。
1011设报童每天批发进n 份报纸时的平均收入为S (n ),若某天需求量x ≤n ,则他售出x 份,退回(n -x )份;若这天需求量x >n ,则n 份报纸全部卖出。
概率统计方法建模PPT课件

第3页/共23页
5.5 随机状态转移模型
状态与状态转移 ➢随机变量Xn:第n年的状态 状态概率 ai (n)
Xn
1, 2,
第n年健康 第n年疾病
ai (n) P(Xn i), i 1, 2, n 0,1,
➢今年处于状态i, 来年处于状态j的概率 pi:j 转移概率
存贮策略是周末库存量为零时订购3架 周末的库存量可 能是0, 1, 2, 3,周初的库存量可能是1, 2, 3。 用马氏链描述不同需求导致的周初库存状态的变化。 动态过程中每周销售量不同,失去销售机会(需求超过 库存)的概率不同。
可按稳态情况(时间充分长以后)计算失去销售机会的 概率和每周的平均销售量。
马氏链的两个重要类型
设状态i是非吸收状态,j是吸收状态,则首达概率f ij (n) 实际上是i经n次转移被j吸收的概率。而
fij = fij (1) + fij(2) + … + fij(n) + …
则是从非吸收状态i出发终将被吸收状态j吸收的概率。 记 F={f ij} 则 F=MR
例如,可以算出前面第二种情况中
第19页/共23页
5. 6 马尔可夫链的应用模型
模型求解 ➢ 估计这种策略下每周的平均销售量
第n周平均售量Rn
需求不超过存 量,销售需求
需求超过存量, 销售存量
3i
Rn [ jP(Dn j, Sn i) iP(Dn i, Sn i)] i1 j 1 3i [ jP(Dn j Sn i) iP(Dn i Sn i)]P(Sn i) i1 j 1
p23 p33
P(Dn k) e1 / k ! (k 0,1, 2 )
概率统计方法建模讲座

f = 94.5863
f = 94.7597
再算一次!!
每天平均的收益的模拟值
随机模拟求最佳订购数量
clear a=2.0;b=1.0;c=0.5;M=[50:150]'; k=length(M);A=zeros(k,1); for J=1:k N=5000;G=zeros(N,1);m=M(J); X=poissrnd(100,N,1); for I=1:N if X(I)>=m G(I)=(a-b)*m; else G(I)=(a-b)*X(I)-(b-c)*(m-X(I)); end end A(J)=mean(G); end [Y,I] = max(A) n=M(I) Y
案例1:如何估计池塘中鱼的数量
一、问题:要估计一个池塘里有多少条鱼,可以采用“标志重捕法”,即:先重池塘中捕出r条,每条鱼都做上记号,经过一段时间后,再从池塘中捕出s条(s>r),统计其中标有记号的鱼的条数t,利用这些信息,估计池塘中鱼的条数N.
需要作哪些假设?
1、实验期间,标记个体不会变化。 2、标记不会对鱼造成伤害。 3、期间没有迁出、迁入、新生和死亡。 4、所有鱼被捕获的概率相等。
随机性存储模型.pptx

Pi
xi s
C2 K C1 C2
(10.34)
成立的最小S值。订货量为Q*=S*—I
【例10.10】假定在例中月初已有5件产品,问还应该订货多少件。 【解】由例10.9解的结果得知,月初存量为零时应订货S*=20件,在引当初存量I=5件时, 还应该订货Q*=20-5=15件。
2024年9月30日星期一
【解】已知C1=100,C2=180,K=100,泊松分布函数为
P(x) x e , x!
x=0,1,2,…
平均需求量为5,则λ=5,由式(10.31)得
C2 K 180 100 0.2857 , C1 C2 100 180
4
Q
P(x)
5 x e 5
xQ
x0 x!
查泊松分布表,当Q=4时,
0.20
0.35
0.25 0.08
报童每售出100张赚15元,如果当天末能售出,每百张赔7元,问报童每天应准备多少 份报纸最佳。
2024年9月30日星期一
第2页/共13页
【解】以每天为一个时期,当订货量大于需求量时,存储费为C1=7,当订货量小于 需求量时,因缺货而失去赚钱的机会,因此缺货费为C2=15,不计订货费,利用式 (10.32)可得
第8页/共13页
二、连续型存储模型 离散型存储策略的分析方法同样适合连续型。设需求量x的概率密
度为
(x), 满足
(x)dx 1,x 0
0
1.期初存量为零的情形 当x≤Q时,总存储费期望值为
当 x>Q时,总缺货费期望值为
订货费为C3+KQ,总费用期望值为
Q
C1 0 (Q x)(x)dx
C2 Q (x Q) (x)dx
数学建模方法之概率统计分析法

Obs
Prin1 Prin2 Prin3 Prin4 Prin5 Prin6 1 -0.38118 -0.32367 -0.04450 0.30363 0.00430 0.06437 2 0.57795 -0.35416 0.49279 0.55119 -0.18726 0.17414 3 0.69219 -0.21588 0.40557 0.40041 -0.10461 0.05393 4 0.22635 -0.39419 0.27521 0.63296 0.13851 -0.06481 5 -0.82981 -0.40293 0.47330 -0.42964 -0.55401 -0.35020 6 -1.19410 -0.40627 -0.36848 0.14000 0.02221 0.01063 7 -1.63568 -0.26394 -0.67179 -0.15189 0.01702 -0.03769 8 0.95195 -0.46156 1.61851 -0.92520 0.08394 0.25530 9 0.46501 -0.14888 0.19070 0.16273 -0.30327 0.20883 10 -1.45693 -0.18670 -0.55658 -0.17088 -0.10267 -0.00922 11 -0.29401 3.71727 -0.02727 -0.02382 -0.06419 0.03517 12 0.08041 0.22542 1.71694 0.12718 0.45539 -0.26668 13 -2.11628 -0.16312 -0.90179 -0.16784 0.14422 -0.03334 14 -0.94513 -0.31477 -0.39513 0.09760 0.11375 -0.03132 15 6.74015 -0.06989 -1.12895 -0.16618 0.04080 -0.11394 16 -0.88090 -0.23673 -1.07853 -0.38025 0.29589 0.10482
数学建模概率模型案例

则传送系统效率为:d=s/n=mp/n
=
m[1(1 1)n]
n
m
mn Dm [1(1nn(n1)) ]1n1
n m 2m 2
2m
D87 .5% 当n=10,m=40
报童的诀窍
问题:报童每天清晨从报社购进报纸零售,晚上 将没有卖掉的报纸退回。设报纸每份的购进价为b, 零售价为a,退回价为c,假设a>b>c。即报童售出
每位被挤掉的乘客获得的赔偿金为常数b。
4 模型建立
先不考虑社会声誉的影响。
公司的经济利益用平均利润(数学期望)S 来衡量
订票的总人数是 m,m有可能超出 n
当有 k个人误机时,
航空公司可能从航班中得到的利润为
s m kg r,
m k n
s n g r (m k n )b , m k n
E ( X )x ip i ( i 1 ,2 , ,n )
连续型随机变量 X 的概率密度函数为 f ( x) 则随机变量 X 的数学期望值为
E(X) xf(x)dx
期望值反映了随机变量取值的“平均”意义!
传送系统的效率
在机械化生产车间里,你可以看到这样的 情景:排列整齐的工作台旁工人们紧张的 生产同一种产品,工作台上方一条传送带 在运转,带上若干个钩子,工人们将产品 挂在经过他上方的钩子上带走,当生产进 入稳态后,请大家构造一个衡量传送系统 效率的指标,并建立模型描述此指标与工 人数量、钩子数量等参数的关系。
mnj1
minPj(m) Pk k0
mJ (a m ) x S r 0 .1 6 n p m 1 b g m k n 0 1 P km n k 1
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
§2 随机存储模型
模型一、销售量为随机的存储模型
报童每天清晨从报社购进报纸零售,晚上将没有卖出的报纸退回。
如果购进报纸太少不够卖,会少赚钱;如果购进太多买不完,将要赔钱。
报童应如何确定每天购进的报纸数量,以求获得最大的收入。
模型假设1、报纸每份购进价b ,零售价a ,退回价c ,且c b a >>
2、市场需求量是随机的,报童已通过经验掌握了需求量r 的随机规律,r 视为连续随
机变量,其概率密度函数)(r p 。
模型建立 记 n —每天购进量,报童每天的收入R 是n 的函数
()()()()()⎩
⎨⎧>----≤-=r n r n c b r b a r n n b a n R ,, 但目标函数不应是报童每天的收入,而应是他长期卖报的日平均收入。
从大数定律的观点看,这相当于每天收入的期望值,即日平均收入:
()()()()[]()()()⎰⎰∞-+----=n n dr r p n b a dr r p r n c b r b a n G 0 ()()()()()()()()⎰⎰∞-+-----=n n dr r p b a n np b a dr r p c b n np b a dn dG 0 ()
()()()⎰⎰∞-+--=n n
dr r p b a dr r p c b 0 令0=dn
dG ,得到 ()()c b b a dr r p dr
r p n n
--=⎰⎰∞
又因为()10=⎰∞
dr r p ,上式又可表示为 ()c
a b a dr r p n
--=⎰0 (1) 使报童平均日收入最大购进量n 由(1)确定
评注 由()()c b b a dr r p dr r p n
n --=⎰⎰∞0,()⎰=n
dr r p p 01是卖不完的概率, ()⎰∞
=n dr r p p 2是卖完的概率。
上式表明,购进的份数应使卖不完与卖完的概率之比等于卖出一份赚的钱b a -与退回一份赔的钱c b -之比。
模型二、到货时间为随机的存储模型
模型假设1、商品订货费1c ,每件商品单位时间的储存费为2c ,缺货费3c ,单位时间需求量为r ;
2、当储存量降至L 时订货,订货量使下周期初的储存量达到固定值Q ;
3、交货时间x 是随机的,如下图中的,...,21x x ,设x 的概率密度函数()x p 。
模型建立
为使总费用最小,选择合适的目标函数建立模型,确定最佳订货点L 。
t
由储存量()t q 的图形可写出一个订货周期内的储存量和缺货量分别为
()()⎪⎪⎩⎪⎪⎨⎧≥<--=r L x r
Q r L x r rx L Q t q ,2,2222 ()()⎪⎪⎩⎪⎪⎨⎧≥-<=r L x r L rx r L x t q ,2,02 于是得到一个订货周期的平均费用为
()()()()()⎰⎰∞⎥⎦⎤⎢⎣⎡-++--+=r
L r L dx x p r L rx c r Q c dx x p r rx L Q c c L c 0
23222221222 目标函数应取为单位时间的平均费用()L S ,由于订货周期的平均长度为
()()x E r L Q L T +-=
这里()()⎰+∞
=0dx x xp x E
所以 ()()()
L T L c L S = 由0=dL
dS ,可以解出最佳订货点*L 满足方程 ()()
()x rE Q L c L c L +='-**
*。