Pythonscipy之希尔伯特变换hilbert.pdfVIP

Pythonscipy之希尔伯特变换hilbert

Python/scipy之希尔伯特变换hilbert

具体来源可参考：

⽂档中说明返回的信号是Analyticsignal解析信号，与Matlab返回信号形式相同。（ps我使⽤转化⼀次时域对称信号，但是返回的值是

实数，还没有深究是什么问题）。在⽂中也是把结果当做希尔伯特变换来看的。

希尔伯特变换获取包络

22

信号f(t)的包络为E(t)=sqrt([H(f(t))]+f(t))

importnumpyasnp

importmatplotlib.pyplotasplt

importscipy.signalassignal

fromscipyimportfftpack

t=np.arange(0,0.3,1/20000.0)

x=np.cos(2*np.pi*1000*t)*(np.cos(2*np.pi*20*t)+np.cos(2*np.pi*8*t)+3.0)

hx=fftpack.hilbert(x)

plt.plot(x,label=uu)

plt.plot(np.sqrt(x**2+hx**2),,linewidth=2,label=uu)

plt.plot(-np.sqrt(x**2+hx**2),,linewidth=2,label=uu)

plt.title(uu)

plt.legend()

plt.show()