Digital System Processing (Audio Processing)
显示 更早的评论
I'm working on project to find the piano notes of a song. From the time domain graph, I was able to find the correct frequency where the note is present. I used: [S,F,T] = spectrogram (y, 512 , [] , 512, fs); to get the frequency of my song. But the size of the matrix S is too big and I cannot find the correcr frequency row. Can someone please help me understand how to get the correct row of frequency knowing the time?
Thanks
采纳的回答
Elimelech Schreiber
2020-12-9
编辑：Elimelech Schreiber
2020-12-9
If you have a specific time, then you should probably find the maximum frequency at that time, in the matrix S.
Each column in the matrix S contains the magnitudes of all the processed frequencies at a certain timestep.
e.g. Lets say t=1.13[sec].
first find the correct time index in T, which index 'i' yields T[i] =~ t (= 1.13) :
[~, i ] = min(abs(T - t))
lets say you've found i = 11,
now find the INDEX of the maximum magnitude in the i-th column of S:
[~, maxIdx] = max(S(:, i) ) % S(:, i)- The i-th column, containing freq. magnitudes at timestep i
finaly convert that index into a frequency, using F:
maxFrequency = F[maxIdx];
5 个评论
Lina Ngoopos
2020-12-9
Elimelech Schreiber
2020-12-9
编辑：Elimelech Schreiber
2020-12-9
no, i is in the interval: 0 <= i <= 2761.
i is the time index you should find, that matches the time you've found in the time domain.
So, back to my example, say you are interested in the note frequency playing at time t=1.13[sec].
You need to find which column of S coresponds to that time. one way of doing this is by lookig for the closest time in the time vector: T.
Throughout this answer i am using S,F,T from your question: [S,F,T] = spectrogram (y, 512 , [] , 512, fs);
Lina Ngoopos
2020-12-9
编辑：Lina Ngoopos
2020-12-9
Elimelech Schreiber
2020-12-9
Yes, replace it if you have already found it.
And skip that step: that line is meant to find the correct index.
Lina Ngoopos
2020-12-9
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Audio Processing Algorithm Design 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
