如何在matlab中实现对延迟几何分布的采样

% 定义延迟几何分布的参数
mu = 1; % 平均延迟
sigma = 0.5; % 延迟的标准差

% 定义延迟几何分布的概率密度函数和累积分布函数
pdf = @(x) (x >= 0) .* (1/(2*mu*sigma)) .* exp(-(abs(x-mu) + abs(x+mu))/(2*sigma^2));
cdf = @(x) (x >= 0) .* (1 - exp(-abs(x-mu)/mu)) .* (1 - exp(-abs(x+mu)/mu));

% 生成随机数
random_number = rand();

% 计算累积概率值
cumulative_probability = cdf(random_number);

% 计算延迟值
delay = fzero(@(x) cdf(x) - cumulative_probability, [0, 10*mu]);

% 打印结果
disp(['采样得到的延迟值为：', num2str(delay)]);