Re: Real and Complex Roots presented in a single plot
- To: mathgroup at smc.vnet.net
- Subject: [mg91852] Re: Real and Complex Roots presented in a single plot
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Wed, 10 Sep 2008 05:09:04 -0400 (EDT)
- Organization: The Open University, Milton Keynes, UK
- References: <ga5koj$r82$1@smc.vnet.net>
Narasimham wrote:
> z[x_] = 1.3 Sin[1.7* x] + 0.6 Sin[4* x] ;
> Plot[z[x], {x, 0, 18}]
>
> In the above plot we can see all the real roots of z very
> approximately at {2.1,3.4, 5.5, 7.6, 9, 11, 13.1, 14.5, 16.5}, as the
> curve crosses x-axis.
>
> We can also recognize and see the real parts of all the
> complex roots of z where the curve is nearest to x -axis at {1.4,
> 4.2, 6.6, 9.6, 12.3, 15.3, 17.7}. They are near to x-values where the
> local maxima/minima occur.But we cannot 'see' their complex parts, as
> they need to be computed.
>
> After computation of all real and complex roots I would like to
> represent all roots, real and complex roots alike, on an ( x - y )
> Argand diagram (in a 2D plot) with the real part on x - axis and
> imaginary part on y - axis. How to do this ?
Assuming I have correctly understood what you are looking for,
visualizing real as well as complex roots can be achieve by formatting
the list of solutions as a list of pairs {Re, Im} that can be displayed
with ListPlot. For instance,
z[x_] = 1.3 Sin[1.7*x] + 0.6 Sin[4*x]
g = Plot[z[x], {x, 0, 18}]
sols = Reduce[Rationalize[z[x], 0] == 0, x] // N
Select[List[sols /. C[1] -> 0 // ToRules], (0 <= Re[x /. #] <= 18 &)]
pts = Through[{Re, Im}[x /. %]] // Transpose
p = ListPlot[pts, PlotStyle -> {Orange, PointSize[Large]}]
Show[g, p]
HTH,
-- Jean-Marc