hardware FFT-based FIR filter, Помогите новичку!!! Опции

[thermit]

Есть только 2 способа секционирования апериодической свертки:
Перекрытие с накоплением и перекрытие с суммированием. Все остальное - лженаука.

Код

clear all;
N=1024;

h=randn(1,N);

x=randn(1,100*N);


H=fft([h zeros(1,N)]);
%Перекрытие с накоплением
y1=[];
for i=1:N:length(x)-N
    X=fft(x(i:i+2*N-1));
    Y=H.*X;
    z=ifft(Y);
    y1=[y1 z(N+1:end)];
end;
%Перекрытие с суммированием
y2=[];
mem=zeros(1,N);
for i=1:N:length(x)-N
    X=fft([x(i:i+N-1) zeros(1,N)]);
    Y=H.*X;
    z=ifft(Y);
    
    y2=[y2 z(1:N)+mem];
    mem=z(N+1:end);
end;

%образец
yy=conv(h,x);
yy=yy(N+1:end);
y2=y2(N+1:end);

grid on;
plot(yy);
hold on;
plot(y1,'r');
plot(y2,'g');

[_Anatoliy]

Спасибо!
Да,overlap-add/overlap-save я моделировал,прекрасно работают.
Подвёл тёзка.