马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 wmjoan 于 2011-6-12 17:27 编辑
{; |1 h1 n* Y! a# t% f" m+ v* L' E, @, \# l$ V6 M: D
最近学习小波,想问一下我是不是可以这样理解,小波的细节信号就是类似于通过高通后的信号,其中d1的频率比d2,d3的高? / i. Z" }) K2 x6 B4 [
程序: 前一段程序是以前看吉伯斯现象的时候写的,高通和小波的程序都是搬的网上的又改的,不知道选的对不对,有没有影响。 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作比较,
$ J6 C6 P1 C7 S; _9 H3 W9 w# ^' `2 b
A; S% R3 w: A; @2 l
放大300点处,
高通滤波,
- j9 I8 x' E+ M5 q
小波滤波,d1,
类似于,
/ S* h" g5 v% t- r& r1 {
问题: 1.是否对于小波变换,突变点就是模极大值点;而对高通滤波来说,原函数的突变点就是滤波后分量的突变点,例如图中300点处? 2.小波变换模极大值点的正负代表什么含义,对于此例,如上图所示,突变点为100,300,500,700,900点 高通滤波器的结果也像是有极性似的,上上下下的。 4 b$ F5 d, D9 D5 t
我现在弄得乱七八糟,没什么是确定的了。 - t8 P2 O# d1 w! ^6 }& E) A6 I+ t2 p1 x4 Y
! y& i9 C9 e0 h* g( s9 t: F
|