%报童问题的计算机仿真
clear all;
close all;
aoverb=[0.1:0.2:4];%设定a/b的变化范围
len=length(aoverb);%得到测试数据的长度
opt=zeros(1,len);%对于每一个a/b值存在一个最优订报量，初始化这个数组
a=zeros(1,len);%初始化a的值
min_loss=ones(1,len)*1000;%初始化每日平均损失值，设为1000
for i = 1 : len%对于每一个a/b值进行仿真
z=1;%初始化z,即订报量
b=1;%退回的价格
tm=300;%一轮实验的预定模拟天数为300天
max=1000;%最大订报量为1000份
a(i)=b.*aoverb(i);%订报的价格
while z %order = 300+400.*rand(1,tm); %产生销售量均匀分布随机数,表示报童当天卖出的报纸数
%order=-500.*log(rand(1,tm));%负指数分布
order=400.*randn(1,tm)+300;%高斯分布
loss=0;%对于一个a/b值的总损失
for j=1:tm%
if z>order(j)%供过于求的情况
loss = loss + (z - order(j)).*b;
else%供不应求的情况
loss = loss + (order(j) - z).*a(i);
end
end
if min_loss(i) >= loss/tm%计算平均损失
min_loss(i) = loss/tm;
opt(i) = z;%最优订报量
end
z = z + 1;%订报数增加1，进行下一次仿真
end
end
plot(aoverb,opt);
title('订报量z与a/b关系的曲线图')
figure;
plot(aoverb,min_loss);
title('最小平均损失与a/b的关系的曲线图')
