clc;
clear;
a = [1 0.4 1];
b = [0.2 0.3 1];
w = logspace(-1, 1);%生成从10的a次方到10的b次方之间按对数等分的n个元素的行向量。n如果省略,则默认值为50。自动挑选200个频率点来计算频率响应h
h = freqs(b,a,w);%根据系数向量计算返回模拟滤波器的复频域响应。freqs 计算在复平面虚轴上的频率响应h,角频率w确定了输入的实向量,因此必须包含至少一个频率点。
mag = abs(h);%注意这里取了绝对值
phase = angle(h);
subplot(2,1,1), loglog(w,mag) %x 轴和y 轴均为对数刻度
grid;
subplot(2,1,2), semilogx(w,phase)%x 轴为对数刻度,y 轴为线性刻度
grid;