赞
踩
多谢董老师,董老师是个好老师!
心情久久不能平静,主要是高频这门课的分析方法实在是让我难以理解,公式也背不过,还是放放吧。
最近厌恶了Matlab臃肿的体积和频繁的读写对我的Mac的损害,所以学习了一下Python这一轻量级的脚本,发现“Python自诞生那天就跟科学计算分不开”这个事实。无聊,写写心得。
配置环境什么的还是弄了几个晚上的。在Mac下用PyCharm还是很好滴,装上NumPy,SciPy等等一众免费的,很不错的Python包,就可以灰了!
1.Hilbert变换及其在单边带(SSB)包络检波的应用
定义神马的,性质神马的自己百度去。我也懒得写公式了,大家将就着看。
先定义个东西,H(t)为Hilbert变换后的时域信号,f(t)为原始时域信号。那么其包络为:
Envelop = sqrt(H^2(t)+f^2(t))。
好了,写代码什么的都简单了。
import numpy as np
import pylab as pl
import scipy.signal as signal
from scipy import fftpack
t = np.arange(0, 0.3, 1/20000.0)
x = np.sin(2*np.pi*1000*t) * (np.sin(2*np.pi*20*t) + np.sin(2*np.pi*8*t) + 3.0)
hx = fftpack.hilbert(x)
pl.subplot(221)
pl.plot(x, label
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。