|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
本帖最后由 wmjoan 于 2011-6-12 17:27 编辑
% ]5 d8 h' S' R, D1 T8 h9 N" o# S$ v, e, a1 ?
最近学习小波,想问一下我是不是可以这样理解,小波的细节信号就是类似于通过高通后的信号,其中d1的频率比d2,d3的高? 7 T* o$ e6 e5 p% f
程序: 前一段程序是以前看吉伯斯现象的时候写的,高通和小波的程序都是搬的网上的又改的,不知道选的对不对,有没有影响。 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作比较,
X/ E0 n( o! B
8 X8 U! D. |5 h3 A# Z5 Z
放大300点处,
高通滤波,
4 ^4 G- E. E+ g' F1 P! n
小波滤波,d1,
类似于,
% H7 d2 Q3 u, A, e3 x
问题: 1.是否对于小波变换,突变点就是模极大值点;而对高通滤波来说,原函数的突变点就是滤波后分量的突变点,例如图中300点处? 2.小波变换模极大值点的正负代表什么含义,对于此例,如上图所示,突变点为100,300,500,700,900点 高通滤波器的结果也像是有极性似的,上上下下的。 5 I4 q: m) x- L& e2 k" r' t: m
我现在弄得乱七八糟,没什么是确定的了。 0 Z' H/ x4 z- o; u" V: }( |
# |8 s# q. V$ D
|