Фильтр Баттерворта

@Arhel · Регистрация: 04.01.2012

Студворк — интернет-сервис помощи студентам

Доброго времени суток!

Подскажите, пожалуйста, как рассчитать коэффициенты полосового фильтра Баттерворта 8-го порядка методом билинейной трансформации? (Частота дискретизации = 3000 Гц, f1 = 450Гц, f2 = 750 Гц)

Пробовал использовать функцию fdatool, но, к сожалению, там нельзя выбрать какой-либо метод - билинейной трансформации или к примеру импульсной инвариации.

Можно, конечно, вручную..но лень =)

Спасибо!

@vital792 · 19.01.2012, 10:02

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
clear all;
close all;
 
Fs = 3000;
t = 0 : 1/Fs : 4095/Fs;
signal = sin(2*pi*100*t)+sin(2*pi*300*t)+sin(2*pi*600*t)+0.5*randn(1, length(t));
 
spec = fft(signal);
figure(1)
plot(abs(spec(1:floor(length(spec)/2))));
 
n = 8;
Wn = [450 750]/Fs;
 
% test
[b_t, a_t] = butter(n, Wn, 'bandpass');
signal1 = filter(b_t, a_t, signal);
spec1 = fft(signal1);
figure(2)
plot(abs(spec1(1:floor(length(spec1)/2))));
 
Wn = 2*pi*Wn;
[b, a] = butter(n, Wn, 'bandpass', 's');
 
% [bz, az] = bilinear(b, a, 2);
[bz, az] = impinvar(b, a, 2);
 
filtredSignal = filter(bz, az, signal);
 
specFiltred = fft(filtredSignal);
figure(1)
hold on;
plot(abs(specFiltred(1:floor(length(specFiltred)/2))), 'r');

на figure1 синий график - спектр исходного сигнала, на figure2 - спектр после фильтра баттерворта сгенерированного функцией butter, на первом красный график - спектр сигнала после фильтра восстановленного из аналогового функцией impinvar

@margadon · 21.12.2015, 11:56

я сделал бы как то так:

скажите, у вас в исходном сигнале частоты 100, 300 и 600 Гц, центральная частота полосы пропускания фильтра - 600 Гц , при этом он, судя по спектру на выходе фильтра, отфильтровывает частоту 300 Гц, т.е. вдвое меньшую.
Поясните пожалуйста, почему так, и наконец - какие единицы у вас в оси абсцис на графиках спектров? они примерно в 1.37 раз больше чем f сигнала. Что это за величина? Спасибо.

Новые блоги и статьи Все статьи Все блоги /
GitHub Copilot в JetBrains: модель по умолчанию, Fix и контроль MCP and_y87 12.10.2026 10 октября 2026 года GitHub обновил Copilot для JetBrains: добавил настройки стартовой модели, исправления из диагностики и управление автозапуском MCP. Поддержка IDE 2025. 1 прекращена — требуется. . .	День 5 из 28: Codex подсказывает следующий запрос and_y87 12.10.2026 9 октября 2026 года OpenAI запустила бету Composer predictions — подсказок следующего сообщения в текущем диалоге Codex. Это помощь с формулировкой запроса, а не автоматическое продолжение задачи. . . .	Запустил конкурс "тем и промптов для текстовых квестов созданных почти чисто ИИ" Adler 06.10.2026 Всем привет! За последние три-четыре дня я создал более 16 текстовых квестовых игр используя преимущественно по одному запросу к ИИ на игру. Мне так понравилось смотреть все ветки/ сцены во всех. . .	ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .
Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .	Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .	Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .

@Arhel 0 / 0 / 0 Регистрация: 04.01.2012 Сообщений: 15

	Фильтр Баттерворта 07.01.2012, 23:17. Показов 20162. Ответов 21 Метки нет (Все метки) Доброго времени суток! Подскажите, пожалуйста, как рассчитать коэффициенты полосового фильтра Баттерворта 8-го порядка методом билинейной трансформации? (Частота дискретизации = 3000 Гц, f1 = 450Гц, f2 = 750 Гц) Пробовал использовать функцию fdatool, но, к сожалению, там нельзя выбрать какой-либо метод - билинейной трансформации или к примеру импульсной инвариации. Можно, конечно, вручную..но лень =) Спасибо! 0