plot HHT in imagesc

dedenlaut

2012 2 29

0 个回答

6 次查看（30 天）

0 个投票

hello, can anyone help my problem,

I use script plot_hht.m from Mr. Alan Tan in Matlab http://www.mathworks.com/matlabcentral/fileexchange/19681-hilbert-huang-transform/content/plot_hht.m

% function plot_hht(x,Ts)
% Plot the HHT.
% plot_hht(x,Ts)
% 
% :: Syntax
%    The array x is the input signal and Ts is the sampling period.
%    Example on use: [x,Fs] = wavread('Hum.wav');
%                    plot_hht(x(1:6000),1/Fs);
% Func : emd
% Get HHT.
N=1024;
t=linspace(1,2,N);
Ts=t(2)-t(1);
fs=1/Ts;
x(1:512)=5*sin(2*pi*100*t(1:512))+5*sin(2*pi*35*t(1:512));
x(513:1024)=5*sin(2*pi*20*t(513:1024));
plot(t,x);
aa=emd(x,0.05,100);
for ii=2:size(aa')
   imf{ii-1}=aa(:,ii)  
end
% imf = emd(x);
for k = 1:length(imf)
   b(k) = sum(imf{k}.*imf{k});
   th   = angle(hilbert(imf{k}));
   d{k} = diff(th)/Ts/(2*pi);
   energy{k} = abs(hilbert(imf{k})).^2;
end
[u,v] = sort(-b);
b     = 1-b/max(b);
% Set time-frequency plots.
N = length(x);
c = linspace(0,(N-2)*Ts,N-1);
figure, 
for k = v(1:11) 
   hold on;
   plot(c,d{k},'k.','Color',b([k k k]),'MarkerSize',3,'linewidth',1.5);
   set(gca,'FontSize',8,'XLim',[0 c(end)],'YLim',[0 1/2/Ts]); xlabel('Time'), ylabel('Frequency');
end

I want the plot result in time-frequency-energy domain (x-y-z).... how I implement that in 'imagesc'

thank you

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Santhana Raj 2012-6-17

hi,

I hav a doubt in ur question..... In the algorithm, u r calc energy for every k value. So u get a scalar value for each time frequency graph. For a 3D plot to exist, you require the third variable (energy) to be present at all 5999 values of the time and frequency plot.

i mean to say tat, u require to know the energy at each and every point in the time frequency graph. But the code above calculates net energy in the plot.

kindly point me out, if am mistaken.

regards,

A.S.Raj

https://sites.google.com/site/santhanarajarunachalam/

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question