1つのグラフに複数プロットしたい

Question

前提・実現したいこと

現在、FFTプログラムを作成しています。
n個(n=1,2,3,.....)のデータを一つのグラフに重ねてプロットしたいと思い、for文を用いて試みているのですが、うまく実行することが出来ません。
現在のソースコードは該当ソースコードのようになっています。
ご教授して頂けると幸いです。

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

TypeError                                 Traceback (most recent call last)
<ipython-input-106-3b91932838c2> in <module>
     61 
     62         plt.figure(figsize=(6,4))
---> 63         plt.plot(fq[:int(N/2)+1],Data[:int(N/2)+1],color=colorlist(range(1,j+1)))
     64         plt.title(str(n)+"回目")
     65         plt.xlabel("frequency[Hz]",fontsize=14)

TypeError: 'int' object is not callable

該当のソースコード

import numpy as np
import pandas as pd
import os
import matplotlib.pyplot as plt

file_name = input("Please input filename:")
Sampling_frequency = int(input("Please input Sampling_frequency: "))
Number_of_iterations = int(input("Please input Number_of_iterations:"))

#読み込んだデータの次元確認
print("配列の次元:\n"+str(data.ndim))
print("配列次元変更後:\n"+str(data2.ndim))

#データ数
N=len(data)
print("サンプル数:\n"+str(N))

#サンプリング周波数
fs=Sampling_frequency
print("サンプリング周波数:\n"+str(fs))

#窓関数
window=np.hamming(N)

#窓関数の補正値
acf=(1/(sum(window)/N))
print("窓関数適用後のスペクトル補正値:\n"+str(acf))

#時間軸のデータ作成
t=np.arange(0,N*(1/fs),1/fs)

#周波数軸のデータ作成
fq=np.linspace(0,fs,N) #周波数軸　linspace(開始,終了,分割数)

#線色のリスト
color_list=["black","gray","lightcoral","brown","darkred","red","tomato","sienna","sandybrown","darkorange","tan","gold","darkkhaki"
            ,"yellowgreen","olivedrab","darkturquoise","lightblue"
            ,"steelblue","cornflowerblue","navy","blue","indigo","darkviolet","purple","magenta","crimson","lightpink"]

colorlist=len(color_list)

for n in range(1,Number_of_iterations+1):
    for j in range(0,colorlist):
        #データの読み込み
        data=pd.read_csv(os.path.join(r"Z:/private/結果/",file_name),header=None,usecols=list(range(1,n+1)),names=[f"{n}回目" for n in range(1,n+1)])
    
        #2次元配列から1次元配列に変更
        data2=data[str(n)+"回目"].values
    
        #窓関数後の信号
        data_w=window*data2
        
        #読み込んだデータをフーリエ変換
        data_w_FFT=np.fft.fft(data_w)
        
        #FFTの複素数結果を絶対値に変換
        data_w_FFT_abs=abs(data_w_FFT)
        
        #窓補正
        Data=(data_w_FFT_abs)*acf
        print(Data)
    
        plt.figure(figsize=(6,4))
        plt.plot(fq[:int(N/2)+1],Data[:int(N/2)+1],color=colorlist(range(1,j+1)))
        plt.title(str(n)+"回目")
        plt.xlabel("frequency[Hz]",fontsize=14)
        plt.ylabel("signal intensity",fontsize=14)
        plt.xlim(0,2000)
        plt.show()

試したこと

補足情報（FW/ツールのバージョンなど）

Comments

[meg_]

> for文を用いて試みているのですが、うまく実行することが出来ません。
どうなっているのでしょうか？

[keisokukunn]

エラーメッセージを載せることを忘れていました。申し訳ございません。
修正したので見て頂けると幸いです。

Answer 1

エラー発生原因は下記でしょう。

color_list=["black","gray","lightcoral","brown","darkred","red","tomato","sienna","sandybrown","darkorange","tan","gold","darkkhaki","yellowgreen","olivedrab","darkturquoise","lightblue","steelblue","cornflowerblue","navy","blue","indigo","darkviolet","purple","magenta","crimson","lightpink"]
colorlist=len(color_list)
j = 0
color=colorlist(range(1,j+1))
# TypeError: 'int' object is not callable

colorlist(range(1,j+1))ですが、まずcolorlistではなくcolor_listが適切でしょう。
color=colorlist(range(1,j+1))で実現したかったことが何かは分かりませんが、単純にcolor=color_list[j]で良いのではないでしょうか？