matlab在信号与图像处理中应用6.docVIP

数字滤波器设计与实现 在上一章中，已经对MATLAB信号处理的一些基本知识作了详细介绍，在此基础上，本章将重点介绍在MATLAB中如何按照指定要求设计并实现一个数字滤波器。这里所谓的滤波器设计是根据特定的任务需求确定一组滤波器系数的过程，而滤波器实现则涉及选取滤波器结构并将结构与滤波器系数相结合的过程。只有完成了这两个步骤，含噪声的信号才能够进行有效的滤波处理。. 数字滤波器可以分为两种类型，即无限长单位冲激响应（IIR）滤波器和有限长单位冲激（FIR）响应滤波器。在本章中，将分别介绍IIR和FIR数字滤波器设计的原理及方法，并运用大量实例说明在MATLAB中如何设计并实现数字滤波器。 熟悉数字滤波器的设计及实现步骤 掌握IIR滤波器的设计方法 掌握FIR滤波器的设计方法 了解IIR和FIR滤波器的区别与优缺点对比 设计两个Butterworth滤波器，并采用直接Ⅱ型结构实现对数据x的滤波，要求分别利用函数dfilt和filter。 设计两个10阶高通Butterworth滤波器，要求截止频率为200Hz和300Hz，给定采样频率1000Hz。 设计一个8阶带通和一个8阶带阻Butterworth滤波器，要求通带频率为200Hz到300Hz，阻带频率为200Hz到300Hz，给定采样频率1000Hz。 设计一个8阶低通ChebyshevⅠ型滤波器，一个8阶高通ChebyshevⅡ型滤波器，要求通带波纹1dB，阻带衰减20dB，截止频率300Hz，定采样频率1000Hz。 设计一个8阶低通椭圆滤波器，要求阻带衰减40dB，通带波纹3dB，截止频率300Hz，给定采样频率1000Hz。 设计一个低通ChebyshevⅠ型滤波器，要求阻带衰减至少50dB，通带波纹小于3dB，通带截止频率200Hz，阻带截止频率300Hz，给定采样频率1000Hz。 采用冲激响应不变法和双线性变换法将一个模拟低通滤波器变换为数字滤波器，采样频率10Hz。 利用冲激响应不变法和双线性变换法各设计一个Butterworth低通数字滤波器，要求通带范围200Hz以下，阻带范围300Hz以上，通带波纹Rp小于3dB，阻带衰减Rs大于20dB，采样频率1kHz。 利用直接设计法设计一个8阶低通滤波器，并绘出期望频响与实际频响的图形。 根据Butterworth滤波器的冲击响应，恢复滤波系数。 设计一个48阶的FIR带通滤波器，要求通带 设计一个25阶FIR低通滤波器，给定频率向量f = [0 0.6 0.6 1]，幅值向量m = [1 1 0 0]，绘出期望响应与实际响应的对比图。 利用频率取样法设计FIR高通滤波器，要求截止频率3/4πrad，采样间隔为π/6，并假设理想频率响应的一个过渡点。 设计一个24阶的反对称分段线性带通FIR滤波器，并将实际频率响应与期望频率响应绘图比较。这里期望频率响应为：f = [0 0.3 0.4 0.6 0.7 0.9]；a = [0 1 0 0 0.5 0.5]。 设计一17阶Parks-McClellan带通滤波器，要求期望频率响应f = [0 0.3 0.4 0.6 0.7 1]；a = [0 0 1 1 0 0]。 一混合正弦信号包含5Hz、15Hz、30Hz的三种频率信号，现要求设计适当的滤波器，保留15Hz频率分量信号,给定采样频率100Hz。 概述 实际中接收到的信号往往掺杂着噪声或干扰，如何从接收信号中消除噪声或干扰是信号处理领域的关键问题之一，即我们这里所要介绍的滤波问题。所谓滤波，就是从接收信号中滤除噪声和干扰从而提取有用信号的过程，相应的系统或装置被称为滤波器。按照滤波器输入输出信号的类型，滤波器可分为模拟滤波器和数字滤波器，这里我们主要研究的是后者。 数字滤波器的分类 数字滤波器的种类多种多样，相应的分类方法也不尽相同，可以从功能上、从结构上、以及从实现方法或设计方法上进行分类，甚至也可以从滤波域进行分类，即根据滤波器是在时间域、频率域还是在空间域实现，对应有时域滤波器、频域滤波器以及空域滤波器。不过，从总体上而言，可以将滤波器分为两大类：经典滤波器和现代滤波器。 经典滤波器 也称为选频滤波器，其特点是假定输入信号中的有用成分和希望滤除的成分各自分布在不同的频带，这样通过选择合适的频率进行滤波，就可达到滤除噪声的目的。可见，经典滤波器都是在频域内实现滤波处理的，因而可将其纳入频域滤波器范畴。 现代滤波器 比经典滤波器更为复杂的滤波器形式，适用于当有用信号和噪声成分的频带相互重叠、经典滤波器不能有效滤除噪声的情况。这类滤波器可以按照随机信号内部的一些统计分布规律，从干扰中最佳地提取有用信号。比较典型的现代滤波器包括：维纳滤波器、卡尔曼滤波器和自适应滤波器等。与经典滤波器只集中于频域处理不同，现代滤波器已