Корректор АЧХ в широкой полосе Опции

[_Anatoliy]

Вошёл в ступор на следующей задаче.Есть АЦП с тактовой Fs,на его вход подаётся сигнал (например синус) на любой частоте в полосе от 0.05*Fs до 0,95*Fs (сигнал может быть как в первой так и во второй зонах Найквиста).В выходном сигнале наблюдается небольшая неравномерность АЧХ,причём она несимметрична относительно частоты Fs/2.Существует ли какой нибудь способ откорректировать АЧХ в заданной полосе частот с помощью фильтрации выходных данных? Если бы сигнал был в полосе хотя бы 0.05*Fs до 0,45*Fs, то задача решилась бы банальным корректирующим КИХ фильтром. А тут как быть?

[andyp]

Ну у авторов на картинках все даже сходится. Может и прокатит. По вычислениям на отсчет - будет дорого. (M+1 ФВЧ кроме всего остального). Да и судя по статье требуется ~1e6 отсчетов чтобы алгоритм сбежался.

[_Anatoliy]

Ну у авторов на картинках все даже сходится. Может и прокатит. По вычислениям на отсчет - будет дорого. (M+1 ФВЧ кроме всего остального). Да и судя по статье требуется ~1e6 отсчетов чтобы алгоритм сбежался.

Угу,спасибо,сяду за модельку.

[_Anatoliy]

Не работает эта штука.Сделал по рисунку 7.
Т.е. как-то работает,если у АЦП есть разброс по К одинаковый во всей полосе,то да - алгоритм сходится и спектр на выходе становится чистым. А вот если разброс по К неодинаковый в полосе то спектр не очищается.И почему-то работает только y0 (y1,y2,y3 практически к нулю стягиваются).

[andyp]

Не работает эта штука.Сделал по рисунку 7.
Т.е. как-то работает,если у АЦП есть разброс по К одинаковый во всей полосе,то да - алгоритм сходится и спектр на выходе становится чистым. А вот если разброс по К неодинаковый в полосе то спектр не очищается.И почему-то работает только y0 (y1,y2,y3 практически к нулю стягиваются).

можно попробовать задать квадратичное рассогласование ЧХ, подставить вместо результатов адаптации известные коэффициенты и посмотреть, уйдет ли помеха. Поможет найти глюки в реализации, если они есть. При беглом просмотре статья производила впечатление настоящей :-). Лучше на многотоновом сигнале смотреть с частотами, не кратными частоте дискретизации. Сразу будет видно, ушел ли алиасинг или остался.

[Anatol']

Насколько я понимаю, можно интерполировать данные с каждого АЦП отдельно фильтром с нужной АЧХ, а потом сложить не ошибившись с разницей во времени. Если разброс характеристик АЦП компенсирован правильно, то настоящий сигнал сложится, а "фальшивый" (Fs - F) сложится в противофазе и исчезнет. При чём независимо от того, в первой он зоне найквиста был, или во второй, всё восстановится правильно. Подозреваю, что для случая >2 АЦП это тоже должно работать.

Только надо не забыть, что компенсировать надо не только АЧХ, но и ФЧХ (вариацию aperture delay - от экземпляра к экземпляру АЦП она может различаться на 20..50..100..чёрт-знает-сколько пс). Любое расхождение даст неполную компенсацию и остатки фантомных частот на Fs-F.

Ещё может вылезти фантомная Fs/2 (для двух АЦП) из-за разницы нуля.

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