Отправлено andy_P 04 июля 2006 г. 12:31 В ответ на: Вроде так. Но тогда получается, что это fft(exp(-j*2*pi*n/(2*N))) n=0..N-1 синк и есть, а я в матлабе поглядел - не очень-то он синк напоминает, скорее совсем не напоминает. отправлено <font color=gray>homekvn</font> 04 июля 2006 г. 11:49

Ну не знаю...

Если удобнее то так

если циклически свернете прореженную нулями ЧХ Вашего фильтра с
fft([ones(1,N/2) zeros(1,N/2)]) - это вроде бы sinc комплексный
то получите
1. Фильтр у которого в четных позициях получаются отсчеты Вашей ЧХ
2. ИХ фильтра будет равна 0 во второй половине отсчетов.

F = fft([ones(1,N/2) zeros(1,N/2)]):
cyмма N/2членов геом. прогрессии
Fk = sum{i=0..N/2-1}(exp(-j*2*pi*i*k/N)) = (1-exp(-j*2*pi*i*k/N)^(N/2))/(1-exp(-j*2*pi*i*k/N))

Что вроде бы должно Вам подойти.

• Там немножко по-другому получается (+) — homekvn   (04.07.2006 13:04 212.185.161.237, 1027 байт, ссылка)
  • Ответ: Можно и так. То что я написал нужно делать для разбавленной нулями длинной выборки. У Вас нечетные (или четные??? запутался) коэффициентя получаются сверткой с H0. Результат я думаю будет аналогичен.+ — andy_P   (04.07.2006 13:37 80.82.63.185, 193 байт)
    • Если просто добавить нулей сзади, то четные отсчеты БПФ от такого сигнала совпадают с отсчетами БПФ исходного сигнала, в который мы нули не добавляли. Нечетные отсчеты и есть моя головная боль. Но то, что Вы написали это интересно (+) — homekvn   (04.07.2006 13:50 212.185.161.237, 502 байт)
      • Я просто предлагаю проинтерполировать ЧХ H0 чтобы получить H. Если выбрать интерполирующий фильтра таким, как я написал, то отсчеты ИХ = ifft(H) будут равны 0 для n>N/2. А четные отсчеты H будут совпадать с отсчетами H0 — andy_P   (04.07.2006 14:15 80.82.63.185, 1 байт)
        • Чего-то не получается (+) — homekvn   (04.07.2006 15:16 212.185.161.237, 1039 байт)
          • Ответ:+ — andy_P   (04.07.2006 15:29 80.82.63.185, 182 байт)
            • Спасибо. Понял. Действительно так работает. Только это верно Вы ранее заметили, что сворачивать с синком в частотной области накладнее будет, чем в лоб сделать иффт, затем добавить нули, затем опять ффт. Это выгоднее будет даже несмотря на прореженность ЧХ, которую мы сворачиваем. (Только что сказанное справедливо для больших длин фильтров) — homekvn   (04.07.2006 15:51 212.185.161.237, пустое)
              • Тут единственный момент в том,что то с чем сворачиваете усекать можно и проверить, может рипплы в ИХ от N/2 до N будут приемлемыми+ — andy_P   (04.07.2006 16:01 80.82.63.185, 276 байт)
                • + — andy_P   (04.07.2006 18:51 80.82.63.185, 109 байт)
                  • Еще не проверял, но, чувствую, в аудиообработке не пройдет - слышно будет. Проще уж фильтр в раза в полтора-два по длине сократить. — homekvn   (04.07.2006 18:56 212.185.161.237, пустое)
                    • Hi-end у Вас немеренный получается :-) В коммуникашке все проще и сложнее - таких требований по разрядности нет, зато полосы широкие шибко :-) — andy_P   (04.07.2006 19:28 80.82.63.185, пустое)
                      • А я-то причем? Мне сказали пилить - я пилю. По мне так вообще КИХ-фильтры для такого приложения, как у меня, не нужны. До селе прекрасный звук и при помощи БИХ-ов в кабине машины получался. А заказчик - увы - непереубеждаем. — homekvn   (04.07.2006 20:32 212.185.161.237, пустое)
                        • Надеюсь, влезете в DSPшку со всем этим хозяйством. Успехов. — andy_P   (04.07.2006 20:57 80.82.63.185, пустое)

E-mail: info@telesys.ru