Отправлено AR 16 января 2005 г. 14:34 В ответ на: Ответ: Конечно :-) Я сразу и предупредил :-) (+) отправлено Vadim Kudryavtsev 16 января 2005 г. 13:13

Частота дискретизации задана жестко, менять низя.
Время наблюдения также.
Как бы не фильтровал (upfirdn, filter, conv и т.д.), что бы не отрезал - в начале, в конце - в фильтрованном сигнале остается просачивание от отфильтрованного сигнала????!!!!
Понятно, когда тем или иным способом частоты поставишь на бины, то проблема снимиться, но это же в реальной жизни невозможно.

clear all
fd=1000; %!!!
N=1024; %!!!
t=0:(1/fd):(N-1)*1/fd;
f1=100;f2=120; %!!!
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));

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

sk=filter(b,1,s); %!!!

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

• Ответ: Передний хвост свертки тоже отрезать надо (+) — Vadim Kudryavtsev   (16.01.2005 15:04, 913 байт)
  • Трудно быть тормознутым! Все получилось, спасибо! Но... — AR   (16.01.2005 18:17, 263 байт)