fitting a normal distribution function to a set of data

Question

cgo 2014-12-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/167123-fitting-a-normal-distribution-function-to-a-set-of-data

回答： Star Strider 2014-12-18

Hi, I have a set of data, in the form of a histogram (with actual data also ready) and I want to fit a normal distribution curve on it. Is their an efficient way to do it?

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2014-12-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/167123-fitting-a-normal-distribution-function-to-a-set-of-data#answer_162762

在 MATLAB Online 中打开

If you have the Statistics Toolbox, use the histfit function.

Otherwise, this works:

d = 2*randn(250,1)+20;                          % Created Data
binrs = min(d):(max(d)-min(d))/25:max(d);       % Bin Definitions
knt = histc(d, binrs);                          % Histogram Counts
s = std(d);                                     % Intial Parameter Estimate
m = mean(d);                                    % Intial Parameter Estimate
% b(1) = mean,  b(2) = std,  b(3) = amplitude
pdfnrm = @(x,b) b(3) * 1./(b(2)*sqrt(2*pi)) .* exp(-((x-b(1)).^2./(2*b(2).^2)));
SSECF = @(b) sum((knt-pdfnrm(binrs,b)').^2);    % Sum-Squared-Error Cost Function
[B,SSE] = fminsearch(SSECF, [m; s; 10]);
figure(1)
bar(binrs,knt,'g')                              % Plot Histogram
hold on
plot(binrs,pdfnrm(binrs,B),'r')                 % Plot Normal Distribution
hold off
grid