matlab freqz函数的使用

freqz函数是MATLAB中用于绘制数字滤波器频率响应的函数。它可以绘制数字滤波器的幅频响应和相频响应。freqz函数的语法如下：

[h,w] = freqz(b,a,n,fs)

其中，b和a表示数字滤波器的系数，n表示要绘制的频率点数，fs表示数字滤波器的采样频率。freqz函数返回滤波器的频率响应h和对应的频率向量w。

下面是一个示例代码，演示如何使用freqz函数绘制数字低通滤波器的频率响应：

% 设计数字低通滤波器
fs = 1000;               % 采样频率
fc = 200;                % 截止频率
[b,a] = butter(4,fc/(fs/2));  % 4阶巴特沃斯低通滤波器

% 绘制频率响应
n = 512;                 % 频率点数
[h,w] = freqz(b,a,n,fs); % 计算频率响应
figure;                  % 创建新的图形窗口
subplot(2,1,1);          % 创建子图1并设置为第1行第1列
plot(w,abs(h));          % 绘制幅频响应
xlabel('Frequency (Hz)');ylabel('Magnitude');
title('Butterworth Lowpass Filter Frequency Response');
grid on;

subplot(2,1,2);          % 创建子图2并设置为第2行第1列
plot(w,angle(h));        % 绘制相频响应
xlabel('Frequency (Hz)');ylabel('Phase (rad)');
grid on;

运行上述代码，即可得到数字低通滤波器的幅频响应和相频响应的图像。