БПФ, поиск максимума спектральной плотности, поиск экстремума (максимума) в отсчетах БПФ - Matlab - Ответ 9463621

[ans2,nans]=max(pyy)

%% Параметры
clear all% Очистка памяти
load sif.mat
N=16384;%Общее количество отсчетов БПФ
Fd=20000;%Частота дискетизации (Гц)
Tm=(N-1)/Fd;% Длина сигнала (с)
FftL=2048;% Количество линий Фурье спектра
 
%% Генерация рабочих массивов
T=0:1/Fd:Tm;% Массив отсчетов времени
%% Спектральное представление сигнала
FftS=abs(fft(sig,FftL));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./FftL;% Нормировка спектра по амплитуде
FftS(1)=FftS(1)/2;% Нормировка постоянной составляющей в спектре
 
%% Построение графиков
subplot(2,1,1);% Выбор области окна для построения
plot(T,sig);% Построение сигнала
title('Сигнал');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда (Попугаи)');% Подпись оси у графика
subplot(2,1,2);% Выбор области окна для построения
F=0:Fd/FftL:Fd/2-Fd/FftL;% Массив частот вычисляемого спектра Фурье
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда (Попугаи)');% Подпись оси у графика