Re: Re: Bode Plots in Mathematica
- To: mathgroup at smc.vnet.net
- Subject: [mg57280] Re: [mg57249] Re: Bode Plots in Mathematica
- From: "David Park" <djmp at earthlink.net>
- Date: Mon, 23 May 2005 02:20:51 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Nasser,
As I said in the earlier posting, the BodePlot routine is in the
DrawGraphics Help under Examples, LogPlots & Bode Plots. It is not a
built-in routine in DrawGraphics, but the code is there in the example and
can be copied and modified if desired.
The BodePlot allows plot options to be passed and these can override the
built-in options. Thus you could pass FrameLabel -> {"frequency
(rad/sec)",...} if you wished.
You could modify the nformat routine for the horizontal ticks, but I don't
think that 10^(-1), 10^0, and 10^1 look that great. Everyone to their own
choice.
It is certainly much more compact to plot phase and magnitude on the same
plot. Why repeat the whole frequency scale? Also it shows the relation
between the phase and magnitude more clearly. Separate plots could be made
by splitting the code and either using GraphicsArray or using two side by
side Rectangles. Much better would be one Rectangle above another. In that
case I would leave the ticks off the individual plots and use XSymbolicScale
to label the frequency values on the lower plot only, and YTickLine to label
the magnitude and phase values just outside the Rectangles. The light grid
lines would carry the frequencies values up through both rectangles. By
leaving the ticks and labels off the internal Rectangle plots you can easily
make them perfectly align.
David Park
djmp at earthlink.net
http://home.earthlink.net/~djmp/
From: Nasser Abbasi [mailto:nma at 12000.org]
To: mathgroup at smc.vnet.net
> David Park
> djmp at earthlink.net
> http://home.earthlink.net/~djmp/
>
>
>
David;
Nice work.
few comments (this is on 5.1):
The x-axis (the frequency axis) usually has the values displayed all
as powers of 10. i.e. as in 10^-2, 10^-1, 10^0, 10^1, etc...
Your output shows these as 10^-2, 0.1, 0, 10, but then uses 10^2 after
that, I think it is better to have all these shown in one format which
is 10^n.
You might want to label the x-axis also as 'frequency (rad/sec)' as
that is the standard instead of 's'.
The 'standard' also is to show the phase and magnitude on separate
plots and not on the same one.
You might want to label the y-axis as 'Phase(degrees)' and
'Magnitude(db)' instead of just 'db' and 'phase'.
Is the BodePlot[] function supposed to be inside one of your packages?
I had to cut/paste the code from the news server posting to define it
even though I just downloaded your packages and loaded
Needs["DrawGraphics`DrawingMaster`"]
Needs["Graphics`Graphics`"]
and thanks for your work.
Nasser