複数のグラフを重ねて表示する方法を教えてください

前提・実現したいこと

複数（今回は、２７個）のデータに対し、FFTを行い、その結果をすべて１つのグラフに表示したいのですが、色々調べているのですが良い方法が見つからず、実現できていません。簡単にできる方法について知りたいです。やはりfor文を用いるのでしょうか。ご教授お願いします。

発生している問題・エラーメッセージ

該当のソースコード

python
1import numpy as np
2import pandas as pd
3import matplotlib.pyplot as plt
4
5#データの読み込み
6df=pd.read_csv("/data/python検証.csv")
7
8#読み込んだデータの次元確認
9print("配列の次元:\n"+str(df.ndim))
10
11#データ数カウント（サンプル数）
12N=len(df)
13print("サンプル数:\n"+str(N))
14
15#サンプリング周波数
16fs=50000
17print("サンプリング周波数:\n"+str(fs))
18
19#窓関数
20window=np.hamming(N)
21
22#窓関数の補正値
23acf=(1/(sum(window)/N))
24
25#時間軸のデータ作成
26t=np.arange(0,N*(1/fs),1/fs)
27
28#周波数軸のデータ作成
29fq=np.linspace(0,fs,N) #周波数軸　linspace(開始,終了,分割数)
30
31#df= np.array(df["200[Hz]"]).reshape(-1)
32#df=np.ravel(df["200[Hz]"])
33for i in range(200,1550,50):
34    
35        #2次元配列から1次元配列に変更
36        data=df[str(i)+"[Hz]"].values
37        #窓関数後の信号
38        data_w=window*data
39        #読み込んだデータをフーリエ変換
40        data_FFT=np.fft.fft(data_w)
41        #FFTの複素数結果を絶対値に変換
42        data_FFT_abs=abs(data_FFT)
43        #窓補正
44        Data=(data_FFT_abs)*acf
45    
46        #print(str(i)+"[Hz]"+str(data_FFT))
47        #print(str(i)+"[Hz]:\n"+str(data_FFT_abs))
48    
49        """
50        #グラフ表示（時間信号）
51        plt.figure(figsize=(6,4))
52        plt.plot(t,data)
53        plt.title(str(i)+"[Hz]")
54        plt.xlabel("time[s]",fontsize=14)
55        plt.ylabel("signal intensity",fontsize=14)
56        plt.show()
57        """
58        #グラフ表示（PA信号）
59        plt.figure(figsize=(6,4))
60        #plt.plot(fq[:int(N/2)+1],data_FFT_abs[:int(N/2)+1])
61        plt.title(str(i)+"[Hz]")
62        plt.xlabel("frequency[Hz]",fontsize=14)
63        plt.ylabel("signal intensity",fontsize=14)
64        plt.xlim(0,2000)
65        plt.show()
66        if i==1500:
67            break
68
69#次元変更後の確認
70print("次元変化後の配列次元:\n"+str(data.ndim))
71