matlab - 如何将双峰分布拟合到函数 matlab

Question

我有一个名为speed包含 200 次步行速度的向量

speed = [normrnd(80,2,100,1); normrnd(120,10,100,1)];

该向量遵循双峰分布。

steps是另一个包含每次行走步数的向量：

a = 8;
b = 100;
steps = (b-a).*rand(200,1) + a;

我创建了根据速度执行的步骤的直方图：

binstep = 1.5;
binranges = (min(speed):binstep:max(speed)+binstep)';
[~, ind] = histc(speed, binranges);
bincounts = accumarray(ind, steps, size(binranges));
hFig = figure(); axh = axes('Parent', hFig); hold(axh, 'all'); grid(axh, 'on');
bar(axh, binranges, bincounts); axis(axh, 'tight');

现在我想

1. 将双峰分布拟合到条形图，
2. 估计双峰分布的参数 u1,u2,sigma1,sigma2,
3. 评估它是否合适（即分布是双峰的）。

请你帮助我好吗？

Answer 1

要将数据拟合到分布，只需使用 fitdist 函数。在您的情况下，使用 Normal 将是：

pd = fitdist(y,'Normal')