a problem while using an anonymous function

1 次查看（过去 30 天）
Eliraz Nahum
Eliraz Nahum 2018-10-10
评论： Eliraz Nahum 2018-10-11
hello,
I can't understand what's wrong. I will appreciate any help. I have a problem while trying to using the integral function with the handle of an anonymous function. I am adding a photo of the output error.
clear all
close all
clc
lamda=5;
x=[1:1:30];
p=zeros(1,length(x));
poiss_PDF=zeros(1,length(x));
poiss_int=zeros(1,length(x));
step=0;
poiss_pdf=@(x) (1/factorial(x))*(lamda^x)*exp(-lamda);
for x_i=x
step=step+1;
p(step)=poiss_pdf(x_i);
poiss_PDF(step)=poisspdf(x_i,lamda);
poiss_int(step)=integral(poiss_pdf,x(1),x(end));
end
plot(x,p,'b')
title (['Poisson PDF (calculated) with the parameter lamda= ',num2str(lamda)])
xlabel('x')
ylabel ('Distribution')
xlim([x(1) x(end)])
hold on
plot (lamda,poiss_pdf(lamda),'dr')
figure
plot(x,poiss_PDF,'g')
title (['Poisson PDF (from Matlab) with the parameter lamda= ',num2str(lamda)])
xlabel('x')
ylabel ('Distribution')
xlim([x(1) x(end)])
hold on
plot (lamda,poisspdf(lamda,lamda),'dr')
figure
plot(x,poiss_int,'k')
title (['Poisson integral (calculated) with the parameter lamda= ',num2str(lamda)])
xlabel('x')
ylabel ('Probability')
xlim([x(1) x(end)])
hold on
plot (lamda,poiss_pdf(lamda),'dr')
  2 个评论
jonas
jonas 2018-10-10
编辑：jonas 2018-10-10
The problem is that factorial only works for integers, so your function is only defined in those points. This is not my field of expertise, but wikipedia tells me you can use a gamma-function instead.
Bruno Luong
Bruno Luong 2018-10-10
Why you are using integral to handle discrete PDF like Poisson? For discrete you need to do a sum.
Or transform the discrete PDF to a sum of weighted dirac, not sur integral can handle though.

请先登录，再进行评论。

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

采纳的回答

the cyclist
the cyclist 2018-10-10
编辑：the cyclist 2018-10-10
Replace
poiss_pdf=@(x) (1/factorial(x))*(lamda^x)*exp(-lamda);
with
poiss_pdf=@(x) (1./gamma(x+1)).*(lamda.^x).*exp(-lamda);
and your code will execute to completion. I did not try to figure out if the output is sensible.
Note that I did two things here:
  • Replaced your factorial with the appropriate gamma function
  • Used element-wise operations in the function, which I'm pretty sure is what you intended
  3 个评论
显示 1更早的评论
the cyclist
the cyclist 2018-10-10
编辑：the cyclist 2018-10-10
Because internally, the integral function creates a vector from x(1) to x(end), and evaluates your function using that vector.

请先登录，再进行评论。

更多回答（1 个）

Jos (10584)
Jos (10584) 2018-10-10
You do not show the whole error message!
But I also assume that the problem is indeed in the use of factorial inside the function poiss_pdf that is being passed to integral. The function factorial(x) is only defined for integer values of x. You can use gamma(x+1) instead. See the documentation of gamma for more details.

类别

在 Help Center 和 File Exchange 中查找有关 Poisson Distribution 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by