如何在matlab中计算出定积分

Question

我已运行以下命令来计算以下函数的数值积分。

但是在Matlab中，而不是数字作为输出，它给了我很大的表达。

gamma = sqrt(alpha^2 - beta^2);
K1=besselk(1,alpha*sqrt(delta^2+(x-mu)^2));
pdf = alpha*delta* K1/(pi*sqrt(delta^2+(x-mu)^2)) * exp(delta*gamma+ beta*(x-mu));

prob = int(pdf , x, 0,0.5);

具有以下参数：

mu = 0.034 ; 


delta = 2.12; 

alpha = 0.05;
beta = -0.001; 

我已经收到了这个结果。

int((53*exp(59679899628370215233/562949953421312000000 - x/1000)*besselk(1, ((x - 17/500)^2 + 2809/625)^(1/2)/20))/(500*pi*((x - 17/500)^2 + 2809/625)^(1/2)), x, 0, 1/2)

任何解决此问题的方法，我将不胜感激。 如何在Matlab中计算此表达式？

Answer 1

您得到的结果尚未评估。 要获得n有效数字的数字答案，请使用，

prob_vpa = vpa(prob, n); 

结果将仍然是符号变量，用于进一步的符号计算。 您也可以将其转换为double ，

prob_double = double(prob);