Pythonでグラフをリアルタイムで描画させる方法

pip3 install matplotlib

from matplotlib import pyplot as plt      # matplotlibをインポート

X = [0,1,2,3,4,5,6,7,8,9,10]
Y = [0,1,1,3,4,4,6,7,8,6,10]

plt.figure(figsize=(12,8))          # グラフのサイズを指定
plt.plot(X, Y)                      # X,Yをプロット
plt.xlim(0, 10)                     # X軸の表示範囲
plt.ylim(0, 10)                     # Y軸の表示範囲
plt.xlabel('X-axis')                # X軸ラベル
plt.ylabel('Y-axis')                # Y軸ラベル
plt.rcParams["font.size"] = 20      # フォントサイズ
plt.grid()                          # グリッドを表示
plt.show()                          # グラフを表示

import time

t0 = time.perf_counter()   #プログラム開始時刻を取得

while True:
   time.sleep(1)
   t1 = time.perf_counter()               #終了時刻を取得
   print('Elapsed time : %.2f'%(t1-t0))   #開始時刻-終了時刻(経過時間)を出力

import RPi.GPIO as GPIO
import time
import math                          # Mathモジュール
import smbus                         # SMBus用モジュール
from matplotlib import pyplot as plt # グラフ描画ライブラリ

class ADCDevice: 
    def __init__(self):           # コンストラクタ。初期化用の関数。
        self.cmd = 0
        self.address = 0
        self.bus=smbus.SMBus(1)   # SMBusの引数に1を指定する。ラズパイのI2Cバスの番号

class ADS7830(ADCDevice):                   # ADCDeviceクラスを継承
    def __init__(self):
        super(ADS7830, self).__init__()
        self.address = 0x48                 # ADS7830の初期アドレス（A1,A0端子を共にGND接続）  
        self.cmd = 0x84                     # コマンドバイト
        
    def analogRead(self):
        value = self.bus.read_byte_data(self.address, self.cmd)
        return value

    def close(self):
        self.bus.close()

adc = ADS7830()
plt.figure(figsize=(12,8))        # グラフのサイズを指定
plt.rcParams["font.size"] = 20    # フォントサイズ
plt.grid()                        # グリッドを表示
X, Y = [], []                     # XとYの空のリスト
start_time = time.perf_counter()  # プログラム開始時の現在の時刻

def setup():
    global adc
    
def loop():
    while True:
        value = adc.analogRead()                              # 指定したピンのAD値を読む
        voltage = value / 255.0 * 3.3                         # AD値を電圧に変換
        Rt = 10 * voltage / (3.3 - voltage)                   # サーミスタの抵抗値を計算
        tempK = 1/(1/(273.15 + 25) + math.log(Rt/10)/3950.0)  # ケルビン温度を計算
        tempC = tempK -273.15                                 # セルシウス温度に変換
        print('Time : %.2f, Temperature : %.2f'%(time.perf_counter() - start_time,tempC))
        
        if len(X) > 10:            # Xに10個以上格納されたら
                del X[0]           # Xの0番目(リストの左端)を削除
                X_max = X[0] + 14  # X軸の最大値設定
                del Y[0]           # Yの0番目(リストの左端)を削除
        else:
                X_max = 14         # Xのリストの要素が10個以下ならX軸の最大値は14とする

        plt.cla()                                    # 画面をクリア
        Y.append(float(tempC))                       # リストの末尾にアイテムを追加
        X.append(time.perf_counter() - start_time)   # プログラム開始からの経過時間をリストの末尾に追加
        plt.plot(X, Y)                               # X,Yをプロット
        plt.ylim(20, 30)                             # Y軸の表示範囲
        plt.xlim(X[0], X_max)                        # X軸の表示範囲
        plt.xlabel('Time [s]')                       # X軸ラベル
        plt.ylabel('Temperature [deg]')              # Y軸ラベル
        plt.pause(1)                                 # 更新時間間隔1秒で描画する

def destroy():
    adc.close()
    GPIO.cleanup()
    
if __name__ == '__main__':
    print ('Program is starting ... ')
    setup()
    try:
        loop()
    except KeyboardInterrupt:
        destroy()

if len(X) > 10:        # Xに10個以上格納されたら
    del X[0]           # Xの0番目(リストの左端)を削除
    X_max = X[0] + 14  # X軸の最大値設定
    del Y[0]           # Yの0番目(リストの左端)を削除
 else:
    X_max = 14         # Xのリストの要素が10個以下ならX軸の最大値は14とする