Как написать функцию, которая включает численное интегрирование?

Question

Мне нужно рассчитать интеграцию во многих точках. Так,

f = f(r,theta,k,phi); 
q =integral2(f,0,1,0,2*pi,'AbsTol',0,'RelTol',1e-10); % integration should be by k and phi 

Я хочу д быть функцией г и тета, что я могу назвать его в любое время, чтобы вычислить интеграцию в данном г и тета точки. Как я могу это сделать? Проблема заключается в том, что я не мог использовать неопределенные функции @ или matlabFunction(), потому что кажется, что интеграция выполняется в первую очередь, и когда Matlab обнаруживает, что он не имеет всех аргументов, он приносит некоторые ошибки.

Answer 1

это все, что вы ищете (я до сих пор не знаю, что f возвращается) ?:

r = ...  % Define 
theta = ... % Define 
g = @(k,phi)f(r,theta,k,phi); % g is now a function of k and phi 
q = integral2(g,0,1,0,2*pi,'AbsTol',0,'RelTol',1e-10); 

Это создает анонимную функцию g где значение r и theta захватываются в качестве параметров и k и theta все еще являются аргументами. Эта концепция известна как closure в области информатики.

Если вы хотите, чтобы превратить все это в зависимости от r и theta, которая возвращает q вы можете можете создать следующую анонимную функцию:

q = @(r,theta)integral2(@(k,phi)f(r,theta,k,phi),0,1,0,2*pi,'AbsTol',0,'RelTol',1e-10); 

, что вы можете позвонить с q(r,theta). Конечно, вы могли бы использовать обычные функции (которые, как правило, быстрее и делают ваш код более понятным для других).