matlab报童的诀窍指数分布代码

matlab报童问题的诀窍指数分布代码如下：

function [Q_star, S_star] = newsboy(lam, c, p)
    % lam: 日销售量期望
    % c: 报纸进价
    % p: 报纸售价
    % Q_star: 最优订货量
    % S_star: 最优期望利润

    syms Q;
    R(Q) = p * min(Q, lam) - c * Q; % 利润函数

    dR = diff(R); % 求导
    Q_star = solve(dR == 0, Q); % 解方程，得到最优订货量

    S_star = double(subs(R, Q, Q_star)); % 计算最优期望利润
end

其中，lam为销售量期望，c为报纸进价，p为报纸售价，Q_star为最优订货量，S_star为最优期望利润。代码中，利润函数R(Q)采用了min函数来保证当销售量小于订货量时，实际售出的报纸数量不会超过订货量。求解最优订货量采用了求导的方法，将一阶导数dR=0的Q解出即为最优订货量。