Not enough Input arguments in integral function

Question

Alix DUMAS 2021-10-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1564591-not-enough-input-arguments-in-integral-function

评论： Alix DUMAS 2021-10-15

采纳的回答： David Hill

Good afternoon,

I am trying to compute the exact Wagner function on MATLAB with the following equation :

, where c(k) is the Theodorsen function. I already created a function computing c(k), but I am having trouble computing the integral. I have defined k = 0.1 and s = linspace(0,100,100).

For now I have tried two options with the integral(fun, xmin, xmax) function of MATLAB :

1) Create a function

in another matlab script and then write W = (2./pi).*integral(fun_W,0, Inf). The result was an error saying :

"Not enough input arguments.

Error in W_integrande (line 2)

f = F_k(k).*sin(k.*s)./k;

Error in HMWK3_DUMAS (line 14)

W = (2./pi).*integral(W_integrande,0, Inf)"

2) On another forum I saw people using

"fun_W = @(k,s) real(c_k(k)).*sin(k.*s)./k;

W = (2./pi).*integral(fun_W,0, Inf)"

So I tried it and this time the error is

"Not enough input arguments.

Error in HMWK3_DUMAS>@(k,s)real(c_k(k)).*sin(k.*s)./k (line 13)

fun_W = @(k,s) real(c_k(k)).*sin(k.*s)./k;

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 83)

[q,errbnd] = vadapt(@AToInfInvTransform,interval);

Error in integral (line 87)

Q = integralCalc(fun,a,b,opstruct);

Error in HMWK3_DUMAS (line 14)

W = (2./pi).*integral(fun_W,0, Inf);"

I don't understand why I am getting the error "not enough input arguments", could someone help me with this ?

Thank you in advance !

2 个评论
显示无隐藏无

David Hill 2021-10-15

It would be helpful if you provided your function for c(k).

Alix DUMAS 2021-10-15

Sure ! My function for c(k) is :

function[c] = c_k(k)

c = besselh(1,2,k)./(besselh(1,2,k)+j.*besselh(0,2,k));

end

请先登录，再进行评论。

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

Answer 1

David Hill 2021-10-15

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1564591-not-enough-input-arguments-in-integral-function#answer_809431

在 MATLAB Online 中打开

s=0:100;
for t=1:numel(s)
    fun_W=@(k)real(c_k(k))./k.*sin(k*s(t));
    W(t)=2/pi*integral(fun_W,0,Inf);
end
function c= c_k(k)
    c = besselh(1,2,k)./(besselh(1,2,k)+1j.*besselh(0,2,k));%need 1j
end