Отправлено Vadim Kudryavtsev 15 января 2005 г. 12:47 В ответ на: C фильтром вроде получилось, но помогите разобраться... отправлено AR 14 января 2005 г. 18:14

Гляньте на приведенный код и подумайте, почему циклическая свертка работает лучше. ФФТ ведь инвариантно к циклической свертке.
А вот как сделать, чтоб при линейной свертке хорошо фильтровало,
подумайте пока сами.
Пока времени нет выдать готовый результат, да и задачу может быть немного подробнее раскрыть не мешало бы. :-)

clear all
fd=1000;
N=1000;
t=0:(1/fd):(N-1)*1/fd;
f1=100;f2=125;
A1=30;A2=1;
s1=A1*sin(2*pi*f1*t)+(A1-A1/4)*sin(2*pi*(f1+5)*t);
s2=A2*sin(2*pi*f2*t);
n=0*rand(1,N);
s=s1+s2+n;

b = fir1(680,[0.18 0.23],'stop',chebwin(681,80));

figure(1)
plot(10*log10(abs(fft(s))))
hold on

plot(40+10*log10(abs(fft(b,length(s)))),'r')

%sk=upfirdn(s,b,1,1);
%sk=sk(length(b):end);
sk=cyclconv(s,b);

plot(3+10*log10(abs(fft(sk))),'g')
grid on
zoom on
hold off

function y=cyclconv(x,z);

y=conv(x,z);
y=y(1:length(x))+[y(length(x)+1:end) zeros(1,2*length(x)-length(y))];

• Да вроде бы подробнее некуда. Обнаружение слабого сигнала находящегося рядом по частоте с сильным... — AR   (15.01.2005 18:14, 400 байт)
  • P.S. Посмотрел Ваш код и пришел к выводу, что Вы батюшка схитрили. Изменили число точек преобразования. В данном случае да, но в жизни такой возможности как правило не бывает :( — AR   (15.01.2005 18:33, пустое)
    • Ответ: Конечно :-) Я сразу и предупредил :-) (+) — Vadim Kudryavtsev   (16.01.2005 13:13, 908 байт)
      • Еще раз - где я сыплюсь? ... — AR   (16.01.2005 14:34, 823 байт)
        • Ответ: Передний хвост свертки тоже отрезать надо (+) — Vadim Kudryavtsev   (16.01.2005 15:04, 913 байт)
          • Трудно быть тормознутым! Все получилось, спасибо! Но... — AR   (16.01.2005 18:17, 263 байт)
• Ответ: В своем исходном примере увеличьте N (к примеру до 60000) и тогда все встанет на свои места. Это оценка короткого преобразования Фурье гадит — Vadim Kudryavtsev   (15.01.2005 16:11, пустое)