Re: Real and Complex Roots presented in a single plot
- To: mathgroup at smc.vnet.net
- Subject: [mg91896] Re: Real and Complex Roots presented in a single plot
- From: magma <maderri2 at gmail.com>
- Date: Fri, 12 Sep 2008 05:26:49 -0400 (EDT)
- References: <ga5koj$r82$1@smc.vnet.net>
>
> 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.
Please note that your claim that the real part of the complex roots is
found at the local minima of the given function as the indipendent
variable moves on the real axis, is not valid in general and wrong in
this particular instance.
For example in the interval {1,2} there are 2 local extrema found
using the derivative:
f[z_] := 1.3 Sin[1.7 z] + 0.6 Sin[4 z]
The derivative:
In[43]:= fp[z_] := D[f[z], z]
In[44]:= fp[z]
Out[44]= 2.21 Cos[1.7 z] + 2.4 Cos[4 z]
In[73]:= extrema = z /. FindRoot[fp[z], {z, #}] & /@ {1.4, 1.6}
Out[73]= {1.33751, 1.69029}
The first is the local minimum (see the plot).
But the complex roots in interval {1,2} are
1.27946 + 0.374308 I and 1.27946 - 0.374308 I
The real part is 1.27946, while the minimum was at 1.33751
This can be visualized on the complex plane using Presentations:
f[z_] := 1.3 Sin[1.7 z] + 0.6 Sin[4 z]
Needs["Presentations`Master`"]
fp[z_] := D[f[z], z]
extrema = z /. FindRoot[fp[z], {z, #}] & /@ {1.4, 1.6}
roots = z /.
FindRoot[
f[z], {z, ToComplex@#}] & /@ {{-0.005546, -0.02102}, {1.285,
0.3771}, {1.279, -0.3678}, {2.081,
0.01109}, {3.462, -0.008173}, {4.194, -0.3485}, {4.213,
0.345}} // Chop;
In the following plot we see the roots (blu) and the extrema (red) in
the interval {1,2}
With[{zmin = -.5 - I, zmax = 5 + I,
contourlist = {0, 0.1, 0.2, .5, .505, .5095, .51, .53, 0.55, 0.57,
0.6, .61, .62, .6205, 1}},
Draw2D[{ComplexCartesianContour[f[z], {z, zmin, zmax}, Abs,
Contours -> contourlist, ColorFunctionScaling -> False,
ColorFunction -> (ContourColors[contourlist,
ColorData["SolarColors"]][#] &), PlotPoints -> {30, 15},
MaxRecursion -> 3, PlotRange -> {0, 11}],
ComplexCirclePoint[#, 3, Black, Blue] & /@ roots,
ComplexCirclePoint[#, 3, Black, Red] & /@ extrema,
ComplexLine[{-5, 5}]}, AspectRatio -> Automatic,
PlotRange -> {{1, 2}, {-.7, .7}}, Frame -> True,
FrameLabel -> {Re, Im}, RotateLabel -> False,
PlotLabel -> Row[{"Modulus of ", f[z]}], ImageSize -> 900,
Background -> Legacy@Linen]]