声音信号希尔伯特黄变换

要求：把去噪后的信号（已有）希尔伯特黄变换得到  经验模态分解的结果 瞬时频率图 希尔伯特谱

clear all

close all

[y,Fs]=audioread('output.wav');

x1=decimate(y(:,1),4);

x2=decimate(y(:,2),4);

x1=x1';

indx=1:1000; 

s=x1(indx);

figure(10)

subplot(2,2,1);

plot(s);

title('原始信号');

[c,l]=wavedec(s,3,'db1');

a3=appcoef(c,l,'db1',3);

d3=detcoef(c,l,3);

d2=detcoef(c,l,2);

d1=detcoef(c,l,1);

dd3=zeros(1,length(d3));

dd2=zeros(1,length(d2));

dd1=zeros(1,length(d1));

c1=[a3 dd3 dd2 dd1];

s1=waverec(c1,l,'db1');

subplot(2,2,2);

plot(s1);grid;

title('强制消噪后的信号');

%% 

figure(20)

spectrogram(s1,100,80,100,Fs,'yaxis'), % 100,80,100是做变化的常数量

view(-27,62)

shading interp

plot_hht(s1,Fs)