Отправлено AlexeyV 30 ноября 2004 г. 09:15

Уважаемые Господа, обьясните убогому суть.(от куда то скачал)
for (j=0;j{
co=cos(pi*j/nsize);si=sin(pi*j/nsize);
a=b=u=uo=0;u1=0;//xx[0];
//' for i=NSIZE-2 to 0 step -1
for (i=0;i{u=2*u1*co-uo+xx[i];
//' ?i;u; xx(i)
//' a=a+u*co-u1+xx(0)
//' b=b+u*si
uo=u1;
u1=u;
}
a=u1*co-uo+xx[0];b=u1*si;a=a/nsize;b=b/nsize;
c=sqrt(a*a+b*b);

cout<<"\nk="<}
Мне надо определить одну j-ю спектральную составляющую (только её модуль)
В моём случае nsize=80, j=20. Внешний цикл не нужен,
закоментированные стоки во внутреннем, особенно 1я и 2я вызывают сомнения.
Что должно измениться, чтобы алгоритм правильно заработал?
Спасибо.