БПФ ОБПФ в звуковизоре Опции

[Acvarif]

Пытаюсь разобраться в математике обработки сигналов звуковизора http://www.cta.ru/cms/f/441553.pdf.
Возникают непонятки. Прошу помочь разобраться.
1. Может у кого имеется скан этой книжки http://urss.ru/cgi-bin/db.pl?lang=Ru&b...ok&id=69052
2. К звуковизору Нептун имеются некоторые матвыкладки но сомневаюсь правильно-ли я понимаю данную математику. Матвыкладки на 2_х страницах.
Все построено на ПЛИС. Имеется антенна 48 элементов. В воде через эти элементы излучается ЛЧМ сигнал длительностью 1мс средней частотой 450 кГц, и на эти-же элементы осуществляется прием отраженных сигналов. Далее принятые сигналы с каждого элемента фильтруются цифровым фильтром (частота выборок 600 кГц) и передаются по Ethernet на ПК. На картинках последовательные шаги обработки принятых сигналов.
Не понятны некоторые шаги данной последовательности обработки.
1. Преобразование Фурье входной выборки - тоесть Фурье по частоте. Зачем?
2. Формирование характеристик направленности антенны для каждой спектральной компоненты. Я так понимаю это должен быть БПФ по пространству. Проще говоря БПФ в котором входними элементами являются данные сразу из 48 элементов антенны. Таким образом можно сформировать ДФУ(диаграммно формирующее устройство). Тут вроде понятно, если это действительно так. Сомневаюсь правильн-ли я понял...
3. Корреляционная обработка сигнала в каждом сформированном угловом канале - это тоже понятно. Все деляется для увеличения разрешаюшей способности по дистанции.
4. Квадратурная обработка - что имеется ввиду? Что нужно проделать дальше с сигналами и с какими?

Если кому не сложно, поясните пожалуйста популярно данную математику... Как на ее базе пулучается изображение. В первую очередь зачем в этой цепочке первый пункт?

[serjj]

Есть N антенных элементов, требуется сформировать M лучей.

Во временной области для формирования задержек используем фильтры, на которых можно с высокой точностью эти задержки выставить. Допустим на такой фильтр требуется Р умножителей. Итого получаем, что на один набор входных отсчётов с N элементов потребуется N*M*P комплексных умножителей (КУ).

Используя формирование ДН в частотной области, потребуется делать N быстрых преобразований Фурье с базой K > P. На N каналов, это N*K*log2(K) КУ. Для формирования ДН в частотной области нужно N*K*M КУ. Чтобы вернуть лучи во временную области потребуется ещё N*K*log2(К) КУ. Итого получаем 2*N*K*log2(K)+N*K*M КУ. Но для частотного формирования обработка выполняется блоками, размер каждого блока равен базе Фурье, если нет перекрытия. Указанное количество КУ - на блок а не на выборку АЦП, тогда на выборку АЦП можно оценить количество КУ: 2*N*log2(K)+N*M. Если использовать перекрытие 50% (распространённый вариант), то количество КУ выборку = 4*N*log2(K) + 2*M*N.

Сравнивая N*M*P и 4*N*log2(K) + 2*M*N нетрудно оценить границы применимости частотного формирования ДН.

Из этих же рассуждений можно увидеть, где возможно сэкономить для частотного метода. Если возвращаться во временную область не требуется, то количество КУ можно сократить до 2*N*log2(K) + 2*M*N. Если нас интересуют не все частоты, то второе слагаемое может быть дополнительно уменьшено, т.к. формировать ДН нужно будет не для K частот, а для K' < K. Если принять M=N, то в частотной области операция будет похожа на DFT, с той разницей, что матрица преобразования W может не быть ортогональной. Существует ряд быстрых алгоритмов для данной ситуации на подобие FFT, тогда M*N -> N*N -> N*log2(N). И ещё одно достоинство частотного метода - если используются алгоритмы адаптации, их реализация в частотной области сводится к обычным sum&shift beamformer'ам, в отличие от адаптации фильтров для временных методов.

Сообщение отредактировал serjj - Nov 4 2015, 07:39

[Acvarif]

В общем и целом процесс формирования лучей в частотной области понятен, кроме некоторой конкретики.
Допустим имеется 2 датчика и на их базе необходимо сформировать 32 луча в частотной области с шагом 0.3 градуса. При этом сектор обзора для этих датчиков всего 10 градусов.
Если применить БПФ для каждого датчика, то сколькиточечный нужен будет БПФ, чтобы обеспечить такой мелкий шаг по углу? Получается 360*3 или ~ 1024 точки. Из них нужно будет выбрать только необходимые 32 луча. Или все же както можно выкрутиться БПФ с меньшим количеством точек?
Это все рассуждения к тому как достигается высокое разрешение по углу в звуковизоре? БПФ на большое количество точек? Если так то тогда и разрешение по дальности будет низкое, поскольку оно будет определяться шагами равными количеству точек БПФ * на частоту дискретизации.

[petrov]

В общем и целом процесс формирования лучей в частотной области понятен, кроме некоторой конкретики.
Допустим имеется 2 датчика и на их базе необходимо сформировать 32 луча в частотной области с шагом 0.3 градуса. При этом сектор обзора для этих датчиков всего 10 градусов.
Если применить БПФ для каждого датчика, то сколькиточечный нужен будет БПФ, чтобы обеспечить такой мелкий шаг по углу? Получается 360*3 или ~ 1024 точки. Из них нужно будет выбрать только необходимые 32 луча. Или все же както можно выкрутиться БПФ с меньшим количеством точек?
Это все рассуждения к тому как достигается высокое разрешение по углу в звуковизоре? БПФ на большое количество точек? Если так то тогда и разрешение по дальности будет низкое, поскольку оно будет определяться шагами равными количеству точек БПФ * на частоту дискретизации.

Разрешение угловое определяется эффективной площадью антенны. А требуемое количество полос тем, насколько диаграмма направленности зависит от частоты, если не зависит, то одна точка БПФ нужна, а разрешение при этом может быть сколь угодно высоким. Банк фильтров можно сделать вообще без децимации, БПФ не вещь в себе, а всего длишь быстрый алгоритм для нужного банка фильтров.

[Acvarif]

Разрешение угловое определяется эффективной площадью антенны. А требуемое количество полос тем, насколько диаграмма направленности зависит от частоты, если не зависит, то одна точка БПФ нужна, а разрешение при этом может быть сколь угодно высоким. Банк фильтров можно сделать вообще без децимации, БПФ не вещь в себе, а всего длишь быстрый алгоритм для нужного банка фильтров.

Стоп. Давайте отвлечемся от антенны. Речь идет о БПФ и формировании диаграммы направленности в частотной области. Предположим что диаграмма зависит от частоты. Значит разрешение по углу определяется какраз количеством точек БПФ. Или я не прав? http://www.eetimes.com/document.asp?doc_id=1278838 Иначе как можно получить 256 лучей при 48 элементах антенны? http://www.cta.ru/cms/f/441553.pdf

[petrov]

Речь идет о БПФ и формировании диаграммы направленности в частотной области. Предположим что диаграмма зависит от частоты. Значит разрешение по углу определяется какраз количеством точек БПФ. Или я не прав? http://www.eetimes.com/document.asp?doc_id=1278838

Разрешение определяется шириной основного лепестка диграммы направленности. Можно использовать FFT для суммирования сигналов с решётки с фазовращающими коэффициентами как по ссылке, количеством точек определяется количество углов, а не разрешение, и это всё для одной частоты частотно-зависимой диграммы направленности.

[Acvarif]

Разрешение определяется шириной основного лепестка диграммы направленности. Можно использовать FFT для суммирования сигналов с решётки с фазовращающими коэффициентами как по ссылке, количеством точек определяется количество углов, а не разрешение, и это всё для одной частоты частотно-зависимой диграммы направленности.

Ну вот. Значит количество точек это количество углов. Теперь если антенна (например линейка, нпример 32 элемента) имеет ширину основного лепестка например 30 градусов. Как с помощью бимформера в частотной области построить в пределах 30 градусов веер из 256 лучей? На сколько точек в этом случае нужет FFT на каждый элемент антенны?

[petrov]

Ну вот. Значит количество точек это количество углов. Теперь если антенна (например линейка, нпример 32 элемента) имеет ширину основного лепестка например 30 градусов. Как с помощью бимформера в частотной области построить в пределах 30 градусов веер из 256 лучей? На сколько точек в этом случае нужет FFT на каждый элемент антенны?

Фазовращающие коэффициенты при суммировании FFT знаете, конструкцию решётки знаете, вот и смотрите набор диаграмм направленности который получаете.