Exponenta Banner

freqs

Частотная характеристика аналоговых фильтров

свернуть все на странице

Синтаксис

h = freqs(b,a,w)
[h,wout] = freqs(b,a,n)
freqs(___)

Описание

пример

h = freqs(b,a,w) возвращает комплексную частотную характеристику аналогового фильтра, заданного векторами коэффициентов b и a, оцененный на угловых частотах w.

[h,wout] = freqs(b,a,n) использование n частота указывает, чтобы вычислить h и возвращает соответствующие угловые частоты в wout.

пример

freqs(___) без выходных аргументов строит величину и фазовые отклики как функции угловой частоты в окне текущей фигуры. Можно использовать этот синтаксис с любым из предыдущих входных синтаксисов.

Примеры

свернуть все

Скрипт Open Live Script

Найдите и изобразите частотную характеристику в виде графика передаточной функции

H(s)=0.2s2+0.3s+1s2+0.4s+1.

a = [1 0.4 1];
b = [0.2 0.3 1];
w = logspace(-1,1);

h = freqs(b,a,w);
mag = abs(h);
phase = angle(h);
phasedeg = phase*180/pi;

subplot(2,1,1)
loglog(w,mag)
grid on
xlabel('Frequency (rad/s)')
ylabel('Magnitude')

subplot(2,1,2)
semilogx(w,phasedeg)
grid on
xlabel('Frequency (rad/s)')
ylabel('Phase (degrees)')

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

Можно также сгенерировать графики путем вызова freqs без выходных аргументов.

figure
freqs(b,a,w)

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

Скрипт Open Live Script

Спроектируйте аналог 5-го порядка lowpass фильтр Бесселя с приблизительно постоянной групповой задержкой до 104 рад/с. Постройте частотную характеристику фильтра с помощью freqs.

[b,a] = besself(5,10000);   % Bessel analog filter design
freqs(b,a)                  % Plot frequency response

Figure contains 2 axes. Axes 1 contains an object of type line. Axes 2 contains an object of type line.

Входные параметры

свернуть все

Коэффициенты передаточной функции в виде векторов.

Пример: [b,a] = butter(5,50,'s') задает пятый порядок Фильтр Баттерворта с частотой среза 50 рад/секунда.

Типы данных: single | double

Угловые частоты в виде положительного вектора действительных чисел описываются в рад/секунда.

Пример: 2*pi*logspace(6,9) задает 50 логарифмически расположенных с интервалами угловых частот от 1 МГц (2π рад/секунда × 106) и 1 ГГц (2π рад/секунда × 109).

Типы данных: single | double

Количество оценки указывает в виде положительного целочисленного скаляра.

Типы данных: single | double

Выходные аргументы

свернуть все

Частотная характеристика, возвращенная как вектор.

Угловые частоты, в который h вычислен, возвращен как вектор.

Алгоритмы

freqs возвращает комплексную частотную характеристику аналогового фильтра, заданного b и a. Функция оценивает отношение полиномов Преобразования Лапласа

H(s)=B(s)A(s)=b(1)sn+b(2)sn−1+⋯+b(n+1)a(1)sm+a(2)sm−1+⋯+a(m+1)

вдоль мнимой оси на частоте указывает s = jω:

s = 1j*w;
h = polyval(b,s)./polyval(a,s);

Смотрите также

| | | | |

Представлено до R2006a

Документация Signal Processing Toolbox

Поддержка