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