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

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

|
Здравствуйте. Подскажите, пожалуйста, есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling, digital up-converter), хочется что-то аналогичное дизерингу при генерации синуса в DDS, когда периодичность ошибки фазы вызывает группировку шума на каких-то частотах спектра, а дизеринг этот шум размазывает, улучшая SFDR. Мне надо что-то такое же. Поясню на картинках.
На этой картинке синим спектр исходного сигнала (ЛЧМ), до передискретизации. Красным - спектр после повышения тактовой и линейной интерполяции. Зелёным - спектр после повышения тактовой и фильтрации обычным КИХ-фильтром.
Эта картинка из временного домена, после повышения тактовой (увеличено вблизи конца временной выборки, когда частота ЛЧМ близка к максимальной). Синим - результат линейной интерполяции, красным - результат интерполяции через КИХ-фильтр (зелёный спектр на предыдущей картинке), а зелёным - разница между синей и красной кривыми. По первой картинке (красная кривая) видно, что присутствуют зеркальные спектры от ЛЧМ. Хочется их маленько ослабить, размазав, т.к. есть запас по частоте в половину диапазона (полоса ЛЧМ - четверть частоты дискретизации). Т.е. аналогично размазывают дизерингом составляющие от DDS. Вот я подумал, нельзя ли что-то подобное сделать для линейной интерполяции? Куда тут только шум влепить... в узловые точки интерполяции?
--------------------
Зная себе цену, нужно ещё и пользоваться спросом...
|
|
|
|
|
 |
Ответов
|
Apr 17 2015, 05:18
|

Гуру
     
Группа: Свой
Сообщений: 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)  Предлагаю поставить несколько циков с разными по частоте полюсами в зеркалах. Спасибо, идея понятна, но по ресурсам фильтры не годятся. Это во много раз сложнее, чем реализованная линейная интерполяция.
--------------------
Зная себе цену, нужно ещё и пользоваться спросом...
|
|
|
|
Сообщений в этой теме
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
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|