|
Здрасте. Задачка по радиолокации.
Хочу реализовать в матлабе обработку сигналов с элементов линейной антенной решетки(АР). В результате нужно определить углы прихода сигналов, т.е. получить пространственный спектр и угловые координаты целей.
Имеется N элементов АР, с них снимаем пространственные отсчеты комплексной огибающей, получаем вектор отсчетов
X = [1, exp(j*2*pi*sin(tetta)*d/lambda), exp(j*2*2*pi*sin(tetta)*d/lambda).......exp(j*(n-1)*2*pi*sin(tetta)*d/lambda)]Т
Где d - расстояние между элементами АР, lambda - длина волны (интересно только их соотношение), tetta - угол прихода сигнала. Пусть источник один.
Для примера все амплитуды примем равными 1.
Далее формируем пространственно-корреляционную матрицу: R=X*XH
Дальше вводим опорный вектор, характеризующий плоскую волну волну с угла альфа:
F0 = [1, exp(j*2*pi*sin(alpha)*d/lambda), exp(j*2*2*pi*sin(alpha)*d/lambda).......exp(j*(n-1)*2*pi*sin(alpha)*d/lambda)]Т
alpha принимает по очереди все значения из заданного диапазона углов.
И производим, к примеру, оценку Бартлета:
PБ(alpha) = F0H*R*F0 (1)
Значение РБ[/b](alpha) должно характеризовать мощность сигнала, распространяющегося в нвправлении tetta.
Я строю получившуюся функцию PБ от alpha, и получаю график с размытым пиком, соответствующем заданному тетта и гиганскую погрешность, И абсолютную неспособность алгоритма разрешить две цели, находящиеся на любом угловом расстоянии. Если я задаю два угла тетта, то максимум функции Р один и соответствует углу, посередине между заданными.
Кто-нибудь знает как реализуются такие методы? И какой физический смысл выражения (1)?
Весь теоретический материал подчерпнут из книги Канащенкова "Защита радиолокационных систем от помех"
Сообщение отредактировал blondex - Jan 9 2010, 10:58
|
Мтоды спектрального оценивания
Jan 9 2010, 10:52
