M概率模型剖析.ppt

M09概率模型.PPT MATLAB 程序 %输入概率向量p, 零售价a,进价b,退价c,输出最佳购报数m及最大期望收入G function [m, G]=maibao(p,a,b,c) G=0; P=0; for m = 0 : inf P=P+p(m+1); %求和, p(m+1)是卖出m 份报纸的概率 D=(a-b)-(a-c)*P; %增加一份(第m+1份)报纸的期望收入 if D=0; return; end G=G+D; %购m+1份报纸的期望收入 end MATLAB求解程序 function m=zhagangex(L,sigma) % L 规定长度, sigma 均方差, m均值. (单位均为cm) lambda=l/sigma; z=fsolve(@(x) zhagang(x,lambda),-1.78); m=sigma*(lambda-z); %函数子程序 function f=zhagang(x,lambda) f=(1-normcdf(x,0,1))/normpdf(x,0,1)+x-lambda; 第9章实验 实验内容: 要设定粗轧钢材的长度m, 初轧钢材长度x服从均值为m, 均方差为σ的正态分布, m大于产品甲的长度L1, 若 L1 x≧L2, 则还可以生产长度为L2的副产品乙. 若xL2则报