多相滤波器的设计及仿真

摘要软件雷达是现代雷达的重要发展方向。

其中数字化中频接收系统是关键技术之一。

本论文以某雷达数字化接收系统为背景，主要研究其中的关键技术——多相滤波器的原理及设计与仿真。

为了更好的理解，本文同时对数字滤波器作了简单的介绍，包括滤波器的定义、分类与实现方式。

对无线电中的两个最基础的概念内插与抽取也做了介绍。

对多相滤波器的原理进行了详细的说明，从公式推导上进行了理论实现的方法与可能。

讨论了多相滤波如何实现信道化。

最后介绍了数字滤波器的设计步骤并实现了多相滤波器的MATLAB仿真。

关键词：多相滤波器、软件无线电、数字滤波器。

AbstractSoftware radar is an important developing direction of modern radar. Which digital intermediate frequency receiving system is the key technical one. This paper taking a radar system for background, digital receiving main research key technology - multiphase the principle of filter and design and simulation. In order to better understanding of digital filter, the paper also makes brief introduction, including the definition, classification and filter implementation. The two most basic to radio the concept interpolation and extract presented also. The principle of multiphase filter for a detailed instructions from the formula, the method to realize the theory with possible. Discusses how to realize the multiphase filtering channelized. At last, the paper introduces the design procedures of the digital filter and realize the multiphase filter MATLAB simulation.Keywords: multiphase filter, software radio, digital filters摘要I第一章问题的提出IV第二章数字滤波器概论IV引言 (IV)2.1、数字滤波器的定义 (V)2.2、数字滤波器的实现方式 (V)2.3、数字滤波器的分类 (V)2.4实际滤波器的设计指标 (VI)第三章多相滤波器的理论原理VII引言 (VII)3.1整数倍抽取 (VII)3.2整数倍内插 (VII)3.3抽取内插器的实时处理结构——多相滤波结构 ........... V III3.4频域抽取 (IX)3.5用加权函数展宽输出滤波器 (XI)3.6改变输出采样速率 (XII)3.7多相滤波器实现信道化 ................................................... X III 第四章多相滤波器的MATLAB仿真XIV引言 (XIV)4.1数字滤波器设计的理论基础 (XIV)4.2FIR窗函数设计法 .............................................................. X V4.3多相滤波器的MATLAB仿真 (XVI)第五章总结XIX参考文献XX致谢 XXI第一章问题的提出随着A/D(analog-to-digital)变换技术、DSP(digital signal processing)技术、FPGA(field programmable gatearray)技术及ASIC(application specific integrated circuit)等技术的发展，宽带数字化接收机正逐渐成为现代雷达、遥测及通信系统中必不可少的重要组成部分。

但不管什么类型的中频数字化接收机，其基本原理框图都可采用中频数字化接收机原理框图如图1所示。

其中多相滤波器是其中的关键技术，多相滤波可以利用抽取因子实现高效滤波，也起到抑制镜像干扰和邻道干扰的作用，因此多相滤波器的设计与研究就显得很重要，本文就在此基础上对多相滤波的原理和实现作了一些简单的讨论。

第二章数字滤波器概论引言出自滤波器与模拟滤波器都是一种选频器件，它对某些频率的信号给予很小的衰减，使具有这些频率分量的信号比较顺利地通过，而对其他不需要的频率分量的信号给予较大幅度衰减，尽可能阻止这些信号通过。

数字滤波器和模拟滤波器具有不同的滤波方法，数字滤波器是通过对输入信号进行数值运算的方法来实现滤波的，而模拟滤波器则用电阻、电容、电感及有源器件等构成电路对信号进行滤波。

因此，数字滤波器具有比模拟滤波器精度高、稳定性强、灵活度大、体积小、重量轻、不要求阻抗匹配及实现模拟滤波器无法实现的特殊滤波功能等优点。

数字滤波器要求输入、输出信号均为数字信号。

本章介绍数字滤波器的定义、分类及实际滤波器的设计指标。

2.1、数字滤波器的定义数字滤波器（Digital Filter ）通常是指一个用有限精度算法实现的离散线性时不变系统。

因此它具有线性时不变系统的所有特性。

通常用的数字滤波器一般属于选频滤波器。

假设数字滤波器的频率响应()j H e ϖ用下式表示：式中，()j H e ϖ称为滤波器幅频响应；()θϖ称为滤波器相频响应。

幅频响应表示信号通过该滤波器后各频率成分的衰减情况，而相频响应反映各频率成分通过滤波器后在时间上的延时情况。

因此，即使两个滤波器幅频响应相同，只要相频响应不同，对应相同的输入，滤波器的输出信号波形也是不一样的。

滤波器的特性最容易通过它的幅频响应的形状来描述。

滤波器在某个频率的幅度增益决定了滤波器对此频率输入的放大因子，增益可任意取值。

增益高的频率范围，信号可以通过，称之为滤波器的通带；相反增益低的频率范围，滤波器对信号有衰减或阻塞作用，称之为滤波器的阻带。

例如低通滤波器使低频成分通过，阻碍高频成分；高通滤波器则相反，使高频成分通过，阻碍低频成分。

理想滤波器的幅频响应是矩形，即通带的增益为1，阻带的增益为0，然而这种理2.2、数字滤波器的实现方式数字滤波器的实现方式一般可以分为两种，即软件实现和硬件实现。

软件实现指的是在通用计算机上执行滤波程序。

这种方法灵活，但一般不能完成实时处理。

硬件实现指的是在单片机、FPGA 或DSP 芯片上实现，由于硬件运算速度快，可以实现实时处理，因此在实际系统中经常用硬件来实现各种数字滤波器。

2.3、数字滤波器的分类数字滤波器按照不同的分类方法，有许多种分类，但总体可以分为两大类。

()()()j j j H e H e e θϖϖϖ=)(|ωj e H 通带阻带一类称为经典滤波器，即一般的线性系统滤波器。

另一类即所谓的现代滤波器。

现代滤波器的理论简历在随机信号处理的理论基础上，它利用了随机信号内部的统计特性对信号进行滤波，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等，在此不做讨论。

经典滤波器的分类可以从滤波功能和实现的网络结构或者单位脉冲响应来划分。

从滤波功能上分类，和模拟滤波器一样，可以分为低通、高通、带通和带阻等滤波器。

需注意的是数字滤波器的频率响应()j H e ϖ都是以2π为周期的，滤波器的低通频处于2π的整数倍处，而高通频带处于π的奇数倍附近，这一点和模拟滤波器是有区别的。

从实现的网络结构或者单位脉冲响应分类，可以分成无限脉冲响应（IIR ）滤波器和有限脉冲响应（FIR ）滤波器。

系统函数如下：()011M r rr Nkk k b z H z a z -=-==-∑∑ ()()10N n n H z h n z --==∑2.4实际滤波器的设计指标理想滤波器的脉冲响应为非因果且无限长序列。

一次它不能通过时移来转变为因果系统。

另外，无限长脉冲响应不能直接转换为非递归差分方程。

简单的方法就是把理想脉冲响应两边响应值很小的采样点截去，将脉冲响应变为有限长，再进行时移得到因果系统，使得脉冲响应所描述的滤波器可用。

截短对滤波器的影响：截短后，滤波器幅频响应曲线不再是理想矩形。

通带不再平坦，有过渡带。

同时阻带衰减不再为零。

脉冲响应保留的采样点越多，即滤波器阶数越高，滤波器形状越接近理Wp 和ws 分别称为通带截止频率和阻带截止频率。

1πp ω11δ+11δ-c ω2δ0sω通带阻带过渡带)(max A )(min A )(s A |)(|ωj e H参数1δ定义了通带波纹，即滤波器通带内偏离单位增益的最大值。

参数2δ定义了阻带波纹，即滤波器阻带内偏离零增益的最大值。

参数t B 定义了过渡带宽度，即阻带下限和通带上限之间的距离，或：t s pB ϖϖ=-。

过渡带一般是单调下降的。

第三章 多相滤波器的理论原理引言多相滤波过程式按照相位均匀划分把数字滤波器的系统函数H(z)分解成若干个具有不同相位的组，形成多个分支，在分支上实现滤波。

采用多相滤波结构，可利用多个阶数较低的滤波来实现原本阶数较高的滤波，而且每个分支滤波器处理的数据率仅为原数据率的I/D ，这为工程上高速实时信号处理提供了实现途径。

本章介绍了多相滤波的相关知识及FIR 数字滤波器设计的一般步骤。

3.1整数倍抽取所谓整数倍抽取是指把原采样序列x(n)每隔（D-1个数据取一个，以形成一个新序列x(M)，即：()()D x m x mD =式中，D 为正整数，抽取过程如图2.16所示，抽取器用符号表示则如图2.17所示。

很显然如果x(n)序列的采样率为fs ，则其无模糊带宽为fs/2.当以D 倍抽取率对x(n)进行抽取后得到的抽取序列x(m)之取样率为fs~D ，其无模糊带宽为fs/(2D)，当x(n)含有大于fs/(2D)的频率分量时，x(m)就必然产生频谱混叠，导致从x(m)中无法恢复x(n)中小于fs/(2D)的频率分量信号。

此处不予证明（证明过程见软件无线电原理与应用2.3.1整数倍抽取。

）由此可以得出一个完整的D 倍抽取器结构如图2.20所示。

途中H()为其带宽小于 的低通滤波器。