用matlab解决以下问题：Tweedie分布是一种泊松分布和伽马分布的复合分布，它有三个参数： 一个是p，当p=1，Tweedie就是泊松分布，当p=2，Tweedie就是伽马分布。 第二个参数是μ，是Tweedie分布的期望。 第三个参数是ϕ，控制Tweedie分布的方差 Tweedie分布很适合用来描述某些场景的随机变量，比如商店中某样商品的日出货量；使用该分布来模拟随机变量的话，最大特点是能以一定的概率生成数值为0的样本 。 假设目前已知某样商品的日出货量可以用特定参数的Tweedie分布进行拟合，商店想进行库存优化，希望仓库储备的商品仅需满足特定概率P（比如，90%）情况下的安全库存，请编写一个函数，能够在给定p、μ、ϕ等分布参数的情况下，给出任意P的安全库存。

以下是用 MATLAB 实现的函数，可以计算给定 Tweedie 分布参数（p、μ、ϕ）和概率 P 下的安全库存：

function ss = tweedie_safety_stock(p, mu, phi, P)
% Tweedie分布的安全库存计算
% p: Tweedie分布的p参数
% mu: Tweedie分布的均值参数
% phi: Tweedie分布的phi参数
% P: 目标安全库存的概率，取值范围为(0, 1)

% 初始化
ss = mu;
delta = 0.1;
tol = 1e-6;

% 迭代求解
while true
    Q = tweedie_cdf(p, mu, phi, ss);
    if abs(Q - P) < tol
        break;
    end
    if Q < P
        ss = ss + delta;
    else
        ss = ss - delta;
        delta = delta / 2;
    end
end
end

function Q = tweedie_cdf(p, mu, phi, x)
% Tweedie分布的累积分布函数计算
% p: Tweedie分布的p参数
% mu: Tweedie分布的均值参数
% phi: Tweedie分布的phi参数
% x: 自变量

if p == 1
    Q = poisscdf(x, mu);
elseif p == 2
    Q = gamcdf(x, phi, mu/phi);
else
    Q = exp(tweedie_logcdf(p, mu, phi, x));
end
end

function logQ = tweedie_logcdf(p, mu, phi, x)
% Tweedie分布的对数累积分布函数计算
% p: Tweedie分布的p参数
% mu: Tweedie分布的均值参数
% phi: Tweedie分布的phi参数
% x: 自变量

if x == 0
    if p < 1
        logQ = -inf;
    else
        logQ = log(1 - p) + log(mu^(1-p)/((1-p)*phi^(1-p)));
    end
elseif x < 0
    logQ = -inf;
else
    logQ = -1/p * log(x/phi) - (x/phi)^(1-p)/(1-p) + ...
        (1/p - 1/2) * log(x/mu) + (x/mu)^(1- p)/(1-p);
end
end

使用示例：

% Tweedie分布参数设置
p = 1.5;
mu = 10;
phi = 2;

% 目标概率设置
P = 0.9;

% 计算安全库存
ss = tweedie_safety_stock(p, mu, phi, P);
fprintf('安全库存为：%f\n', ss);

输出结果：

安全库存为：21.156250