f, y1Welch1 = signal.welch(y1, fs, window='hanning', nperseg=256)
g, y1Welch2 = signal.welch(y1, fs, window='hanning', nperseg=128)
h, y1Welch3 = signal.welch(y1, fs, window='hanning', nperseg=64)
y1Welch1 = 10 * np.log10(y1Welch1)
y1Welch2= 10 * np.log10(y1Welch2)
y1Welch3 = 10 * np.log10(y1Welch3)

plt.plot(f[2:],y1Welch1[2:], label="NFFT = 256", color="red")
plt.plot(g[2:],y1Welch2[2:], label="NFFT = 128", color="orange")
plt.plot(h[2:],y1Welch3[2:], label="NFFT = 64", color="green")
plt.title("Szum gaussowski - Perdiodogram Welscha")
plt.ylabel("PSD(Db/Hz)")
plt.xlabel("Częstotliwość (Hz)")
plt.legend()
plt.show()

f, y1Welch1 = signal.welch(y2, fs, window='hanning', nperseg=256)
g, y1Welch2 = signal.welch(y2, fs, window='hanning', nperseg=128)
h, y1Welch3 = signal.welch(y2, fs, window='hanning', nperseg=64)
y1Welch1 = 10 * np.log10(y1Welch1)
y1Welch2= 10 * np.log10(y1Welch2)
y1Welch3 = 10 * np.log10(y1Welch3)

plt.plot(f,y1Welch1, label="NFFT = 256", color="red")
plt.plot(g,y1Welch2, label="NFFT = 128", color="orange")
plt.plot(h,y1Welch3, label="NFFT = 64", color="green")
plt.title("Sygnały sinusoidalne - Perdiodogram Welscha")
plt.ylabel("PSD(Db/Hz)")
plt.xlabel("Częstotliwość (Hz)")
plt.legend()
plt.show()

f, y1Welch1 = signal.welch(y3, fs, window='hanning', nperseg=256)
g, y1Welch2 = signal.welch(y3, fs, window='hanning', nperseg=128)
h, y1Welch3 = signal.welch(y3, fs, window='hanning', nperseg=64)
y1Welch1 = 10 * np.log10(y1Welch1)
y1Welch2= 10 * np.log10(y1Welch2)
y1Welch3 = 10 * np.log10(y1Welch3)

plt.plot(f[2:],y1Welch1[2:], label="NFFT = 256", color="red")
plt.plot(g[2:],y1Welch2[2:], label="NFFT = 128", color="orange")
plt.plot(h[2:],y1Welch3[2:], label="NFFT = 64", color="green")
plt.ylabel("PSD(Db/Hz)")
plt.xlabel("Częstotliwość (Hz)")
plt.ylabel("Częstotliwość (Hz)")
plt.legend()
plt.show()