《基于FPGA的FIR数字滤波器的设计.doc》由会员分享，可在线阅读，更多相关《基于FPGA的FIR数字滤波器的设计.doc（51页珍藏版）》请在三一办公上搜索。
1、基于FPGA的FIR数字滤波器的设计The Design Of FIR Digital Filter Based On FPGA系（院）名称：电子信息与电气工程系 QQ 号： 309810851 目 录摘要IAbstractII引言.1第一章 FIR数字滤波器设计.21.1 FIR数字滤波器.21.1.1数字滤波器简介.21.1.2 FIR数字滤波器.31.1.3 FIR数字滤波器的结构41.2 FIR数字滤波器设计方法.61.2.1时窗函数法61.2.2频率采样法81.2.3等波纹最佳逼近法101.2.4三种设计方法的比较.11第二章 基于MATLAB的FIR数字滤波器的设计与仿真122.1
2、 MATLAB软件122.2使用MATLAB函数进行FIR数字滤波器的设计.122.2.1时窗函数设计法.122.2.2频率采样设计法.132.2.3等波纹最佳逼近设计法142.3使用FDATool设计FIR数字滤波器.152.4使用Simulink进行仿真.16第三章 EDA技术与分布式算法.193.1 EDA技术.193.1.1 EDA技术简介.193.1.2 QUARTUS II软件介绍.193.1.3 FPGA介绍.203.2 FPGA设计方法.223.3分布式算法.233.3.1分布式算法基础.233.3.2并行的分布式算法243.3.3串行的分布式算法24第四章 15阶并行FIR数
3、字滤波器的设计与仿真.264.1顶层文件的设计.264.2并行FIR各个模块的设计284.2.1输入数据预处理模块.284.2.3预相加模块.304.2.4查找表模块.324.2.5加法器模块.354.2.6移位累加器模块.374.2.7输出处理模块.384.3顶层综合与仿真.384.4数据误差分析.41结论.43致谢.44参考文献.45基于FPGA的FIR数字滤波器的设计摘要：线性相位是FIR数字滤波器设计中的一个重要条件。线性相位FIR滤波器的设计方法主要包括时窗函数法、频率采样法、等波纹最佳逼近法。本文基于MATLAB软件，分别用这三种方法设计FIR低通滤波器，并进行比较分析。为了验证F
4、IR滤波器的性能，本文使用MATLAB自带的FDATool工具设计等波纹法15阶低通滤波器，在Simulink里进行仿真，最后得到预期的结果。设计好15阶FIR低通滤波器的系数后，本文就FIR低通滤波器的FPGA实现进行了研究。EDA技术的发展和大规模可编程逻辑器件的应用，使得数字信号处理借助于FPGA实现变为可能。为了实现这一目标，本文讨论了FPGA器件的结构、特点和分布式算法实现FIR滤波器的可行性，然后给出了基于并行分布式算法的FIR滤波器的VHDL描述，最后借助于EDA软件QUARTUSII进行了综合和仿真，并取得了成功。关键词：分布式算法 FIR数字滤波器 线性相位 The Desi
5、gn Of FIR Digital Filter Based On FPGAAbstract：Linear phase is one of the important conditions of FIR digital filter design. The design method of the linear phase FIR filter mainly includes window function method, frequency sampling method and the optimal approximation corrugated . This paper based
6、on the MATLAB software with the three methods designs FIR low-pass filters and analysis the results. In order to verify the performance of FIR filter, this paper use FDATool owned by MATLAB software to design low-pass filter of 15 order with corrugated method and convent it into a model for the simu
7、lation later carried out. The final simulation gets the expected results. Finishing design 15 order FIR low-pass filter coefficients, this paper will study FIR low-pass filter on FPGA realizing. With the development of EDA technology and the wide application of large-scale programmable logic devices
8、, the application of digital signal processing on FPGA realizing become possible. In order to achieve this goal, this paper discusses the structure and characteristics of FPGA device and the feasibility of FIR filter design based on distributed algorithm, and propose the VHDL code based on the paral
9、lel distributed FIR filter. Finally we use EDA software called QUARTUSII for the synthesizing and simulation , and achieve a success.Key words：Distributed Algorithm; FIR Digital Filter; Liner Phrase引 言近年来，随着数字信号处理技术的发展和数字集成电路的广泛应用，越来越多的工程领域在将模拟信号转为数字信号后，面临的首要任务即是滤除干扰信号，滤波的好坏直接决定着系统能不能得到有用的信号，因此滤波在数字
10、信号处理中处于基础性的地位。同时，计算机技术的发展为数字滤波器的设计与仿真提供了便利。例如，使用MATLAB软件即可在很短的时间内就能设计性能很高的数字滤波器，并且还可以进行模拟仿真，这使得我们可以在硬件实现以前对系统的性能可以有个初步的感性认识。此外，大规模可编程逻辑器件的应用，为数字滤波器的硬件实现提供了新的方式，在设计合理的情况下，能提高数字滤波器的工作速度，并且还能节省资源和功耗。基于此，本文将充分利用便捷的EDA技术，采用MATLAB进行FIR数字滤波器的设计和仿真，并研究使用VHDL进行硬件描述，然后在FPGA上开发的可行性。通过本文，可以感受到计算机技术带给数字滤波系统设计的快捷
11、与便利。第一章 FIR数字滤波器设计1.1 FIR数字滤波器数字滤波器在数字信号处理中属于预处理的部分，因而起着基础性的作用，数字滤波器包括IIR和FIR数字滤波器。数字滤波器具有精度高、稳定性好、灵活性强、不要求阻抗匹配，易于修改等特点。下面将首先介绍一下数字滤波器，然后重点讨论FIR数字滤波器的设计原理和结构。1.1.1 数字滤波器简介一个简单的数字滤波系统如图1-1所示。图中，x(t)为模拟信号，经过A/D转换器后变为一个有着先后顺序的数字序列x(n)。然后x(n)通过数字滤波系统H(z)，即得到数字滤波器的输出y(n)。H(z)为该数字滤波系统的单位脉冲响应h(n)的Z变换，即： (1
12、-1)若h(n)为无限长序列，则得到的数字滤波器为IIR数字滤波器，又称递归滤波器；反之，若h(n)为有限长序列，则得到的数字滤波器为FIR滤波器，也称非递归滤波器。一个线形时不变因果滤波器可表示为： (1-2)其中N为h(n)的长度，即滤波器的长度。滤波器的阶数为N-1。图1-1 数字滤波系统 IIR滤波器主要是基于对模拟滤波器如巴特沃斯滤波器、椭圆滤波器等的幅频响应进行逼近，而其相频响应是非线性的。与IIR滤波器不同，FIR滤波器可以把相位特性设计成线性。这使得FIR数字滤波器在信号无失真传输、数据通信、图像传输与处理、语音信号处理等有线性相位要求的领域应用广泛。FIR滤波器的优点是软硬件
13、实现结构简单，不用考虑系统的稳定性问题；缺点是实现较高性能的频率响应需要较高的滤波器阶数。1.1.2 FIR数字滤波器由于线性相位FIR滤波器的广泛应用，本文将首先讨论一下如何实现线性相位FIR滤波器。由式（1-2）将代入，得到： (1-3)为复数，因此可以表示成如下形式： (1-4)式中，为实函数，称为幅度特性函数，为相位特性函数。 当满足: (1-5)为初始相位，为常数。该滤波器具有广义的线性相位，将（1-5）式变换为： (1-6)由（1-6）式可知， 不为常量。这样，当不同频率的信号通过该滤波器时，便会产生相位的失真。不过通常较小，相比来说小得多，由引起的相位失真可忽略不计。 在（1-5
14、）式中，当时，为常量，这时滤波器有严格的线性相位，即对于不同频率的信号，通过该滤波器都有恒定的延迟，而不产生相位的失真。 将式（1-3），（1-4），（1-5）联立，可得 (1-7a) (1-7b)将（1-7a）式除以(1-7b）式，消去 ，得到： =0 (1-8)当时，式（1-8）变为：=0 (1-9)观察可知，若关于求和区间中心奇对称，则（1-9）式成立。由于关于奇对称，令关于偶对称，则满足关于求和区间中心奇对称的要求，即（1-9）式成立。 当时，式（1-8）变为: (1-10)若关于奇对称，关于偶对称，则（1-10）式成立。 基于的对称不同和长度N的奇偶区别，线性相位FIR数字滤波器的幅
15、度频率特性有所不同，因而所实现的滤波器的功能有所不同，具体如表1-1所示。第三列为能够实现的滤波器的性能。表1-1 四种类型的线性相位滤波器I型N为奇数低通、带通、高通、带阻II型N为偶数低通、带通III型N为奇数带通IV型N为偶数带通、高通1.1.3 FIR数字滤波器的结构 根据FIR数字滤波器实现算法的不同，可以把FIR滤波器的结构划分为直接型、级联型、频率采样型和快速卷积型四种基本形式。本文主要讨论前两种结构。A 直接型结构由式子（1-2）可直接画出FIR数字滤波器的直接型结构，如图1-2所示。对于直接型结构来说，一个长度为N的FIR滤波器，每产生一个输出数据，要经过N次乘法，N-1次加
16、法。对于使用FPGA开发FIR数字滤波器，这样的结果显然不令人满意。图1-2 直接型结构于是本文做了一下改进，这种改进是基于线性相位的FIR数字滤波器的。以严格线性相位，N为偶数的FIR滤波器为例，如图1-3所示。图1-3 直接型的改进由于关于对称，我们可以将经过延时环节的位置关于对称的数据预先相加，然后可以再乘以相应的滤波器系数进行累加得到最终的输出结果。这样，每产生一个输出，经过次乘法，次加法，比原来减少次乘法。B 级联型结构对式子（1-2）进行因式分解，并将零点共轭的因式放在一起，这样产生了若干个一阶子式和二阶子式，将一阶子式看作二阶子式的一个特例，则系统函数可以表示为： (1-11)F
17、IR数字滤波器的级联型结构如图1-4所示。图1-4 级联型结构从图中可以看出，级联型结构每产生一个输出，需要次乘法，次加法。级联型的最大特点是可以分别独立调整每个子系统的零点之值，当需要精确控制滤波器的零点位置时，往往采用这种结构。1.2 FIR数字滤波器设计方法 FIR数字滤波器的设计方法主要有三种：时窗函数法、频率采样法和等波纹设计法。下面将分别讨论一下这三种方法的设计原理，然后给出设计步骤，最后在比较一下这三种方法的优缺点。1.2.1 时窗函数法 一个理想低通滤波器的幅频特性如图1-5所示。的表达式为 （1-11）图1-5 理想低通滤波器频谱图对进行离散时间傅里叶逆变换，得 （1-12）
18、是一个关于对称的无限长序列，是系统的单位脉冲响应，这说明系统是非因果的，在物理上不可实现。为了能够得到一个可实现的系统，可以将截取一部分，并顺序右移，使之成为一个因果的有限长序列。这种方法就好比在时域打开一个窗口一样，因而称为时窗法。可以预见，使用这种方法截取的序列越长，对理想的幅频特性逼近的越好。常见的时窗函数有矩形窗、汉宁窗、汉明窗、布莱克曼窗、凯塞窗等。矩形窗的表达式为：，为整数。当用矩形窗设计FIR滤波器时，所产生的频率响应，与理想滤波器的频率特性相比，在理想特性不连续点处附近形成了过渡带，并在过渡带两侧形成持续时间很长，逐渐衰减的波纹，即通带内产生了波动，而阻带内产生了余振，这种现象
19、称之为吉布斯(Gibbs)效应。吉布斯效应直接影响到滤波器的性能，因为通带内的波动会影响到滤波器的平稳性，阻带内的波动则影响阻带最小衰减。为了克服吉布斯效应，可以从两个方面着手：一是在频域，避免理想滤波器频谱中出现的跃变现象，把它改造成一条连续光滑的曲线，所用的方法是镶边法；二是在时间域，对截尾函数进行改造，也即设计出好的时窗函数。一个好的时窗函数的要求：A.主瓣宽度尽可能地小。B.旁瓣水平（振幅或能量）相对于主瓣来说也尽可能地小。但是这两个标准之间彼此是有矛盾的，即主瓣宽度越大，旁瓣水平才可能越低。因此实际上，我们只能在这两个标准之间作一权衡，针对具体问题，找出一个适当的时窗函数。五种窗函数
20、的主要参数如表1-2所示。应当说明的是，用时窗函数法设计的FIR滤波器通带波纹幅度近似等于阻带波纹幅度。一般阻带最小衰减达到40dB以上，则通带最大衰减就小于0.1dB。最后，总结一下利用时窗函数设计FIR数字滤波器的步骤：A. 根据所需设计的数字滤波器类型（低通、高通、带通、带阻），确定线性相位数字滤波器类型。B. 根据滤波器阻带衰减，选择窗函数的类型，根据过渡带宽度确定时窗函数的长度,并根据线性相位条件进行修正。C 确定理想数字滤波器的频率响应函数，其中为幅度特性函数，为相位特性函数。D 计算理想滤波器的单位脉冲响应，即。E 加窗得到设计结果，即。表1-2 五种窗函数参数表窗函数主瓣宽度/
21、rad旁瓣水平/db过渡带宽度/rad阻带最小衰减/db矩形窗-1321汉宁窗-3144汉明窗-4153布莱克曼窗-5774凯塞窗（=7.865）-57801.2.2 频率采样法 频率采样法是一种基于频率域抽样来逼近所要设计的FIR滤波器的频率特性的一种方法。频率设计法的基本设计流程如图1-6所示。 图1-6 频率设计法流程为了保证具有线性相位条件，其单位采样响应是实序列，且满足条件。对于偶对称的情况来说，线性相位条件如下： (1-13) (1-14)令代入式（1-13）、（1-14），得到 (1-15)式中， ，N为奇数 ，N为偶数。 频率响应函数在各采样点上等于，而采样点之间的值是由各采样
22、值之间的内插函数叠加形成，因而有一定的逼近误差，误差的大小取决于理想频率响应的曲线形状。采样点之间的理想频率特性变化越陡，则内插值与理想值之间的误差就越大，因而在理想频率响应特性的不连续点附近，就会产生肩峰和波纹；反之，理想频率响应特性变化越平缓，则内插值越接近理想值，逼近误差小。因此，可以在理想频率下响应特性的通带与阻带之间设置过渡带，从而减小逼近误差。过渡带的采样点个数m与所设计的滤波器的阻带最小衰减有关，具体如表2-3所示。此外过渡带的带宽与采样点数m之间的关系为： (1-16)表1-3 过渡带采样点数与阻带最小衰减的关系m0123/dB20445465758595综上，频率采样法的设计
23、步骤为：A. 根据阻带最小衰减，确定过渡带的样点数m。B. 根据过渡带宽度的要求，估算滤波器的长度。C. 构造希望逼近的频率响应函数，一般为具有理想频率特性，并满足线性相位的要求。D. 对进行频域等间隔N点采样，得到。E. 对进行N点IDFT，得到。1.2.3 等波纹最佳逼近法等波纹最佳逼近法是基于最大误差最小化的设计原则。最大误差最小化能使幅度误差在整个逼近频段上均匀分布，即所设计的FIR数字滤波器的幅度特性在通带和阻带范围内是等波纹的，且可以分别独立控制通带和阻带的波纹幅度，并且在滤波器长度给定的情况下，加权的幅度波纹误差最小。 定义加权幅度误差函数为 (1-17)，式中为幅度误差加权函数
24、，用来控制不同频带的幅度逼近误差。一般地，在要求逼近精度高的频带，取值大，而要求逼近误差精度低的频带，取值小。设计过程中是由设计者根据通带最大衰减和阻带最小衰减的指标要求取定的已知函数。对于FIR数字低通滤波器，常取: (1-18)式中，和分别为滤波器设计指标中通带和阻带的振荡波纹幅度，k为正的系数。 滤波器的通带最大衰减和阻带最小衰减与通带和阻带的振荡波纹幅度和的换算关系为： (1-19) (1-20)等波纹最佳逼近法的设计，在于找到滤波器的系数向量，使得在通带和阻带频带内的最大绝对值幅度误差为最小。 帕克斯-麦克莱伦采用基于交替定理的雷米兹交替算法，通过逐次迭代逼近的运算求得滤波器的系数向
25、量，从而实现等波纹最佳逼近法。 一般情况下，估计滤波器长度N的凯塞经验公式为： (1-21) 综上，用等波纹最佳逼近法设计FIR数字滤波器的步骤为：A. 根据滤波器的设计指标要求：边界频率、通带最大衰减、阻带最小衰减等，估计滤波器的长度，并确定幅度误差加权函数。B. 采用雷米兹交替算法，获得所设计的滤波器的单位脉冲响应。1.2.4 三种设计方法的比较时窗函数法和频率采样法，设计简单方便，易于实现，缺点是不能够精确指定滤波器的边界频率和，并且不能够分别控制通带和阻带波纹。等波纹法能够分别独立控制通带和阻带波纹，确定边界频率，并且所设计的滤波器的阶数与前两种方法相比要少。第二章 基于MATLAB的
26、FIR数字滤波器的设计与仿真 2.1 MATLAB软件 MATLAB是由MATRIX和LABORATORY两个英文单词的前三个字母组合而成，意为矩阵实验室。它有着强大的数值分析、符号运算、信号处理和图形显示的功能，能够满足绝大部分的计算要求。MATLAB程序语言语法结构和C语言相似，简单易用，可移植性好，易于开发。MATLAB内嵌了许多科学研究和工程应用领域的工具箱（Tool box），诸如优化算法、神经网络、信号处理、图像处理、嵌入式系统开发、DSP与通讯、电力系统仿真等，使用户可以直接使用工具箱进行学习、应用和研究而不需要自己编写代码。 MATLAB内嵌了一个优秀的仿真软件Simulink
27、，主要用来实现对工程问题的建模和动态仿真。Simulink应用了模块化设计和系统级仿真的思想，使得建模仿真如同搭积木一样简单。Simulink可以应用于动力系统、信号控制、通信系统设计、金融财会及生物医学等多个领域。Simulink还可以与MATLAB中的DSP工具箱、信号处理工具箱以及通信工具箱等联合使用，进而实现软硬件的接口，从而成为实用的控制软件。2.2 使用MATLAB函数进行FIR数字滤波器的设计设计指标：归一化通带截止频率：0.2归一化阻带截止频率：0.4通带最大衰减：1db阻带最小衰减：47d类型：FIR低通2.2.1 时窗函数设计法根据阻带最小衰减，查看表1-2，选择汉明窗，过
28、渡带带宽为。MATLAB程序如下：clear all;Wp=0.2*pi;Ws=0.4*pi;B=Ws-Wp;N=ceil(6.6*pi/B)+1;n=0:1:N-1;a=(N-1)/2;Wc=(Ws+Wp)/2;hd=sin(Wc*(n-a+eps)./(pi*(n-a+eps);w_ham=(hamming(N);hn=hd.*w_ham;freqz(hn,1);求得滤波器的长度为34，可知该滤波器类型为II型低通。频率响应如图2-1所示，实际的归一化通带截止频率为0.247，阻带最小衰减为50db，滤波器线性相位。图2-1 使用窗函数法设计的FIR低通的频率响应2.2.2 频率采样设计法
29、根据阻带最小衰减，可确定m=1。过渡带带宽B=Ws-Wp=0.2,由N=20，取N=21。令过渡带样值为0.3，Wp位于n=3处。MATLAB程序如下：N=21;T1=0.3;a=(N-1)/2;n=0:(N-1);wl=(2*pi*n/N);Hrs=1,1,T1,zeros(1,15),T1,1,1;Hdr=1,1,0,0;Wd=0,0.2,0.2,1;k1=0:floor(N-1)/2);k2=(floor(N-1)/2+1):(N-1);angH=-a*2*pi*k1/N,a*2*pi*k2/N;H=Hrs.*exp(j*angH);h=real(ifft(H,N);freqz(h,1)
30、;频率响应如图2-2所示。实际的归一化通带截止频率为0.15，阻带最小衰减为50db。滤波器类型为I型低通。图2-2使用频率采样法设计FIR低通滤波器的频率响应2.2.3 等波纹最佳逼近设计法等波纹最佳逼近法的MATLAB程序如下：wp=0.2*pi;ws=0.4*pi;Rp=1;As=47;a1=(10(Rp/20)-1)/(10(Rp/20)+1);a2=10(-As/20);rip=a1,a2;f=wp/pi,ws/pi;m=1,0;M,f0,m0,w=remezord(f,m,rip);hn=remez(M,f0,m0,w);freqz(hn,1);使用等波纹最佳逼近法设计的FIR滤波
31、器的滤波器的长度为16，归一化通带截止频率为0.20，阻带最小衰减为43db。频率响应如图2-3所示。图2-3 使用等波纹法设计FIR低通滤波器的频率响应2.3使用FDATool设计FIR数字滤波器 MATLAB提供了功能齐全的滤波器设计，与信号处理交互式图形用户界面（Interactive Graphicla User Interface），主要包括FDATool和SPTool两种交互工具。其中FDATool主要用于数字滤波器的设计和分析，而SPTool不仅可以设计，分析滤波器，而且可以对信号进行时域和频域进行分析。本次设计的硬件实现就是采用FDATool工具设计的滤波器。在MATLAB命令
32、窗口中，执行“FDATool”命令，就会启动FDATool滤波器设计与分析工具，其图形界面如图2-4所示。 图2-4 FDATool界面Response type选择低通，选中FIR类型，再选择等波纹。然后选择归一化频率，通带截止频率选择0.2，阻带截止频率选择0.4，通带最小衰减选择1db，阻带最小衰减选择47db，点击DESIGN FILTER按键，即可设计FIR滤波器了。滤波器的频率响应和图2-3中相同。滤波器的系数如式(2-1)所示。 （2-1）2.4 使用Simulink进行仿真 仿真原理：一个50Hz的正弦波和450Hz的正弦波相加得到受干扰的信号，然后通过所设计的FIR低通滤波器
33、滤出。 仿真过程：打开Simulink，新建一个mdl文件，然后将3.3节中设计的FIR滤波器转换成一个Simulink模块，并拖入到新建的mdl文件中。然后在Simulink Library 里找到DSP工具箱，拖出两个正选波模块，一个设置为：幅度为1，频率为50Hz，初始相位为0，采样时间为0.001s,每帧采样点数为100，另一个除了频率为450Hz，初始相位为以外，其它均与第一个相同。再分别拖入频谱示波器和矢量示波器。频谱示波器缓冲区设为512，缓冲输出设为256，矢量示波器帧数设为5帧。最后在拖入相加模块，便可以搭建模型了。连接方式如图2-5。由于篇幅限制，本文只给出了输出信号波形及
34、频谱，以及混杂有噪音的信号波形。具体如图2-6，2-7，2-8所示。可以看到，450Hz的噪音已经衰减到45db以上，所设计的滤波器较好的实现了滤波功能。图2-5 Simulink仿真图2-6 混有噪音的信号波形图2-7 输出的信号波形图2-8 输出的信号的频谱第三章 EDA技术与分布式算法3.1 EDA技术3.1.1 EDA技术简介EDA是Electronic Design Automation 的缩写，意思是电子设计自动化，即利用计算机自动完成电子系统的设计。EDA技术是以计算机和微电子技术为先导，汇集了计算机图形学、拓扑、逻辑学、微电子工艺与结构学和计算数学等多种计算机应用学科最新成果的
35、先进技术。EDA技术作为现代电子设计技术的核心，通过功能强大的计算机和EDA工具软件平台，对以硬件描述语言HDL为系统逻辑描述手段完成的设计文件，自动地完成逻辑编译、逻辑化简、逻辑分割、逻辑综合，以及逻辑优化和仿真测试，以至实现既定的电子线路系统功能。3.1.2 QUARTUS II软件介绍 QUARTUS II是ALTERA公司在21世纪初推出的FPGA/CPLD开发环境，是ALTERA前一代FPGA/CPLD集成开发环境MAX+PLUS II的更新换代产品，其功能强大，界面友好，使用简便。QUARTUS II软件集成了ALTERA的FPGA/CPLD开发流程中所涉及的所有工具和第三方软件接
36、口。图3-1给出了QUARTUS II软件的交互界面。图3-1 QUARTUS II软件的交互界面 QUARTUS II具有一下特点： A支持多时钟定时分析、Logic Lock基于块的设计、SOPC（可编程的片上系统)、内嵌SignalTapII逻辑分析仪和功率估计器等高级工具。 B易于引脚分配和时钟约束。 C强大的HDL综合能力。 D支持的器件种类多。 利用QUARTUS II软件进行FPGA/CPLD开发的流程如图3-2所示。图3-2 FPGA/CPLD开发流程图3.1.3 FPGA介绍可编程逻辑器件PLD(Programmable Logic Devices)是ASIC(Applica
37、tion Specific Integrated Circuits)的一个重要分支。FPGA属于复杂高密度的PLD器件。FPGA使用了可编程的查找表(Look Up Table,LUT)结构,其中LUT是可编程的最小的逻辑构成单元。LUT示意图如图3-3所示。图3-3 查找表结构图由于设计人员可以将存储在片外的EPPROM或者计算机的配置数据控制加载到FPGA器件中进而实现在现场修改器件的逻辑功能，FPGA得到了普遍的应用。 下面以CYLONEII系列FPGA器件为例来介绍一下FPGA器件的结构。CYLONEII系列FPGA器件是由美国Altera公司生产的中端产品。CYLONEII系列FPG
38、A器件采用了90nm工艺，片内逻辑单元数量最多可达68416个逻辑单元，片内存储器容量最多可达1.1Mbts，用户可用引脚最多有622个。CYLONEII系列FPGA器件嵌有乘法器，这些乘法器可用于完成高速乘法操作，使得CYLONEII系列FPGA器件的数字信号处理能力得到增强。CYLONEII系列FPGA器件的速度等级有三个：-6，-7和-8，其中-6的速度最快。 CYLONEII系列FPGA器件的内部资源是按行、列的方式呈二维分布，如图3-4所示。这些资源主要包括逻辑阵列、M4k存储器块、乘法器等。这些资源模块通过FPGA内部的各种连接通路连接起来。图3-4 EP2C20资源分布图 逻辑单
39、元(Logic Element,LE)是FPGA内部用于完成用户资源的最小单元。一个逻辑阵列包含16个逻辑单元以及一些其他资源。一个逻辑单元主要有以下部件组成：一个4输入的查找表，一个可编程的寄存器，一条进位链和一条寄存器级联链。逻辑单元结构图如图3-5所示。 图3-5 逻辑单元结构图CYLONEII系列FPGA器件内部的存储器是以M4k存储器块的形式按列排列的，每个M4k存储器块的大小为4608bit。M4k存储器块包括输入/输出寄存器，输入寄存器用于同步输入信号，输出寄存器在设计中增加一级流水线，使用输入/输出寄存器可以改善电路的性能。M4k存储器除了可以用作标准的存储器使用外，还可以被配置为移位寄存器、先入先出寄存器、只读存储器等。CYLONEII系列FPGA器件内部嵌有硬件乘法器，可以完成高速乘法运算操作。实现许多数字信号处理运算如滤波、快速傅里叶变换、卷积、解相关等。CYLONEII系列FPGA器件关于时钟控制的部分主要是全局时钟网络和锁相环（PLL）。全局时钟网络负责把时钟分配到器件内部的各个单元，控制器件内部的所有资源。锁相环可以完成分频、倍频、移相等关于时钟的基本操作。3.2 FPGA设计方法 传统的数字系统设计一般是从选择具体的元器件开始，通过这些元器件进行逻辑电路设计，完成系统中各独立功能模块的设计，再把这些模块连接起来，组装成整个硬件系统。上述过程
