在MATLAB中，可以使用freqz函數來計算和繪制數字濾波器的頻率響應。
freqz函數的基本用法如下：
freqz(b,a,n,fs)
其中，b和a是數字濾波器的系數，n是頻率響應的點數，fs是采樣頻率。
具體步驟如下：
定義數字濾波器的系數b和a。
調用freqz函數來計算濾波器的頻率響應。可以指定繪制頻率響應的點數n和采樣頻率fs。
freqz函數返回兩個輸出參數，第一個是頻率響應的幅度響應，第二個是頻率響應的相位響應。
可以使用plot函數來繪制頻率響應。
下面是一個例子：
% 定義數字濾波器的系數
b = [0.1 0.2 0.3 0.2 0.1];
a = 1;
% 計算頻率響應
n = 1024; % 頻率響應的點數
fs = 1000; % 采樣頻率
[h, w] = freqz(b, a, n, fs);
% 繪制頻率響應
figure;
subplot(2,1,1);
plot(w, abs(h));
title('Amplitude Response');
xlabel('Frequency (Hz)');
ylabel('Magnitude');
subplot(2,1,2);
plot(w, angle(h));
title('Phase Response');
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
在上面的例子中，我們定義了一個5階的數字濾波器，然后使用freqz函數計算并繪制了該濾波器的頻率響應。注意，我們使用subplot函數來將幅度響應和相位響應分別繪制在兩個子圖中。
