Не понял. Зачем её вести? нужно просто получить сигнал на выходе без пост. составляющей. IIR первого порядка поможет?

Спасибо, Сергей, что не устаёте отвечать на вопросы. непонятно зачем понижать частоту до 2Гц? Ведь меня интересует полоса 0,5-1000Гц. Требования(идеализированые): минимальный алиасинг и линейная фаза. Если искать компромис то допустимо снизить требования к алиасингу, но не к линейности фазы. Если ст...

Необходимо осуществить фильтрацию сигнала. Сигнал оцифрован с частотой 128кГц. Интересуемый диапазон частот лежит в пределах 1000Гц. Именно до 1000Гц - это интеесующий диапазон частот, в котором не должно быть помех (помехи должны быть минимизированы максимально) Так же необходимо произвести фильтра...

у меня получается такая петрушка: функция farrow3 - сам фильтр выдаёт правильный результат (в сравнении с "плавучей" версией). Но сам результат на выходе отличается достаточно сильно (если ошибка в спектре для плавучки составляет 0,00038, то для фиксированой точки 0,011 - это просто пример, чтобы да...

Чистый С не предоставлю. Да думаю тут сложного нет. Цикл и переменные. Truc отбрасывает дробную часть остаётся целая. Диапазон Х тут -1...0 вот, вся загвоздка в том, что как можно представить -1..0 в целых числах? Напрмиер числа -0,5 прсото нет. Вместо него может быть -16384 (потому что перём число...

В вашем коде при определенных условиях X выходит за пределы -1<=X<0 или j ошибочно вычесляется не помню точно. Вот мой код. Где X всегда в рамках -1<=X<0 и j правильная. procedure TForm1.Resempling; var I, J: Integer; X, T, Steep: Extended; begin // Размер исходного массива N = Trunc((N1 - 1) * f0 ...

Чуть более шустрый алгоритм Герцеля http://www.dsplib.ru/forum/viewtopic.php?f=5&t=111&p=665#p665 всего 3 умножения на один семпл. Вот код procedure TForm1.Hercel; var I: Integer; S, V0, V1, V2, a, Re, Im, fRe, fIm: Double; begin Re := Cos(2 * Pi / 32); Im := Sin(2 * Pi / 32); a := Re * 2; V2 := 0;...

Вот так должно работать #include <stdio.h> #include <stdlib.h> #include <math.h> #define K 32 short FIX_MPY(short a, short b) { /* shift right one less bit (i.e. 15-1) */ int c = ((int)a * (int)b) >> 14; /* last bit shifted out = rounding-bit */ b = c & 0x01; /* last shift + rounding bit */ a = (c >...

Взялся переписать подпрограмму фильтрации Фильтра Фарроу 3-го порядка, позаимствованную с этого сайта. До недавнего времени, она успешно функционировала как есть. Но сейчас появилась задача переписать код для использования только с целыми числами. Сам фильтр прекрано функционирует (проверено в сверк...

Я использовал вашу бесценную статью о фильтрах Фарроу и ваш код. Результаты приятно удивили. Теперь дело за железными экспериментами - надо как-то упихать это дело в мой микроконтроллер, предварительно переведя всё на fixed point Большое вам спасибо! Выбор мой такой: фарроу 3-го порядка и БПФ на 32 ...