Integrating a truncated normpdf

2017 4 11
0 个回答
更新时间：2017 4 12
2 次查看（30 天）

0 个投票

I have a function multiplied by a normal distribution that I can integrate in the following manner:
f2 = @(theta) f1(theta)*normpdf(theta)
ans = integral(f2,-k,k)
and everything is okay and great. However, I need normpdf(theta) to be truncated from -k to k. Truncate() only seems to work for makedist() which would get rid of the anonymous function I need for the integral. The closest solution I can find is here: https://www.mathworks.com/matlabcentral/newsreader/view_thread/100921 but the link to the tutorial is broken and I'm not sure if mle is the answer to my problem...
Thank you in advance for any advice

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

Christina
Christina 2017-4-11
f1 =@(theta)(1 - 1/(5+2*x))*((100+theta+s)^(0.5)) / (0.5) + (1/(5+*x))*((100+theta+s-5)^(0.5))/(0.5) where x and s are syms found using fsolve(); the integral is part of the equations being passed into fsolve
Torsten
Torsten 2017-4-12
Why don't you insert the explicit formula for the truncated normal from -k to k ?
https://en.wikipedia.org/wiki/Truncated_normal_distribution
Best wishes
Torsten.
Christina
Christina 2017-4-12
编辑：Christina 2017-4-12
Thank you Torsten for your reply (and you too Star Strider). That link really helped me understand what was going on in that question I linked to. I did find a work-around that I'll share.
pd = truncate(makedist('Normal'),-k,k);
integral(@(theta)pdf(pd,theta),-k,k)

请先登录，再进行评论。

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

回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

提问:

Christina
Christina
2017-4-11

编辑：

Christina
Christina
2017-4-12

Community Treasure Hunt

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

Start Hunting!

Translated by