sparse discrete fourier transform

Hi all, I am trying to use sparse discrete fourier transform on a time series in matlab. I could not find any sample code written in matlab in that regard. Would you please help me resolving that?!

I came across this website but the code in there is written in c.

Thanks,

Answer 1

[Matt J]

How about something like this,

 K=length(yourSignal);

 %Reduced DFT matrix
 M=exp( j*2*pi/K * bsxfun(@times,(0:K-1).', find(yourSignal).' - 1) );

 result=M*nonzeros(yourSignal),

[tafteh]

Thanks Matt, My question is that whether the sparse algorithm would be a different technique in power spectrum estimation? or it is just faster way of estimating using DFT method?

[Matt J]

What I presented would be a faster DFT for sparse input signals. If that's not the goal, you should clarify what is.

Answer 2

[dongyang gao]

I also want to this code maded in matlab, could you help me ?