А почему с синком сворачивать надо? Там же хрен знает сворачивать надо, вот с ней: fft(exp(-j*2pi*n/(2*N))) при n=0..N-1. А она ж вроде как не очень синк?... Да и правы Вы, будь она даже синком, укорачивать бы ее я не стал - опасно бо зело. («Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

[an error occurred while processing this directive] А почему с синком сворачивать надо? Там же хрен знает сворачивать надо, вот с ней: fft(exp(-j*2pi*n/(2*N))) при n=0..N-1. А она ж вроде как не очень синк?... Да и правы Вы, будь она даже синком, укорачивать бы ее я не стал - опасно бо зело.

(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

Отправлено homekvn 03 июля 2006 г. 22:05
В ответ на: Либо свертка с синком в частотной области, чтобы окно на ИХ во временной наложить. Синк подусечь можно, но тогда алиасить будет. Кстати, если синк недлинный взять, то вычислений немного потребуется. Хотя о чем я :-) Ведь FFT c немеренной разрядностью у Вас в этом проекте? отправлено andy_P 03 июля 2006 г. 21:56

Составить ответ ||| Конференция ||| Архив

Ответы

С синком прореженную нулями ЧХ сворачивать - эквивалентно наложению прямоугольного окна во временной области. — andy_P (04.07.2006 10:47 80.82.63.185, пустое)
- Вроде так. Но тогда получается, что это fft(exp(-j*2*pi*n/(2*N))) n=0..N-1 синк и есть, а я в матлабе поглядел - не очень-то он синк напоминает, скорее совсем не напоминает. — homekvn (04.07.2006 11:49 212.185.161.237, пустое)
  - Ответ:+ — andy_P (04.07.2006 12:31 80.82.63.185, 522 байт)
    - Там немножко по-другому получается (+) — 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, пустое)
  - Думаю, что так будет выглядеть "периодический" синк. — homekvn (04.07.2006 12:13 212.185.161.237, пустое)
Ну а fft у меня на плавучке будет. — homekvn (03.07.2006 22:06 212.185.161.237, пустое)

Перейти к списку ответов ||| Конференция ||| Архив ||| Главная страница ||| Содержание

E-mail: info@telesys.ru

Ответы

Отправка ответа