MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Plotter for complex polynomials (complex coefficients)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg124896] Re: Plotter for complex polynomials (complex coefficients)
  • From: Chris Young <cy56 at comcast.net>
  • Date: Sat, 11 Feb 2012 06:36:33 -0500 (EST)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • References: <jgauk3$eu3$1@smc.vnet.net> <jgdmtm$qo5$1@smc.vnet.net>

On 2012-02-02 10:01:26 +0000, DrMajorBob said:

> This seems a little faster:
> 
> Manipulate[
>   Module[{f, z},(*Convert 2D point to complex point*)
>    f[z_] = a.{1, I} z^3 + b.{1, I} z^2 + c.{1, I} z + d.{1, I};
>    Plot3D[Abs@f[x + y I], {x, -6, 6}, {y, -6, 6}, PlotPoints -> 100,
>     MaxRecursion -> 2, Mesh -> 11,
>     MeshStyle -> Directive[Gray, AbsoluteThickness[0.01]],
>     MeshFunctions -> ({x,
>         y} \[Function] (\[Pi] - Abs@Arg[f[x + I y]])/\[Pi]),
>     ColorFunctionScaling -> False,
>     ColorFunction -> ({x, y} \[Function]
>        Hue[0.425 \[LeftFloor]12 (\[Pi] -
>                Abs@Arg@f[x + I y])/\[Pi]\[RightFloor]/12, sat, bri]),
>     PlotStyle -> Opacity[opac],
>     AxesLabel -> {"x", "i y",
>       "|f(x + iy)|"}]],(*Item["The complex coefficients"],*){a, {-2, \
> -2}, {2, 2}}, {b, {-2, -2}, {2, 2}}, {c, {-2, -2}, {2,
>     2}}, {d, {-2, -2}, {2, 2}}, {{opac, 0.75, "Opacity"}, 0,
>    1}, {{sat, 0.75, "Saturation"}, 0, 1}, {{bri, 1, "Brightness"}, 0,
>    1}, ControlPlacement -> {Left, Left, Left, Left, Bottom, Bottom,
>     Bottom}]
> 
> Bobby

Thanks, I incorporated this into the latest plotter.

http://home.comcast.net/~cy56/Mma/ComplexCoeffPlotter.nb

http://home.comcast.net/~cy56/Mma/ComplexCoeffPlotterPic.png

Took me a long time to figure out that the two options

   PlotRange -> {{-2, 2}, {-2, 2}, {6, 0}},
   PlotRangePadding -> None

have to be in the same "option set" in order to really cut out the 
padding, so that the grid will go right through the dots.

Chris
cy56 at comcast.net


Manipulate[
 Module[
  {\[ScriptCapitalC], \[ScriptCapitalH],
   f, z, 
   arg},
  
  \[ScriptCapitalC] =  ( \[ScriptCapitalP] \[Function] \
\[ScriptCapitalP] . {1, I});                 (* 
  Convert 2D vector to complex point *)
  \[ScriptCapitalH] = ({\
\[ScriptCapitalP], h} \[Function] Append[\[ScriptCapitalP], h]);   (* 
  Convert 2D point \[ScriptCapitalP] to 3D point with height h *)
  
  (* The complex polynomial, with complex coefficients: *)
  
  f[z_] = (\[ScriptCapitalC] @ a) z^3 + (\[ScriptCapitalC] @ 
       b) z^2 + (\[ScriptCapitalC] @ c) z + (\[ScriptCapitalC] @ d);
  
  (* The "complex argument" (i.e., angle), 
  modified to run from 0 at -\[Pi] to 0.4 at 0 and back again, 
  in 12 steps: *)
  
  arg[z_] = 1/12 (0.4) Round[12  (\[Pi] - Abs[ Arg[z]])/\[Pi]]; 
  
  Show[
   Graphics3D[
    {
     Red,        Sphere[\[ScriptCapitalH][a, 0], dotR],
     Yellow, Sphere[\[ScriptCapitalH][b, 0], dotR],
     Green,    Sphere[\[ScriptCapitalH][c, 0], dotR],
     Blue,      Sphere[\[ScriptCapitalH][d, 0], dotR]
     }
    ],
   
   Plot3D[
    Abs[f[x + y I]], {x, -3, 3}, {y, -3, 3}, 
    
    PlotPoints -> plotPts,
    MaxRecursion -> maxRecurs, 
    Mesh -> {12, 5},
    MeshStyle -> Directive[Gray, AbsoluteThickness[0.01]],
    
    ColorFunctionScaling -> False,
    
    MeshFunctions -> {
      ({x, y} \[Function] arg[f[x + y I]]),
      ({x, y, z} \[Function] z)
      },
    ColorFunction -> ({x, y} \[Function] Hue[arg[f[x + y I]], sat, bri]),
    PlotStyle -> Opacity[opac],
    ClippingStyle -> None
    ],
   
   Lighting -> "Neutral",
   Axes -> True,
   PlotRange -> {{-2, 2}, {-2, 2}, {6, 0}},
   PlotRangePadding -> None,
   BoxRatios -> {4, 4, 6},
   AxesLabel -> {"x", "y \[ImaginaryI]", "|f(x + iy)|"},
   FaceGrids -> {{0, 0, 1}, {0, 0, -1}}
   ]
  ],
 
 {{a, {1, 0}}, {-2, -2}, {2, 2}, 0.5},
 {{b, {0, 1}}, {-2, -2}, {2, 2}, 0.5},
 {{c, {-1, 0}}, {-2, -2}, {2, 2}, 0.5},
 {{d, {0, -1}}, {-2, -2}, {2, 2}, 0.5},
 
 {{opac, 0.95, "Opacity"}, 0, 1, 0.125},
 {{sat, 0.75, "Saturation"}, 0, 1, 0.125},
 {{bri, 1, "Brightness"}, 0, 1, 0.125},
 {{plotPts, 100}, 10, 200, 10},
 {{maxRecurs, 2}, 1, 4, 1},
 
 {{dotR, 0.05}, 0.01, 0.2, 0.01},
 
   ControlPlacement -> {Left, Left, Left, Left, 
   Bottom, Bottom, Bottom, Bottom, Bottom, Bottom}
 ]



  • Prev by Date: Re: Why does the order of down values come back?
  • Next by Date: Repackaging function arguments
  • Previous by thread: Re: Plotter for complex polynomials (complex coefficients)
  • Next by thread: Re: Trying to close some loose plotting