What's the best snipet for get the first batch of fft abs greater than half the maximum?

1 次查看（过去 30 天）
Georges Theodosiou
Georges Theodosiou 2022-10-21
编辑： Georges Theodosiou 2022-10-31
Ran in:
Dear sirs, pls let me ask your forgive for my previous wrong question. I expect catch first harmonic (FFT's pitch) by catcing the first batch of fft abs greater than half the maximum. I suppose, in human's voice and in violin's sound first harmonic's fft abs is greater than half the maximum, 2.3610e+08 in code. Then by max fucnion catch first harmonic's. What's the best snipet for catch first batch greater than half the "FFTabs"?
SampFreq = 16000;
Segm = 1:2048;
Pitch = 45;
FirstHarmAngles = Pitch*2*pi/SampFreq*Segm+1.9*pi;
SinFirstHarmAngles = sin(FirstHarmAngles);
SecondHarmAngles = Pitch*2*2*pi/SampFreq*Segm+2.9*pi;
SinSecondHarmAngles = sin(SecondHarmAngles);
ThirdHarmAngles = Pitch*3*2*pi/SampFreq*Segm+0.3*pi;
SinThirdHarmAngles = sin(ThirdHarmAngles);
Xn = 170000*SinFirstHarmAngles+220000*SinSecondHarmAngles+150000*...
SinThirdHarmAngles;
FFTXn = fft(Xn,16384);
FFTabs = abs(FFTXn);
plot(FFTabs(1:200));

请先登录，再回答此问题。

回答（1 个）

Georges Theodosiou
Georges Theodosiou 2022-10-26
编辑：Georges Theodosiou 2022-10-31
Dear sirs let me answer my qustion.
maxabs = max(FFTabs);
firstBatch = zeros(1,2048);
for a = 1:2048
firstBatch(a) = FFTabs(a);
if FFTabs(a)>maxabs/2 && FFTabs(a+1)<maxabs/2
break
end
end
[maxInFirstBatch, FFTpitch] = max(firstBatch)

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by