matlab 火柴人_小波分析检测信号奇异点matlab代码

【实例简介】

【实例截图】

【核心代码】

%信号奇异点检测

clear,clc,close;

load StartEndPressure;

signal = StartEndPressure';

dataNum = length(signal);

figure(1)

plot(signal)

xlabel('时间');ylabel('幅值');

title('频率突变信号');

figure(2)

[c,l]=wavedec(signal,5,'db6');

cfd=zeros(5,dataNum);

for k=1:5

d=detcoef(c,l,k);

d1=d(ones(1,2^k),:);

d2 = d1(:)';

cfd(k,:)= wkeep(d2,dataNum)

end

cfd=cfd(:);

%cfd1 = cfd;

I=find(abs(cfd)

%cfd2 = cfd;

%cfd2(I)=zeros(size(I));

cfd=reshape(cfd,5,dataNum);

colormap(pink(64));

img=image(flipud(wcodemat(cfd,64,'row')));

set(get(img,'parent'),'YtickLabel',[]);

title('离散小波变换后系数的绝对值')

ylabel('层数');

figure(3)

ccfs=cwt(signal,1:32,'haar','plot');

title('连续小波变换系数的绝对值')

colormap(pink(64));

ylabel('尺度')

xlabel('时间(或者空间)')