Question

我遇到了Matlab集成的问题。它告诉我有关尺寸的错误。但是他们确实同意并且使用向量化运算符（。^。* etc ...）正确地完成了操作。这是简单的代码，但我被卡住了

A = 1:10;
B = 1:10;

K_fun = @(x) (x ./ sqrt((x + A.^2 ) .* (x + B.^2) .* (x + B.^2)) );
K = integral( K_fun, 0,Inf );

这里是命令窗口中的错误消息：

Error using  + 
Matrix dimensions must agree.

Error in @(x)(x./sqrt((x+A.^2).*(x+B.^2).*(x+B.^2)))

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

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

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

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

Error in PROVA_New_Drag3 (line 21)
K = integral( K_fun , 0 , Inf);

提前谢谢

Answer 1

它会评估向量中的积分限制，因此函数中的x实际上是[0 Inf]，这就是长度不同意的原因。

您可以设置＆＃39; ArrayValued＆＃39;在积分调用中标记为true

K = integral( K_fun, 0,Inf,'ArrayValued', true );

让它分别评估它们。

我收到了关于奇点的另一个警告，但这很可能是因为你使用了0，无穷大和除法，所以它与你的函数的关系比积分调用更多。在eps中散布一些K_fun可能会有所帮助。

更新：请参阅特洛伊对奇点警告的解释，并注意eps实际上对此没有帮助。

数值积分，尺寸不同意，尽管使用。^和。*

1 个答案: