如何绘制涉及Matlab中集成的函数？

Question

我有一个可用的特征函数CF(t, param1, param2, ...) （这里param s是某种分布的参数，被认为是固定的）。 我想使用Gil-Pelaez公式来获取CDF，所以我写了以下内容：

function [ F ] = CDF( x, param1, param2, ... )

integrand = @(v) imag(exp(-1i.*v.*x) .* CF(t, param1, param2, ...)) ./ v;
F = 0.5 - (1./pi) .* integral(integrand, 0, 100);

end

这适用于单个值，例如CDF(0.1, param1, param2, ...)给我想要的结果。 现在我想针对x的范围绘制CDF ，所以我做了

x = linspace(-1,1,100);
y = CDF(x, param1, param2, ...)
plot(x, y)

我最终遇到了以下错误：

Not enough input arguments.

Error in integralCalc/iterateScalarValued (line 314)
                fx = FUN(t);

Error in integralCalc/vadapt (line 132)
            [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);

但是y = CDF(x, param1, param2, ...)确实起作用，所以罪魁祸首似乎是被积数的exp(-1i.*v.*x)部分，因为v和x的大小不起作用比赛。 但我不确定如何解决此问题。

Answer 1

您已经猜对了，v和x的大小不匹配会导致错误。 integral将向量值传递到函数句柄，以通过一次计算多个点的函数值来加快计算速度。

但是在这种情况下，您应该在CDF函数中对积分的调用中将参数'ArrayValued'为true 。 这将迫使积分将v的标量值传递给函数句柄，然后CDF将能够为向量x返回向量y的向量。