There was one mistake and its symptoms in my 3rd extension to lennon310's answer. 4th extension to lennon310's answer
I run
h = hamming(64);
h2 = hamming(38);
h = conv(h, h2);
B = 0; T = 0; F = 0;
data1 = filter(data(1 : 64),1,h); [B,T,F] = tfrwv(data1, 1:length(data1), length(data1));
for i=1:133
data1 = filter(data( 1 + i*37 : 64 + i*37 ),1,h); [b,t,f] = tfrwv(data1, 1:length(data1), length(data1));
B = [B b']; T = [T t]; F = [F; f];
end
data1 = filter(data(4959 : 5001),1,h); [b,t,f] = tfrwv(data1, 1:length(data1), length(data1));
B = [B b']; T = [T t]; F = [F; f];
T = 49.8899*T; % dummy constant to get appropriate time interval
and get pictures like this
I have not managed to show all thin peaks in one picture. A new question about it here.
I am plotting this by
t = 1/360; % 360 samples per second
fs = 360.5;
imagesc(T*t, F*fs, abs(B))
The algorithm is accumulating the points to the right dimension. I am not sure if multiplying by the dummy constant is the right way to go earlier.