Отправлено Mega 16 февраля 2005 г. 17:35

Есть сигнал при входной частоте дискретизации. Есть мысль скинуть его на более низкую частоту дискретизации. В качестве модели делаю так:

clear all
fd=200;
f1=2;f2=40;
N=1024;
t=0:1/fd:(N-1)/fd;
s1=0.5*sin(2*pi*f1*t);
s2=sin(2*pi*f2*t);
s=(1+s1).*s2;

sp=fft(s);
n1=150;n2=250;
c1=sp(n1:n2);
c2=sp((N-n2+2):(N-n1+2));
c=[c1 c2];
figure,plot(abs(sp))
figure,plot(abs(c))

cv=ifft(c);
figure,plot(real(cv))

На выходе что-то похожее на желаемый результат, и одновременно совсем не то что хотелось бы получить.

• Просю пардону. Не матлаблю. — -=ВН=-   (16.02.2005 18:27, 666 байт)
  • Ответ: (+) — Mega   (16.02.2005 18:54, 483 байт)
    • Ответ: — -=ВН=-   (16.02.2005 19:35, 1867 байт)
      • Окна FFT идут стык-в-стык. Если обрезается часть спектра - то при iFFT могут появится краевые эффекты. — AntZ   (16.02.2005 19:54, пустое)
        • Так я и писал, что наложения будут. — -=ВН=-   (16.02.2005 19:57, пустое)
• А нельзя ли просто отбросить каждые N-1 из N самплов, предварительно пропустив сигнал через антиалиасинговый low-pass фильтр? — AntZ   (16.02.2005 17:40, пустое)
  • Дело не в данном конкретном случае, а в идее - как можно строить DDC на основе Фурье — Mega   (16.02.2005 17:49, пустое)
    • Строить можно, но гораздо проще использовать полифазную декомпозицию с последующей децимацией и фильрацией на низкой. — AntZ   (16.02.2005 17:54, 65 байт)
      • Я же и говорю, что это просто пример примитивный, чтоб легче было разобраться, а так задача стоит в многоканальном DDC, где каналы не обязательно кратны бинам Фурье. — Mega   (16.02.2005 17:59, пустое)
        • Можно я задам несколько вопросов — AntZ   (16.02.2005 18:13, 551 байт)
          • Ответ: (+) — Mega   (16.02.2005 18:29, 515 байт)