Отправлено AR 14 января 2005 г. 18:14

Есть две синусоиды имитирующие сильный мешающий сигнал(сигнал s1).
Есть полезный сигнал (s2). Когда смотришь спект s1+s2, то понятно - сильный узкополосный сигнал
маскирует слабый за счет эффектов просачивания и т.д. Но если я этот мешающий сигнал отфильтровал,
то почему он продолжает просачиваться и маскировать слабый?! Как можно вытащить слабый сигнал
на фоне сильной помехи, о которой кроме того, что можно получить из спектра (типа центральной, полосы, уровня)
ничего не известно?

clear all
fd=1000;
N=1024;
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));

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);
plot(3+10*log10(abs(fft(sk))),'g')
hold off

• Я тут кое-что изменил в угоду теории, но ушел от практики (+) — Vadim Kudryavtsev   (15.01.2005 12:47, 954 байт)
  • Да вроде бы подробнее некуда. Обнаружение слабого сигнала находящегося рядом по частоте с сильным... — 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, пустое)