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

 
 
> есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling), хочется что-то аналогичное дизерингу при генерации синуса
Krys
сообщение Apr 16 2015, 04:10
Сообщение #1


Гуру
******

Группа: Свой
Сообщений: 2 002
Регистрация: 17-01-06
Из: Томск, Россия
Пользователь №: 13 271



Здравствуйте. Подскажите, пожалуйста, есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling, digital up-converter), хочется что-то аналогичное дизерингу при генерации синуса в DDS, когда периодичность ошибки фазы вызывает группировку шума на каких-то частотах спектра, а дизеринг этот шум размазывает, улучшая SFDR. Мне надо что-то такое же.
Поясню на картинках.
Прикрепленное изображение

На этой картинке синим спектр исходного сигнала (ЛЧМ), до передискретизации. Красным - спектр после повышения тактовой и линейной интерполяции. Зелёным - спектр после повышения тактовой и фильтрации обычным КИХ-фильтром.

Прикрепленное изображение

Эта картинка из временного домена, после повышения тактовой (увеличено вблизи конца временной выборки, когда частота ЛЧМ близка к максимальной). Синим - результат линейной интерполяции, красным - результат интерполяции через КИХ-фильтр (зелёный спектр на предыдущей картинке), а зелёным - разница между синей и красной кривыми.

По первой картинке (красная кривая) видно, что присутствуют зеркальные спектры от ЛЧМ. Хочется их маленько ослабить, размазав, т.к. есть запас по частоте в половину диапазона (полоса ЛЧМ - четверть частоты дискретизации). Т.е. аналогично размазывают дизерингом составляющие от DDS. Вот я подумал, нельзя ли что-то подобное сделать для линейной интерполяции? Куда тут только шум влепить... в узловые точки интерполяции?


--------------------
Зная себе цену, нужно ещё и пользоваться спросом...
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
ViKo
сообщение Apr 17 2015, 04:33
Сообщение #2


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Зачем обижаться? И почему вы не принимаете во внимание очевидное (что линейная интерполяция не подходит)?
Самая первая картинка в теме, зеленый график, полоса сигнала зарезана - разве это интерполяция?
Go to the top of the page
 
+Quote Post
Krys
сообщение Apr 17 2015, 05:18
Сообщение #3


Гуру
******

Группа: Свой
Сообщений: 2 002
Регистрация: 17-01-06
Из: Томск, Россия
Пользователь №: 13 271



Цитата(ViKo @ Apr 17 2015, 11:33) *
Зачем обижаться?
за сказанное другими словами "спасибо, Кэп" ))))


Цитата(ViKo @ Apr 17 2015, 11:33) *
И почему вы не принимаете во внимание очевидное (что линейная интерполяция не подходит)?
Потому что все в этой теме пытаются подменить вопрос. Мне не нужен ответ на вопрос, какой тип интерполяции здесь подойдёт. Мне нужен ответ, как размазать спектр. Потому что этот вопрос может быть актуален и для других типов интерполяции. Пусть используем более качественную интерполяцию, получим лучше результат. Но пусть он нас по SFDR не устраивает. Зато мы видим, что спектр ошибки неравномерен, откуда возникает желание его размазать, как это сделано в DDS. Неважно какая интерполяция, т.к. нет предела совершенству.


Цитата(ViKo @ Apr 17 2015, 11:33) *
Самая первая картинка в теме, зеленый график, полоса сигнала зарезана - разве это интерполяция?
Ну там подписано. Это классическая интерполяция через КИХ-фильтр. Полоса зарезана - что Вы имеете в виду? Я использовал ЛЧМ с полосой в половину возможной, я об этом где-то выше упомянул.

Цитата(Hose @ Apr 17 2015, 11:57) *
Далее из картинки очевидно, что при первичном повышении частоты дискретизации было копирование отчетов вместо добивания нулями, что уже дало искажения спектра в основной полосе.
Копирования не было, была линейная интерполяция. Это показывает картинка во временной области. Там и подписи есть, что где.



Цитата(Hose @ Apr 17 2015, 11:57) *
Красная картинка это правильное повышение частоты + сик фильтр с полюсом в центре зеркала. Зато эффект от цика виден...
Это обычный КИХ через fdatool. Он больше нужен, чтобы его считать эталонным, а потом вычесть его из линейно-интерполированного, чтобы увидеть ошибку.


Цитата(Hose @ Apr 17 2015, 11:57) *
Предлагаю поставить несколько циков с разными по частоте полюсами в зеркалах.
Спасибо, идея понятна, но по ресурсам фильтры не годятся. Это во много раз сложнее, чем реализованная линейная интерполяция.


--------------------
Зная себе цену, нужно ещё и пользоваться спросом...
Go to the top of the page
 
+Quote Post
ViKo
сообщение Apr 17 2015, 06:52
Сообщение #4


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(Krys @ Apr 17 2015, 08:18) *
за сказанное другими словами "спасибо, Кэп" ))))

Я не спрашивал, за что. rolleyes.gif
Цитата
Ну там подписано. Это классическая через КИХ-фильтр. Полоса зарезана - что Вы имеете в виду? Я использовал ЛЧМ с полосой в половину возможной, я об этом где-то выше упомянул.

А, у вас там частота относительная, и полоса синего сигнала равна полосе зеленого. Это же очевидно!
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Krys   есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling)   Apr 16 2015, 04:10
- - ViKo   Может, лучше сделать кубическую интерполяцию? Вмес...   Apr 16 2015, 04:31
- - Krys   Конечно лучше )) Тот же Фарроу )) Но это надо умно...   Apr 16 2015, 05:15
|- - rloc   Есть CIC фильтры только со сложением, можно и FIR ...   Apr 16 2015, 05:33
|- - rloc   ..   Apr 16 2015, 05:33
- - Krys   CIC-фильтры из другой темы - это совсем другая тем...   Apr 16 2015, 05:53
|- - rloc   В DDS другая история, там уровень составляющих ниж...   Apr 16 2015, 06:20
|- - Maverick   кусочно-линейная интерполяция тоже не подходит, т....   Apr 16 2015, 06:31
- - Krys   Цитата(rloc @ Apr 16 2015, 13:20) В DDS д...   Apr 16 2015, 06:42
|- - Maverick   Цитата(Krys @ Apr 16 2015, 09:42) Не поня...   Apr 16 2015, 06:59
- - serjj   Добавление на вход интерполятора шума ничего не да...   Apr 16 2015, 07:37
|- - Krys   Цитата(serjj @ Apr 16 2015, 14:37) Добавл...   Apr 16 2015, 07:59
|- - petrov   Цитата(Krys @ Apr 16 2015, 10:59) Чаще ко...   Apr 16 2015, 08:11
- - Krys   Maverick, Понял Вас. Но ресурсов тоже сожрёт кучу....   Apr 16 2015, 07:39
- - ViKo   Если вопрос только в том, куда лепить шум - лепите...   Apr 16 2015, 07:40
- - serjj   Что вы понимаете под сигналом ошибки в данном случ...   Apr 16 2015, 07:42
- - petrov   Может лучше ЛЧМ сразу на высокой частоте сгенериро...   Apr 16 2015, 07:46
- - ViKo   Покажите картинку с добавленным шумом.   Apr 16 2015, 08:15
|- - Fat Robot   "у меня нет ресурсов на то, чтобы реализовать...   Apr 16 2015, 08:31
|- - Krys   Цитата(ViKo @ Apr 16 2015, 15:15) Покажит...   Apr 16 2015, 08:58
- - serjj   ЦитатаМне чудеса и не нужны. Мне не опустить лепес...   Apr 16 2015, 08:31
- - serjj   Цитата"у меня нет ресурсов на то, чтобы реали...   Apr 16 2015, 08:37
|- - Fat Robot   Положим Есть также быстрые схемы деления, наприме...   Apr 16 2015, 08:49
- - blackfin   Цитата(Fat Robot @ Apr 16 2015, 11:49) По...   Apr 16 2015, 08:55
|- - Fat Robot   10<{16,32,64}<100 а имея делители на и на ...   Apr 16 2015, 09:00
- - serjj   Цитата10<{16,32,64}<100 а имея делители на ...   Apr 16 2015, 09:08
- - Krys   Цитата(serjj @ Apr 16 2015, 15:31) Почему...   Apr 16 2015, 09:13
|- - Fat Robot   Процитирую одного интернет-пользователя с утонченн...   Apr 16 2015, 09:24
|- - blackfin   Цитата(ViKo @ Apr 16 2015, 13:28) Подскаж...   Apr 16 2015, 10:33
|- - ViKo   Цитата(blackfin @ Apr 16 2015, 13:33) В у...   Apr 16 2015, 10:42
|- - blackfin   Цитата(ViKo @ Apr 16 2015, 13:42) А на 10...   Apr 16 2015, 10:45
- - serjj   ЦитатаНе понял, откуда неравномерность? У меня по ...   Apr 16 2015, 09:21
- - Krys   Цитата(serjj @ Apr 16 2015, 16:21) Неравн...   Apr 16 2015, 09:30
|- - Fat Robot   Когда у нас будет нужна консультация о том, что мн...   Apr 16 2015, 09:37
|- - petrov   Цитата(Krys @ Apr 16 2015, 12:30) Так что...   Apr 16 2015, 20:04
- - ViKo   Krys, вы добавили шум ко входному сигналу. А что б...   Apr 16 2015, 10:28
|- - Fat Robot   Здесь, например. Цитата(ViKo @ Apr 16 2015, ...   Apr 16 2015, 10:50
- - Krys   petrov, Фарроу не годится в любом случае по ресурс...   Apr 17 2015, 02:54
|- - petrov   Цитата(Krys @ Apr 17 2015, 05:54) petrov,...   Apr 17 2015, 09:57
- - Hose   Первое что видно из картинок: первичный спектр име...   Apr 17 2015, 04:57
- - serjj   ЦитатаМне нужен ответ, как размазать спектр. Я нав...   Apr 17 2015, 06:55
|- - Fat Robot   А как же amplitude dithering? Он, правда, по очев...   Apr 17 2015, 07:09
|- - Krys   Цитата(serjj @ Apr 17 2015, 13:55) В LE/п...   Apr 17 2015, 07:48
|- - ViKo   Цитата(Krys @ Apr 17 2015, 10:48) да заче...   Apr 17 2015, 08:18
- - ViKo   Добавили промежуточных выборок, добавьте к ним шум...   Apr 17 2015, 07:25
- - serjj   ЦитатаНу я где-то выше уже расписывал. Типа на лин...   Apr 17 2015, 07:51
|- - Krys   Цитата(serjj @ Apr 17 2015, 14:51) А если...   Apr 17 2015, 08:35
|- - ViKo   Цитата(Krys @ Apr 17 2015, 11:35) Спасибо...   Apr 17 2015, 08:47
- - Krys   Ну раз уж фильтр ставить, то зачем же ещё шумов до...   Apr 17 2015, 09:19
|- - ViKo   Цитата(Krys @ Apr 17 2015, 12:19) Ну раз ...   Apr 17 2015, 09:28
|- - Krys   Цитата(ViKo @ Apr 17 2015, 16:28) Мне пре...   Apr 17 2015, 09:52
- - Krys   Дайте пожалуйста ссылки на описанный Вами простой ...   Apr 17 2015, 10:00
|- - petrov   Цитата(Krys @ Apr 17 2015, 13:00) Дайте п...   Apr 17 2015, 11:21
- - blackfin   Цитата(ViKo @ Apr 16 2015, 13:42) А на 10...   Apr 17 2015, 13:30
- - Hose   В ддс спурс это одна частотная составляющая. Она з...   Apr 17 2015, 14:07
- - Hose   Если это курсовой или дирломная, то качающийся пол...   Apr 17 2015, 15:20
|- - Maverick   Цитата(Hose @ Apr 17 2015, 18:20) Если эт...   Apr 17 2015, 17:11
- - Krys   Всем спасибо за подсказки. ЗЫ: это не курсовая и ...   Apr 18 2015, 04:22
|- - shf_05   простая последовательность рассуждений приведет ва...   Apr 20 2015, 05:31
|- - Krys   Цитата(shf_05 @ Apr 20 2015, 12:31) прост...   Apr 20 2015, 08:10
|- - blackfin   Цитата(Krys @ Apr 20 2015, 11:10) Спасибо...   Apr 20 2015, 11:51
|- - Krys   Цитата(blackfin @ Apr 20 2015, 18:51) Из ...   Apr 21 2015, 02:05
- - ViKo   Krys, в вашем скрипте используются функции plot_re...   Apr 20 2015, 09:39
- - serjj   ЦитатаТ.е. линейная интерполяция ничего в полосе н...   Apr 21 2015, 07:00
|- - Krys   Цитата(serjj @ Apr 21 2015, 14:00) Смотря...   Apr 21 2015, 07:32
- - serjj   ЦитатаЯ имел в виду под портить наличие интермодул...   Apr 21 2015, 08:04
|- - Krys   Цитата(serjj @ Apr 21 2015, 15:04) Тоже п...   Apr 21 2015, 08:36
|- - petrov   Цитата(Krys @ Apr 21 2015, 11:36) Вот это...   Apr 21 2015, 09:56
- - Krys   Ой, нецелое число раз - это уже совсем высший пило...   Apr 21 2015, 11:00
- - mihalevski   Если Вы реализовали частотный детектор, используя ...   May 9 2015, 14:00
- - Krys   да нет, у меня не частотный детектор. У меня прост...   May 18 2015, 06:18


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

 


RSS Текстовая версия Сейчас: 27th August 2025 - 00:06
Рейтинг@Mail.ru


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