Tony Reina

This software works well and can handle large signals but could use some improved documentation. However if you are very familiar with the Jones and Baraniuk paper you can piece together what it is doing and especially the form of the output matrix. I am a fan of the AOK-TFR method for adaptive time-frequency estimation.

There seems to a minor glitches in the figure produced in main_AOK.m:
The plot it produces is inverted so that the positive frequency components are shown as negative frequencies and vice versa. Also the final plot is not trimmed so that the time axis does not correspond to the time interval of the original signal and the time units on this final plot are in observation points not seconds (or other original time units of the signal)

Consequently the time axis is of questionable value, however this is easily fixed by using something like on the output matrix ofp:

contour(1/samp_freq:1/samp_freq:(length(real(sig_in_tmp)))/samp_freq,...
f_axis,...
(abs(ofp((nraf/2+1):(length(real(sig_in_tmp))+nraf/2),:)')));
xlabel('Time');
ylabel('Frequency (Hz)');