小波变换

小波学习 2019.7.15《matlab》

%利用meyer小波基画连续小波变换系数分析图

t=0:0.2:2;%11

y=sin(2*t);

len=length(y);

cw1=cwt(y,1:32,‘meyr’,‘plot’);%矩阵=(尺度范围a)32*(信号长度b)11

title(‘系数’);

xlabel(‘信号长度’);

ylabel(‘尺度’);

【问题：怎么看图？？】

例2 ?mesh图

load vonkoch

vonkoch=vonkoch(1:510);

coefs=cwt(vonkoch,1:32,‘db9’);

colormap([0 0 0])

mesh(abs(coefs));

例3

x=zeros(1000,1);

x(500)=1;%点500处为1

cw1=cwt(y,1:128,‘haar’,‘plot’);%1281000-----纵坐标横坐标

title(‘系数’);

xlabel(‘信号长度’);

ylabel(‘尺度a’);

colormap jet

colorbar

%比较尺度系数10,50,90?的小波系数（共128）

subplot(311)

plot(cw1(10,:));%1*1000

title(‘尺度10’)

subplot(312)

plot(cw1(50,:));

title(‘尺度50’)

subplot(313)

plot(cw1(90,:));

title(‘尺度90’)

m=cw1(10,:);%1*1000

%%结论：尺度a越大,视野越宽（时域宽）–低频信号–大概