Скомпилируйте, с помощью codegen, функционального EstimateDOA.m. Эта функция оценивает направления прибытия (DOA) двух сигналов с добавленным шумом, которые получены стандартным Универсальным массивом строки (ULA) с 10 элементами. Антенна рабочая частота составляет 150 МГц и элементы массива, является распределенной половиной длины волны независимо. Фактическое направление прибытия первого сигнала является азимутом на 10 °, повышением на 20 °. Направление прибытия второго сигнала является азимутом на 45 °, повышением на 60 °. Сигналы и шум сгенерированы с помощью функции sensorsig.
function [az] = EstimateDOA() % Example: % Estimate the DOAs of two signals received by a standard % 10-element ULA with element spacing one half-wavelength apart. % The antenna operating frequency is 150 MHz. % The actual direction of the first signal is 10 degrees in % azimuth and 20 degrees in elevation. The direction of the % second signal is 45 degrees in azimuth and 60 degrees in % elevation. c = physconst('LightSpeed'); fc = 150e6; lambda = c/fc; fs = 8000; nsamp = 8000; sigma = 0.1; ang = [10 20; 45 60]'; antenna = phased.IsotropicAntennaElement('FrequencyRange',[100e6,300e6]); array = phased.ULA('Element',antenna,'NumElements',10, ... 'ElementSpacing',lambda/2); pos = getElementPosition(array)/lambda; sig = sensorsig(pos,nsamp,ang,sigma^2); estimator = phased.RootMUSICEstimator('SensorArray',array,... 'OperatingFrequency',fc,... 'NumSignalsSource','Property','NumSignals',2); doas = estimator(sig); az = broadside2az(sort(doas),[20,60]); end
Запустите codegen в командной строке, чтобы сгенерировать MEX-функцию, EstimateDOA_mex, и затем запустить MEX-функцию:
codegen EstimateDOA.m
EstimateDOA_mexПредполагаемые углы падения:
az = 10.0036 45.0030
Программа содержит фиксированное значение для шумового отклонения. Если вы хотели снова использовать тот же код для различного уровня шума, можно передать шумовое отклонение в качестве аргумента в функцию. Это сделано в функциональном EstimateDOA1.m, показанном здесь, который имеет входной параметр sigma.
function [az] = EstimateDOA1(sigma) % Example: % Estimate the DOAs of two signals received by a standard % 10-element ULA with element spacing one half-wavelength apart. % The antenna operating frequency is 150 MHz. % The actual direction of the first signal is 10 degrees in % azimuth and 20 degrees in elevation. The direction of the % second signal is 45 degrees in azimuth and 60 degrees in % elevation. c = physconst('LightSpeed'); fc = 150e6; lambda = c/fc; fs = 8000; nsamp = 8000; ang = [10 20; 45 60]'; antenna = phased.IsotropicAntennaElement('FrequencyRange',[100e6,300e6]); array = phased.ULA('Element',antenna,'NumElements',10, ... 'ElementSpacing',lambda/2); pos = getElementPosition(array)/lambda; sig = sensorsig(pos,nsamp,ang,sigma^2); estimator = phased.RootMUSICEstimator('SensorArray',sULA,... 'OperatingFrequency',fc,... 'NumSignalsSource','Property','NumSignals',2); doas = estimator(sig); az = broadside2az(sort(doas),[20,60]);
Запустите codegen в командной строке, чтобы сгенерировать MEX-функцию, EstimateDOA1_mex, с помощью опции -args, чтобы задать тип входного параметра. Затем запустите MEX-функцию с несколькими различными входными параметрами:
codegen EstimateDOA1.m -args {1} EstimateDOA1_mex(1) az = 10.0130 45.0613 EstimateDOA1_mex(10) az = 10.1882 44.3327 EstimateDOA1_mex(15) az = 8.1620 46.2440
Увеличение значения sigma ухудшает оценки углов азимута.