compare.m

希尔伯特黄变换的程序源码

clear
clc
%加载数据
load trans_nodamg_var2to5.txt
load trans_e607_damg5perc_var2to5.txt
%------------------------
l_time=size(trans_nodamg_var2to5(:,1),1);
time=trans_nodamg_var2to5(:,1);
node2034_e607_nodamg=trans_nodamg_var2to5(:,3);
node2034_e607_damg5perc=trans_e607_damg5perc_var2to5(:,3);
%--------------------------
% for i=1:size(node2034_e607_nodamg)
%     damg5perc_noise(i,1)=node2034_e607_nodamg(i,1)+rand(1)*0.00003*abs(node2034_e607_nodamg(i,1)); %to add the noise
% end
%------------------------------------------
shuju=[time,node2034_e607_nodamg];
% shuju=[time,node2034_e607_damg5perc];
%--------------------------
shuju1=shuju(:,1);%时间
f=shuju(:,2);
r=f;
shuju2=shuju(:,2);%位移
%------------------------------------
%------------------------------------
%emd分解
imf = emd(shuju2);
[A,f,tt] = hhspectrum(imf,shuju1,1,1);
%-----------------------------------
%画出各阶imf分量
% for i=1:size(imf,1)
%     subplot(size(imf,1),1,i)
%     plot(imf(i,:))
% end
emd_visu(shuju2,shuju1,imf)
pause
close