数学建模实验答案-概率模型.docVIP

PAGE

PAGE4

实验10概率模型（2学时）

（第9章概率模型）

1.（验证）报童的诀窍p302~304,323（习题2）

关于每天报纸购进量的优化模型：

已知b为每份报纸的购进价，a为零售价，c为退回价（abc），每天报纸的需求量为r份的概率是f(r)(r=0,1,2,…)。

求每天购进量n份，使日平均收入，即

达到最大。

视r为连续变量，f(r)转化为概率密度函数p(r)，则所求n*满足

已知b=0.75,a=1,c=0.6，r服从均值=500（份），均方差=50（份）的正态分布。报童每天应购进多少份报纸才能使平均收入最高，这个最高收入是多少？

[提示：normpdf,normcdf]

(i)计算正态变量的概率密度函数的调用形式为：Y=normpdf(X,mu,sigma)

正态变量的概率密度函数为

其中：

X是x的一组值，Y对应一组函数值。

mu为，sigma为。

当=0，=1时，为标准正态变量的概率密度函数。

(ii)计算正态变量的分布函数的调用形式为：P=normcdf(X,mu,sigma)

正态变量的分布函数为

且

标准正态变量的概率密度函数对应标准正态变量的分布函数。

要求：

(1)在同一图形窗口内绘制和的图形，观察其交点。

程序：

n=500:530;

mu=500;sigma=50;

y1=normcdf(n,mu,sigma)-normcdf(0,mu,sigma);

a=1;b=0.75;c=0.6;

y2=(a-b)/(a-c)*ones(size(n));

plot(n,[y1;y2]);

gridon;

[提示]，

☆(1)运行程序并给出结果：

(2)求方程的根n*（四舍五入取整），并求G(n*)。

程序：

functiony=fun(n)

mu=500;sigma=50;

a=1;b=0.75;c=0.6;

y=normcdf(n,mu,sigma)-normcdf(0,mu,sigma)-(a-b)/(a-c);

(3)在同一图形窗口内绘制和的图形，观察它们的交点。（参考题1的(1)）

★(3)给出程序及运行结果（比较[309]图2）：

z=-2:0.1:2;

y1=(1-normcdf(z,0,1))./normpdf(z,0,1);

l=2;sigma=0.2;

y2=l/sigma-z;

plot(z,[y1;y2]);

gridon;

(4)求方程的根z*，并求m=l-σz*。（参考题1的(2)）

提示：由(3)得到的图形可观察到z*的大概位置。

★(4)给出程序及运行结果（比较[310]）：

functiony=fun(z)%方程

l=2;sigma=0.2;

y=l/sigma-z-(1-normcdf(z,0,1))./normpdf(z,0,1);

3.（验证）航空公司的预订票策略p313~316

模型如下：

给定λ,n,p,b/g，求m使单位费用获得的平均利润J(m)最大。

约束条件为

其中：

m 预订票数量的限额。

λ(1) 利润调节因子。

n 飞机容量。

p 每位乘客不按时前来登机的概率，q=1–p。

b 每位被挤掉者获得的赔偿金。

g 机票价格。

b/g 赔偿金占机票价格的比例。

不按时前来登机的乘客数K服从二项分布，其概率为

被挤掉的乘客数超过j人的概率为

（等价于m位预订票的乘客中不按时前来登机的不超过m–n–j–1人）

该模型无法解析地求解，我们设定几组数据，用程序作数值计算。

[提示：binopdf,binocdf]

(i)二项分布的概率密度函数：Y=binopdf(X,N,P)

计算X中每个X(i)的概率密度函数，其中，N中对应的N(i)为试验数，P中对应的P(i)为每次试验成功的概率。Y,N,和P的大小类型相同，可以是向量、矩阵或多维数组。输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。

N中的值为正整数，P中的值从[0,1]取。

已知x和参数n,p，二项分布概率密度函数为

q=1–p。y为n次独立试验中成功x次的概率，其中，每次试验成功的概率为p。x=0,1,...,n。

(ii)二项式累积分布函数：Y=binocdf(X,N,P)

计算X中每个X(i)的二项式累积分布函数，其中，N中对应的N(i)为试验数，P中对应的P(i)为每次试验成功的概率。Y,N,和P的大小类型相同，可以是向量、矩阵或多维数组。输入的标量将扩展成一个数组，使其大小类型与其它输入相一致。

N中的值为正整数；X中的值从[0,N]取；P中的值从[0,1]取。