Отправлено andy_P 03 июля 2006 г. 20:43 В ответ на: Хм-м! Про четные и нечетные отсчеты я в основном посте написал. Про свертку там же тоже упомянул, но проблема была в том, с чем сворачивать (+) отправлено <font color=gray>homekvn</font> 03 июля 2006 г. 20:02

Вроде бы все так. У меня тоже циклическая свертка в конце получилась

Интересно вот что:

real(fft(exp(-j*pi/N*(0:N/2-1))))) = ones(1,N/2);

соответственно циклическая свертка с единицами будет...

• Но из этого можно что-то выжать только для определенного вида ИХ фильтра. Вобщем зря я это... Вы правы, ничего хорошего из этого не получится. — andy_P   (03.07.2006 21:09 80.82.63.185, пустое)
  • Жалко, что так. А как заманчиво выглядит обертка! Синтезировал фильтр в частотной области (прямо на лету) и тут же пропустил через него сигнал все в той же частотной области... Ан нет! Накаси-выкуси! Что ж буду брать иффт, брать буду ее действительную часть, затем снабжать все это добро нулями и обять ффтить... :-( — homekvn   (03.07.2006 21:27 212.185.161.237, пустое)
    • Либо свертка с синком в частотной области, чтобы окно на ИХ во временной наложить. Синк подусечь можно, но тогда алиасить будет. Кстати, если синк недлинный взять, то вычислений немного потребуется. Хотя о чем я :-) Ведь FFT c немеренной разрядностью у Вас в этом проекте? — andy_P   (03.07.2006 21:56 80.82.63.185, пустое)
      • А почему с синком сворачивать надо? Там же хрен знает сворачивать надо, вот с ней: fft(exp(-j*2pi*n/(2*N))) при n=0..N-1. А она ж вроде как не очень синк?... Да и правы Вы, будь она даже синком, укорачивать бы ее я не стал - опасно бо зело. — homekvn   (03.07.2006 22:05 212.185.161.237, пустое)
        • С синком прореженную нулями ЧХ сворачивать - эквивалентно наложению прямоугольного окна во временной области. — 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, пустое)
    • Тк ИХ у Вас действительная, то FFT половинной длины можно использовать. Все легче. — andy_P   (03.07.2006 21:48 80.82.63.185, пустое)
      • Да, эт-точно, можно. А можно и два канала в одном флаконе обработать, сделав один комплексный сигнал из двух действительных каналов. — homekvn   (03.07.2006 22:00 212.185.161.237, пустое)
• тьфу, real(fft(exp(-j*2*pi/N*(0:N/2-1))))) = ones(1,N/2). — andy_P   (03.07.2006 20:46 80.82.63.185, пустое)
  • Да, результат действительно интересный. Раньше такого не замечал. Просто не приходило в голову рассматривать действительную часть от fft. Обычно меня только модуль интересует. Только вот поможет ли это?... — homekvn   (03.07.2006 21:06 212.185.161.237, пустое)

E-mail: info@telesys.ru