离散系统频响特性的计算

离散系统的频响特性是理解和分析数字信号处理系统性能的关键。在MATLAB中，可以使用freqz函数来计算离散系统的幅频和相频响应。

要使用freqz函数进行计算，需要提供以下参数：

- b: 分子的系数（多项式）数组；

- a: 分母的系数（多项式）数组；

- w: 抽样的频率点，可以是向量或标量。为了获得准确的频响特性，通常至少需要选择两个以上的频率点。

freqz函数返回的是系统的幅频响应和相频响应，分别用绝对值abs和角度angle表示。例如，下面的代码演示了如何使用freqz计算一个离散系统的频响特性：

b = [1 -0.5 0];  % 分子的系数数组
a = [1 -2 3];    % 分母的系数数组
w = linspace(-pi, pi, 512); % 抽样的频率点，从-π到+π，共512个点
H = freqz(b, a, w);  % 计算频响特性
figure; plot(w/pi*100, abs(H)); title(\"幅频响应\"); xlabel(\"频率（Hz）\"); ylabel(\"幅度（dB）\"); grid on;
clear H
figure; plot(w/pi*100, angle(H)*180/pi); title(\"相频响应\"); xlabel(\"频率（Hz）\"); ylabel(\"相位（度）\"); grid on;
clear w b a