The Math Forum

Search All of the Math Forum:

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

Math Forum » Discussions » Software » comp.soft-sys.matlab

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

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Thomas Kragh

Posts: 8
Registered: 12/7/04
Trying to Vectorize Double Integrals
Posted: Aug 5, 1996 2:13 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

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.:

x = dblquad('function_xy',xlow,xhigh,ylow,yhigh)

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...


Thomas Kragh

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.