马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 wmjoan 于 2011-6-12 17:27 编辑 3 V5 g, B) W$ o* [& K P
_8 c: Z' @0 S. u$ Q. V( o
最近学习小波,想问一下我是不是可以这样理解,小波的细节信号就是类似于通过高通后的信号,其中d1的频率比d2,d3的高?
' K( `! [3 U' p/ ?
程序: 前一段程序是以前看吉伯斯现象的时候写的,高通和小波的程序都是搬的网上的又改的,不知道选的对不对,有没有影响。 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作比较,
" O2 L# b l, |( }5 r c, z f- G
# W; |0 m3 `. w( g. |
放大300点处,
高通滤波,
' {7 |4 ~7 C) M& e) b/ X
小波滤波,d1,
类似于,
2 x9 h8 e- f" a2 v9 W' W
问题: 1.是否对于小波变换,突变点就是模极大值点;而对高通滤波来说,原函数的突变点就是滤波后分量的突变点,例如图中300点处? 2.小波变换模极大值点的正负代表什么含义,对于此例,如上图所示,突变点为100,300,500,700,900点 高通滤波器的结果也像是有极性似的,上上下下的。
7 F- r! _6 c0 e7 y我现在弄得乱七八糟,没什么是确定的了。
/ ?: d/ `+ q5 }3 ^3 B! O
" d ?0 x& a: ^$ ]
|