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

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

40563 просмотров, 21 ответов

Егм 3 декабря 2021 г. в 06:33#1
Попалась интересная статья в ЦОС 2 01 Дворкович .
Но не фига не понял как считаются коэффициенты оптимизированной оконной функции .
Кто силён в математике , объясните на пальцах как считать .
Игорь 2 3 декабря 2021 г. в 11:49#2
ra0ahc 4 декабря 2021 г. в 03:18#3
Не, не хватит 24 бит. Там смеситель свернётся от таких напруг. Нужен атт перед смесителем. Но при этом если поставить атт то не хватит чуйки. Вообщем 24 бита это круто, у меня он сразу после смесителя стоит и вытягивает чуйку достойную. Во всяком случаи, шум ограничен входным смесителем и первыми каскадами. Если не надо было на нч переходить тогда можно было говорить об одном ацп после кф и все, но таких по моему нет . Там вроде 18 бит до 24 Мгц есть ацп и все.
Оконные функции я просто использовал не задумываясь. Панораму через оф пропускал…. Что ещё … не помню уже , больше года прошло. Как перешёл на адау так и забыл про все это. Адау - легко и красиво все фильтрует и это - очень круто, так как в процессоре не надо заморачиваться фильтрацией и занимаешься только общей логикой аппарата.
ra0ahc 4 декабря 2021 г. в 08:48#4
Не совсем понял вопрос
Здесь я где-то тему создал про оконные функции.
Вам пример на С показать ?
Там реально одна функция для нас подходит.
А так, перед скоростной обработкой массивов с замерами ацп, создаётся матрица с расчетами этих функций (чтобы проц синусами не грузить) и дальше просто подставляются индекс массива и оттуда берётся предрассчитанный коэфф и умножается на замер ацп… и тд
Егм 4 декабря 2021 г. в 10:56#5
У ОФ(оконной функции) есть два взаимосвязанных параметра ширина главного лепестка(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 .
Но как их рассчитать пока вопрос .
Нужна помощь разобраться как Дворкович ОФ считал .
Егм 13 декабря 2021 г. в 10:44#6
Разобрался как считают ОФ.
Формула №5 . Для расчета на лету очень сложно получается .
6Ж2П 13 декабря 2021 г. в 11:00#7
Поясните, какая связь формулы №5 с расчётами на лету?
Егм 13 декабря 2021 г. в 11:12#8
Коэффициенты рассчитываются из условия минимальной разницы бесконечного интеграла и интеграла с пределами
-С...С(ширина главного лепестка) . Подбор коэффициентов а по минимуму разности получился довольно сложным .
Простого способа подбора я не знаю.
GenaSPB 14 декабря 2021 г. в 11:45#9
Для расчета на лету очень сложно получается
Расчёт на лету он же при смене фильтра нужен, при смене режима. Не на каждый же блок обрабатываемых семплов...
6Ж2П 14 декабря 2021 г. в 11:50#10
Наверное все расчёты для применяемых режимов вообще сделаны заранее и сведены в таблицы коэффициентов. Потому не понимаю, о каких расчётах на лету и Вы пишете?
Егм 15 декабря 2021 г. в 10:11#11
Расчёт на лету он же при смене фильтра нужен, при смене режима. Не на каждый же блок обрабатываемых семплов...
Верно . Но переход с режима на режим за несколько минут - как то долговато .

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

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

Порой не хватает элементарных знаний.
GenaSPB 17 декабря 2021 г. в 01:06#16
Егм 17 декабря 2021 г. в 09:22#17
Спасибо !
Попробую разобраться .
GenaSPB 17 декабря 2021 г. в 11:19#18
Разбираться проще посмотрев в этом файле откуда вызывается эта функция
Егм 17 декабря 2021 г. в 11:56#19
Посмотрел листинг , пока не разобрался.
Нужно время чтобы переварить информацию .
GenaSPB 21 декабря 2021 г. в 06:19#20
Какой лисинг Вы смотрели? Есть ли вопросы?
Егм 22 декабря 2021 г. в 01:33#21
Какой лисинг Вы смотрели? Есть ли вопросы?
Ваша ссылка Ответ #15 .
Для вопросов нужно разобраться что зачем и почему , пока анализом занимаюсь урывками ...
Как вопросы появятся задам здесь или в личку постучу .