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

clear all% Очистка памяти
load 'D:\Matlab\p3030.mat'% Загрузка файла
%% Параметры и переменные
N=16384;%Общее количество отсчетов БПФ
i=(1:N);%Отсчеты БПФ
Fd=16394;%Частота дискретизации (Гц)
%Сигнал в виде массива в файле данных
k=kaiser(N);%Весовая функция Кайзера-Бесселя
sig=p3030;
 
%% Построение временной диаграммы         
figure(1)%Временная диаграмма
plot(sig(i)),grid;%Построение графика с сеткой
title ('Временная диаграмма');%Название диаграммы
xlabel('Отсчеты БПФ(N)');%Метка диаграммы по оси Ox
ylabel('Амплитуда сигнала');%Метка диаграммы по оси Oy
 
%% Выполнение БПФ и вычисление спектральной плотности  
y=fft(sig,N);%БПФ сигнала (сигнал в комплексном виде)
spa = y.*conj(y)/N;%Вычисление спектральной плотности амплитуд(сигнал в числовом виде)
Udc = spa(1);%Постоянная составляющая сигнала
spa(1)=0;%вычитание постоянной составляющей из сигнала (первая ячейка массива=0)
 
%% Построение частотной диаграммы
figure(2)%Частотная диаграмма
f = Fd*(0:8191)/N;%Нормировка частотной шкалы (исключение зеркальной составляющей спектра)
plot(f, spa(1:N/2)), grid;%Построение графика с сеткой
title ('Спектр сигнала');%Название диаграммы
xlabel('Частота сигнала(Гц)');%Метка диаграммы по оси Ox
ylabel('Амплитуда сигнала');%Метка диаграммы по оси Oy
 
%% Поиск максимальной спектральной составляющей в области главного лепестка
[ans2,nans]=max(spa);