Форум Радиолюбителей

Оконные функции оптимизация

40573 просмотров, 21 ответов — стр. 1 из 2

Егм
Сообщений: 912
3 декабря 2021 г. в 06:33#1
Попалась интересная статья в ЦОС 2 01 Дворкович .
Но не фига не понял как считаются коэффициенты оптимизированной оконной функции .
Кто силён в математике , объясните на пальцах как считать .

Игорь 2
Сообщений: 21288
3 декабря 2021 г. в 11:49#2
ra0ahcra0ahc
Сообщений: 4868
4 декабря 2021 г. в 08:48#3
Не совсем понял вопрос
Здесь я где-то тему создал про оконные функции.
Вам пример на С показать ?
Там реально одна функция для нас подходит.
А так, перед скоростной обработкой массивов с замерами ацп, создаётся матрица с расчетами этих функций (чтобы проц синусами не грузить) и дальше просто подставляются индекс массива и оттуда берётся предрассчитанный коэфф и умножается на замер ацп… и тд
Егм
Сообщений: 912
4 декабря 2021 г. в 10:56#4
У ОФ(оконной функции) есть два взаимосвязанных параметра ширина главного лепестка(Fш) и уровень подавления
боковых лепестков (Wmax).

Прямоугольное окно Fш=1 бин Wmax=13 дБ , динамически диапазон всего 13 дБ .
АЦП =3 разряда (DR = 20lg(2^N), где N — разрядность квантования).

ОФ Хэмминга Fш=1.36 бин Wmax=43 дБ , динамически диапазон 43 дБ .
АЦП =7 разрядов . Чтобы сохранить частотное разрешение массив замеров нужно увеличить минимум в 1.36 раза .

У ЗК разрядность 24 ДД=144 дБ . Простые окна не подходят , только окно Кайзера-Бесселя может обеспечить
Wmax>80 дБ но и Fш>2 , массив замеров ещё больше.

Для полного использования АЦП нужно согласовать ДД(ЗК) и Wmax(ОФ) и при этом получить минимальный Fш .
Для реального времени обработки размер массива замеров имеет существенное значение .
Оптимизированное окно в принципе может сгладить противоречие Fш и Wmax .
Но как их рассчитать пока вопрос .
Нужна помощь разобраться как Дворкович ОФ считал .
ra0ahcra0ahc
Сообщений: 4868
4 декабря 2021 г. в 03:18#5
Не, не хватит 24 бит. Там смеситель свернётся от таких напруг. Нужен атт перед смесителем. Но при этом если поставить атт то не хватит чуйки. Вообщем 24 бита это круто, у меня он сразу после смесителя стоит и вытягивает чуйку достойную. Во всяком случаи, шум ограничен входным смесителем и первыми каскадами. Если не надо было на нч переходить тогда можно было говорить об одном ацп после кф и все, но таких по моему нет . Там вроде 18 бит до 24 Мгц есть ацп и все.
Оконные функции я просто использовал не задумываясь. Панораму через оф пропускал…. Что ещё … не помню уже , больше года прошло. Как перешёл на адау так и забыл про все это. Адау - легко и красиво все фильтрует и это - очень круто, так как в процессоре не надо заморачиваться фильтрацией и занимаешься только общей логикой аппарата.
Егм
Сообщений: 912
13 декабря 2021 г. в 10:44#6
Разобрался как считают ОФ.
Формула №5 . Для расчета на лету очень сложно получается .
6Ж2П6Ж2П
Сообщений: 1065
13 декабря 2021 г. в 11:00#7
Поясните, какая связь формулы №5 с расчётами на лету?
Егм
Сообщений: 912
13 декабря 2021 г. в 11:12#8
Коэффициенты рассчитываются из условия минимальной разницы бесконечного интеграла и интеграла с пределами
-С...С(ширина главного лепестка) . Подбор коэффициентов а по минимуму разности получился довольно сложным .
Простого способа подбора я не знаю.
GenaSPB
Сообщений: 74
14 декабря 2021 г. в 11:45#9
Цитата: Егм
Для расчета на лету очень сложно получается

Расчёт на лету он же при смене фильтра нужен, при смене режима. Не на каждый же блок обрабатываемых семплов...
6Ж2П6Ж2П
Сообщений: 1065
14 декабря 2021 г. в 11:50#10
Наверное все расчёты для применяемых режимов вообще сделаны заранее и сведены в таблицы коэффициентов. Потому не понимаю, о каких расчётах на лету и Вы пишете?
Егм
Сообщений: 912
15 декабря 2021 г. в 10:11#11
Цитата: GenaSPB
Расчёт на лету он же при смене фильтра нужен, при смене режима. Не на каждый же блок обрабатываемых семплов...

Верно . Но переход с режима на режим за несколько минут - как то долговато .

Цитата: 6Ж2П
Наверное все расчёты для применяемых режимов вообще сделаны заранее и сведены в таблицы коэффициентов.

Таскать за программой десятки мегабайт таблиц мне не нравится . На все случаи таблиц не напасёшся.
ra0ahcra0ahc
Сообщений: 4868
15 декабря 2021 г. в 11:24#12
Какая ваша цель? Что вы делаете? Если вам нужен фильтр - так скажите. Если просто Лаба - супер , надо чём-то заниматься.
Если что у меня есть модуль на С с окном и переменной полосой пропускания. Все коэффициенты рассчитываются на лету при смены полосы. Надо?
Егм
Сообщений: 912
15 декабря 2021 г. в 11:44#13
Спектральный анализ.
Нужно подгонять уровни подавления боковых лепестков под шумы ЗК .
И не нужно забывать про ширину главного лепестка .
В дальнейших расчётах это важные величины .

Цитата: ra0ahc
Если что у меня есть модуль на С с окном и переменной полосой пропускания. Все коэффициенты рассчитываются на лету при смены полосы. Надо?

Пока нет . Ещё куча общих вопросов в которых пока не разобрался . Обычный дефицит времени ...
GenaSPB
Сообщений: 74
16 декабря 2021 г. в 08:45#14
Какие минуты? На 1024 бина расчёт полинома из пяти членов... какие минуты тут могут быть?
Егм
Сообщений: 912
16 декабря 2021 г. в 10:17#15
Цитата: GenaSPB
Какие минуты? На 1024 бина расчёт полинома из пяти членов...

Если не проблемно киньте математику .

Порой не хватает элементарных знаний.