Integrating The Normal Probability Density Function

2017 10 15
1 个回答
回答已采纳
更新时间：2017 10 16
12 次查看（30 天）

0 个投票

Hello community
I was wondering if the community could help me with integration with Matlab.
I am new to matlab so my apologies for the "basic" question.
I have a set of data from which I have calculated the:-
• Mean = 1853.910 • Standard Deviation = 1.829
I want to integrate the normpdf function in Matlab from x = 1855.739 to x = 1852.081 however I am getting silly results, I should be getting an approx value of 68.26%.
I have searched the Matlab directory and found the normpdf M file and amended it so it now looks like:-
y = exp(-0.5 * ((x - 1853.910)./1.829).^2) ./ (sqrt(2*pi) .* 1.829);
Where I have replaced the variables MU and SIGMA with 1853.910 and 1.829.
I cannot seem to integrate this function from x = 1855.739 to x = 1852.081.
Can anyone help?
Its driving me nuts.
Thank you.

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

 采纳的回答

Star Strider
Star Strider 2017-10-15
编辑：Star Strider 2017-10-15

2 个投票

The normcdf function will do the integration for you:
Mean = 1853.910;
Standard_Deviation = 1.829;
lims = [1855.739 1852.081];
cp = normcdf(lims, Mean, Standard_Deviation);
Prob = cp(1) - cp(2)
Prob =
0.6826895
EDIT —
If you want to do the integration yourself, this works:
p = @(x,m,s) exp(-((x-m).^2)/(2*s.^2)) / (s*sqrt(2*pi));
c = integral(@(x) p(x, Mean, Standard_Deviation), 1852.081, 1855.739);
The result is the same as ‘Prob’.

6 个评论
显示 4更早的评论 隐藏 4更早的评论

John hope
John hope 2017-10-15
Star Strider thank you for your response.
After looking at the normcdf function I have realised that it does not give me what I am looking for.
I actually need to integrate the normpdf function because I created a graph by passing in all of my X values (1851.400, 1851.700, 1852.600, ...) along with my mean (1853.910) and standard deviation (1.829) using the normpdf function. It is this normpdf function I am trying to integrate between.
I took your steps but replaced normcdf with normpdf, I got the following:-
>> Mean = 1853.910;
>> Standard_Deviation = 1.829;
>> lims = [1855.739 1852.081];
>> cp = normpdf(lims, Mean, Standard_Deviation);
>> Prob = cp(1) - cp(2)
Prob =
1.643130076445232e-014
As you can see I get a very wired result.
Can you help?
I have attached my X values along with the corresponding Y values from the normpdf function for reference.
Any ideas?
Thanks
Star Strider
Star Strider 2017-10-15
You’re using normpdf, not normcdf as I used in my code. That makes all the difference!
Change to my code, or use my additional code after the ‘EDIT’ text that uses an anonymous function implementation of normpdf), and you will get the correct result.
John hope
John hope 2017-10-16
Hi Star Strider
I know I am being really stupid so my apologies and thanks for putting up with me.
But in your example you have used the function normcdf :-
cp = normcdf(lims, Mean, Standard_Deviation);
But I want to integrate normpdf.
I am not sure what I am missing.
Would it be possbile to explain the line of code and why you have used this:-
cp = normcdf(lims, Mean, Standard_Deviation);
I have taken your additional code:-
p = @(x,m,s) exp(-((x-m).^2)/(2*s.^2)) / (s*sqrt(2*pi));
and altered it to look like:-
p= @(x)exp(-((x-1853.910).^2)/(2*1.829.^2))/(1.829*sqrt(2*pi))
I then used the following code to integrate it:-
>> (int (p,x,1852.081,1855.739))*100
ans =
68.268949213710298306169235885366
So with your help the problem has been solved so thank you but I would still be eager to know why you have used the normcdf function just to learning purposes.
Thanks again.
Star Strider
Star Strider 2017-10-16
My pleasure.
The normcdf function (the cumulative distribution function) is defined as the integral of the normpdf (the probability density function).
To illustrate —
cum_prob = integral(@(x) normpdf(x, Mean, Standard_Deviation), 1852.081, 1855.739)
cp =
0.6826895
Integrating the normpdf function gives the same result as using the normcdf function.
The probability density function is just that — the probability at a specific value of the independent variable. Use the integral of that — the cumulative distribution function — to get the probability that a specific event has occurred between the bounds of the integration. See any decent basic statistics text for a full explanation.
If my Answer helped you solve your problem, please Accept it!
John hope
John hope 2017-10-16
Much appreciated answer accepted.
Thanks again.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 MATLAB Support for MinGW-w64 C/C++ Compiler 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by