python signal matplotlob パワースペクトル密度

Question

pythonのビギナーです．

下記のリンク先のコードを書くと
https://org-technology.com/posts/power-spectral-density.html
リンク先の通りのグラフが出力されます．

ただ，ひとつ疑問のですが，
freq1, P1 = signal.periodogram(y, fs)
というコードの表記は
freq1とP1は同じものであるという事でしょうか？

P1とはたぶんPower（freq1の2乗）だと思うのですが
この表記の意味がよくわかりません．
どなたかご教授ください．

import numpy as np
from scipy import signal
import matplotlib.pyplot as plt

n = 1024
dt = 0.001
fs = 1/dt
f1 = 120
f2 = 150
t = np.linspace(1, n, n)*dt-dt
y = np.sin(2*np.pi*f1*t)+2*np.sin(2*np.pi*f2*t)+0.1*np.random.randn(t.size)

freq1, P1 = signal.periodogram(y, fs)


plt.figure()
plt.plot(freq1, 10*np.log10(P1), "b", label="periodogram")
plt.ylim(-60, 0)
plt.legend(loc="upper right")
plt.xlabel("Frequency[Hz]")
plt.ylabel("Power/frequency[dB/Hz]")
plt.show()

Answer 1

freq1とP1は同じものであるという事でしょうか？

いいえ、違います。
これはsignal.periodogramが返す２つの戻り値をそれぞれfreq1とP1に代入している記述です。
関数の仕様を見るとf(周波数[ndarray型]), Pxx(パワースペクトル[ndarray型])の２つの戻り値を返すようです。
なので、freq1には周波数、P1にはパワースペクトルが入ります。

scipy.signal.periodogramの仕様
https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.periodogram.html