В Matlab интегрируйте функцию, которая вызывает другую функцию - PullRequest
Новая
       78

В Matlab интегрируйте функцию, которая вызывает другую функцию

1 голос
/ 21 июня 2020

Я численно оцениваю двойной интеграл с помощью команды matlab integ2. Мой код выглядит следующим образом: 

l1=0; m1=0; 
funn =  @(theta,phi)(Y_hsph_new(l1,m1,theta,phi));
q = integral2(funn,0,pi/2,0,2*pi);

где 'Y_hsph_new' - это моя собственная функция MatLab, написанная как: 

function [ YYY ] = Y_hsph_new(l1,m1,theta,phi)
% OUTPUT: (I x (N+1)^2) matrix
I = length(theta);
Y=zeros(I,(l1+1)^2);
for i=1:I
    k=1;
    for n=0:l1
        P=legendre(n,2*cos(theta(i))-1); % (2*x-1) for transformation
        for m=-n:n   % degree
            if m>0
               Y(i,k)=(sqrt(((2*n+1)*factorial(n-abs(m)))/((2*pi)*factorial(n+abs(m)))))*P(abs(m)+1)*sqrt(2)*cos(m*phi(i));
            elseif m == 0
               Y(i,k)=(sqrt(((2*n+1)*factorial(n))/((2*pi)*factorial(n))))*P(1);
            else
               Y(i,k)=(sqrt(((2*n+1)*factorial(n-abs(m)))/((2*pi)*factorial(n+abs(m)))))*P(abs(m)+1)*sqrt(2)*sin(abs(m)*phi(i));
            end
            k=k+1;
         end
     end
 end
 YYY = Y(:,m1+l1+1+(l1)^2);

Я использовал функцию Y_hsph_new в других моих обычных кодах, и это работает нормально. Я считаю, что, возможно, код функции в порядке, но что-то, связанное с его конфигурацией, которое будет использоваться при интеграции, вызывает проблему. Ошибка, которую я получаю при запуске: - 

Error using integral2Calc>integral2t/tensor (line 241)
Integrand output size does not match the input size.

Error in integral2Calc>integral2t (line 55)
[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);

Error in integral2Calc (line 9)
[q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);

Error in integral2 (line 106)
Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);

Я серьезно дергаю за волосы после многих попыток. Было бы очень полезно, если бы кто-нибудь мог мне помочь в этом отношении. Спасибо.

1 Ответ

1 голос
/ 21 июня 2020

Способ определения funn несовместим с интегралом2. Например, если я ввожу массив значений theta и phi в funn, я ожидаю, что на выходе будет тот же размер массива, но это не то, что дает funn. например: 

[th ph]=meshgrid(linspace(0,pi/2),linspace(0,2*pi))

>> size(th)

ans =

100   100

, тогда как: 

>> size(funn(th,ph))

ans =

100     1

Итак, funn не работает должным образом для использования в интегральном2. Я думаю, что первое место, с которого нужно начать в вашей функции, - это строка I = length(theta);, поскольку длина теты не подходит, если тета является 2-мерным массивом.

...