Создайте функциональный periodogram_data.m
это возвращает модифицированную оценку степени спектральной плотности (PSD) периодограммы входного сигнала с помощью окна. Функция задает много дискретных точек преобразования Фурье, равных длине входного сигнала.
type periodogram_data
Используйте codegen
сгенерировать MEX-функцию.
%#codegen
директива в функции указывает, что код MATLAB® предназначается для генерации кода.
-args
опция задает аргументы в качестве примера, которые задают размер, класс и сложность входных параметров к файлу MEX. В данном примере задайте inputData
как 1024 1 двойная точность случайный вектор и window
как Окно Хэмминга длины 1024. В последующих вызовах MEX-функции используйте входные сигналы с 1024 выборками и окна.
Если вы хотите, чтобы MEX-функция имела другое имя, используйте -o
опция.
Если вы хотите просмотреть отчет генерации кода, добавьте -report
опция в конце codegen
команда.
codegen periodogram_data -args {randn(1024,1),hamming(1024)}
Вычислите оценку PSD шумной синусоиды с 1024 выборками с помощью функции периодограммы и MEX-функции, которую вы сгенерировали. Укажите, что синусоида нормировала частоту рад/выборка и окно Hann. Постройте две оценки, чтобы проверить, что они совпадают.
N = 1024; x = 2*cos(2*pi/5*(0:N-1)') + randn(N,1); periodogram(x,hann(N)) [pxMex,fMex] = periodogram_data(x,hann(N)); hold on plot(fMex/pi,pow2db(pxMex),':','Color',[0 0.4 0]) hold off grid on legend('periodogram','MEX function')