Отправлено Vadim Kudryavtsev 16 января 2005 г. 13:13 В ответ на: P.S. Посмотрел Ваш код и пришел к выводу, что Вы батюшка схитрили. Изменили число точек преобразования. В данном случае да, но в жизни такой возможности как правило не бывает :( отправлено AR 15 января 2005 г. 18:33

Я смоделировал ситуацию (может и гипотетическую) где все получается, и намекнул, что отсюда плясать и надо (то есть найти тот шаг от моего кода к Ващему, где эта прелесть теряется).
Во-первых, в первом моем ответе хорошо то, что частоты попадают в сетку ДПФ.
Во-вторых, во втором моем ответе, рекомендующем увеличить длину наблюдения опять хороший результат.
Почему?
Я тогда сам еще не знал.

Внимательно посмотрите на выходной сигнал - sk - во временной области
Функция upfirdn дает на выходе сигнал аналогичный функции conv, те
последние length(b)-1 отсчеты расчитаны при помощи добавлени виртуальных нулей к сигналу s. (Спектр сигнала такого сигнала заметно шире)
Вот Вы и получили растекание спектра на выходе.

Или нужно отрезать эти отсчеты
или вместо строчки
sk=upfirdn(s,b,1,1);

использовать filter
sk=filter(b,1,s);

Она не расчитывает этот часто вредный хвост

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