|
|
| Author |
Message |
Bart science forum beginner
Joined: 07 Jul 2005
Posts: 39
|
Posted: Fri Jun 23, 2006 12:32 pm Post subject:
weight functions for (numerical) integration
|
|
|
Could anybody comment on the following:
Suppose you want to calculate the integral
I = \int_{0}^{1} f(x)dx
over the unit interval or unit hypercube. Then using Monte Carlo
methods, this can be done by simply taking the average of the
function values at certain randomly chosen x_k:
I_N = 1/N * sum(f(x_k), k=1..N)
Now suppose you introduce a weightfunction w(x), and as an
alternative, you approximate the integral by
I_N = 1/N * sum(w(k/N)*f(x_k), k=1..N)
Then my question is: in how far is it necessary for the weights
to be normalized such that the integral still converges to I?
That is, in how far must we have that
1/N * sum(w(k/N), k=1..N) = 1
or
lim_{N->infty} 1/N * sum(w(k/N), k=1..N) = 1
and what happens if e.g.
1/N * sum(w(k/N), k=1..N) = C
or
lim_{N->infty} 1/N * sum(w(k/N), k=1..N) = C
with C a certain constant.
Can anybody comment on this or give me some usefull references
that are related to this question?
Thanks,
Bart
-- |
|
| Back to top |
|
 |
Peter Spellucci science forum Guru
Joined: 29 Apr 2005
Posts: 702
|
Posted: Fri Jun 23, 2006 3:17 pm Post subject:
Re: weight functions for (numerical) integration
|
|
|
In article <1151065959.64787@seven.kulnet.kuleuven.ac.be>,
Bart <Reply.In@This.Group> writes:
| Quote: | Could anybody comment on the following:
Suppose you want to calculate the integral
I = \int_{0}^{1} f(x)dx
over the unit interval or unit hypercube. Then using Monte Carlo
methods, this can be done by simply taking the average of the
function values at certain randomly chosen x_k:
I_N = 1/N * sum(f(x_k), k=1..N)
Now suppose you introduce a weightfunction w(x), and as an
alternative, you approximate the integral by
I_N = 1/N * sum(w(k/N)*f(x_k), k=1..N)
Then my question is: in how far is it necessary for the weights
to be normalized such that the integral still converges to I?
That is, in how far must we have that
1/N * sum(w(k/N), k=1..N) = 1
or
lim_{N->infty} 1/N * sum(w(k/N), k=1..N) = 1
and what happens if e.g.
1/N * sum(w(k/N), k=1..N) = C
or
lim_{N->infty} 1/N * sum(w(k/N), k=1..N) = C
with C a certain constant.
Can anybody comment on this or give me some usefull references
that are related to this question?
Thanks,
Bart
--
|
a standard theorem about numerical quadrature:
theorem:
a rule
1/N * sum(w(k/N)*f(x_k), k=1..N) converges to \int_{0 to 1 } f(x)dx
as N to infinity for any f continuous in [0,1]
if and only if
1) limsup (1/N sum_{k=1 to N} ABS(w(k/N) ) N-> infinity is finite
2) the rule converges for any polynomial
(theorem of Polya and Stekloff)
hence, taking f(x)=1 for all x
you see that
lim {N to infinity } 1/N sum_{k=1 to N} w(k/N) = 1
is a must.
see for example: H. Engels, Numerical qyadrature and cubature, Acda. Press 1980,
page 183
hth
peter |
|
| Back to top |
|
 |
Bart science forum beginner
Joined: 07 Jul 2005
Posts: 39
|
Posted: Fri Jun 23, 2006 3:53 pm Post subject:
Re: weight functions for (numerical) integration
|
|
|
On 2006-06-23, Peter Spellucci <spellucci@fb04373.mathematik.tu-darmstadt.de> wrote:
| Quote: |
a standard theorem about numerical quadrature:
theorem:
a rule
1/N * sum(w(k/N)*f(x_k), k=1..N) converges to \int_{0 to 1 } f(x)dx
as N to infinity for any f continuous in [0,1]
if and only if
1) limsup (1/N sum_{k=1 to N} ABS(w(k/N) ) N-> infinity is finite
2) the rule converges for any polynomial
(theorem of Polya and Stekloff)
hence, taking f(x)=1 for all x
you see that
lim {N to infinity } 1/N sum_{k=1 to N} w(k/N) = 1
is a must.
|
Thanks for confirming what i was thinking intuitively. :-)
| Quote: | see for example: H. Engels, Numerical qyadrature and cubature, Acda. Press 1980,
page 183
|
Thanks. I see we have this book in our university library. I
will definitely look it up!
Best wishes,
Bart
-- |
|
| Back to top |
|
 |
Google
|
|
| Back to top |
|
 |
|
|
The time now is Sat Jan 10, 2009 3:41 am | All times are GMT
|
|
Guitar Lessons | Loans | Hotels in Rome | Bad Credit Mortgages | Free Credit Scores
|
|
Copyright © 2004-2005 DeniX Solutions SRL
|
|
Other DeniX Solutions sites:
Electronics forum |
Medicine forum |
Unix/Linux blog |
Unix/Linux documentation |
Unix/Linux forums
|
Powered by phpBB © 2001, 2005 phpBB Group
|
|