Area : DSP

Date : Wed Sep 30, 11:06                                                       
From : Vladimir L Vasilevsky                               2:5020/211.60
To   : All                                 
Subj : Ламеры                                                                
────────────────────────────────────────────────────────────────────────────────

 Здравствуйте, сабжевики (или как там у вас принято говорить). Как поживаете ?
Отличать БПФ от ДПФ уже научились ? Hу надо же. Это же вам не какой-нибудь
мастдайный пень под крутую кваку глючно разгонять.

  Вот вам пример: Hello Filter называется. По секрету скажу, что, например, из
таких фильтров можно даже многополосный эквалайзер собрать. Hо что-то я весьма
сомневаюсь, что у кого-либо из вас хватит на это толку.


 Vladimir

 P.S. Какие же вы все-таки сирые и убогие,,,


/* Расчет коэффициентов для полосового рекурсивного фильтра 2-го порядка */
/* Сам фильтр:  Y[n]=X[n]+A*Y[n-1]+B*Y[n-2] */

#include <stdio.h>
#include <math.h>
#if !defined M_PI
#define M_PI 3.14159265
#endif

double T=0.05;       /* Постоянная времени затухания */
double Fdis=8000.0;  /* Частота дискретизации */
double F=1234.5678;  /* Центральная частота   */

void main(void)
{
double dt,A,B,C1,C2;
printf("\n Hello Filter \n");
dt=1.0/Fdis;
C1=exp(-dt/T)*cos(2*M_PI*F*dt);
C2=exp(-2*dt/T)*cos(4*M_PI*F*dt);
A=(C2-C1)/(C1-1);
B=C1-A;
printf("\n A=%f B=%f",A,B);
}

--- Yet another GoldED 2.41+
 * Origin:  Hет такого незнания, которым бы не гордились (2:5020/211.60)