Как использовать функцию quadgk в Matlab в качестве дескриптора функции?

Я пытаюсь решить двойной интеграл, где я решаю как внутренний, так и внешний интеграл с функцией quadgk.

% The integrand is of course a function of both x and y
integrand = @(x,y) (phi(j,y,X) - phi(j,x,X))*(phi(i,y,X) - phi(i,x,X))/abs(y-x)^(2*s+1)

% The inner integral is a function of x, and integrates over y
inner = @(x) quadgk(@(y)integrand(x,y), x-lambda, x+lambda)

% The inner integral is integrated over x to yield the value of the double integral 
dblIntegral = quadgk(inner, -(1+lambda), 1+lambda)

и я получаю следующую ошибку:

integrand = @(x,y)(phi(j,y,X)-phi(j,x,X))*(phi(i,y,X)-phi(i,x,X))/abs(y-x)^(2*s+1)
inner = @(x)quadgk(@(y)integrand(x,y),x-lambda,x+lambda)

??? Error using ==> quadgk at 108
A and B must be scalar floats.

Error in ==> @(x)quadgk(@(y)integrand(x,y),x-lambda,x+lambda)
Error in ==> quadgk>evalFun at 344
        fx = FUN(x);

Error in ==> quadgk>f1 at 362
    [y,too_close] = evalFun(tt);

Error in ==> quadgk>vadapt at 258
        [fx,too_close] = f(x);

Error in ==> quadgk at 197
[q,errbnd] = vadapt(@f1,interval);

Error in ==> frational_laplacian at 29
dblIntegral = quadgk(inner, -(1+lambda), 1+lambda)

matlab numerical-methods
person Vikas Gautam    schedule 22.10.2013    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Вы получаете эту проблему, потому что quadgk внутренне оценивает переданную функцию, используя вектор значений. Это означает, что в строке

dblIntegral = quadgk(inner, -(1+lambda), 1+lambda)

функция inner вызывается с вектором x. quadgk требует, чтобы a и b (в данном случае x±лямбда) были скалярами, поэтому вы получаете ошибку, которую вы видели. Если вы хотите убедиться, что inner всегда получает скаляр в качестве входных данных, вы можете использовать

dblIntegral = quadgk(@(x)arrayfun(inner,x), -(1+lambda), 1+lambda)
person MrAzzaman    schedule 22.10.2013