Русскоязычный форум закрыт из-за отсутствия активности (доступен только для чтения).
Пожалуйста, пользуйтесь англоязычной его версией. Приносим извинения за неудобства
Добро пожаловать, Гость! Чтобы использовать все возможности Вход. Новые регистрации запрещены.

Уведомление

Icon
Error

Вход


Опции
К последнему сообщению К первому непрочитанному
Offline уни  
#1 Оставлено : 6 ноября 2010 г. 23:45:18(UTC)
уни


Статус: 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)  | Причина: Не указана

Россия навсегда!
Вячеслав Мезенцев
thanks 1 пользователь поблагодарил уни за этот пост.
smath оставлено 09.11.2010(UTC)

Wanna join the discussion?! Login to your Форум проекта SMath forum account. Новые регистрации запрещены.

Offline уни  
#2 Оставлено : 8 ноября 2010 г. 18:09:31(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)  | Причина: Не указана

Россия навсегда!
Вячеслав Мезенцев
Пользователи, просматривающие эту тему
Быстрый переход  
Вы не можете создавать новые темы в этом форуме.
Вы не можете отвечать в этом форуме.
Вы не можете удалять Ваши сообщения в этом форуме.
Вы не можете редактировать Ваши сообщения в этом форуме.
Вы не можете создавать опросы в этом форуме.
Вы не можете голосовать в этом форуме.