| Topic: |
Science > Physics |
| User: |
"Luna Moon" |
| Date: |
26 Jul 2007 06:55:52 AM |
| Object: |
question about complex line (Bromwich) integral |
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
.
|
|
| User: "Robert Israel" |
|
| Title: Re: question about complex line (Bromwich) integral |
26 Jul 2007 01:01:39 PM |
|
|
"Luna Moon" <lunamoonmoon@gmail.com> writes:
Hi all,
I want to use the numerical integration to do the Bromwich type integral,
as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Under the circumstances, I wouldn't trust any of them. It certainly looks
like there's something wrong here. Can you tell us what your function is?
Or better yet, find a simpler function that still gives you the same trouble.
--
Robert Israel
Department of Mathematics http://www.math.ubc.ca/~israel
University of British Columbia Vancouver, BC, Canada
.
|
|
|
| User: "Luna Moon" |
|
| Title: Re: question about complex line (Bromwich) integral |
26 Jul 2007 11:21:38 PM |
|
|
"Robert Israel" <> wrote in message
news:rbisrael.20070726175607$14f3@news.ks.uiuc.edu...
"Luna Moon" <lunamoonmoon@gmail.com> writes:
Hi all,
I want to use the numerical integration to do the Bromwich type integral,
as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along
the
vertical line x=c in the complex plane such that c is greater than the
real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values
go
from low to high.
Here is my headache: with no closed-form solution exist for such
integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Under the circumstances, I wouldn't trust any of them. It certainly looks
like there's something wrong here. Can you tell us what your function is?
Or better yet, find a simpler function that still gives you the same
trouble.
--
Robert Israel
Department of Mathematics http://www.math.ubc.ca/~israel
University of British Columbia Vancouver, BC, Canada
I thought they were all close to the correct value, the only problem is that
I don't know which one is the most correct.
Okay! Let me try to figure out how to put a huge formulae here, or a
simplified one, which pin-points the problem. I don't want people to lose
interest when seeing a huge function...
.
|
|
|
|
|
| User: "Nico Temme" |
|
| Title: Re: question about complex line (Bromwich) integral |
30 Jul 2007 02:39:19 PM |
|
|
On Thu, 26 Jul 2007, Luna Moon wrote:
Date: Thu, 26 Jul 2007 07:55:52 -0400
From: Luna Moon <lunamoonmoon@gmail.com>
Newsgroups: comp.soft-sys.math.maple, comp.soft-sys.matlab,
sci.math.num-analysis, sci.physics, sci.math
Subject: question about complex line (Bromwich) integral
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
You should take the value c where F(s) exp(st) has a minimum on the
real positive axis
(and of course c>max(RealPartOf(All Singularities))).
In general, at this minimum the analytic integrand has a saddle point.
Nico.
.
|
|
|
| User: "Luna Moon" |
|
| Title: Re: question about complex line (Bromwich) integral |
02 Aug 2007 11:58:26 AM |
|
|
On Jul 30, 3:39 pm, Nico Temme <Nico.Te...@cwi.nl> wrote:
On Thu, 26 Jul 2007, Luna Moon wrote:
Date: Thu, 26 Jul 2007 07:55:52 -0400
From: Luna Moon <lunamoonm...@gmail.com>
Newsgroups: comp.soft-sys.math.maple, comp.soft-sys.matlab,
sci.math.num-analysis, sci.physics, sci.math
Subject: question about complex line (Bromwich) integral
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
You should take the value c where F(s) exp(st) has a minimum on the
real positive axis
(and of course c>max(RealPartOf(All Singularities))).
In general, at this minimum the analytic integrand has a saddle point.
Nico.
Min or Max?
I plotted the F(c)*exp(c*t),
sometimes there is a min, sometimes there is no min(min is at the left
bound or right bound of the inteval), but there is a max...
.
|
|
|
| User: "Luna Moon" |
|
| Title: Re: question about complex line (Bromwich) integral |
02 Aug 2007 03:07:47 PM |
|
|
On Aug 2, 12:58 pm, Luna Moon <lunamoonm...@gmail.com> wrote:
On Jul 30, 3:39 pm, Nico Temme <Nico.Te...@cwi.nl> wrote:
On Thu, 26 Jul 2007, Luna Moon wrote:
Date: Thu, 26 Jul 2007 07:55:52 -0400
From: Luna Moon <lunamoonm...@gmail.com>
Newsgroups: comp.soft-sys.math.maple, comp.soft-sys.matlab,
sci.math.num-analysis, sci.physics, sci.math
Subject: question about complex line (Bromwich) integral
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
You should take the value c where F(s) exp(st) has a minimum on the
real positive axis
(and of course c>max(RealPartOf(All Singularities))).
In general, at this minimum the analytic integrand has a saddle point.
Nico.
Min or Max?
I plotted the F(c)*exp(c*t),
sometimes there is a min, sometimes there is no min(min is at the left
bound or right bound of the inteval), but there is a max...
If when I did the integral using the "optimal c" and I got a value
that is different than the theoretical integral value, whereas using a
c that is slightly different than the "optimal c", I got a value that
agrees with the theoretical integral value.
What might be the problem?
Can anybody show how is the "optimal c" obtained correctly? Thanks
.
|
|
|
|
|
| User: "Luna Moon" |
|
| Title: Re: question about complex line (Bromwich) integral |
02 Aug 2007 07:35:04 AM |
|
|
On Jul 30, 3:39 pm, Nico Temme <Nico.Te...@cwi.nl> wrote:
On Thu, 26 Jul 2007, Luna Moon wrote:
Date: Thu, 26 Jul 2007 07:55:52 -0400
From: Luna Moon <lunamoonm...@gmail.com>
Newsgroups: comp.soft-sys.math.maple, comp.soft-sys.matlab,
sci.math.num-analysis, sci.physics, sci.math
Subject: question about complex line (Bromwich) integral
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
You should take the value c where F(s) exp(st) has a minimum on the
real positive axis
(and of course c>max(RealPartOf(All Singularities))).
In general, at this minimum the analytic integrand has a saddle point.
Nico.
Sometimes the root of that saddle point is too costly to find. And
sometimes the root-finder even cannot find a root. Suppose I don't
find that root as the optimal "c", but why should that matter? Theory
tells us as long as "c" is away from any singularities, it should
always give back the same integration results... right?
.
|
|
|
|
| User: "Luna Moon" |
|
| Title: Re: question about complex line (Bromwich) integral |
02 Aug 2007 10:45:56 AM |
|
|
On Jul 30, 3:39 pm, Nico Temme <Nico.Te...@cwi.nl> wrote:
On Thu, 26 Jul 2007, Luna Moon wrote:
Date: Thu, 26 Jul 2007 07:55:52 -0400
From: Luna Moon <lunamoonm...@gmail.com>
Newsgroups: comp.soft-sys.math.maple, comp.soft-sys.matlab,
sci.math.num-analysis, sci.physics, sci.math
Subject: question about complex line (Bromwich) integral
Hi all,
I want to use the numerical integration to do the Bromwich type integral, as
shown in the URL below:
http://en.wikipedia.org/wiki/Bromwich_integral
The key is to find the line "x=c", where the integration is done along the
vertical line x=c in the complex plane such that c is greater than the real
part of all singularities of F(s).
Suppose I have identified all the singularities, and determined that in a
large range of
c>max(RealPartOf(All Singularities)),
the location of "c" should not matter and I should be able to freely vary
"c" as long as it doesn't violate our rules above, and the integration
should yield precisely the same values.
Now I plot the numerical integratal results along with varying "c" in the
safe region.
I found that the results differ, by about 10%. The change is very smooth.
From one end of "c_min" to the other end of "c_max", the integral values go
from low to high.
Here is my headache: with no closed-form solution exist for such integral,
only numerical solutions exist.
Which value of the numerical integral results shall I trust? (I used
Matlab's quadl, which is Gaussian Lobatto).
Secondly, the speed of integration varies a lot with the change of c. Does
anybody have some suggestions/pointers on what is a good theoretical
guideline of choosing "c" when both accuracy and speed are needed?
Thanks!
You should take the value c where F(s) exp(st) has a minimum on the
real positive axis
(and of course c>max(RealPartOf(All Singularities))).
In general, at this minimum the analytic integrand has a saddle point.
Nico.
What's the procedure of finding such "c"? Do you mean taking
derivative of F(c)*exp(c*t) for c as a real number?
.
|
|
|
|
|
|
Related Articles |
|
|
