用MATLAB实现报童模型,求大神帮我看下程序,万分感谢!!
我在做“需求信息对零售商订货策略影响”的论文,用报童模型,用MATLAB编了出来,但是验算出来的和答案不同。求大神帮我看看,很急,万分谢谢。function Untitled2
clear;clc
w1 = 35;
p = 60;
w2 = 45;
r = 10;
miu = 2000;
sigma = 600;
for w1 = 35:45
v1 = norminv((w2-w1)/(w2-r),miu,sigma); % 正态分布反函数
v1 = normpdf(v1,miu,sigma); % 正态分布概率密度函数
fun1 = @(x) (x-v1).*normpdf(x,0,1);
L1 = quadl(fun1,v1,10000); % 无穷积分,上限用10000
E1 = (p-w1)*miu-(w1-r)*v1*sigma-(p-r)*L1*sigma; % 期望
v2 = norminv((p-w2)/(p-r),miu,sigma);
v2 = normpdf(v2,miu,sigma);
fun2 = @(x) (x-v2).*normpdf(x,0,1);
L2 = quadl(fun2,v2,10000);
E2 = (p-w2)*miu-(w2-r)*v2*sigma-(p-r)*L2*sigma;
plot(w1,E1,'r*',w1,E2,'o')
hold on
end
附上参考:http://www.