马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 wmjoan 于 2011-6-12 17:27 编辑
/ S( d* q8 o$ M, r U2 t" b! u4 @1 [
最近学习小波,想问一下我是不是可以这样理解,小波的细节信号就是类似于通过高通后的信号,其中d1的频率比d2,d3的高?
) v3 d/ { J# d7 o0 p! d' m
程序: 前一段程序是以前看吉伯斯现象的时候写的,高通和小波的程序都是搬的网上的又改的,不知道选的对不对,有没有影响。 clear a0=1; n=5000; f=0; t=0:0.01:10; for i=1:n f=f+2*sin(i*pi*2/4)/(i*pi)*cos(i*2*pi/4*t); end fs=10000;%高通滤波器 fc=10; [b,a]=butter(4,0.6,'high'); [h,w]=freqz(b,a); x=f; y=filter(b,a,x); [d,a]=wavedec(f,6,'db6');%db6族小波 a6=wrcoef('a',d,a,'db6',6); d6=wrcoef('d',d,a,'db6',6); d5=wrcoef('d',d,a,'db6',5); d4=wrcoef('d',d,a,'db6',4); d3=wrcoef('d',d,a,'db6',3); d2=wrcoef('d',d,a,'db6',2); d1=wrcoef('d',d,a,'db6',1); x是构造出的很多频率叠加的矩形波,y是用高通滤波器滤波后的结果,d1是小波变换后尺度1。 将x与y,d1作比较,
" I9 i3 L& T* P' M9 c) G
" F$ j7 b2 [+ F1 E8 ]0 e3 z$ _9 b
放大300点处,
高通滤波,
1 l q/ P% { C/ i2 Y6 a% f% C2 C
小波滤波,d1,
类似于,
9 S! I9 T3 ~; L: R9 n
问题: 1.是否对于小波变换,突变点就是模极大值点;而对高通滤波来说,原函数的突变点就是滤波后分量的突变点,例如图中300点处? 2.小波变换模极大值点的正负代表什么含义,对于此例,如上图所示,突变点为100,300,500,700,900点 高通滤波器的结果也像是有极性似的,上上下下的。
. k. @; f$ N- x我现在弄得乱七八糟,没什么是确定的了。 0 W! H1 r# E9 j/ }4 L! Y2 d: v" Y# Z
4 g4 Q1 t* ]: q* e
|