征费航班分布直方图 Matlab答案

【问题标题】：Levy Flight Distribution Histogram Matlab征费航班分布直方图 Matlab
【发布时间】：2012-11-06 18:32:53
【问题描述】：

我一直在四处寻找将征费分布拟合到直方图的方法，但无济于事。我希望测试 Levy Flight 分布在数据上的表现，无论它是否真的适合我们的数据类型。由于我对 PDF 相当陌生，并且除了 matlab 中的 distfittool GUI 之外还安装了我自己的 PDF，所以我有点不知道我需要做什么才能正确执行此操作。

所以目前，我的数据是一个 208x1 的向量，208 个点代表 208 个不同对象的不同速度。速度仅通过每次总距离计算。

现在，目前我使用了描述 Levy 航班的函数： http://reference.wolfram.com/mathematica/ref/LevyDistribution.html（输出[3]）

我用下面的代码试了一下：

load('Speeds.mat')
modelFun = @(p,x) (exp(-p(1)./(2.*(x-p(2)))).*(p(1)./(x-p(2))).^3/2)./(sqrt(2.*pi).*p(1));
startingVals = [1 1];
coefEsts = nlinfit(LBSpeed,modelFun,startingVals);

我完全清楚，我对 Levy 航班分布不熟悉是为什么我什至不确定这是否是用于分布的正确函数，也不是我需要正确传递给它的参数的根本原因做这个。如果有人能给我更多的见解，我将不胜感激。

【问题讨论】：

您能否更详细地了解您在尝试时究竟出了什么问题？

标签： matlab histogram data-fitting

【解决方案1】：

我已经看到与此问题类似的问题没有答案，所以在得到同事的帮助后，我想发布解决方案

与我最初的问题不同的另一件事是，它变得分段以更好地满足我自己发布的征费飞行方程的参考。我选择的起始值是任意的。

load('Speeds.mat')
[N,X] = hist(Speed,20);
Y = N/(sum(N))/diff(X(1:2));

%Get best parameters
modelFun = @(p,x) (x>p(2)).*(exp(-p(1)./(2.*(x-p(2)))).*(p(1)./(x-p(2))).^(3/2))./(sqrt(2.*pi).*p(1));
startingVals = [1,1];
coefEsts = nlinfit(X,Y,modelFun, startingVals);

%Visualize fit
bar(X,Y);
hold on;
model_eval = modelFun(coefEsts,X);
plot(X,model_eval,'r','LineWidth',2);

我一开始不知道如何拟合直方图，所以希望这对新手有所帮助！

【讨论】：

【解决方案2】：

您可以在本网站上找到matlab代码中的拟合例程

http://math.bu.edu/people/mveillet/html/alphastablepub.html

【讨论】：