carpet of frequencies over time - ideas or hints?

[Erhy]

Hello,

 for a music choir I want to generate

a 3D plot with spectral frequencies over time
and traces of the "should be" keynotes.

Thank you for discussion

Erhy

PS. My first steps of ftt over a piece of music:

stacksize('max')
w=loadwave("H:\temp\abc.wav")
y=w(2,:)  // take only one channel
clear w
// because mono, no channel select y=w(1,:);
Fs = 44100;                   // Sampling frequency
L = length(y);                 // Number of samples
NFFT = L + modulo(L,2);         // uprounded even number of samples
Y = fft(y)/L;                  //fft with array of L samples delivers array
of L complex values
f = Fs/2*linspace(0,1,NFFT/2+1);
p = 2*abs( Y(1:NFFT/2+1) );
plot(f,p)  // over all few
[m,i]=maxi(p(:));  // find index of the max.
scf(2);
plot(f((i-1000):(i+1000)),p((i-1000):(i+1000)))  // detail few