baotong.m

本次上传的代码为经典报童问题的matlab算法

%报童问题的计算机仿真
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<max
    %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的关系的曲线图')