FAQFAQ   SearchSearch   MemberlistMemberlist   UsergroupsUsergroups   RegisterRegister 
 ProfileProfile   PreferencesPreferences   Log in to check your private messagesLog in to check your private messages   Log inLog in 
Forum index » Science and Technology » Math » num-analysis
weight functions for (numerical) integration
Post new topic   Reply to topic Page 1 of 1 [3 Posts] View previous topic :: View next topic
Author Message
Bart
science forum beginner


Joined: 07 Jul 2005
Posts: 39

PostPosted: Fri Jun 23, 2006 12:32 pm    Post subject: weight functions for (numerical) integration Reply with 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

--
Back to top
Peter Spellucci
science forum Guru


Joined: 29 Apr 2005
Posts: 702

PostPosted: Fri Jun 23, 2006 3:17 pm    Post subject: Re: weight functions for (numerical) integration Reply with quote

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

PostPosted: Fri Jun 23, 2006 3:53 pm    Post subject: Re: weight functions for (numerical) integration Reply with quote

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
Display posts from previous:   
Post new topic   Reply to topic Page 1 of 1 [3 Posts] View previous topic :: View next topic
The time now is Sat Jan 10, 2009 3:41 am | All times are GMT
Forum index » Science and Technology » Math » num-analysis
Jump to:  

Similar Topics
Topic Author Forum Replies Last Post
No new posts Generating function for Mathieu funct... cosmicstring@gmail.com Math 1 Fri Jul 21, 2006 8:39 am
No new posts Entire functions, polynomial bounds david petry Math 2 Thu Jul 20, 2006 11:09 pm
No new posts Help in identifying a numerical method Don num-analysis 2 Thu Jul 20, 2006 8:56 pm
No new posts Are Bessel Functions Transcendental? John Schutkeker Math 28 Tue Jul 18, 2006 2:24 am
No new posts scalar functions on non-diagonalizabl... xyz91234@yahoo.com Math 6 Mon Jul 17, 2006 5:37 pm

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
[ Time: 0.6080s ][ Queries: 16 (0.1568s) ][ GZIP on - Debug on ]