这里使用Scipy模块中的fft实现时域信号的FFT变换,如下:
时域信号:该信号为带有噪声的正弦信号经过小波去噪后的图像
转换结果:
T=wpdec(y,5,'db40')
%信号y进行波包解层数5T波树plot看
a10=wprcoef(T,[1,0])
%a10节点[1,0]进行重构信号貌似没层重构说吧能某层某节点进行重构节点编号波树
%以下为滤波程序(主要调节参数c的大小)
c=10
wn=0.1
fs=50000%采样频率;
b=fir1(c,wn/(fs/2),hamming(c+1))
y1=filtfilt(b,1,y)%对y滤波。