Matlab中的DSP应用中的freqs函数

注意这是个S域分析…

返回量H则包含了离散系统频响在0——pi范围内N个频率等分点的值，向量w则包含范围内N个频率等分点。调用中若N默认，默认值为512。

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;