由于MATLAB的广泛使用和功能的不断更新,基于MATLAB的滤波器设计方法以其方便快捷的特点,受到了设计者的欢迎。下面将举例说明基于MATLAB的FIR滤波器的设计。
1 基于FDATool的FIR滤波器设计
使用FDATool设计FIR滤波器的具体步骤如下: 1.1 滤波器指标
假设需要设计一个16阶的FIR滤波器〔h(0)=0〕,给定的参数如下: (1) 低通滤波器
(2) 采样频率FS为48kHz,滤波器FC为10.8kHz (3) 输入序列位宽为9位〔最高位为符号位〕 在此利用MATLAB来完成FIR滤波器系数确实定。 1.2 打开MATLAB的FDATool
MATLAB集成了一套功能强大的滤波器设计工具FDATool〔Filter Design & Analysis Tool〕,可以完成多种滤波器的设计、分析和性能评估。
单击MATLAB主窗口下方的“Start”按钮,如图B.1所示,选择菜单“ToolBox” →“Filter Design” →“Filter Design & Analysis Tool〔FDATool〕”命令,打开FDATool,如图B.2所示。
图B.1 FDATool的启动
图B.2 FDATool的主界面
另外,在MATLAB主命令窗口内键入“fdatool”,同样可打开FDATool程序界面。 1.3 选择Design Filter
FDATool界面左下侧排列了一组工具按钮,其功能分别如下所述: ● 滤波器转换〔TransForm Filter〕
● 设置量化参数〔Set Quantization Parameters〕 ● 实现模型〔Realize Model〕 ● 导入滤波器〔Import Filter〕 ● 多速率滤波器〔Multirate Filter〕 ● 零极点编辑器〔Pole-zero Editor〕 ● 设计滤波器〔Design Filter〕
选择其中的按钮,进入设计滤波器界面,进行以下选择,如图B.3所示。
图B.3 FDATool设计FIR滤波器
● 滤波器类型〔Filer Type〕为低通〔Low Pass〕
● 设计方法〔Design Method〕为FIR,采用窗函数法〔Window〕 ● 滤波器阶数〔Filter order〕定制为15 ● 窗口类型为Kaiser,Beta为0.5 ● FS为48kHz,FC为10.8kHz
最后单击Design Filter图标,让MATLAB计算FIR滤波器系数并作相关分析。 其系统函数H(z)可用下式来表示:
H(z)=bkzk
k116显然上式可以写成:
H(z)=z1bzkk015k
1即可以看成是一个15阶的FIR滤波器的输出结果经过了一个单位延时单元z,所以在FDATool中,把它看成15阶FIR滤波器来计算参数。 1.4 滤波器分析
计算完FIR滤波器系数以后,往往需要对设计好的FIR滤波器进行相关的性能分析,以便了解该滤波器是否满足设计要求。分析操作步骤如下:
选择FDATool的菜单“Analysis”→“Magnitude Response”,启动幅频响应分析如图B.4所示,x轴为频率,y轴为幅度值〔单位为dB〕。
图B.4 FIR滤波器幅频响应
在图的左侧列出了当前滤波器的相关信息:
● 滤波器类型为Direct Form FIR〔直接I型FIR滤波器〕 ● 滤波器阶数为15
选择菜单“Analysis”→“Phase Response”,启动相频响应分析,如图B.5所示。由该图可以看到设计的FIR滤波器在通带内其相位响应为线性的,即该滤波器是一个线性相位的滤波器。
图B.5 滤波器相频响应
图B.6显示了滤波器幅频特性与相频特性的比较,这可以通过菜单“Analysis”→“Magnitude and Phase Response”来启动分析。
图B.6 滤波器幅频和相频响应
选择菜单“Analysis”→“Group Delay Response”,启动群时延分析。 FDATool还提供了以下几种分析工具: ● 群时延响应分析。
● 冲激响应分析〔Impulse Response〕,如图B.7所示。 ● 阶跃响应分析〔Step Response〕,如图B.8所示。 ● 零极点图分析〔Pole/Zero Plot〕,如图B.9所示。
图B.7 冲激响应
图B.8 阶跃响应
图B.9 零极点图
求出的FIR滤波器的系数可以通过选择菜单“Analysis”→“Filter Coefficients”来观察。如图B.10所示,图中列出了FDATool计算的15阶直接I型FIR滤波器的部分系数。
图B.10 滤波器系数
1.5 量化
可以看到,FDATool计算出的值是一个有符号的小数,如果建立的FIR滤波器模型需要一个整数作为滤波器系数,就必须进行量化,并对得到的系数进行归一化。为此,单击FDATool左下侧的工具按钮进行量化参数设置。量化参数有三种方式:双精度、单精度和定点。在使用定点量化前,必须确保MATLAB中已经安装定点工具箱并有相应的授权。 1.6 导出滤波器系数
为导出设计好的滤波器系数,选择FDATool菜单的“File”→“Export”命令,打开Export〔导出〕对话框,如图B.11所示。
图B.11 滤波器系数Export对话框
在该窗口中,选择导出到工作区〔Workplace〕。这时滤波器系数就存入到一个一维变量Num中了。不过这时Num中的元素是以小数形式出现的:
Num=
Columns 1 through 9
-0.0369 0.0109 0.0558 0.0054 -0.0873 -0.0484 0.1805 0.4133 0.4133 Columns 10 through 16
0.1805 -0.0484 -0.0873 0.0054 0.0558 0.0109 -0.0369
2 基于MATLAB内建函数的FIR设计
在Matlab中已经内建有各种滤波器的设计函数,可以直接在程序中调用,这里介绍其中几个函数。 2.1 fir1函数
功能:设计标准频率响应的基于窗函数的FIR滤波器。 语法:b=fir1(n,Wn);
b=fir1(n,Wn,‘ftytpe’); b=fir1(n,Wn,Window);
b=fir1(n,Wn,‘ftype’,Window);
说明:fir1函数可以实现加窗线形相位FIR数字滤波器设计,它可以设计出标准的低通、高通、带通和带阻滤波器。
b=fir1(n,Wn)可得到n阶低通,截至频率为Wn的汉明加窗线形相位FIR滤波器,0≤Wn≤1,Wn=1相当于0.5fs。滤波器系数包含在b中,可表示为
b(z)b(1)b(2)z1b(n1)Zn
当Wn=[W1 W2]时,fir1函数可得到带通滤波器,其通带为W1<w<W2。 当ftype=high时,设计高通FIR滤波器;当ftype=stop时,设计带阻滤波器。
在设计高通和带阻滤波器时,由于对奇次阶的滤波器,其在Nyquist频率处的频率响应为零,不
适合构成高通和带阻滤波器。因此fir1函数总是使用阶数为偶数的滤波器,当输入的阶数为奇数时,fir1函数会自动将阶数加1。
b=fir1(n,Wn,Window)利用参数Window来指定滤波器采用的窗函数类型。其默认值为汉明窗。 b=fir1(n,Wn,‘ftype’,Window)可利用ftype和Window参数,设计各种滤波器。 2.2 fir2函数
功能:设计任意频率响应的基于窗函数的FIR滤波器。 语法:b=fir2〔n,f,m〕;
b=fir2〔n,f,m,Window〕; b=fir2〔n,f,m,npt〕;
b=fir2〔n,f,m,npt,window〕; b=fir2〔n,f,m,npt,lap〕;
b=fir2〔n,f,m,npt,lap,Window〕;
说明:fir2函数可以用于设计有任意频率响应的加窗FIR滤波器,对标准的低通、带通、高通和带阻滤波器的设计可使用fir1函数。
b=fir2〔n,f,m〕可设计出一个n阶的FIR滤波器,其滤波器的频率特性由参数f和m决定。参数f为频率点矢量,且f∈[0,1],f=1对应于0.5fs。矢量f按升序排列,且第一个元素必须是0,最后一个必须为1,并可以包含重复的频率点。矢量m中包含了与f相对应的期望得到的滤波器的幅度。
b=fir2〔n,f,m,Window〕中用参数Window来指定使用的窗函数类型,默认值为汉明窗。 b=fir2〔n,f,m,npt〕中用参数npt来指定fir2函数对频率响应进行内插的点数。 b=fir2〔n,f,m,npt,lap〕中用参数lap来指定fir2在重复频率点附近插入的区域大小。
3 基于FDATool的HDL代码产生
在MATLAB 7 中,对数字滤波器的设计提供了与假设干种现实方案的接口。此类接口提供MATLAB到设计工具的无缝连接,即MATLAB根据设计工具的文件格式,将包含滤波器设计参数的文件输出。设计工具导入该文件,并作为设计模块的一部分。在此类接口中包括与Xilinx公司和TI公司的接口,还包括C头文件以及HDL代码。
与FDATool的启动类似,单击MATLAB主窗口下方的“Start”按钮,选择“ToolBox”→“Filter Design HDL coder”→“Filter Design & Analysis Tool(FDATool)”,打开FDATool。
根据上节相同的设计和分析步骤,对FIR滤波器进行分析和设计,在设计完毕之后,可以得到滤波器系数。此时就可以应用设计工具接口。 3.1 C语言头文件的产生
选择FDATool菜单的“Targets”→“Generate C Header”命令,打开产生C语言头文件的窗口,
如图B.12所示。
图B.12 产生C语言头文件
可以看到输出的头文件中,变量名和变量长度名可以自定义,变量输出的格式也有很多种可以选择。根据在FDATool量化时选用的量化方式,窗口中会显示推荐使用的输出格式。
产生的头文件内容如下:
/*
* Filter Coefficients (C Source) generated by the Filter Design and Analysis Tool * Generated by MATLAB(R) 7.0 * Generate on:22-Dec-2005 11:42:24 */ /*
* Discrete-Time FIR Filter (real) *------------------------------------
* Filter Structure :Direct-Form FIR * Filter Order :15 * Stable :Yes * Linear Phase :Yes (Type 2) */
/* General type conversion for MATLAB generated C-code */ # include“tmwtypes.h” /*
* Expected path to tmwtypes.h
* C:\\MATLAB7\\extern\\include\mwtype.h
*/
const int BL = 16; const real64_T B[16]={
-0.03687003131181, 0.01091744268631, 0.0558306521771, 0.005429393216792, -0.08726921427845, -0.04839711653448, 0.1804973650249, 0.4133400743259 0.4133400743259, 0.1804973650249, -0.04839711653448, -0.08726921427845, 0.005429393216792, 0.0558306521771, 0.01091744268631, -0.03687003131181 };
3.2 Xilinx系数文件的产生
选择FDATool菜单的“Targets”→“Xilinx Coefficients〔.COE〕file”命令,MATLAB直接提示文件的保存位置,保存完毕之后另开一个窗口显示该文件的内容。典型的文件内容如下所示。
;
; XILINX CORE Generator(tm) Distributed Arithmetic FIR filter coefficient (.COE) File ; Generated by MATLAB (R) 7.0 and the Filter Design Toolbox 3.0. ;
; Generated on:22-Dec-2005 12:03:10 ;
Radix = 16 ;
Coefficient_Width = 16 ; CoefData = 0000, 079c, 05a2, f90b, f166, 0000, 2272, 4000, 4000, 2272, 0000, f166, f90b, 05a2,
079c, 0000;
这些系数已经运用二进制补码对其进行了运算,所以可以表示成十六进制。 3.3 CCS接口
CCS〔Code Composer Studio〕是TI公司为其数字信号处理器开发的集成开发环境〔IDE〕。在HDL coder中,可以向CCS输出C语言头文件或者直接写入处理器的存储空间中,然后由CCS本身作进一步的处理。在此接口中,可以选择目标板和目标处理器的型号,在这里目标板的型号一般只能是TI公司提供的专用EVM板。上述型号必须手动输入。
选择FDATool菜单的“Targets”→“Code Composer Studio(R)IDE”命令,打开CCS接口。此接口的窗口如图B.13所示。
图B.13 CCS接口窗口
3.4 HDL语言的生成
选择FDATool菜单的“Targets”→“Generate HDL”命令,打开生成HDL源代码的窗口,该窗口如图B.14所示。
图B.14 生成HDL源代码的窗口
在生成HDL文件时,可以选择生成VHDL或是Verilog HDL。可以对时钟信号、复位信号、进程、构造体和实体的定义等等进行设置。最后可以选择文件输出的文件夹和文件名。
文件一般有如下内容和格式:
-- ---------------------------------------------------------- -- Module : filter
-- Generated by MATLAB(R) 7.0 and the Filter Design HDL Coder 1.0 -- Generated on : 2005-12-22 12:33:35 -- HDL Code
-- 此处显示滤波器特性以及输出HDL代码时的设置 ------------------------------------------------------------------ LIBRARY IEEE;
USE IEEE.std_logic_1164.all; USE IEEE.numeric_std.ALL;
ENTITY filter IS
PORT ( clk : IN std_logic; clk_enable : IN std_logic; reset : IN std_logic;
filter_in : IN std_logic_vector (15 DOWNTO 0) ; -- sfix16_En15 filter_out : OUT std_logic_vector (15 DOWNTO 0) --sfix16_En9
);
END filter;
--------------------------------------------------------------------- --Module Architecture : filter
--------------------------------------------------------------------- ARCHITECTURE rt1 OF filter IS --局部函数 --类定义
TYPE delay_pipeline_type IS ARRAY (NATURAL range<>) OF signed (15 DOWNTO 0); --sfix16_En15 --定义常数
CONSTANT coeff1:signed(15 DOWNTO 0):=to_signed(0,16);--sfix16_En14 CONSTANT coeff2:signed(15 DOWNTO 0):=to_signed(1948,16);--sfix16_En14 CONSTANT coeff3:signed(15 DOWNTO 0):=to_signed(1442,16);--sfix16_En14 CONSTANT coeff4:signed(15 DOWNTO 0):=to_signed(-1781,16);--sfix16_En14 CONSTANT coeff1:signed(15 DOWNTO 0):=to_signed(0,16);--sfix16_En14 CONSTANT coeff2:signed(15 DOWNTO 0):=to_signed(1948,16);--sfix16_En14 CONSTANT coeff3:signed(15 DOWNTO 0):=to_signed(1442,16);--sfix16_En14 CONSTANT coeff4:signed(15 DOWNTO 0):=to_signed(-1781,16);--sfix16_En14 CONSTANT coeff5:signed(15 DOWNTO 0):=to_signed(-3738,16);--sfix16_En14 CONSTANT coeff6:signed(15 DOWNTO 0):=to_signed(0,16);--sfix16_En14 CONSTANT coeff7:signed(15 DOWNTO 0):=to_signed(8818,16);--sfix16_En14 CONSTANT coeff8:signed(15 DOWNTO 0):=to_signed(16384,16);--sfix16_En14 CONSTANT coeff9:signed(15 DOWNTO 0):=to_signed(16384,16);--sfix16_En14 CONSTANT coeff10:signed(15 DOWNTO 0):=to_signed(8818,16);--sfix16_En14 CONSTANT coeff11:signed(15 DOWNTO 0):=to_signed(0,16);--sfix16_En14 CONSTANT coeff12:signed(15DOWNTO 0):=to_signed(-3738,16);--sfix16_En14 CONSTANT coeff13:signed(15DOWNTO 0):=to_signed(-1781,16);--sfix16_En14 CONSTANT coeff14 :signed(15DOWNTO 0):=to_signed(1442,16);--sfix16_En14 CONSTANT coeff15:signed(15 DOWNTO 0):=to_signed(1948,16);--sfix16_En14 CONSTANT coeff16:signed(15 DOWNTO 0):=to_signed(0,16);--sfix16_En14 --信号
SIGNAL delay_pipeline :delay_pipeline_type (0 To 15);sfix16_En15 SIGNAL product15 :signed(31 DOWNTO 0);--Sfix32_En29
SIGNAL product14 :signed(31 DOWNTO 0);--Sfix32_En29 SIGNAL product13 :signed(31 DOWNTO 0);--Sfix32_En29 SIGNAL product12 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product10 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product9 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product8 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product7 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product5 :signed(31 DOWNTO 0);-- Sfix32_En29 SIGNAL product4 SIGNAL product3 SIGNAL product2 SIGNAL sum1 SIGNAL add_temp SIGNAL sum2 SIGNAL add_temp_1 SIGNAL sum3 SIGNAL add_temp_2 SIGNAL sum4 SIGNAL add_temp_3 SIGNAL sum5 SIGNAL add_temp_4 SIGNAL sum6 SIGNAL add_temp_5 SIGNAL sum7 SIGNAL add_temp_6 SIGNAL sum8 SIGNAL add_temp_7 SIGNAL sum9 SIGNAL add_temp_8 SIGNAL sum10 SIGNAL add_temp_9 SIGNAL sum11 SIGNAL add_temp_10 :signed(31 DOWNTO 0);-- Sfix32_En29 :signed(31 DOWNTO 0);-- Sfix32_En29 :signed(31 DOWNTO 0);-- Sfix32_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(32 DOWNTO 0);--sfix33_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29 :signed(39 DOWNTO 0);--sfix40_En29 :signed(40 DOWNTO 0);--sfix41_En29
SIGNAL output_typeconvert :signed (15 DOWNTO 0); -- Sfix16_EN9 SIGNAL output_register :signed (15 DOWNTO 0); -- Sfix16_EN9 BEGIN --块声明
Delay_pipeline_process : PROCESS (clk , reset ) BEGIN
IF reset = ‘1’ THEN
delay_pipeline(0 to 15 )<=(OTHERS=>(OTHERS=>’0’)); ELSIF clk’event AND clk = ‘1’ THEN
IF clk_enable=’1’ THEN
delay_pipeline(0)<=signed(fliter_in);
delay_pipeline(1 to 15 )<=delay_pipeline(0 to 14) ; END IF ; END IF ;
END PROCESS Delay_pipeline_process; Product15 <= delay_pipeline(14) * coeff15 ; Product14 <= delay_pipeline(13) * coeff14 ; Product13 <= delay_pipeline(12) * coeff13 ; Product12 <= delay_pipeline(11) * coeff12 ; Product10 <= delay_pipeline(9) * coeff10;
Product9 <= resize(delay_pipeline(8)(15 DOWNTO 0) & “ 00000000000000”,32); Product8 <= resize(delay_pipeline(7)(15 DOWNTO 0) & “ 00000000000000”,32); Product7 <= delay_pipeline(6) * coeff7; Product5 <= delay_pipeline(4) * coeff5 ; Product4<= delay_pipeline(3) * coeff4 ; Product3 <= delay_pipeline(2) * coeff3 ; Product2 <= delay_pipeline(1) * coeff2 ;
add_temp<=resize(product2 ,33) +resize(product3,33); Sum1<=resize(add_temp,40);
add_temp_1<=resize(sum1 ,41) +resize(product4,41); Sum2<= add_temp_1 (39 DOWNTO 0);
add_temp_2<=resize(sum2 ,41) +resize(product5,41); Sum3<= add_temp_2 (39 DOWNTO 0);
add_temp_3<=resize(sum3 ,41) +resize(product7,41); Sum4<= add_temp_3 (39 DOWNTO 0);
add_temp_4<=resize(sum4 ,41) +resize(product8,41); Sum5<= add_temp_4 (39 DOWNTO 0);
add_temp_5<=resize(sum5 ,41) +resize(product9,41); Sum6<= add_temp_5 (39 DOWNTO 0);
add_temp_6<=resize(sum6 ,41) +resize(product10,41); Sum7<= add_temp_6 (39 DOWNTO 0);
add_temp_7<=resize(sum7 ,41) +resize(product12,41); Sum8<= add_temp_7(39 DOWNTO 0);
add_temp_8<=resize(sum8 ,41) +resize(product13,41); Sum9<= add_temp_8 (39 DOWNTO 0);
add_temp_9<=resize(sum9 ,41) +resize(product14,41); Sum10<= add_temp_9 (39 DOWNTO 0);
add_temp_10<=resize(sum10 ,41) +resize(product15,41); Sum11<= add_temp_10 (39 DOWNTO 0);
output_typeconvert<=resize(shift_right(sum11(35 DOWNTO 0)+(“0”& (sum11 (20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20) & NOT sum11(20))),20 ) ,16);
Output_Register_process : PROCESS (clk ,reset) BEGIN
IF reset =’1’ THEN
Output_register<=(OTHERS=>’0’); ELSIF clk’event AND clk =‘1’THEN IF clk_enable = ‘1’ THEN
Output_register<=output_typeconvert; END IF; END IF;
END PROCESS Output_Register_process; --任务声明
filter_out <=std_logic_vector(output_register);
END rtl;
与此同时,还会生成一个所谓test bench文件,产生用于检验滤波器响应的测试向量。
实验 3: 设计 FIR〔有限冲激响应〕滤波器
针对 Spartan-3E 开发套件
介绍
在这个实验里, 将向你展示通过系统发生器的FIR 和 FDATool模块来指定、模拟和实现FIR滤波器的方法.
FDATool 模块被用来定义滤波器的阶数和系数, FIR模块被用作 Simulink 模拟以及在 FPGA中用 Xilinx ISE来实现设计.你也可以通过实际硬件来运行它以验证这个设计的功能.
注意: 在 c:\\xup\\dsp_flow\\labs\\labsolutions\\lab3\\ 目录下有完整的例子.
目标
在完成这个试验后, 你将能够:
• 用 FDATool 模块输入你的滤波器指数,在设计中使用产生的系统或把它储存在工作区
• 使用具有 FDATool 模块产生的系数的FIR模块, 并在Simulink中运行 bit-true 仿真 • 产生设计并用 Resource Estimator 模块和post-map 报告估计资源利用量
设计描述
你是Cyberdyne系统的一个dsp设计者. 你的公司正在调查用数字滤波器代替安防检测器中的模拟滤波器,以尝试提高性能和降低整个系统的成本. 这将使贵公司可以进一步渗入日益增长的安防市场. 一个单频取样滤波器设计如下:
• Sampling Frequency (Fs)〔采样频率〕 = 1.5 MHz • Fstop 1 = 270 kHz • Fpass 1 = 300 kHz • Fpass 2 = 450 khz • Fstop 2 = 480 kHz
• Attenuation on both sides of the passband〔双边通频带衰减〕 = 54 dB • Pass band ripple〔通频带脉动〕 = 1
因为灵活性和上市时间的原因,Cyberdyne已经选择FPGA来实现它. 你的HDL设计经验是有限的. 因为你对MathWorks的产品比较熟悉,所以System Generator for DSP会是一个在FPGA中实现滤波器的优秀解决方案.
你的经理 Miles Booth 已经要求你创建一个在即将完成的 Spartan-3E™ 原型板上实现的滤波器的原型. 这个原型必须尽可能快地完成,这是因为Aggressive Security会议即将来临, 这个会议是业界今年最大的会议,我们不能错过它.
你的经理已经提供了具有输入源和输出接收端的最初模型. 你的设计必须用随机输入源和来自DSP Blockset的chirp信号来仿真. 为分析滤波器的输出, 输入和输出信号在频谱示波器中显示. 频谱示波器用来比较在FPGA中实现的定点 FIR 滤波器的录放频谱响应.
两个不同的输入源用来仿真这个滤波器:
• chirp 信号模块, 其扫描频率为 6 KHz 到 10 KHz ,不管其瞬时输出频率
• 随机信号源发生器, 它输出范围在-1.9到1.9的均匀分布的随机信号. 因为均匀信号是有界的,它是驱动定点
滤波器的较好的选择.
设计流程
这个试验由六个主要步骤组成. 在步骤1,你将用系统发生器的FDATool模块来产生指定FIR滤波器的系数. 在步骤2, 你将把这个系数与 FIR filter 模块联系起来. 步骤3要求你用两个已提供的输入源仿真这个设计并分析系数和输入信号宽度变化的影响. 在步骤4, 你将要添加一个转换模块以调整输出宽度,添加一个延时模块以提高效率,添加一个资源估计模块以估计资源占用量,从而完成整个设计. 在下一步将实现这个设计. 在每一步总的指示下面, 你会发现手把手的说明和图表,指
出了实现总的指示的更多细节. 如果你对个别指导有信心, 可以略过手把手的说明,跳到下一步的总指示.
注意: 如果你当时不能完成这个实验, 你可以从Xilinx大学计划网站 ://university.xilinx 下载实验文件
产生FIR滤波器所需的系数
步骤 1
实验总流程: 步骤 1: 产生 滤波器Step 5: 系数 实现设计 步骤 2: 把系数与滤波器相关联 Step 6: 估计资源 占用量
步骤 3: 仿真 滤波器 Step 7: 执行 HW-in-the-L
oop 步骤 4: 完成整个 设计
从 Xilinx Blockset DSP blockset中添加FDATool模块到包含一个 DA FIR 滤波器的设计中. 按照以下要求用FDATool模块产生FIR滤波器的系数
• Sampling Frequency (Fs) 〔采样频率〕 = 1.5 MHz • Fstop 1 = 270 KHz • Fpass 1 = 300 KHz • Fpass 2 = 450 Khz • Fstop 2 = 480 KHz
• Attenuation on both sides of the passband〔双边通频带衰减〕 = 54 dB • Pass band ripple〔通频带脉动〕 = 1
在 Matlab中, 改变目录为 c:/xup/dsp_flow/labs/lab3/: 在命令行窗口键入
cd c:/xup/dsp_flow/labs/lab3/ .
从MATLAB控制台窗口打开 bandpass_filter.mdl 模块 从 Xilinx Blockset DSP 添加FDATool模块到这个设计 在 FDATool Design Filter 窗口 (图 3-1)输入以下滤波器参数
• Response Type: Bandpass • Units: KHz
• Sampling Frequency (Fs) = 1.5 MHz
• Fstop 1 = 270 KHz • Fpass 1 = 300 KHz • Fpass 2 = 450 Khz • Fstop 2 = 480 KHz
• Attenuation on both sides of the passband = 54 dB (Astop1 and Astop2 parameters) • Pass band ripple = 1 (Apass)
图 3-1. 在 FDATool中设计一个滤波器. 点击 Design Filter 按钮以确定滤波器指令
频谱窗口将被更新并显示如图 3-2
图 3-2. 所设计滤波器的幅频响应.
1. 基于所定义的技术要求, 最小的滤波器阶数是多少?
用 File Save Session来储存fda格式文件coefficients.fda中的系数
注意: 这是可选的步骤. 这些系数对这个设计仍然可以利用. 如果你把参数储存在 fda-file中,你可以通过 FDATool 模块参数对话框来加载它们.
用File Export导出工作区中的系数,其Numerator variable name为Num〔图3-3〕
注意: 这将在你的 MATLAB 工作区添加Num变量. 对于一个 FIR 滤波器, Num 代表用在滤波器中的系数. 因为通过 FDATool 模块仍能利用这些系数,这也是一可选的步骤
图 3-3. 在工作区中导出系数.
在 MATLAB 控制台窗口键入Num 来查看系数列表
在 MATLAB 控制台窗口键入max(Num) 以确定指出系数宽度和二进制小数点的最大系数值
2. 填入与系数相关的以下信息 最大值: 最小值:
____________________
把系数与 FIR 滤波器相关联
实验总流程:
Step 1: Generate
Filter Step 5: Coefficients Implement Design
步骤 2
Step 2: Associate Coefficients Step 6: to Filter Estimate Resources Step 3: Simulate the Filter Step 7: Perform HW-in-the-L
oop Step 4: Complete the Design
从 Xilinx Blockset DSP library添加FIR 滤波器模块并与产生的系数关联
从 Xilinx Blockset DSP library 添加FIR (DA FIR v9_0)滤波器到设计中 双击 Xilinx FIR Filter 模块 并在模块参数窗口 (图 3-4)中输入以下参数
• • • • •
Coefficients : xlfda_numerator(‘FDATool’)
Coefficient Structure : Inferred from Coefficients Number of bits per Coefficients : 12 Binary Point for Coefficients : 11 Provide Valid Ports : unchecked
图 3-4. FIR 滤波器模块参数.
点击 OK 以接受设置
如图 3-5所示连接设计的各个模块
图 3-5. 设计用来仿真的FIR 滤波器模块.
在Simulink中仿真 FIR 滤波器
步骤 3
实验总流程: Step 4: Step 2: Step 3: Step 1:
Complete the Associate Simulate the Generate
Design Coefficients Filter Filter Step 5:
Step 6: Step 7: to Filter Coefficients Implement
Estimate Perform Design
Resources HW-in-the-L oop
设置样值输入为 FIX_8_6 ,输入采样周期为 1/1500000. 利用频谱示波器, 研究 chirp 和 noise 信号的输出.
双击 Gateway In 模块并设置格式为 FIX_8_6 ,采样周期为 1/1500000 选择 Chirp 信号为输入源并开始仿真
如果你收到以下错误, 根据消息更新FIR的等待时间, 并再次仿真.
调出示波器并验证一直衰减的FIR滤波器的输出信号,这些信号应如图 3-7 和 图 3-8所示.
图 3-7. 无衰减通频带 (频谱示波器).
图 3-8. 禁带衰减 (频谱示波器).
停止仿真
选择 Random Source 并运行仿真
图 3-10. 随机信号输入源 (频谱示波器).
停止仿真
完成 FIR 滤波器设计
实验总流程: Step 1: Generate
Filter Step 5: Coefficients Implement Design
步骤 4
Step 2: Associate Coefficients Step 6: to Filter Estimate Resources
Step 3: Simulate the Filter Step 7: Perform HW-in-the-L
oop Step 4: Complete the Design 加入 convert 模块以得到 FIX_8_6 格式的输出,从而降低显示所需的动态量程. 加入 delay 部件以提高性能. 利用 Resource Estimator 估计设计所需的资源.
从 Xilinx Blockset Basic Elements 库中添加一个Convert模块到 FIR 滤波器的输出 以使输出符合
FIX_8_6 格式并改变 quantization 为Truncate , Overflow 为 Wrap
从 Xilinx Blockset Basic Elements 库中添加一个延时部件到输出以实现流水线技术 并提高性能,流水
线技术会在输出衰减器中得到实现
确信 FIR size 被设置为 FIX_12_11 , Gateway In size 被设置为 FIX_8_6 从 Xilinx Blockset Index 库中添加Resource Estimator模块到设计中
注意: 注意你的设计应如图 3-11所示.
图 3-11. 完成的FIR 滤波器设计.
实现 FIR 滤波器
实验总流程: Step 1: Generate Filter Step 5: Coefficients Implement Design
步骤 5
Step 2: Associate Coefficients Step 6: to Filter Estimate Resources
Step 3: Simulate the Filter Step 7: Perform HW-in-the-L
oop Step 4: Complete the Design
设置 FIR 硬件过采样率为 9, FIR 核等待时间为14, 然后运行仿真. 这将更新采样率为 7.407e-008. 按照以下技术要求用系统发生器产生代码. 在工程向导里打开 bandpass_filter.ise 工程, 综合并执行它.
• Input : Width = FIX_8_6, Quantization = Truncate, Overflow = Wrap • Output Width: FIX_8_6 • FIR Core Latency: 14
• FIR Hardware Over-Sampling Rate: 9 • FIR Coefficients: FIX_12_12 • Compilation: HDL Netlist
• Part: Spartan3e xc3s500e-4fg320 • Synthesis Tool: XST
• Target Directory: c:/xup/dsp_flow/labs/lab3/ise • Create Testbench: Unchecked
• Simulink System Period (sec): 7.407e-008 • FPGA System Clock Period (ns): 20
双击 FIR 模块并设置 Hardware Over-Sampling Rate 为 9 , latency 为 14
3. Hardware Over-Sampling Rate 设置为 9 后,滤波器的实现结果有何变化? 为何要把它设为 9 而不是 8?
运行仿真. 如果采样率更新窗口出现, 接受采样率为 7.407e-008 并再次运行仿真 双击系统发生器符号并设置以下参数
• • • • • • • Compilation: HDL Netlist
Part: Spartan3e xc3s500e-4fg320 Synthesis Tool: XST
Target Directory: c:/xup/dsp_flow/labs/lab3/ise (or ./ise) Create Testbench: Unchecked
Simulink System Period (sec): 7.407e-008 FPGA System Clock Period (ns): 20
双击 Generate 按钮以产生设计
在资源管理器中双击 bandpass_filter_cw.ise 以打开Project Navigator 高亮突出顶层文件并双击执行
4. 利用各种报告, 答复以下问题
Slices的数目: 有无定时约束? 实际时钟周期:
5. 在选择连续执行的情况下, 滤波器的实际时钟是多少?
用 Resource Estimator估计资源利用率
步骤 6
实验总流程: Step 1: Generate
Filter Step 5: Coefficients Implement Design
Step 2: Associate Coefficients Step 6: to Filter Estimate Resources Step 3: Simulate the Filter Step 7: Perform HW-in-the-L
oop Step 4: Complete the Design
用 Resource Estimator 模块和 post-map 报告, 估计此设计所用的资源.
双击 resource estimator 模块
从下拉框中选择 Post-Map 并点击 Estimate 按钮.
ISE 将自动通过综合和映射来运行并产生映射报告.
6. 资源估计器的报告有什么内容?
Slices的数目: FFs的数目: LUTs数目:
实现 Hardware-in-the-Loop 验证 步骤 7
实验总流程: Step 1: Generate
Filter Step 5: Coefficients Implement Design
Step 2: Associate Coefficients Step 6: to Filter Estimate Resources
Step 3: Simulate the Filter Step 7: Perform HW-in-the-Loop Step 4: Complete the Design
利用系统发生器产生硬件电路并通过硬件电路板来验证设计. 通过Simulink仿真这个设计.
储存模型为 bandpass_filter_hwcosim.mdl 双击系统发生器符号并设置以下参数
• • • • • Compilation: Hardware Co-simulation xup sp3e_starter_kit Synthesis Tool: XST
Target Directory: c:/xup/dsp_flow/labs/lab3/sp3e (or ./sp3e) Create Testbench: Unchecked
Simulink System Period (sec): 7.407e-008
点击 Generate 按钮,会出现一个对话框,显示编译处理的过程,如图 3-12
图 3–12. 在命令窗口的编译过程.
当发生器成功完成后, 一个新的 Simulink 库窗口将打开,并出现一个具有适当输入输出数的模块.
图 3–13. 在新的 Simulink 窗口中打开的被编译模块. 复制被编译模块到设计中并连结它如图 3-14所示
图 3–14. 在循环仿真中完成硬件的设计准备工作.
连接硬件板并通过Simulink对设计仿真. 把电源线与硬件板相连. 板子的LEDs 将闪亮. 在板子与PC间连接下载电缆
在硬件协同仿真模块上双击并设置下载电缆为USB. 在Simulink窗口中点击 run 按钮 (
3-15)
) 以运行仿真. 配置文件将被下载,仿真将开始运行
仿真结果将在输出示波器中显示,其中 Simulink 仿真器的输出显示在顶部,硬件输出显示在底部(Figure
图 3–15. 仿真结果,Simulink 仿真器的输出显示在顶部. 硬件输出显示在底部 完成后关闭电源
储存模块并关闭 MATLAB
结论
在这个试验里, 你学会了如何用 FDATool 创建滤波器系数并从FIR滤波器模块使用它们. 学会使用Resource Estimator模块估计资源占用量.
答案
1. 基于定义的技术要求, 最小的滤波器阶数是多少?
92
2. 完成以下与系数相连的信息
最大值: 最小值:
0.1610 -0.1541
3. Hardware Over-Sampling Rate 设置为 9 后,滤波器的实现结果有何变化? 为何要把它
设为 9 而不是 8?
它指出滤波器可以连续实现. 每次输出会占用9个时钟周期
. 设为9而不设为8是因为要实现系数的对称. 在对称的条件下, 可以在通过串行乘法器之前实现采样. 这实现了8位数相加,结果为9位数的执行结果 (位增长).
4. 利用各种报告, 答复以下问题
Slices的数目: 278 有无时间约束? Yes 实际的时钟周期: ~146 MHz
5. 在你选择连续执行的情况下, 滤波器实际的时钟周期是多少? 1.5MHz x 9 = 135 MHz
6. 资源估计器的报告内容?
Slices的数目: FFs的数目: LUTs的数目: 282 481 410
因篇幅问题不能全部显示,请点此查看更多更全内容