Matlab单变量数据拟合Porter-Thomas分布及mle报错解决方法

报错核心原因 #

• 自定义PDF的参数顺序不符合Matlab mle 函数要求：mle 要求自定义pdf的输入参数必须是第一个位置为样本数据，后面依次为待估参数，你现有代码中 f=@(v,y) 的顺序是参数在前、数据在后，调用时传入样本会导致计算逻辑完全错误，产出无效的负/零概率值。
• 额外数值隐患：
  1. 若数据Y包含0或极小值，y.^(v/2 -1) 当v<2时会出现负数次幂运算，生成非有限值
  2. 参数下界设为0会导致 gamma(v/2) 在v趋近于0时发散

修正后的完整实现代码 #

% 1. 预处理数据：过滤非正样本（Porter-Thomas分布仅支持y>0）
Y = Y(Y > eps);

% 2. 修正参数顺序的PDF定义，添加数值稳定性处理
f = @(y, v) ((v./(2*0.99)).^(v/2)) .* (max(y, eps).^(v/2 - 1)) .* exp(-v.*y./(2*0.99)) ./ gamma(v/2);

% 3. 修正MLE调用，调整参数下界避免数值异常
V = mle(Y, 'pdf', f, 'Start', 0.5, 'LowerBound', 1e-3, 'UpperBound', 1);

% 4. 绘制匹配密度尺度的拟合结果
histogram(Y, 50, 'Normalization', 'pdf');
hold on
bins = linspace(min(Y), max(Y), 105);
plot(bins, f(bins, V), 'r--', 'LineWidth', 2);
hold off

额外说明 #

如果拟合结果偏差较大，可根据数据直方图的形态调整 Start 初始值的取值，也可以自定义对数似然函数用优化工具求解，避免mle内置优化陷入局部最优。

内容的提问来源于stack exchange，提问作者amir1122