버터 워스 저역 통과 필터의 주파수 범위 전체에 대해 위상 플롯이 생성되지 않습니다.

Question

버터 워스 저역 통과 필터를 설계하고 크기 및 위상 응답을 확인하고 싶습니다. 다음 코드를 시도했다.

n = 8;  % order of filter 
fs = 1000; % sampling frequency 
fc = 20;  % cutoff frequency 
[b a]= butter(n,fc/(fs/2),'low'); % design low pass filter 
freqz(b,a,1024,fs); 

나는 다음과 같은 플롯을 얻고있다 :

바와 같이, 모든 주파수 범위에 대한 위상 플롯을 제공하지 않습니다 볼 수 있습니다

(300 ~ 500 사이를 Hz의 주파수에 대한 위상 플롯이 없다). 이 문제의 원인을 이해하고이를 수정하는 방법을 알려주십시오.

Answer 1

예, 버그 인 것 같습니다. 하지만 '손으로'동일한 음모를 만들 수 있습니다.

[H,w] = freqz(b,a,1024,fs); 
plot(w, unwrap(angle(H)));