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