Как установить ось y в радианах на графике Python?

Question

Я хотел бы написать радианы единиц осей, как пропорциональных \pi: что-то вроде

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

вместо

0.785, 1.5707 ... 

Есть ли стандартный способ? В качестве примера, что я должен добавить к следующему коду?

from pylab import * 

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

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

Я нашел этот пример http://matplotlib.sourceforge.net/examples/units/radian_demo.html, но он не работает, потому что у меня нет модуля basic_units.

Спасибо!

Answer 1

жесткий код их фракции или принимают плавающие числа

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() 

результат

Answer 2

я написал небольшую функцию, которая дает Вам список меток:

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

Вы можете скачать файл basic_units.py here:

После того, как должен работать так:

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() 

В качестве альтернативы вы могли бы реализовать функции арктангенса, как они реализованы потому функция в basic_units.py