Силовая спектральная плотность синусоида в питоне с использованием scipy ⇐ Python

[Anonymous]

У меня есть следующий код, чтобы генерировать синусоидальный сигнал с некоторым полом шума: < /p>
Kmax = 1000000
fref = 125e6
fin = fref/5
t = 10*np.linspace(0,Kmax,Kmax)/fref

tdelay = (1+np.sin(2*np.pi*fin*t)) + np.random.uniform(0,0.2,Kmax)
plt.plot(t,tdelay)
plt.grid()
plt.show()
f, psd_td = signal.welch((tdelay - np.mean(tdelay)), fref, nperseg = len(tdelay) / 8)
plt.semilogx(f,10*np.log10(psd_td), label = 'PSD Signal')
plt.grid()
plt.show()
< /code>
Форма волны доменного домена выглядит ниже:
< /p>
Но спектральная плотность мощности не показывает меня. src = "https://i.sstatic.net/edjbm2wz.png"/>
Может ли кто -нибудь объяснить, что я делаю не так? И почему я не соблюдаю тон сигнала.

Подробнее здесь: https://stackoverflow.com/questions/796 ... sing-scipy