Deriving Data and Plotting

asked 2012-08-14 16:15:56 +0100

gjm
27 ●4 ●4 ●8

updated 2012-08-14 21:08:26 +0100

calc314
4191 ●22 ●48 ●113

In the following x means multiply.

Let

z1=2 x pi x 650 x 10^6
p1=2 x pi x 1.9 x 10^9
p2=2 x pi x 5 x 10^9
adc=.667
deltaF=.2 x 10^9

In the following j is the imaginary representation Now let N=(adc x p1 x p2)/z1 and

M=(-2 x j x pi x freq x z1)/((-2 x j x pi x freq+p1) x (-2 x j x pi x freq+p2))

The base formula is abs(20 x log(N x M,10)) or the absolute value of 20 log(N x M) where the log base is base 10. Within the base formula is freq (the value of the frequency) which will be swept from 10^8 to 10^11 in steps of deltaF (.2 x 10^9)

The vertical range is from -25 to 5 (linear scale) and the horizontal range is from 10^8 to 10^11 (log scale)

The horizontal axis as stated is a Log scale as opposed to Linear.

How do I set this up in Sage to plot per the above? Thanks.

edit retag flag offensive close merge delete

add a comment

3 Answers

Sort by » oldest newest most voted

answered 2012-08-14 22:50:51 +0100

calc314
4191 ●22 ●48 ●113

I'm not getting the vertical range that you are specifying, but here is what I've got. Note that the log plot only works in the latest release (5.2) of Sage.

z1=2 * pi * 650 * 10^6
p1=2 * pi * 1.9 * 10^9
p2=2 * pi * 5 * 10^9
adc=.667
deltaF=.2 * 10^9

N=(adc * p1 * p2)/z1
M(freq)=(-2 * i * pi * freq * z1)/((-2 * i * pi * freq+p1) * (-2 * i * pi * freq+p2))

g(frq)=20*abs(log(N*M(frq),10))

pts=[(frq,g(frq).n()) for frq in srange(10^8,10^11,deltaF)]

list_plot(pts, scale='semilogx')

edit flag offensive delete link

add a comment

answered 2012-08-15 15:34:28 +0100

gjm
27 ●4 ●4 ●8

It's okay, I got it to work with the abs defined where I wanted it and had one minor change to the equation to get the correct range. Now I need to find out how to get Sage 5.2.