[求助]使用Morlet小波实现时变带通滤波的困扰！
本帖最后由 loverpipi 于 2010-11-11 12:16 编辑最近想用Morlet小波做一个滤波器，但是滤波器原理的基础不是很好，所以想用卷积方式实现。具体代码如下：
clc;clear;
fs=1024;N=1024;%采样频率1024，采样点数1024
t=;
Sig=exp(j*2*pi*70*t)+exp(j*2*pi*100*t)+exp(j*2*pi*130*t);%信号为70，100，130HZ平稳信号
fc=100;B=200;%设置小波的中心频率和带宽
for i=1:N
fb=4/B/B;
dt=t(i); %时间延迟
Morl=(pi*fb).^(-0.5).*exp(i*2*pi*fc*(t-dt)-(t-dt).^2/fb);
fy1=conv(Sig,Morl)*2/length(Sig);%卷积滤波
yy(i)=fy1(i); %yy为滤波结果，但是无法实现。使用传统的滤波器可以实现
end
我的想法是对于每个点处使用小波滤波，然后将所有点连起来实现时变得滤波器功能。使用平稳信号时为了方便比较，现在是无法实现这个功能。请高手指点一下，看看哪里出问题了，或者有其它做法也行(主要是要用Morlet小波实现)
页:
[1]