[an error occurred while processing this directive]
Но из этого можно что-то выжать только для определенного вида ИХ фильтра. Вобщем зря я это... Вы правы, ничего хорошего из этого не получится.
(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)
Отправлено
andy_P
03 июля 2006 г. 21:09
В ответ на:
Ответ+
отправлено andy_P 03 июля 2006 г. 20:43
Составить ответ
|||
Конференция
|||
Архив
Ответы
Жалко, что так. А как заманчиво выглядит обертка! Синтезировал фильтр в частотной области (прямо на лету) и тут же пропустил через него сигнал все в той же частотной области... Ан нет! Накаси-выкуси! Что ж буду брать иффт, брать буду ее действительную часть, затем снабжать все это добро нулями и обять ффтить... :-(
—
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
,
пустое
)
Отправка ответа
Имя (обязательно):
Пароль:
E-mail:
NoIX ключ
:
Запомнить
Тема (обязательно):
Сообщение:
Ссылка на URL:
Название ссылки:
URL изображения:
Перейти к списку ответов
|||
Конференция
|||
Архив
|||
Главная страница
|||
Содержание
E-mail:
info@telesys.ru