08 频率归一化与转换幻灯片.ppt

数字信号处理 8 频率归一化与滤波器转换 模拟滤波器频率转换 巴特沃斯滤波器设计 数字滤波器设计 实例分析 回顾 回顾 回顾 N又分两种情况： 偶数 奇数 回顾 回顾 N分两种情况： 偶数 奇数 回顾 设计任务 选择有限长 ，得到系统函数 ，使得： 幅频特性 满足技术指标要求； 相频特性 达到线性相位。 回顾 —— fir1 函数 MATLAB提供了fir1函数，以实现线性相位FIR滤波器。阅读P213 调用格式如下： hn=fir1(N, wc, ‘ftype’, window) 参数： N：阶数 wc：归一化的数字频率，0≤wc≤1。 ftype：滤波器类型，如高通、带阻等。当ftype=high时，设计高通FIR滤波器；当ftype=stop时，设计带阻FIR滤波器。注意，在设计高通和带阻滤波器时，阶数N只能取偶数(h(n)长度N+1为奇数)。不过，当用户将N设置为奇数时，fir1会自动对N加1。 window：应用的窗函数类型，默认为hamming窗。 回顾 —— fir1 函数 例如： hn=fir1(N,wc,bartlett(N+1)) %使用Bartlett窗设计 hn=fir1(N,wc,chebwin(N+1,R)) %使用Chebyshev窗设计 下面的实例代码： wc=29/66; N=32; %N=h(n)的长度减1 hn=fir1(N,wc,high) 例 2 —— 用函数设计滤波器 用窗函数法设计线性相位的FIR低通滤波器，实现对模拟信号采样后进行数字低通滤波，对模拟信号的滤波要求如下： 通带截止频率： =2 kHz 阻带截止频率： =3 kHz 阻带最小衰减： 40 dB 采样频率： =10 kHz 选择合适的窗函数及其长度，求出 ，并画出幅频特性衰减曲线和相频特性曲线。 例 2 —— 用函数设计滤波器 【解】 （1）根据题意确定相应的数字滤波器指标： 通带截止频率： 阻带截止频率： 阻带最小衰减： 过渡带宽度： 例 2 —— 用函数设计滤波器 （2）选择窗函数w(n)，计算窗函数长度N 由阻带衰减为40 dB知道，hanning和hamming窗函数都满足要求，这两种窗函数过渡带宽度相同。这里选hamming窗： hamming窗函数的过渡带宽度为： 所以应当满足： 解之得到： 例 2 —— 用函数设计滤波器 （3）确定希望逼近的理想低通滤波器截止频率 我们知道， 频率点 (该点为过渡带中心)处滤波器频响幅度为0.5(等价于衰减6dB)， 所以： 例 2 —— 用函数设计滤波器 hn=fir1(N-1,wc); % 求h(n) M=1024; % 以下为计算频率相应函数和绘图部分 hk=fft(hn,M); % 求h(n)的FFT变换 n=0:N-1; figure(1);subplot(1,2,1); stem(n,hn,.); % 绘制h(n)的火柴杆图 xlabel(n); ylabel( h(n)); grid on k=1:M/2+1; w=2*(0:M/2)/M; subplot(1,2,2); plot(w,20*log10(abs(hk(k)))); % 绘制h(n)的角频谱图(rad/s) axis([0,1,-100,5]); xlabel(ω/π); ylabel(20lg|Hg(ω)|); grid on 例 2 —— 用函数设计滤波器 所设计的滤波器如下： 例 2 —— 用函数设计滤波器 绘制频谱（Hz）： figure(2); [HK,F]=freqz(hn，[1]，1024，Fs); % 没有分母，相当于其系数为1 plot(F，20*log10(abs(HK))); xlabel(频率(Hz)); ylabel(幅度); 例 3 —— 信号提取 已知某信号由三个子信号叠加组成，表达式如下： 用窗函数法设计一个FIR滤波器，提取出第1个信号； 绘制滤波前、后的信号波形与频谱。 例 3 —— 信号提取 第1步：初始化 clear all; A1=3; % 信号1的幅度 A2=1.5; % 信号2的幅度 A3=5; % 信号3的幅度 f1=100; % 信号1的频率(Hz) f2=250; % 信号2的频率(Hz) f3=270; % 信号3的频率(Hz) Fs=1024; % 采样频率(Hz) P1=-30; % 信号1的相位(度) P2=90;