Статус: Advanced Member
Группы: Registered
Зарегистрирован: 02.06.2009(UTC)
Сообщений: 346
Сказал «Спасибо»: 50 раз Поблагодарили: 156 раз в 105 постах
Последняя редакция плагина :
FIRFilterDesign.dll Список функций: remez(vg, vr, vw, n) . Документ: КИХ фильтр.sm ( alglib.dll :remez(), FIRFilterDesign.dll :al_convr1d(), al_fftc1d() ) Это аналог маткадовской функции remez(vg, vr, vw, n). Работает аналогично.
Документ:
КИХ фильтр.sm Дополнение (плагин):
FIRFilterDesign.dll П.С. В функции не обрабатываются ошибки, так что использовать нужно с умом, т.е. знать для чего это нужно и как работает.
Пришлось поднапрячь мозг, чтобы без документации написать численную функцию в качестве дополнения.
На картинке ошибка. Частота Найквиста - это половина частоты дискретизации, т.е. 0,5 * Fд. Поэтому на оси абсцисс граничная частота и есть частота Найквиста, т.е. 0,5 - она нормирована по отношению к частоте дискретизации.
Можно кстати профильтровать суперпозицию из нескольких синусоид и посмотреть как работает дискретный фильтр.
Отредактировано пользователем 10 июля 2011 г. 5:42:43(UTC)
| Причина: Не указана
Россия навсегда! Вячеслав Мезенцев
1 пользователь поблагодарил уни за этот пост.
smath оставлено 09.11.2010(UTC)
Статус: Advanced Member
Группы: Registered
Зарегистрирован: 02.06.2009(UTC)
Сообщений: 346
Сказал «Спасибо»: 50 раз Поблагодарили: 156 раз в 105 постах
Лабораторная работа по цифровой фильтрации. Построен фильтр НЧ 27 порядка. Работа фильтра показана на примере фильтрации суперпозиции 10 синусоид одинаковой амплитуды (одной там не видно, т.к. её "скушало" преобразование: 10-я гармоника совпадает с верхней частотой спекта).
Правда спроектированный фильтр не обеспечивает подавления -40 дБ в полосе задерживания. Для этого нужен фильтр большего порядка. Оценка производилась по эмпирической формуле (не знаю откуда взята, судить не могу).
Неравномерность в полосе пропускания, пожалуй, большевата, но для вычислительного эксперимента пойдёт.
Здесь используется вещественная свёртка - функция al_convr1d(). Она входит в состав плагина AlgLib. С её помощью считается реакция фильтра. Пока прикладывать документ и плагин не буду, т.к. последний ещё не доделан.
На последних двух картинках частоты приведены на оси в Гц.
Для наглядности показан спектр сигнала на входе и выходе фильтра. Спектр выходного сигнала совмещен с АЧХ фильтра.
Отредактировано пользователем 8 ноября 2010 г. 18:14:33(UTC)
| Причина: Не указана
Россия навсегда! Вячеслав Мезенцев
Быстрый переход
Форум проекта SMath
SMath Studio
- Скачать SMath Studio
--- Архив
--- Ночные сборки
- Примеры использования
- Вопросы по использованию
- Предложения
- Ошибки
Математика
- Документы
- Алгоритмы
- Математические задачи
Общение
- О проекте
- Общение
Вы не можете создавать новые темы в этом форуме.
Вы не можете отвечать в этом форуме.
Вы не можете удалять Ваши сообщения в этом форуме.
Вы не можете редактировать Ваши сообщения в этом форуме.
Вы не можете создавать опросы в этом форуме.
Вы не можете голосовать в этом форуме.