таки съелось - алгоритм - повтор («Телесистемы»: Конференция «Микроконтроллеры и их применение»)

---

Отправлено VVT 14 июля 2002 г. 11:29 В ответ на: Близкий к нижнему вопрос - анализ спектра и БПФ отправлено VVT 14 июля 2002 г. 11:08

автокорреляция
int16  i,k;
for (i=0; i<=P; i++)
{
  R[i]=0;
  for(k=0; k<(window-1); k++)
    R[i]+=xw[k]*xw[k+i];
}
Сам алгоритм

float alfa;  
int16 i,k;
float E,B,PARCOR;
float b[64];
beta[0]=1;            
E=R[0];
for (i=1; i<=P; i++) 
{ 
  B = R[i];
  for (k=1; k<(i); k++)  B +=R[i-k]*beta[k];
  PARCOR=-B/E;
  for (k=1; k<(i); k++) b[k]=beta[i-k];
  for (k=1; k<(i); k++) beta[k] += PARCOR*b[k];
  beta[i]=PARCOR;
  E *= (1 - beta[i]*beta[i]);
}
alfa=E;