реклама на сайте
подробности

 
 
> Мтоды спектрального оценивания
blondex
сообщение Jan 9 2010, 10:52
Сообщение #1





Группа: Участник
Сообщений: 8
Регистрация: 9-01-10
Пользователь №: 54 682



Здрасте. Задачка по радиолокации.
Хочу реализовать в матлабе обработку сигналов с элементов линейной антенной решетки(АР). В результате нужно определить углы прихода сигналов, т.е. получить пространственный спектр и угловые координаты целей.
Имеется 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
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
bahurin
сообщение Jan 9 2010, 18:40
Сообщение #2


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



чтото вы не так написали по ходу.
Если есть массив

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)]^T,

то думается что аргумент любой из экспонент даст arg = (n-1)*2*pi*sin(tetta)*d/lambda). Поскольку n, d и лямбда известны то достаточно:

sin(tetta) = arg*lambda/((n-1)*2*pi*d);

ну и собсна tetta = asin от всего этого безобразия.

Однако мне лично после фразы
Цитата
Значение РБ(alpha) должно характеризовать мощность сигнала, распространяющегося в направлении tetta.


как то страшно за войска ПВО становиться cranky.gif biggrin.gif
Go to the top of the page
 
+Quote Post
blondex
сообщение Jan 10 2010, 09:23
Сообщение #3





Группа: Участник
Сообщений: 8
Регистрация: 9-01-10
Пользователь №: 54 682



Цитата(bahurin @ Jan 9 2010, 21:40) *
чтото вы не так написали по ходу.

Дело в том, что то, что я написала в векторе Х - это некая модель входных отсчетов. То есть таким образом я задаю (моделирую) входные сигналы на элементах АР. Тут нет шумов, и источник излучения всего один. Если в выражение добавить все остальное, предположить, что число источников равно К, все они разной мощности, добавить шумы, некоррелированные в разных каналах, то вектор входных отсчетов примет следующий вид:
X = [A1*exp(fi1)*(exp(j*2*pi*sin(tetta1)*d/lambda)+exp(j*2*pi*sin(tetta2)*d/lambda)+.....+exp(j*2*pi*sin(tettaК)*d/lambda));
A2*exp(fi2)*(exp(j*2*pi*sin(tetta1)*d/lambda)+exp(j*2*pi*sin(tetta2)*d/lambda);.......]

Короче, выражение усложняется, число неизвестных растет, да и на приемной стороне этого всего неизвестно, ни числа источников, ни отношения шумов, ни случайных фаз fii, которые генерируются randomом.

Кстати, со своей проблемой я наполовину разобралась, если увеличить число элементов АР до астрономических значений, то алгоритм все же работает и разрешает сигналы с близких направлений, вот только число элементов у меня 100 - многовато.

Может кто знаком с реализацией одного из методов, таких как MUSIC, метод максимального Правдоподобия, Прони и тд?
Цитата
как то страшно за войска ПВО становиться

А че с ними будет? Не за тех боитесь)
Go to the top of the page
 
+Quote Post
Bill G8s
сообщение Jan 10 2010, 13:21
Сообщение #4





Группа: Участник
Сообщений: 4
Регистрация: 20-02-07
Из: Microsoft
Пользователь №: 25 529



Цитата(blondex @ Jan 10 2010, 12:23) *
Может кто знаком с реализацией одного из методов, таких как MUSIC, метод максимального Правдоподобия, Прони и тд?


Может кто и знаком
http://dsp-book.narod.ru/marpl.htm
Go to the top of the page
 
+Quote Post
blondex
сообщение Jan 10 2010, 13:38
Сообщение #5





Группа: Участник
Сообщений: 8
Регистрация: 9-01-10
Пользователь №: 54 682



Цитата(Bill G8s @ Jan 10 2010, 16:21) *
Может кто и знаком
http://dsp-book.narod.ru/marpl.htm

Я читала марпла. Кстати, на этом сайте ссылки на все главы битые.
Все методы там описаны подробно,но книжку не попросишь объяснить, поделиться опытом и не спросишь почему программа не работает

Сообщение отредактировал blondex - Jan 10 2010, 13:39
Go to the top of the page
 
+Quote Post
Oldring
сообщение Jan 12 2010, 09:31
Сообщение #6


Гуру
******

Группа: Свой
Сообщений: 3 041
Регистрация: 10-01-05
Из: Москва
Пользователь №: 1 874



Цитата(blondex @ Jan 10 2010, 16:38) *
Все методы там описаны подробно,но книжку не попросишь объяснить, поделиться опытом и не спросишь почему программа не работает


Это непростая область. Нужно понимать используемую математику.
Во временной области ваши сигналы от целей коггерентны или нет?
Это один из ключевых вопросов для выбора методов решения задачи. Пока не поняли почему - Вы не понимаете ничего.


--------------------
Пишите в личку.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- blondex   Мтоды спектрального оценивания   Jan 9 2010, 10:52
|- - Bill G8s   Цитата(blondex @ Jan 10 2010, 16:38) Я чи...   Jan 10 2010, 15:11
|- - AndeyP   Может поможет книжка "Narrowband Direction of...   Jan 10 2010, 18:32
||- - blondex   AndeyP, ЦитатаМожет поможет книжка "Narrowba...   Jan 16 2010, 09:07
||- - Oldring   Цитата(blondex @ Jan 16 2010, 12:07) Сигн...   Jan 18 2010, 07:54
||- - vadimuzzz   Цитата(blondex @ Jan 16 2010, 15:07) А во...   Jan 19 2010, 01:48
||- - AndeyP   ЦитатаМожет поможет книжка "Narrowband Direct...   Jan 19 2010, 20:03
- - hobgoblin   Цитата(blondex @ Jan 9 2010, 13:52) Хочу ...   Jan 12 2010, 08:20
- - blondex   Наверное соболезновали вы мне рано, и вообще зря, ...   Jan 18 2010, 18:28
- - Oldring   Цитата(blondex @ Jan 18 2010, 21:28) Наве...   Jan 18 2010, 18:57
- - blondex   Цитата(Oldring @ Jan 18 2010, 21:57) Толь...   Jan 24 2010, 10:03
- - AndrewN   Цитата(blondex @ Jan 24 2010, 14:03) сигн...   Jan 24 2010, 13:55
- - Oldring   Цитата(blondex @ Jan 24 2010, 13:03) Пожа...   Jan 25 2010, 08:47
- - blondex   Цитата(Oldring @ Jan 25 2010, 11:47) Но в...   Jan 26 2010, 17:24
- - Oldring   Цитата(blondex @ Jan 26 2010, 20:24) Я мо...   Jan 26 2010, 21:10
|- - blondex   Цитата(Oldring @ Jan 27 2010, 00:10) Точн...   Jan 28 2010, 16:51
|- - Oldring   Цитата(blondex @ Jan 28 2010, 19:51) Прин...   Jan 28 2010, 17:26
- - AndrewN   Цитата(blondex @ Jan 26 2010, 20:24) Не с...   Jan 28 2010, 17:52


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 23rd July 2025 - 22:07
Рейтинг@Mail.ru


Страница сгенерированна за 0.01499 секунд с 7
ELECTRONIX ©2004-2016