Jak ustawić oś y w radianach na wykresie Pythona?

Question

chciałbym napisać Radian jednostek osi jako proporcjonalnych do \pi: coś

$\frac{\pi}{4}$, $\frac{\pi}{2}$, ... 

zamiast

0.785, 1.5707 ... 

Czy istnieje standardowy sposób? Jako przykład, co należy dodać do następującego kodu?

from pylab import * 

x=arange(-10.0,10.0,0.1) 
y= arctan(x) 

plot(x,y,'b.') 
show() 

Znalazłem ten przykład http://matplotlib.sourceforge.net/examples/units/radian_demo.html ale to nie działa, ponieważ nie mam moduł basic_units.

Dziękujemy!

Answer 1

ciężko kod ich frakcji lub akceptuje pływających numery

import matplotlib.pyplot as plt 
import numpy as np 

x=np.arange(-10.0,10.0,0.1) 
y=np.arctan(x) 

fig = plt.figure() 
ax = fig.add_subplot(111) 

ax.plot(x,y,'b.') 

y_pi = y/np.pi 
unit = 0.25 
y_tick = np.arange(-0.5, 0.5+unit, unit) 

y_label = [r"$-\frac{\pi}{2}$", r"$-\frac{\pi}{4}$", r"$0$", r"$+\frac{\pi}{4}$", r"$+\frac{\pi}{2}$"] 
ax.set_yticks(y_tick*np.pi) 
ax.set_yticklabels(y_label, fontsize=20) 

y_label2 = [r"$" + format(r, ".2g")+ r"\pi$" for r in y_tick] 
ax2 = ax.twinx() 
ax2.set_yticks(y_tick*np.pi) 
ax2.set_yticklabels(y_label2, fontsize=20) 

plt.show() 

wynik jest

Answer 2

napisałem małą funkcję, która daje z powrotem listę etykiet:

import numpy as np 
from fractions import Fraction 

def create_pi_labels(a, b, step): 

    max_denominator = int(1/step) 
    # i added this line and the .limit_denominator to solve an 
    # issue with floating point precision 
    # because of floataing point precision Fraction(1/3) would be 
    # Fraction(6004799503160661, 18014398509481984) 

    values = np.arange(a, b+step/10, step) 
    fracs = [Fraction(x).limit_denominator(max_denominator) for x in values] 
    ticks = values*np.pi 

    labels = [] 

    for frac in fracs: 
     if frac.numerator==0: 
      labels.append(r"$0$") 
     elif frac.numerator<0: 
      if frac.denominator==1 and abs(frac.numerator)==1: 
       labels.append(r"$-\pi$") 
      elif frac.denominator==1: 
       labels.append(r"$-{}\pi$".format(abs(frac.numerator))) 
      else: 
       labels.append(r"$-\frac{{{}}}{{{}}} \pi$".format(abs(frac.numerator), frac.denominator)) 
     else: 
      if frac.denominator==1 and frac.numerator==1: 
       labels.append(r"$\pi$") 
      elif frac.denominator==1: 
       labels.append(r"${}\pi$".format(frac.numerator)) 
      else: 
       labels.append(r"$\frac{{{}}}{{{}}} \pi$".format(frac.numerator, frac.denominator)) 

    return ticks, labels 

https://github.com/MaxNoe/python-plotting/blob/master/source/create_pi_labels.py

Answer 3

Możesz pobrać plik basic_units.py here:

Po to powinno działać tak:

from pylab import * 
from basic_units import radians 

x = arange(-10.0,10.0,0.1) 
y = map(lambda y: y*radians,arctan(x)) 
x = map(lambda x: x*radians,x) 

plot(x,y,'b.',xunits=radians,yunits=radians) 
show() 

Alternatywnie można zaimplementować funkcję arctan sposób ich wykonania cos funkcjonować w basic_units.py