Частотная характеристика аналоговых фильтров
freqs(___) без выходных аргументов строит график величины и фазовых откликов как функций угловой частоты в окне текущего рисунка. Этот синтаксис можно использовать с любым из предыдущих входных синтаксисов.
Найти и построить график частотной характеристики передаточной функции
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)')
Можно также создать графики, вызвав freqs без выходных аргументов.
figure freqs(b,a,w)
Разработать аналоговый фильтр Бесселя 5-го порядка с приблизительно постоянной групповой задержкой до рад/с. Постройте график частотной характеристики фильтра с помощью freqs.
[b,a] = besself(5,10000); % Bessel analog filter design freqs(b,a) % Plot frequency response
freqs возвращает комплексную частотную характеристику аналогового фильтра, указанного b и a. Функция вычисляет отношение многочленов преобразования Лапласа
sm + a (2) sm−1+⋯+a (m + 1)
вдоль воображаемой оси в частотных точках s = jλ:
s = 1j*w; h = polyval(b,s)./polyval(a,s);