MATLAB的报童模型代码

以下是一个简单的MATLAB代码示例，用于实现报童模型：

% 参数设置
demand_probability = [0.1, 0.15, 0.2, 0.25, 0.3]; % 需求量的概率分布
selling_price = [5, 6, 7, 8, 9]; % 销售价格
purchase_cost = 4; % 进货成本
stock = 10; % 初始库存

% 计算期望利润
expected_profit = zeros(1, length(demand_probability));
for i = 1:length(demand_probability)
    demand = i;
    if demand <= stock
        profit = (demand * selling_price(i) - purchase_cost) * demand_probability(i);
    else
        profit = (stock * selling_price(i) - purchase_cost) * demand_probability(i);
    end
    expected_profit(i) = profit;
end

% 选择最大期望利润对应的需求量
[max_profit, max_profit_index] = max(expected_profit);
optimal_demand = max_profit_index;

% 输出结果
fprintf('最大期望利润为：%f\n', max_profit);
fprintf('最优订货量为：%d\n', optimal_demand);

这段代码中，我使用了一个简单的需求量概率分布和销售价格来计算报童模型中的期望利润。通过比较不同需求量下的期望利润，找到最大利润对应的需求量作为最优订货量。