trofim2 писал(а): А Вы не перепутали к-ты А и B?

В обозначениях, принятых на dsplib

Соответственно:
Y[n]=1/a_0{b_0X[n]+b_1X[n-1+b_2X[n-2]- a_1Y[n-1]-a_2Y[n-2]}

Коэффициенты b:
9.985882797899660000e-001
-1.997176559578946400e+000
9.985882797899660000e-001
Коэффициенты a:
1.000000000000000000e+000
-1.997176559578946400e+000
9.971765595799320100e-001
Но как из них сделать целые я не знаю На 4096 умножить?

Вот такой фильтр попробуйте:
y[n]=0,99505*(X[n]-X[n-1])+0,9901 Y[n-1]
Его только в целые коэффициенты надо перевести....

Исправил.

trofim2 писал(а): Мне бы еще простой ФВЧ с частотой среза не более 1 Гц при частоте выборок 3500 Гц.

А режектор не на 50 Гц, а на 0 Гц не поможет?

[ref]trofim2[/ref], в статье "Цифровые рекурсивные фильтры. Часть 4." http://we.easyelectronics.ru/Theory/cif ... ast-4.html дал ссылку на Вашу реализацию.

8592 А что такое Изи? Сокращение от easy electronics :D 8592 Режекция работает замечательно, но частотная характеристика в области нижних частот от частоты 16 Гц и ниже имеет подъем +4дБ (1,6 раз). В чем у меня может быть ошибка. Я думаю ошибка произошла при "округлении" коэффициентов... Попробуйте...

Не, это ты написал, что h(1)=1
Для преобразования Гильберта у меня это вызвало некоторое недоумение.
Всем известно, что при преобразовании Гильберта постоянная составляющая сигнала не сохраняется.
А h(1) - это как раз сомножитель F(1) - или 0-частота (постоянная составляющая)

8562 Скажите, а если умножить двойку на частоту Найквиста и на нулевую частоту, что произойдет? Ну умножь сам и посмотри. А вообще за этот вопрос - тебе ещё одно предупреждение . Для начала посмотри, что такое БПФ. (Я имею ввиду стандартные программы). Э.Р. Канасевич "Анализ временных последователь...

8559 Я просто не пойму почему такой массив h получается, я думал массив h должен содержать половину нулей и половину двоек и все, объясните а то я не понимаю, тупой наверно Да так оно в принципе и есть (это я про массив h) :D Исключение составляют всего две точки - с 0-частотой и частотой Найквиста...