NaN (Not a Number) problem in integration

1 次查看（过去 30 天）

Serhat 2012-9-19

Hi all. I am trying to calculate this integration, it gives NaN ,instead of pi/2. How can I solve the problem?

mz=zeros(4,1); h=zeros(4); I=eye(4);

f=@(x)(arrayfun(@(x)(exp(-mz'*(1./(h+((sin(x).^2).*I)))*mz)),x));

out = quad(f,0,pi/2)

Thanks

José-Luis 2012-9-19

编辑：José-Luis 2012-9-19

No, it should not give you pi/2. The way you have set up your matrices, the expression

h+((sin(x).^2).*I)))*mz

will always be equal to zero. Division by 0 = Inf.

Serhat 2012-9-19

Now it works, you are right. I put parentheses at wrong places. Thank you very much.

f = @(theta)(arrayfun(@(x)(det(sigmaZ./sin(x).^2+eye(4,4)).^-1)*exp(-1*mZ'*1/(sigmaZ+sin(x).^2.*eye(4))*mZ),theta));

José-Luis 2012-9-19

My pleasure.

MATLAB Language Fundamentals Data Types Numeric Types NaNs

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

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

Translated by