Тактовая синхронизация (ресамплинг), реализация на фильтре Фарроу Опции

[bigo]

Здраствуйте. Помогите разобраться с проблемой. Хочу использовать ресамплер в демодуляторе для подстройке по тактам. Ресамплер хочу реализовать на базе модифицированного фильтра Фарроу 3-го порядка. В схеме ресамплер стоит после КП (комплексного перемножителя) и ФНЧ. Соответственно ФД в ФАПЧ вычисляет ошибку по ресемплированным отсчётам. Реализацию ресамплера взял с http://electronix.ru/forum/index.php?showt...=73192&st=0
Теперь сама проблема.
Есть:
K=700 000; количество отсчётов
Fs0 = 294600 Гц;
Fs1 = 292000 Гц;
Ресамплер работает с массивом y0[i], который заполняется отсчётами после прохождения КП и ФНЧ.
Теперь если, например, при перерасчёте индекса исходного сигнала (взято из вышеуказанного форума - " int j = int((double)i*Fs0/Fs1)-2;") скажем 100 000-го получается j= 100 000*294600/292000 =100 890. Теперь получается что при расчёте коефициентов а3,а1,а2 нужны отсчёты которые находятся в массиве за адресом "y1[i] = farrow3(y0+j, x); // фильтр Фарроу" y0+100 890, а там соответственно на 100 000 отсчёте ещё нули :-( .

Вышеуказанный код проверял отдельно при полностью заполненном массиве, всё работает. А как сделать в моём случае, по отсчётно не знаю, потому прошу Вашей помощи в решении этой проблемы. Заранее спасибо.

[bigo]

Буду очень благодарен Вам за помощь bahurin
Возможно проблема в не правильном переводе с С++ на Delphi '
Вот код:
begin
for n:=0 to K-1 do
begin
t1[n]:=n/Fs1;
j:=(Round(n*Fs0/Fs1))-2;
if j<0 then j:=0;
x:=(((t1[n]-t0[j])*Fs0-2));
a3:=(y0[j+3]-y0[j])/6.0 +(y0[j+1]-y0[j+2])/2.0;
a1:=(y0[j+3]-y0[j+1])/2.0-a3;
a2:=y0[j+3]-a3-a1-y0[j+2];

y1:=x*(x*(x*a3+a2)+a1)+y0[j+2];
form1.Series3.AddXY(n,y1, '', clGreen);
end;

Я просто в С++ полный 0 :-)