《数字信号处理》课件_第6章.ppt

综上所述，调用remez和remezord函数设计线性相位FIR滤波器，关键是根据设计指标求出remezord函数的调用参数f、m、rip和Fs。

下面给出由给定的各种滤波器设计指标确定remezord调用参数f、m和rip的公式，编程时直接套用即可。(1)低通滤波器设计指标。

逼近通带为［0，ωp］，通带最大衰减为αp(dB);逼近阻带为［ωs，π］，阻带最小衰减为αs(dB)。

remezord调用参数：，m=［1，0］，rip=［δ1，δ2］其中，f向量省去了起点频率0和终点频率1；δ1和δ2分别为通带和阻带波纹幅度，由式(6.5.2)计算可得，下面相同。(2)高通滤波器设计指标。

逼近通带为［ωp，π］，通带最大衰减为αp(dB);逼近阻带为［0，ωs］，阻带最小衰减为αs(dB)。

remezord调用参数：，m=［0，1］，rip=［δ2，δ1］(3)带通滤波器设计指标。

逼近通带为［ωpl，ωpu］，通带最大衰减为αp(dB);逼近阻带为［0，ωsl］、［ωsu，π］，阻带最小衰减为αs(dB)。

remezord调用参数：(4)带阻滤波器设计指标。

逼近通带为［0，ωpl］、［ωpu，π］，通带最大衰减为αp(dB);逼近阻带为［ωsl，ωsu］，阻带最小衰减为αs(dB)。

remezord调用参数：【例6.5.1】利用remez交替算法，设计一个线性相位低通FIR数字滤波器，其指标为：通带边界频率fp=800Hz，阻带边界fs=1000Hz，通带最大衰减αp=0.5dB，阻带最小衰

减αs=40dB，采样频率Fs=4000Hz。解在Matlab中可以用remezord和remez两个函数设计，程序如下：

clc;clearall;

fp=800;fs=1000；%设置滤波器参数

Fs=4000; %设置采样频率

alphap=0.5;alphas=40; %设置滤波器参数

f=［fp，fs］; 

m=［10］;delta1=(10^(alphap/20)-1)/(10^(alphap/20)+1);

%求通带波纹幅度delta1

delta2=10^(-alphas/20); %求阻带波纹幅度delta2

delta=［delta1，delta2］;

［L，fpts，mag，wt］=remezord(f，m，delta，Fs); %用remezord求参数

hn=remez(L，fpts，mag，wt);

%用remez求滤波器的脉冲响应h(n)

［H，w］=freqz(hn，1，512);

%用freqz求频率响应H

plot(w*2000/pi，20*log10(abs(H)));grid;

xlabel(′Frequency/Hz′)；

ylabel(′Amplitude/dB′)运行程序后，结果如图6.5.3所示，滤波器阶数为32，FIR数字低通滤波器在fp=800Hz处衰减为0.31dB左右，在fs=1000Hz处衰减为38.3dB左右，基本满足设计指标。

【例6.5.2】利用等波纹最佳逼近法设计带通滤波器，其指标为：①ωpl=0.32π，ωpu=0.6π;②ωsl=0.28π，ωsu=0.66π;③αp≤1dB，αs≥50dB。图6.5.3remez交替算法设计低通滤波器举例解逼近通带为［0.32π，0.6π］，通带最大衰减αp=1dB，逼近阻带为［0，0.28π］、［0.66π，π］,

阻带最小衰减αs=50dB。调用remezord和remez函数求解，

Matlab参考程序如下:

clc;clearall;

f=［0.28，0.32，0.6，0.66］；

m=［0，1，0］;

alphap=1;

alphas=50;delta1=10^(alphap/20-1)/(alphap/20+1);

delta2=10^(-alphas/20);

rip=［delta2，delta1，delta2］;

［L，fo，mo，w］=r