Those 101 signals in "hs" are in different frequency ranges defined by variable f. Choose those indices in hs that are within the required range specified in f then add those signals together and it should give you what you want. For example if f(2) to f(5) are bigger than 0.5 Hz and less than 5 Hz then your final time domain signal should be hs(2,:) + hs(3,:) +... hs(5,:)

Otherwise just use a normal IIR or FIR filter.

Hope this help.

Mani

Hi mani entezami

Thank you for giving reply.

[hs,f,t] = hht(imf,fs);

f(2) to f(11) values are in between 0.5 to 5.

I added hs(2,:)+hs(3,:)+. . . +hs(11,:)

But I didn't get a time-domain signal of 0.5 to 5Hz frequency range. The resultant signal is similar to the original signal of the frequency range in between 0 to 50Hz

I also compared the resultant time domain signal with FIR based results. FIR gives believable results. But I want to do such extraction using EMD.

A = full(hs); % Convert sparse matrix to full storage

resq=zeros(1,100);

for i=2:11

resq=resq+A(i,:);

end