数字滤波器的MATLAB设计实现资料.docVIP

数字滤波器的MATLAB设计实现 （1）设计题目 （2）主要指标和要求 （3）方案设计 （4）MATLAB程序 （5）仿真结果 （6）设计结论 （7）收获、体会和改进设计的建议 设计题目； 用FFT变换对信号Y(t)=sin(50*2*pi*t)+随机噪声(0=t=0.5)进行频谱分析。 用滤波器完全设计函数，设计IIR滤波器将信号中的正弦波信号从噪声中提取出来。 用窗函数法，设计FIR滤波器完成同样的信号提取。 主要指标和要求 数字信号分析 数字信号处理是用数字和符号的序列来表示信号，用数字计算机来处理这些序列，如对信号的滤波，增强信号的有用分量，削弱无用分量，估计信号的某些特征参数等。 离散傅里叶变换（DFT） 快速傅里叶变换（FFT） 用DFT对模拟信号进行谱分析 对模拟信号进行以T为周期的采样 对采样后的序列信号进行截断处理，使成为包含N个采样点的有限长序列 对连续的频域函数进行采样，在频域的一个周期内取N个样点 应用数字计算机对模拟信号进行傅里叶分析时，将时间函数和频率函数都修改成周期序列，用N个时间样本值和N个频率样本值分别表示时域波形和频域波形的一个周期。 快速傅里叶变换（FFT）是DFT的快速算法，特别是当点数N比较大时，在计算量上有明显优势。MATLAB中提供了专门的函数来进行FFT计算： Y=fft(x) Y=fft(x,n) Y=fft(x,[],dim) Y=fft(x,n,dim) 数字滤波器 数字信号处理过程中，所处理的信号往往含有噪声，从接收到的信号中消除或削弱噪声，提取有用信号的过程就是滤波，实现滤波功能的系统就是滤波器。按照时域特性来分，数字滤波器还可以分为无限冲激响应滤波器（IIR）和有限冲激响应滤波器（FIR） 在很多实际应用中，例如话音或音频信号处理中，数字滤波器常常被用来实现选频操作，因此，通常把频域幅度和相位响应作为性能指标，而在通带中，常常希望系统具有线性相位响应。通过FIR滤波器可以得到准确的线性相位。在IIR滤波器中，通带中的线性相位是不可能得到的，因此IIR滤波器的设计一般只考虑幅度指标。幅度指标可以按两种方式给出。第一种是绝对指标，它提供对幅度响应函数的要求，这些指标一般应用于FIR滤波器的设计。IIR滤波器设计则以一种不同的方式给出指标。第二种指标是相对指标，它以分贝（dB）值的形式提出要求， 其值定义： 方案设计 无限冲激响应滤波器（IIR） IIR滤波器有以下几个特点： 1．IIR数字滤波器的系统函数可以写成封闭函数的形式。 2．IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。 3．IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。 4．IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。 巴特沃斯IIR滤波器的设计 在MATLAB下，设计巴特沃斯IIR滤波器可使用butter函数。 Butter函数可设计低通、高通、带通和带阻的数字和模拟IIR滤波器，其特性为使通带内的幅度响应最大限度地平坦，但同时损失截止频率处的下降斜度。在期望通带平滑的情况下，可使用butter函数。 butter函数的用法为： [b,a]=butter(n,Wn,/ftype/) 其中n代表滤波器阶数，Wn代表滤波器的截止频率，这两个参数可使用buttord函数来确定。buttord函数可在给定滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n，同时给出对应的截止频率Wn。buttord函数的用法为： [n,Wn]= buttord(Wp,Ws,Rp,Rs) 其中Wp和Ws分别是通带和阻带的拐角频率（截止频率），其取值范围为0至1之间。当其值为1时代表采样频率的一半。Rp和Rs分别是通带和阻带区的波纹系数。 不同类型（高通、低通、带通和带阻）滤波器对应的Wp和Ws值遵循以下规则： 1．高通滤波器：Wp和Ws为一元矢量且WpWs； 2．低通滤波器：Wp和Ws为一元矢量且WpWs； 3．带通滤波器：Wp和Ws为二元矢量且WpWs，如Wp=[0.2,0.7],Ws=[0.1,0.8]; 4．带阻滤波器：Wp和Ws为二元矢量且WpWs，如Wp=[0.1,