freqz
MATLAB提供了專門用于求離散系統(tǒng)頻響特性的函數(shù)freqz()顶霞，調(diào)用freqz()的格式有以下兩種：
1:[H,w]=freqz(B,A,N)
B和A分別為離散系統(tǒng)的系統(tǒng)函數(shù)分子肄程、分母多項(xiàng)式的系數(shù)向量，N為正整數(shù)选浑，返回量H則包含了離散系統(tǒng)頻響 在 0——pi范圍內(nèi)N個(gè)頻率等分點(diǎn)的值蓝厌，向量w則包含 范圍內(nèi)N個(gè)頻率等分點(diǎn)。調(diào)用中若N默認(rèn)古徒，默認(rèn)值為512拓提。
2:[H,w]=freqz(B,A,N,’whole’)
該調(diào)用格式將計(jì)算離散系統(tǒng)在0—pi范內(nèi)的N個(gè)頻率等分店的頻率響應(yīng)的值。因此隧膘，可以先調(diào)用freqz()函數(shù)計(jì)算系統(tǒng)的頻率響應(yīng)代态，然后利用abs()和angle()函數(shù)及plot()函數(shù)，即可繪制出系統(tǒng)在 或 范圍內(nèi)的頻響曲線疹吃。
例：繪制如下系統(tǒng)的頻響曲線
H(z)=(z-0.5)/z
MATLAB命令如下：
B=[1 -0.5];
A =[1 0];
[H,w]=freqz(B,A,400,'whole');
H是頻率響應(yīng)的幅度蹦疑，w是0---pi內(nèi)的400個(gè)點(diǎn)
Hf=abs(H);
Hx=angle(H);
clf
figure(1)
plot(w,Hf)
title('離散系統(tǒng)幅頻特性曲線')
figure(2)
plot(w,Hx)
title('離散系統(tǒng)相頻特性曲線')
幅頻特性.PNG
相頻特性.PNG
這樣畫出來的是線性的，要想獲得db格式的幅度萨驶，需要轉(zhuǎn)換 20*log10（Hf）
之后再畫就是db格式的
也可以直接用freqz(b,a,w)這樣就會(huì)畫出幅頻響應(yīng)和相頻響應(yīng)歉摧，幅頻響應(yīng)直接是db格式的幅度。
原文地址:http://zfblog.xyz/archives/230.html
