Здраствуйте. Помогите разобраться с проблемой. Хочу использовать ресамплер в демодуляторе для подстройке по тактам. Ресамплер хочу реализовать на базе модифицированного фильтра Фарроу 3-го порядка. В схеме ресамплер стоит после КП (комплексного перемножителя) и ФНЧ. Соответственно ФД в ФАПЧ вычисляет ошибку по ресемплированным отсчётам. Реализацию ресамплера делал на базе кода который есть на этом сайте.
Теперь сама проблема.
Есть:
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 888. Теперь получается что при расчёте коефициентов а3,а1,а2 нужны отсчёты которые находятся в массиве за адресом 100 888, 100 889, 100 890, 100 891, а там соответственно на 100 000 отсчёте в массиве y0[i] ещё нули .
Как мне правильно сделать в моём случае. Или возможно я где то допустил ошибку при переводе из С++ в Дельфи.
Вот код Дельфи:
сам ресамплинг
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);

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