baotong.m

报童问题的计算机仿真 ％tm一轮实验的预定模拟天数 ％t一轮实验的仿真天数累积值 ％z订报量 ％z 最优订报量 ％g订报量z之上界 ％s1损失值之累计值 ％s最小损失值值 ％r按概率分

function [superz,supers]=baotong(tm,g,a,b)
%报童问题的计算机仿真
%tm一轮实验的预定模拟天数
%t一轮实验的仿真天数累积值
%z订报量
%z'最优订报量
%g订报量z之上界
%s1损失值之累计值
%s最小损失值值
%r按概率分布产生随机售报量样本
z=1;
supers=1000;
while z<g
    %r=poisondis(5,1,tm);
    r=round(2*randn(1,tm)+5);
    %r产生均匀分布随机数,表示报童当天卖出的报纸数
    %randn(m,n)是一个m*n随机数矩阵
    %round(X)使矩阵X的元素四舍五入到最接近它的整数
    t=1;
    s=0;
    dv=z>r;
    s=sum(((z-r)*b).*dv);
    s=s+sum(((r-z)*a).*(1-dv));
    aver_s=s/tm;
    if supers>=aver_s;
        supers=aver_s;
        superz=z;
    end;
    z=z+1;
end;