Страница 3 из 4 Первая 1234 Последняя
Показано с 41 по 60 из 62

Тема: О погрешности конвольвера

  1. #1 Показать/скрыть первое сообщение.
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию О погрешности конвольвера

    Заинтересовала меня недавняя дискуссия в теме «О фазолинейности. Продолжение или революция?!» о погрешности, вносимой конвольвером. После слов ViktKors
    Цитата Сообщение от ViktKors Посмотреть сообщение
    Увы, если я слышу грязь, я интерпертирую это как грязь.
    Причем, что характрено, если в наушниках я слушаю после конвольвера _искривленную ФЧХ_ , и именно на искривление можно попытаться списать эффект, то в АС ФЧХ ведь выравнивается, а эффект точно так-же плохой.
    я решил проверить, какая же «грязь» может возникать после применения фазового FIR-фильтра.

    Пусть у нас есть дискретные отсчёты сигнала с частотой выборки 48 кГц. Построим all-pass фильтр с коротким окном на 16 отсчётов. Постараемся сделать так, чтобы коэффициент передачи по амплитуде был равен 1, а фаза менялась по закону косинуса от частоты. При 0 Гц фаза равна 0°, при 12 кГц – минус 90°, а при 24 кГц достигает минус 180°. Для вычисления 16 искомых коэффициентов конвольвера нам нужно решить систему из 16 линейных. Я зафиксировал 7 точек на шкале частот 3,6,9,12,15,18,21 кГц. Для них по условиям задана фаза и амплитуда, что даёт 14 уравнений. Ещё два уравнения даёт условие для коэффиициента передачи +1 на частоте 0 Гц и минус 1 на частоте 24 кГц. Эти частоты стоят особняком, так как на них фаза может быть либо 0° либо ±180°.

    Решаем систему и получаем 16 коэффициентов конвольвера:
    Код:
      0.0000000000000002
      0.9722096263235436
     -0.0000000000000001
     -0.1668931209924746
      0.0000000000000001
     -0.0025951830757338
     -0.0000000000000001
     -0.0016313244799198
      0.0000000000000001
      0.0015642123336106
     -0.0000000000000001
      0.0070858939842292
      0.0000000000000000
      0.0288213444185797
     -0.0000000000000002
      0.1614385514881652
    Проверяем, строим график ФЧХ и радуемся:
    Нажмите на изображение для увеличения. 

Название:	PH16.png 
Просмотров:	310 
Размер:	6.5 Кб 
ID:	339288
    Целевая фаза показана зелёным цветом, а точки аппроксимации синие кружочки.
    Но есть один нюанс. Исходный сигнал содержит, конечно, больше частот, чем наши 9 штук прибитых гвоздями с шагом 3 кГц. Эти частоты тоже буду преобразовываться при свёртке с нашим окном из 16 коэффициентов. Их ФЧХ показана красной кривой. Между узлами интерполяции меняется даже знак производной, т.е. часть частот с положительным ГВЗ будут отставать, а другие с отрицательным ГВЗ – наоборот опережать соседей.

    На графике АЧХ картина столь же ужасна.
    Нажмите на изображение для увеличения. 

Название:	DB16.png 
Просмотров:	249 
Размер:	8.3 Кб 
ID:	339289
    Хотя точки с шагом 3 кГц точно лежат на целевой зелёной полке 0 дБ, но зато между ними мы огребаем погрешность от -4 до +2 дБ.

    Первая мысль была: надо взять больше точек и всё улучшится, но расчёты говорят об обратном.
    Увеличил в 16 раз длину конвольвера до 256 точек.
    Нажмите на изображение для увеличения. 

Название:	PH256.png 
Просмотров:	225 
Размер:	9.4 Кб 
ID:	339290
    Нажмите на изображение для увеличения. 

Название:	DB256.png 
Просмотров:	241 
Размер:	9.4 Кб 
ID:	339291
    Смотрим и видим, что хотя в самих узлах аппроксимации АФЧХ всё волшебно, но между ними творится кошмар. Одни частоты усиливаются, другие ослабляются, а фаза у всех разъехалась в разные стороны от целевой кривой. Причём амплитуда погрешности и не думает уменьшаться с ростом окна.

    Есть у меня подозрения, что именно "лохматость" АФЧХ конвольвера воспринимается на слух как "грязь".
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  2. #41
    самый главный Аватар для Игорь Гапонов
    Регистрация
    03.03.2010
    Адрес
    Одесса
    Возраст
    12
    Сообщений
    3,156

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    /// может эта симметрия и вредит восприятию../
    Ниже обоснованное ИМХО. Но готов измениться, побриться, бросить курить и посыпать голову пеплом из всех моих грязных пепельниц.

    Чисто по "мысленному эксперименту" (по заветам Эйнштейна) при безошибочной коррекции или при реально неслышимой ошибке (хоть каузальной, хоть некаузальной) получается, что объективный стимул ощущений может не соответствовать неискажённому только в начале и конце всей фонограммы как массива данных, т.к. реально за "пределами" этого массива данные не определены, а не равны каким-то конкретным значениям (обычно принимается, что все "остальные семплы"=0). Т.е. на интервалах времени самой фонограммы всё должно быть чики пики. Собственно, этот феномен связан с усечением данных (числа семплов) больше/глубже, чем с округлением.

    Однако....

    1. Реальные фонограммы создаются не в идеальных условиях и с ошибками, не учитывающими или учитывающими не в должной мере реальное усечение данных как "подфайлов" составляющих результирующую фонограмму/продукт при микшировании. Т.е. сложение и обработка происходят с какой-то ошибкой, которая может вылезти на системе с линейной коррекцией, приближенной к идеальной.

    2. Согласно п.1 к неидеальным условиям создания относится и студийный участок "АС-ухИ звукорежа" при контроле и редактировании "файнал кат". Т.е. зависят от решений звукорежа или ещё кого не зависимо объективные или субъективные или ещё какие причины лежат в их основе.

    Поэтому несовпадений мнений/ощущений слушателя с мнениями/ощущениями звукорежа могут быть источником ошибок при слуховой оценке систем вторичного поля с линейной коррекцией, приближенной к идеальной. Вероятность этих ошибок оценки может быть существенно снижена при симуляции условий студии. Однако эти симуляции должны быть выполнены с достаточной точностью на интервале времени, соответствующему интервалу фонограммы. В другие моменты времени оценка должна игнорироваться.
    Электричество дисциплинирует

  3. #42
    Завсегдатай Аватар для E.Sokol
    Регистрация
    17.01.2015
    Адрес
    Урал
    Сообщений
    1,290

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    реальное измерение в комнате
    Кстати. Если немного трансформировать импульс, то легко можно увидеть комнатные отражения:
    Нажмите на изображение для увеличения. 

Название:	Комната из точки прослушивания - peaks (3000 samples, 96kHz).png 
Просмотров:	164 
Размер:	28.9 Кб 
ID:	339553
    Вот главный враг хорошего звука.

  4. #43
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Если немного трансформировать импульс, то легко можно увидеть комнатные отражения:
    Можно и для правого канала посмотреть . По ним можно ориентироваться на уровень амплитуды относительно основного тона?
    Я могу накидать внеосевых характеристик АС без комнаты для каждого отражения, если получится разобраться где кто, но думается мне что самый большой пик - это отражение от дивана, и где-то в этой же области сидят отражения от пола/потолка. Останется потом только грамотно расставить относительно единичного пика, например по совпадению общей ИХ, поставив вместо единицы характеристику АС.
    Миниатюры Миниатюры Нажмите на изображение для увеличения. 

Название:	Right.jpg 
Просмотров:	115 
Размер:	228.1 Кб 
ID:	339567   Нажмите на изображение для увеличения. 

Название:	Right Rectify.jpg 
Просмотров:	112 
Размер:	221.7 Кб 
ID:	339568  

  5. #44
    Завсегдатай Аватар для E.Sokol
    Регистрация
    17.01.2015
    Адрес
    Урал
    Сообщений
    1,290

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    Можно и для правого канала посмотреть
    Это вопрос или ответ? Вообще, разделить прямой сигнал и отражённый - задача нетривиальная и я не нашёл пока для неё удовлетворительного решения.

  6. #45
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Это вопрос или ответ?
    Ответ, но импульс могу приложить .
    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Вообще, разделить прямой сигнал и отражённый - задача нетривиальная и я не нашёл пока для неё удовлетворительного решения.
    Я полагаю, что единственный вариант - симулятор, в котором можно двигать импульсы друг относительно друга и менять их амплитуду, набор внеосевых характеристик АС, учет примерной геометрии и вперед поехал до тех пор пока картинки измеренного не совпадут с симуляцией. Потом меняешь основной импульс на единичку такой же амплитуды, инвертируешь все остальные и фильтр готов.
    Вложения Вложения

  7. #46
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Лабораторная работа №2.

    В качестве подопытного кролика возьмём воображаемую АС, в которой два идеальных широкополосных динамика сшиты на частоте 1,5 кГц при помощи фильтра LR2. Считаем, что в воздухе их излучение точно складываются, в результате чего получается фазовый all pass фильтр.
    Нажмите на изображение для увеличения. 

Название:	LR2.png 
Просмотров:	131 
Размер:	3.5 Кб 
ID:	339949

    Для ровной АЧХ полярность включения НЧ динамика должна быть противоположной полярности ВЧ полосы, что стандартно для двухполосок с фильтрами 2-го порядка. В результате фаза на НЧ и ВЧ отличается на 180 °, а Step response получается знакопеременным:
    Нажмите на изображение для увеличения. 

Название:	StepLR2.png 
Просмотров:	136 
Размер:	2.4 Кб 
ID:	339950

    Несколько модифицировав метод из первого поста, я синтезировал цифровой фильтр, компенсирующий фазовые искажения LR2 и сравнил его точность при разной длине конвольвера N.
    Нажмите на изображение для увеличения. 

Название:	Ampl.png 
Просмотров:	129 
Размер:	12.4 Кб 
ID:	339951Нажмите на изображение для увеличения. 

Название:	Phase.png 
Просмотров:	136 
Размер:	12.2 Кб 
ID:	339952

    Зададим N=81, что обеспечивает ошибку менее 0.03 ° и 0.005 dB.

    кэффициенты конвольвера для любопытных

    -3.394067535792349e-007
    -1.073035276353023e-006
    3.171145920027496e-007
    -1.039798762437149e-006
    4.154506651811293e-007
    -1.051967121840412e-006
    4.695300470071582e-007
    -1.083746586623230e-006
    5.164196740916020e-007
    -1.131102445137578e-006
    5.674010526800140e-007
    -1.195298118183046e-006
    6.297735171748399e-007
    -1.280615784801981e-006
    7.113029231869829e-007
    -1.394591527784939e-006
    8.225663430505983e-007
    -1.549525824524923e-006
    9.797048233622516e-007
    -1.765632497385991e-006
    1.209058990439488e-006
    -2.077251138648720e-006
    1.556104455466599e-006
    -2.545585813800315e-006
    2.104181316946129e-006
    -3.286633899912586e-006
    3.017187605609536e-006
    -4.538005656585458e-006
    4.647166137858660e-006
    -6.837716932521536e-006
    7.842500215335270e-006
    -1.157858664053566e-005
    1.499672373085892e-005
    -2.311673627838549e-005
    3.461745344518448e-005
    -5.966467079234168e-005
    1.108191374659712e-004
    -2.489745155857576e-004
    7.371729915147321e-004
    -4.010096770829227e-003
    8.230326994255250e-001
    -3.198694093481699e-001
    -2.681289324755042e-001
    -2.193144365089487e-001
    -1.805726730114869e-001
    -1.482166774982609e-001
    -1.218824204511451e-001
    -1.001008394638855e-001
    -8.228968401990872e-002
    -6.759662173872803e-002
    -5.556225294570145e-002
    -4.564539388703905e-002
    -3.751687690704045e-002
    -3.082228218599461e-002
    -2.533268941908416e-002
    -2.081297072538952e-002
    -1.710580646911148e-002
    -1.405431940007817e-002
    -1.155091085834554e-002
    -9.490713784821548e-003
    -7.800201742695296e-003
    -6.409327997704577e-003
    -5.267751250107707e-003
    -4.328831047002429e-003
    -3.557948698402691e-003
    -2.924221555286931e-003
    -2.403694414092892e-003
    -1.976071934949767e-003
    -1.624679035567794e-003
    -1.336263894377172e-003
    -1.099233919074448e-003
    -9.048739956342524e-004
    -7.453757625790847e-004
    -6.146261659716933e-004
    -5.081656919406395e-004
    -4.205833601467008e-004
    -3.518071106375566e-004
    -2.939682952048770e-004
    -2.583791209939207e-004
    -2.231776728119541e-004
    -3.223662747824697e-004
    [свернуть]
    Графически импульсная характеристика фильтра выглядит так:
    Нажмите на изображение для увеличения. 

Название:	Impls.png 
Просмотров:	128 
Размер:	2.9 Кб 
ID:	339953
    Коэффициенты слева от центральной точки (отсчёты из прошлого) в среднем в 1000 раз меньше, чем коэффициенты справа от неё, т.е. сигнал на выходе корректирующего фильтра на 99.9% зависит от значений входных семплов в ближайшем будущем. Так получилось потому что у исходного LR2 фильтра ситуация противоположная. Он, как всякий IIR фильтр, зависит только от значений входного сигнала в прошлом, но ничего не знает о будущем сигнала.

    Теперь проверим точность свёртки. В качестве входного сигнала зададим точный аналитический Step response LR2: f(t) = 0 для t < t0, f(t) = 2*exp((t0-t)/(RC))-1 для t > t0, где t0 – момент времени единичного скачка на входе. Перемножаем, строим график и видим симпатичную восстановленную ступеньку:
    Нажмите на изображение для увеличения. 

Название:	Res.png 
Просмотров:	214 
Размер:	4.4 Кб 
ID:	339954

    И последний шаг: вычтем из результатов свёртки единичный степ, т.е. выделим отклонение от точного решения. Одновременно мы избавимся от колебаний вызванных явлением Гиббса вокруг ступеньки, которые сильно мешали бы восприятию, начни мы строить интерполяцию сигнала.
    Нажмите на изображение для увеличения. 

Название:	Delta.png 
Просмотров:	149 
Размер:	3.0 Кб 
ID:	339955
    Видно, что максимальная погрешность наблюдается в окрестности скачка, но даже там она составляет 0.4 %. Ошибка в минус 48 дБ далеко не рекордная величина, но для скромной длины конвольвера очень даже неплохо.

    Последнее замечание. В расчётах я опирался на точно известную фазовую характеристику LR2. Понятно, что при измерении АФЧХ реальных АС, точность 0.03 ° и 0.005 dB не достижима, что автоматически ограничивает точность синтеза фильтра коррекции. При исходных данных с ошибками качество восстановленной ступеньки будет естественно хуже, а самое неприятное состоит в том, что мы не получим протяжённых горизонтальных полочек как на последней картинке. ИМХО там неизбежно возникнут колебания, которые часто называют звоном FIR.
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  8. #47
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    _Сам_, а если дополнить срезом снизу и посмотреть результат на НЧ? Например, добавить входной срез 20Гц, а затем сравнить ступеньки без и с конвольвером.
    Относительно недавно пробовал похожий фильтр (корректировка раздела СЧ-ВЧ), не трогая минимальную фазу на НЧ, басы хуже звучат, чем с IIR. Применительно 22 поста корректировка провала фазы с 800Гц до 30кГц.

  9. #48
    самый главный Аватар для Игорь Гапонов
    Регистрация
    03.03.2010
    Адрес
    Одесса
    Возраст
    12
    Сообщений
    3,156

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от _Сам_ Посмотреть сообщение
    явлением Гиббса вокруг ступеньки, которые сильно мешали бы восприятию, начни мы строить интерполяцию сигнала.
    ..... При исходных данных с ошибками качество восстановленной ступеньки будет естественно хуже, а самое неприятное состоит в том, что мы не получим протяжённых горизонтальных полочек как на последней картинке. ИМХО там неизбежно возникнут колебания, которые часто называют звоном FIR.
    Если мы слышим после какой-то обработки какого-то сигнала "звон" и/или "явление Гиббса", то частота дискретизации недостаточна или (что то же самое) дискретное время/частота анализа не адекватны биологии/физиологии/психике.

    Вот вам/нам готовая стратегия хитрых психоакустических опытов, позволяющих вырваться из инерционных цепей КРГ/"спектров"

    С другой стороны, мы должны чётко определить физический интервал времени оценки тест сигнала, игнорируя оценку за пределами этого интервала.

    Интересно, что в дискретном домене времени обычное суммирование растягивается во времени, т.е. происходит совсем "не мгновенно". Поэтому простая "склейка файлов", как той магнитофонной ленты, содержит ошибку суммирования (к Вашему анализу это не относится, ясень пень). И все слова, что "PCM легче обрабатывать, чем SD" становятся отфонарными.
    Последний раз редактировалось Игорь Гапонов; 18.02.2019 в 20:11.
    Электричество дисциплинирует

  10. #49
    Завсегдатай Аватар для E.Sokol
    Регистрация
    17.01.2015
    Адрес
    Урал
    Сообщений
    1,290

    По умолчанию Re: О погрешности конвольвера

    Эффект Гиббса - это явление из математики, а не реального мира. Если сдвигать (задерживать во времени) единичный импульс на дробную часть семпла, то он будет то появляться, то исчезать. Но это не значит, что мы можем услышать разницу, обеспечивая эту задержку изменением физического расстояния до источника звука.

    [щёлкнуть для анимации]
    Нажмите на изображение для увеличения. 

Название:	fractional sample delay.gif 
Просмотров:	152 
Размер:	167.4 Кб 
ID:	339984

  11. #50
    самый главный Аватар для Игорь Гапонов
    Регистрация
    03.03.2010
    Адрес
    Одесса
    Возраст
    12
    Сообщений
    3,156

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Эффект Гиббса - это явление из математики, а не реального мира. Если сдвигать (задерживать во времени) единичный импульс на дробную часть семпла, то он будет то появляться, то исчезать. Но это не значит, что мы можем услышать разницу, обеспечивая эту задержку изменением физического расстояния до источника звука.

    [щёлкнуть для анимации]
    Нажмите на изображение для увеличения. 

Название:	fractional sample delay.gif 
Просмотров:	152 
Размер:	167.4 Кб 
ID:	339984
    Нет, Гиббс - это как раз реальный мир с непрерывным однонаправленным временем (подчёркивает однонаправленность времени в реалиях- при стремлении дискрета времени к нулю меняется только "период" волн, но их энергия остаётся неизменной, что требует реальный мир). А вот твоя анимация нет. - невозможно засинхронизировать моменты отсчёта с событием, которое ещё не произошло (каузальность). Но вот, когда мы вычислим условия слышимости "по феномену Гиббса", тогда таки да. Кстати, первые опыты уже проводились, доказывающие, что анализ в ухе вовсе не спектральный - превышение предела Габора для спектральной границы 20кГц. Точнее - точно нелинейный, если граница устоит .
    Последний раз редактировалось Игорь Гапонов; 18.02.2019 в 21:00.
    Электричество дисциплинирует

  12. #51
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    _Сам_, а если дополнить срезом снизу и посмотреть результат на НЧ? Например, добавить входной срез 20Гц, а затем сравнить ступеньки без и с конвольвером.
    Форма ступеньки что-то скажет нам про НЧ? Степ правильной ступенчатой формы лишь свидетельствует о фазолинейности тракта и близости его АЧХ к const.
    Цитата Сообщение от mr-marlen Посмотреть сообщение
    басы хуже звучат, чем с IIR.
    эх, знать бы чему математическому соответствует "хуже / лучше звучит"
    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Эффект Гиббса - это явление из математики, а не реального мира.
    гиббсовские колебания абсолютно реальны. Достаточно взять осциллограф и посмотреть на меандр с выхода ЦАП.
    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Если сдвигать (задерживать во времени) единичный импульс на дробную часть семпла, то он будет то появляться, то исчезать.
    у ЦАП с корректно реализованным восстанавливающим фильтром вы увидим, как на дробную часть шага выборки сдвигается всплеск вида f(x)=sin(x)/x. Исчезать ничего не должно.
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  13. #52
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от _Сам_ Посмотреть сообщение
    Форма ступеньки что-то скажет нам про НЧ? Степ правильной ступенчатой формы лишь свидетельствует о фазолинейности тракта и близости его АЧХ к const.
    Да вот я тоже хз. Можно конечно списать на кривую реализацию свертки в miniSHARC, но думаю вряд-ли.
    Цитата Сообщение от _Сам_ Посмотреть сообщение
    эх, знать бы чему математическому соответствует "хуже / лучше звучит"
    Фактически я тупо подгружаю в каждую полосу FIR длиной 564тапс ЧД 96кГц и должен как-бы слышать небольшие улучшения по сцене и более ясную подачу в диапазоне 800Гц - 30кГц. По факту там может и есть небольшие улучшения, но звучание басов с таким фильтром сразу же заставляет вернутся на конфиг без этой коррекции.
    Кто-то мне сказал, что это проблема от недостатка разрешения на НЧ из-за короткой длины FIR и нужно использовать мультирэйт фир (децимация, фильтр с существенно большим кол-вом тапас, ресемплер), но по-идее там вообще ничего не должно меняться.
    Прицепил то, что подгружал в FIR.
    Вложения Вложения

  14. #53
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    Фактически я тупо подгружаю в каждую полосу FIR длиной 564тапс
    у меня получается АЧХ = const, а ФЧХ вашего фильтра вот такая
    Нажмите на изображение для увеличения. 

Название:	564.png 
Просмотров:	124 
Размер:	3.0 Кб 
ID:	340052
    Наверное это больше похоже на фазу самой АС, а у фильтра фаза должна быть с обратным знаком?
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  15. #54
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    _Сам_, у меня цифровые IIR с определенными задержками между полосами. Картинку АЧХ, ФЧХ моей АС с IIR фильтрами и без FIR приложил.
    ВЧ часть появляется немного раньше НЧ-СЧ части. Если делать идеально по фазе, то внеосевые будут хуже (ниже по частоте не дают КНИ пищалки).
    Цитата Сообщение от _Сам_ Посмотреть сообщение
    Наверное это больше похоже на фазу самой АС, а у фильтра фаза должна быть с обратным знаком?
    Да, на твоей картинке фаза АС. Прикрепил как выглядит сравнение АС и фильтра.
    Но в целом, как мне кажется, оно не должно влиять на НЧ.
    Миниатюры Миниатюры Нажмите на изображение для увеличения. 

Название:	phase.jpg 
Просмотров:	121 
Размер:	47.6 Кб 
ID:	340056   Нажмите на изображение для увеличения. 

Название:	AC_Filter.jpg 
Просмотров:	122 
Размер:	152.9 Кб 
ID:	340062  
    Последний раз редактировалось mr-marlen; 19.02.2019 в 13:01.

  16. #55
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    В моём случае наоборот.
    ок, поменял порядок коэффициентов на обратный, чтобы получить зеркальную фчх
    Нажмите на изображение для увеличения. 

Название:	фчх.png 
Просмотров:	113 
Размер:	2.9 Кб 
ID:	340058

    степ респонс после FIR выглядит как-то так:
    Нажмите на изображение для увеличения. 

Название:	step.png 
Просмотров:	129 
Размер:	1.7 Кб 
ID:	340060
    только что из этого можно извлечь полезного я не знаю. Разве только наложить сверху степ АС.

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    Но в целом, как мне кажется, оно не должно влиять на НЧ.
    так оно и не влияет на НЧ (в полосе 0-500 Гц фаза меняется менее чем на 0.5 градуса, а амплитуда просто постоянна). Дело наверное в том, что тот инструмент, который на слух вы считаете НЧ, на самом деле содержит обертоны в области 3-15 кГц. А они подвергаются трансформации.
    Последний раз редактировалось _Сам_; 19.02.2019 в 13:21.
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  17. #56
    Завсегдатай Аватар для E.Sokol
    Регистрация
    17.01.2015
    Адрес
    Урал
    Сообщений
    1,290

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от _Сам_ Посмотреть сообщение
    гиббсовские колебания абсолютно реальны. Достаточно взять осциллограф и посмотреть на меандр с выхода ЦАП.
    К сожалению, мы слышим не глазами. Ну и цитата:

    «In mathematics, the Gibbs phenomenon, discovered by Henry Wilbraham (1848) and rediscovered by J. Willard Gibbs (1899), is the peculiar manner in which the Fourier series of a piecewise continuously differentiable periodic function behaves at a jump discontinuity.»
    Никто не слышал эффект Гиббса, пока не увидел его в результате разложения ступеньки в ряд Фурье.

    Цитата Сообщение от _Сам_ Посмотреть сообщение
    у ЦАП с корректно реализованным восстанавливающим фильтром вы увидим, как на дробную часть шага выборки сдвигается всплеск вида f(x)=sin(x)/x. Исчезать ничего не должно.
    Не увидим мы sinc в ЦАПе - он же не умеет в прошлое, ещё и бесконечное.

    ---------- Сообщение добавлено 22:03 ---------- Предыдущее сообщение было 19:38 ----------

    Собственно вот.

    В архиве 6 файлов:
    1. оригинальный аудио-отрывок
    2. пропущенный через фазовращатель с сильным эффектом Гиббса
    3. пропущенный через фазовращатель с 3-мя частотами 4-го порядка
    4. пропущенный через тот же фазовращатель с 3-мя частотами 4-го порядка, но реверсированный
    5, 6. сами фазовращатели
    Нажмите на изображение для увеличения. 

Название:	Hilbert Transform (512 samples, 96kHz).png 
Просмотров:	139 
Размер:	22.6 Кб 
ID:	340138Нажмите на изображение для увеличения. 

Название:	allpass (512 samples, 96kHz).png 
Просмотров:	138 
Размер:	20.4 Кб 
ID:	340137

    Нужно на слух определить, что где.

  18. #57
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от E.Sokol Посмотреть сообщение
    Не увидим мы sinc в ЦАПе - он же не умеет в прошлое, ещё и бесконечное.
    бесконечно не умеет, но конечный sinc-подобной всплеск получится вроде фильтра Ланцоша
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  19. #58
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    Прицепил то, что подгружал в FIR.
    прицепил свою версию фильтра на 255 тапов на основе сглаживания АЧХ ваших АС
    Fltr256.wav.Rar
    отличия от сгенерированного рефейзом состоит в том, что я аппроксимировал ФЧХ в области 7 кГц двумя звеньями 2-го порядка. В результате у меня ГВЗ в окрестности 0 Гц стремится не к нулю, как у рефейза, а к постоянной величине.

    Ниже изображено сравнение ФЧХ моего (зелёный) и вашего (красный) фильтров
    Нажмите на изображение для увеличения. 

Название:	ФЧХ.png 
Просмотров:	126 
Размер:	4.7 Кб 
ID:	340715
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

  20. #59
    Завсегдатай Аватар для mr-marlen
    Регистрация
    03.05.2009
    Адрес
    Москва, Бун.Луга
    Возраст
    36
    Сообщений
    4,016

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от _Сам_ Посмотреть сообщение
    прицепил свою версию фильтра на 255 тапов на основе сглаживания АЧХ ваших АС
    Fltr256.wav.Rar
    Послушал. По впечатлениям результат аналогичный. IIR версия нравится больше. Приложил импульсы для одной АС примерно из точки прослушивания в формате Arta.
    Вложения Вложения

  21. #60
    котег
    Автор темы
    Аватар для _Сам_
    Регистрация
    04.05.2010
    Адрес
    ☂☂☂Питер☂☂☂
    Сообщений
    2,305

    По умолчанию Re: О погрешности конвольвера

    Цитата Сообщение от mr-marlen Посмотреть сообщение
    Послушал. По впечатлениям результат аналогичный.
    В том смысле, что коррекция фазы FIRом на 7 кГц существенно изменяет восприятие НЧ, я правильно понимаю?
    «Cобрать стадо из баранов легко, трудно собрать стадо из кошек» — Сергей Капица

Страница 3 из 4 Первая 1234 Последняя

Социальные закладки

Социальные закладки

Ваши права

  • Вы не можете создавать новые темы
  • Вы не можете отвечать в темах
  • Вы не можете прикреплять вложения
  • Вы не можете редактировать свои сообщения
  •