akrabad

Я получаю G1 и G2 вот так for iter = 1:1:Nfreq-2 %Цикл по всем сигналам Temp = fft(ch_1, Nfft); %fft первой АЦП G1(iter) = Temp(iter+1); %G1 - 1x34 % Temp - 1x35 Temp = fft(ch_2, Nfft); %fft второй АЦП G2(iter) = Temp(iter+1); end F0 = conj(G1) ./ conj(G2);% По формуле нам надо только(-ws/2 + w) Правильно ли? Под нормированием вы подразумеваете только приведение суммы коэфф. к 1? После ifft ИХ выглядит так. "Сдвигали" ли вы ее?

Что вы подразумеваете под эталонными частотами? А для чего вы использовали FFT в алгоритме получения фильтра? Подал 14 частот на АЦП; По этим сигналам получил 14 значений для разностных АЧХ и ФЧХ для фильтра F2; Произвел интерполяцию; Получил 1025 значений АЧХ и ФЧХ(теперь буду меньше); Задал комплексный коэффициент передачи; Выполнил IFFT и получил 1025 значений фильтра.

Количество тапов = 1025. F1 - задерживающий фильтр. F2 - фильтр, полученный обратным преобразование фурье по необходимой АЧХ и ФЧХ + окно Ханна.

Если я корректирую канал фильтром от 0 до Fs/4, где Fs - результирующая частота дискретизации двух АЦП, то сигналы, поступающие от 0 до Fs/4 корректируются. Но при использовании апсамплинга и фильтра от 0 до Fs/2 сигнал не корректируется. Они не корректируются потому что при апсамплинге возникает копия спектра. Если задавать АЧХ и ФЧХ по значениям для конкретной частоты, то мы можем корректировать эту частоту.

При корректирующей фильтрации канала, интерполированного нулями в два раза, происходит корректировка не только полезного сигнала, но и его "копии". Т.к. мы эту "копию" откорректировали, то она появилась на выходе. Как это надо учесть при расчете фильтров(а)? По формулам из статьи, предоставленной пользователем _Anatoliy, я так до конца и не уяснил метод корректировки, а точнее, его реализацию.

При моделировании фильтров F1 и F2, по формулам из статьи, с частотой дискретизации = 206438400, тау1 = 1.001*10^(-9), тау2 = 1.003*10^(-9), где тау1 и тау2 - постоянные времени RC-цепочек, получил график АЧХ, похожий на график АЧХ в статье. Синий - F1. Зеленый - F2. Эти фильтры зависят, по-сути, только от этих тау. Как узнать их для реального устройства? Как Вы их узнавали? Каким-то детектором? И еще один вопрос. Для замены двух фильтров одним необходимо поделить F2(w) на F1(w)?

Хорошо, попробую. Спасибо за статью! Я ее до этого видел, вроде бы, но из-за моего незнания отбросил.

Посмотрел по-быстрому статью. Схема корректировки содержит два фильтра, которые фильтруют "эффекты", вносимые SHA, и уничтожают(eliminate) мнимые компоненты. Результаты корректировки в статье отличные. Пики на уровне - 90dB, что мне и нужно. Но я не могу использовать SHA. Его попросту нет. У меня есть две АЦП в устройстве и ПЛИС, которая и будет выполнять цифровую корректировку на заднем плане (digital background calibration) Если я сказал что-то не так по поводу "вашей" статьи, то, пожалуйста, поправьте. В первом сообщение я спрашивал о "физике" процесса корректировки фильтром. Фильтр sinc компенсирует фазовое отклонение сигнала. Пускай отсчеты идеального сигнала приходят раньше, чем надо. Фильтром sinc мы их пересчитываем на те значение, которые должны были быть. Создаем задержку. Для идеального сигнала мы ее знаем. Тогда почему после Fs/4 такой косяк со знаком возникает? http://www.labbookpages.co.uk/audio/beamfo...ionalDelay.html Пример, после которого я решил заменить фильтр из статьи ( https://drive.google.com/file/d/0B-a0qo3gvn...1eTA/view?pli=1 ) на фильтр sinc.

https://drive.google.com/file/d/0B-a0qo3gvn...iew?usp=sharing https://drive.google.com/file/d/0B-a0qo3gvn...iew?usp=sharing Автор первой статьи продвигает немного измененную формулу для расчета импульсной характеристики фильтра, которая не подходит в моем случае - для разных частот нужны разные иксы. Из второй статьи беру детектор, который как раз и сходится к тому x = 0.001 для частот меньших Fs/4. Что и нужно.

Здравствуйте. Для повышения частоты дискретизации с наименьшими затратами используют простой способ - параллельное включение 2-x АЦП или более, что дает увеличение частоты дискретизации в два раза. Минус такого способа состоит в том, что возникают ошибки усиления, смещения и фазы. Это прекрасно видно в спектре сигнала. К примеру, для сигнала с частотой f0, пропущенного через такую систему АЦП, будут присутствовать "ложные" пики. Для ошибки смещения - на частоте Fs/2. Для гэйна и фазовой ошибок на частоте (Fs/2 - fo). Существуют варианты корректировок таких сигналов. Самый простой из них - с помощью FIR фильтра. Этот фильтр, по сути, является задержкой во втором канале относительно первого. Он пересчитывает отсчеты сигнала так, чтобы они соответствовали своему "реальному положению". Импульсная характеристика такого фильтра есть не что иное, как кардинальный синус( sinc(Pi *x ) ), где x - рассчитанная задержка. Пускай x = 0.001; (Для Fs = 2064384 и задержки во втором канале в 1 нс). Сигнал корректируется. С определенными результатами, но корректируется. После Fs/4 сигнал корректируется только если x = - 0.001, хотя мне надо, чтобы x был равен 0.001. Объясните пожалуйста, в чем проблема. И как это можно исправить. Информацию о методах корректировок беру из статей на английском, ибо наши об этом почти не пишут.