freqs

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

Описание

пример

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

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

пример

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

Примеры

свернуть все

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

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.

Разработайте аналоговый lowpass фильтр Бесселя 5-го порядка с приблизительно постоянной групповой задержкой до 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)sn1++b(n+1)a(1)sm+a(2)sm1++a(m+1)

вдоль мнимой оси в частотных точках s = :

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

См. также

| | | | |

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте