function harmonic_3
%本程序用来带通滤波;
clc
clear all
fs=4096;%采样频率;
N=1000;%采样点数;
f=fs*(0:N-1)./N;%频率轴;
dt=1/fs;
t=0:dt:(N-1)*dt;%时间轴;
t=t/max(t);
s=40*sin(2*pi*12*t)+0.05*sin(2*pi*150*t);
fs0=fft(s,N);%信号的傅里叶变换;
fN=fs/2;%分析频率;
m=11;%实际带宽下限;
n=13;%实际带宽上限;
Wd=hann(N);%汉宁窗;
k=1/2*(n+m)/N;
fsai=(exp(2*n*pi*(t-k/(n-m))*i)-exp(2*pi*m*(t-k/(n-m))*i))./(2*pi*(n-m)*(t-k/(n-m))*i).*Wd';%谐波小波基;
fsai=fsai./sum(abs(fsai));
f1=conj(fft(fsai,N));%谐波小波基的傅里叶变换的共轭;
W=f1.*fs0;
w1=ifft(W);
plot(t,real(w1));
评论6