[Date Index]
[Thread Index]
[Author Index]
Re: Limit and Root Objects
*To*: mathgroup at smc.vnet.net
*Subject*: [mg72620] Re: Limit and Root Objects
*From*: Paul Abbott <paul at physics.uwa.edu.au>
*Date*: Fri, 12 Jan 2007 06:05:49 -0500 (EST)
*Organization*: The University of Western Australia
*References*: <em606c$2o1$1@smc.vnet.net> <4586C045.2020805@metrohm.ch> <em8jfr$pfv$1@smc.vnet.net>
In article <em8jfr$pfv$1 at smc.vnet.net>,
Andrzej Kozlowski <andrzej at akikoz.net> wrote:
> What you describe, including the fact that the numbering or roots
> changes is inevitable and none of it is not a bug. There cannot exist
> an ordering of complex roots that does not suffer from this problem.
> What happens is this.
> Real root objects are ordered in the natural way. A cubic can have
> either three real roots or one real root and two conjugate complex
> ones. Let's assume we have the latter situation. Then the real root
> will be counted as being earlier then the complex ones. Now suppose
> you start changing the coefficients continuously. The roots will
> start "moving in the complex plane", with the real root remaining on
> the real line the two complex roots always remaining conjugate
> (symmetric with respect to the real axis). Eventually they may
> collide and form a double real root. If this double real root is now
> smaller then the the "original real root" (actually than the root to
> which the original real root moved due the the changing of the
> parameter), there will be a jump in the ordering; the former root
> number 1 becoming number 3.
> This is completely unavoidable, not any kind of bug, and I am not
> complaining about it. It takes only elementary topology of
> configuration spaces to prove that this must always be so.
But is there a continuous root numbering if the roots are not ordered?
What I mean is that if you compute the roots of a polynomial, which is a
function of a parameter, then if you assign a number to each root, can
you follow that root continuously as the parameter changes? Two examples
are presented below.
Here is some code to animate numbered roots using the standard root
ordering, displaying the root numbering:
rootplot[r_] := Table[ListPlot[
Transpose[{Re[x /. r[a]], Im[x /. r[a]]}],
PlotStyle -> AbsolutePointSize[10],
PlotRange -> {{-3, 3}, {-3, 3}},
AspectRatio -> Automatic,
PlotLabel -> StringJoin["a=", ToString[PaddedForm[Chop[a], {2, 1}]]],
Epilog -> {GrayLevel[1],
MapIndexed[Text[#2[[1]], {Re[#1], Im[#1]}] & , x /. r[a]]}],
{a, -6, 10, 0.5}]
First, we have a polynomial with real coefficients:
r1[a_] = Solve[x^5 - a x - 1 == 0, x]
Animating the trajectories of the roots using
rootplot[r1]
we observe that, as you mention above, when the complex conjugate roots
2 and 3 coalesce, they become real roots 1 and 2 and root 1 becomes root
3. But, ignoring root ordering, why isn't it possible for these roots to
maintain their identity (I realise that at coelescence, there is an
arbitrariness)?
Second, we have a polynomial with a complex coefficient:
r2[a_] = Solve[x^5 + (1+I) x^4 - a x - 1 == 0, x]
Animating the trajectories of the roots using
rootplot[r2]
we observe that, even though the trajectories of the roots are
continuous, the numbering switches:
2 -> 3 -> 4
5 -> 4 -> 3
3 -> 4 -> 5
4 -> 3 -> 2
and only root 1 remains invariant. Again, ignoring root ordering, why
isn't it possible for all these roots to maintain their identity and so
be continuous functions of the parameter? And wouldn't such continuity
be nicer than enforcing root ordering?
Cheers,
Paul
_______________________________________________________________________
Paul Abbott Phone: 61 8 6488 2734
School of Physics, M013 Fax: +61 8 6488 1014
The University of Western Australia (CRICOS Provider No 00126G)
AUSTRALIA http://physics.uwa.edu.au/~paul
Prev by Date:
**Re: Re: ListDimension function**
Next by Date:
**Re: speed of multiplying polynomials**
Previous by thread:
**16bit image writing**
Next by thread:
**Re: Re: Limit and Root Objects**
| |