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

[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]

Это получает так - x:=(((t1[n]-t0[j])*Fs0-2))+ErrT;, где ErrT - это и есть временная ошибка, так будет правильно ?

[svalery]

Это получает так - x:=(((t1[n]-t0[j])*Fs0-2))+ErrT;, где ErrT - это и есть временная ошибка, так будет правильно ?

имхо не лучший способ разбираться с такими вещами начиная с кода, вот м-скрипт примерно того что у вас, если я правильно понял, что вам нужно

Код

clc;
clear;

N = 1; % кол-во периодов синуса
spp = 20; % кол-во семплов на период
ii = 1:spp*N;

y0 = sin((ii/spp)*2*pi);

plot(y0,'--*');

x = -0.3; % ошибка по времени, изменяется в пределах [-1...0]

for jj=1:length(ii)-3
    a3 = (y0(jj+3)-y0(jj))/6.0 +(y0(jj+1)-y0(jj+2))/2.0;
    a1 = (y0(jj+3)-y0(jj+1))/2.0-a3;
    a2 = y0(jj+3)-a3-a1-y0(jj+2);
    y1(jj) = x*(x*(x*a3+a2)+a1)+y0(jj+2);
end;

figure;
plot(y1,'--*');

При изменении x (на рисунке mu), на 2й картинке будет видно как выборка "плывёт" по синусу относительно выборок синуса 1й картинки