Перейти к содержанию
    

Есть только 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');

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Спасибо!

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

Подвёл тёзка.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Присоединяйтесь к обсуждению

Вы можете написать сейчас и зарегистрироваться позже. Если у вас есть аккаунт, авторизуйтесь, чтобы опубликовать от имени своего аккаунта.

Гость
Ответить в этой теме...

×   Вставлено с форматированием.   Вставить как обычный текст

  Разрешено использовать не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отображать как обычную ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставлять изображения напрямую. Загружайте или вставляйте изображения по ссылке.

×
×
  • Создать...