1、解差分方程y=filter(b,a,x) b,a 分别为输出、输入各阶差分的系数向量,x为输入序列。
2、绘制零极点图zplane(b,a) b,a分别为系统函数的分子分母系数向量。
3、绘制频率函数曲线
(1)绘制数字系统H(z)的幅频特性和相频特性曲线(只绘图)
freqz(b,a) b,a分别为系统函数的分子分母多项式系数向量。
(2)计算[0, π)区间512个频率点的频率函数值,并将计算结果存放在H中,512个频率点存放在w中:
[H,w]=freqz(b,a)
Plot(w,abs(H)) %绘制幅频特性曲线
Plot(w,20*log10(abs(H))) %绘制衰减响应曲线
Plot(w,angle(H)) %绘制相频特性曲线
4、IIR滤波器设计
(1)巴特沃斯型
[N,wc]=buttord(wp,ws,Rp,As,’s’)
[B,A]=butter(N,wc,’ftype’,’s’)
(2)椭圆型
[N,wpo]=ellipord(wp,ws,Rp,As,’s’)
[B,A]=ellip(N,Rp,As,wpo,’ftype’,’s’)
①
wp,ws,Rp,As分别是通带截止频率,阻带截止频率,通带允许的最大衰减,阻带允许的最小衰减。
②
ftype 缺省,wp,ws为一元向量,设计低通
ftype 缺省,wp,ws为二元向量,设计带通
ftype =high,wp,ws为一元向量,设计高通
ftype =stop,wp,ws为二元向量,设计带阻
③
设计数字滤波器,去掉’s’选项,且wp,ws为归一化的频率。
(3)模拟滤波器转换为数字滤波器
①脉冲响应不变法[b,a]=impinvar(B,A,fs) %B,A为模拟滤波器系统函数分子分母多项式系数向量,b,a为数字滤波器系统函数分子分母多项式系数向量,fs为采样频率。
②双线性变换法[b,a]=bilinear(B,A,fs) %B,A为模拟滤波器系统函数分子分母多项式系数向量,b,a为数字滤波器系统函数分子分母多项式系数向量,fs为采样频率。
5、FIR滤波器设计
h=fir1(M,wc,’ftype’,window)
(1)M为滤波器的阶数,wc为通带截止频率,取值与选项ftype有关;
(2)window为选定的窗函数名称,缺省时,选用‘hamming’。其他窗函数名称为boxcar(矩形),Bartlett(三角),hann(汉宁),blackman(布莱克曼)。
(3)输出参数h长度为M+1,表示用选定窗函数截断后的待设计滤波器的单位脉冲响应的M+1个值。
(4)ftype选项缺省,wc为标量,设计低通滤波器;
ftype=high,wc为标量,设计高通滤波器;
ftype选项缺省,wc=[wcl,wcu],设计带通滤波器;