Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Trying to Vectorize Double Integrals
Replies: 4   Last Post: Aug 9, 1996 9:34 AM

 Messages: [ Previous | Next ]
 Thomas Kragh Posts: 8 Registered: 12/7/04
Trying to Vectorize Double Integrals
Posted: Aug 5, 1996 2:13 PM

Has anyone worked out a way to allow numerical 2D integration, with vector inputs for the
upper/lower bounds in the two dimensions? The mathworks FAQ page describes an m-file called
"dblquad", which calls "quad.m" twice in a nested fashion to perform the integration over the
two dimensions, but only returns a single value, i.e.:

where x, xlow, xhigh, ylow, yhigh are scaler values (1x1).

I'm modeling a focal plane array in Matlab as an MxN array of pixels, with a continuous
point-spread function. I want to evaluate the total energy over each pixel by integrating the
2D point-spread function over each pixel, like in the following pseudo-code...

rows = 1:1:M
cols = 1:1:N
Energy_on_FPA = vector_dblquad('psf', (rows - 0.5), (rows + 0.5), ...
(cols - 0.5), (cols + 0.5))

Where Energy_on_FPA is a MxN matrix, with each element being the integral of the 2D
point-spread function over a given pixel, centered at coordinates (row,col).

Any suggestions? I already looked in the numerical integration toolbox in the user-supplied
function section of the ftp site, but to no avail...

Thanks,

Thomas Kragh
tjkragh@naa.rockwell.com
tjkragh@chelsea.ios.com

Date Subject Author
8/5/96 Thomas Kragh
8/6/96 Patrick C Hew
8/6/96 Allan White
8/6/96 Cleve Moler
8/9/96 Joe Deasy