《数字信号处理》2.doc

1、数字信号处理实 验 指 导 书周春 编 写适用专业： 电子信息工程 通信工程 安 徽 建 筑 工 业 学 院 电 子 与 信 息 工 程 学 院2008 年 09 月1前 言数字信号处理是一门理论和实践密切结合的课程，为深入掌握课程内容， 应当在学习理论的同时上机实验。上机实验不仅可以帮助学生深入的理解和消化基本理论，而且能锻炼初学者独立解决问题的能力。本书根据 课程重点编写了四个实验。采样定理，时域离散系统的时域和频域分析以及系统对输入信号的响应是数字信号处理的基本内容，在学习完该部分内容后，进行实验一，以帮助读者掌握以上内容。同时加深学生对“ 数字信号处理是通过对输入信号的一种运算达到处理

2、目的的一种运算”的理解DFT、FFT 是数字信号处理的重要数学工具，它具有广泛的使用内容。学习完 该部分内容后，进行实验二“用 FFT 对信号进行谱分析”，通过该实验加深理解 DFT 的基本概念、基本性质，及其快速算法FFT。DFT 和 FFT 是学习数字信号处理必须 掌握的内容。数字滤波器的基本理论和设计方法是数字信号处理技术的重要内容。学习这 一部分时，应重点掌握 IIR 和 FIR 两种滤波器的基本设计方法。IIR 滤波器的单位脉冲响应是无限长的， 设计方法是先设计模拟滤波器，然后通过 S-Z 平面转换，求出响应滤波器的系统函数。这里的平面转换有两种方法，即脉冲响应不变法和双线性变换法。

3、后者没有频率混叠的缺点，且转换简单，是一种普遍应用的方法。FIR 滤波器的单位脉冲响应是有限长的，设计滤波器的目的即是求出符合要求的单位脉冲响应。窗函数法是一种基本的，也是重要的设计方法。数字信号处理实验适用于电子信息工程、通信工程等专业学生2的实践性教学环节，通过该课程的学习，使学生能够熟练掌握数字信号处理的一般方法，熟练使用 Matlab 作为编程工具，实现在课堂上学习过的各种数字信号处理方法与滤波器设计方法，在此基础上养成严谨认真的科学工作作风，为今后的工作打下良好的基础。在本课程的学习中，学生应达到如下要求：1、提前预习课本、 实验指导书和参考书籍，巩固和加深对实验中要用到基本理论知识

4、的理解；2、在实验过程中认真观察分析实验结果,运用所学知识解决实验中发生的问题；3、能熟练使用编程工具 Matlab；4、正确处理实验数据，并完成符合要求的实验报告。本实验指导书共包括四个实验，分别是：信号、系统及系统响应、用 FFT 作谱分析、用双 线性变换法设计 IIR 数字滤波器、用窗函数法设计 FIR 数字 滤波器，均 为必做实验。主 讲教师可以根据实际教授的专业和课时情况适当增、减部分实验内容。学生在做完每个实验后，必须独立写出详细的实验报告，包括实验方法、实验过 程（步骤）及结果、实验结果分析等等。3实验简要介绍序号 实验项目名称 内 容 提 要主要仪器设备及台套数实验类型实验要求

5、计划学时每组人数1信号、系统及系统响应1、用 matlab 画出指定的连续信号，并对其采样，观察不同采样频率的采样结果。2、用 matlab 观察信号的时域、频域响应，验证卷积定理计算机Matlab 验证 必做 2 402用 FFT 作谱分析1、用 matlab 显示指定信号并做谱分析2、观察不同点数 FFT 结果计算机Matlab验证 必做 2 403用双线性变换法设计 IIR 数字滤波器用双线性变换法设计巴特沃斯低通 IIR 数字滤波器，显示其幅频率响应，并用其对心电信号滤波计算机Matlab验证 必做 2 404用窗函数法设计 FIR 数字滤波器编写程序实现矩形窗、升余弦窗、改进升余弦窗

6、和二阶升余弦窗的 FIR 数字滤波器计算机Matlab验证 必做 2 40合计 84目 录前 言 1实验一：信号、系统及系统响 应 5实验二：用 FFT 作谱分析 .11实验三：用双线性变换法设计 IIR 数字 滤波器 16实验四：用窗函数法设计 FIR 数字滤波器 .195实验一：信号、系统及系统响应 实验学时：2实验类型：验证 实验要求：必修一、实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对采样定理的理解。2、熟悉时域离散系统的时域特性。3、利用卷积方法观察分析系统的时域特性。4、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分

7、析。二、实验内容1、认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容。2、编制实验用主程序及相应子程序，包括（1）信号产生子程序，用于产生实验中要用到的下列信号序列：a、采样信号序列：对下面的连续信号：（1.1）)(sin)(0tuAetxata进行采样，可得到采样序列（1.2）50),(sin()()( 0 nuTenTxaTaa其中 A 为幅度因子，a 为衰减因子， 0 是模拟角频率， T 为采样间隔。这些参数都要在实验过程中由键盘输入，产生不同的 xa(t)和 xa(n)。b、单位脉冲序列： )(nxbc、矩形序列： 10,NRNc6（2）系统单位脉冲响应序列

8、产生子程序。本实验要用到两种 FIR 系统a、 )(10nRhb、 )3()2(5.)1(5.2n（3）有限长序列线性卷积子程序，用于完成两个给定长度的序列的卷积。可以直接调用 MATLAB 语言中卷积函数 conv。conv 用于两个有限长度序列的卷积，它假定两个序列都从 n=0 开始。调用格式如下：y = conv(x,h)其中参数 x 和 y 是两个已赋值的行向量序列。在完成编制上述子程序的基础上，编制本实验主程序。图 1.1 给出主程序流程框图，供参考。开始调用信号产生子程序调用序列傅氏变换数值计算子程序 ， 求 )(kjeX两次调用绘图子程序 ， 分别绘制 x ( n ) ， 的图形

9、kj改变信号序列否 ？YN调用系统单位脉冲响应序列产生子程序 ， 求 h ( n )调用傅氏变换数值计算子程序 ， 求kjeH两次调用绘图子程序 ， 分别绘制 h ( n ) ， 的图形kj改变 h ( n ) 否 ？YN调用卷积子程序 ， 求 y ( n ) = x ( n ) * h ( n )调用傅氏变换数值计算子程序 ， 求kjeY两次调用绘图子程序 ， 分别绘制 y ( n ) ， 的图形kj结束图 1.17三、实验原理、方法和手段采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域发生的变化以及信号信息不丢失的条件，而且可以加深对傅里叶变换、Z

10、 变换和序列傅里叶变换之间关系式的理解。我们知道，对一个连续信号 xa(t)进行理想采样的过程可以用下式表示（1.3）)(tpa其中 为 的理想采样， p(t)为周期冲激脉冲，即)(txa)(tanT（1.4）的傅里叶变换 为)(txa )(jXa（1.5）msajT)(1上式表明， 为 的周期延拓，其延拓周期为采样角频率)(ja)(ja（ s=2/T） 。只有满足采样定理时，才不会发生频率混叠失真。在计算机上用高级语言编程，直接按照 1.5 式计算 的频谱 很不方)(txa)(jXa便，下面导出用序列的傅里叶变换来计算 的公式。)(jXa将 1.4 式代入 1.3 式，并进行傅里叶变换，（1

11、.6）nTjnatjtjnaaexdtexjX)()()(式中的 xa(nT)就是采样后的序列 x(n)，即：x(n)= xa(nT)x(n)的傅里叶变换为 （1.7）njj eeX)(比较 1.7 和 1.6 可知 （1.8）Tjaj)(说明两者之间只在频率度量上差一个常数因子 T。实验过程中应注意这一差8别离散信号和系统在时域均可以用序列来表示，序列图形给人以形象直观的印象，它可以加深我们对信号和系统的时域特征的理解。本实验还将观察分析几种信号及系统的时域特性。为了在数字计算机上观察分析各种序列的频域特性，通常对 X(ej )在0，2 上进行 M 点采样来观察分析。对长度为 N 的有限长序

12、列 x(n)，有（1.9）10)()(Nnnjkj kemxeX其中 ,.,2kk通常 M 应取的大一些，以便观察谱的细节变化。取模 可绘出幅)(kjeX频特性曲线。一个时域离散线性非移变系统的输入/输出关系为（1.10）mnhxnhxny )()(*)(这里 y(n)为系统的输出序列，x(n)为输出序列。h(n )、x (n)可以是无限长，也可以是有限长，为了计算机绘图观察方便，主要讨论有限长情况。如果 h(n)和x(n)的长度分别为 N 和 M，则 y(n)的长度为 L=N+M-1。这样（1.10）式所描述的卷积运算就是序列移位、相乘和累加的过程，所以编程十分简单。上述卷积运算也可以在频域

13、实现)()(jjj eHXeY（1.11）1.11 式右边的相乘是在个频点 k上的频谱值相乘。四、实验组织运行要求根据本实验的特点、要求和具体条件，采用集中授课形式为主的方式进行。五、实验条件1、计算机92、matlab六、实验步骤1、分析采样序列的特性，产生采样信号序列xa(n)，使 。25025018.4， aA图 1.2 给出了连续信号 xa(t)的幅频特性曲线。由此图可以确定对 xa(t)应采用的采样频率。a、取采样频率 fs=1kHz，即 T=1ms。观察所得采样 xa(n)的幅频特性 和图 1.2 中的)(jeX)(jXa在折叠频率附近有无明显差别。应当注意，实验中所得频率是用序列

14、 xa(n)的傅里叶变换公式求得的，所以在频率度量存在关系： ， 为数字频率， 为T模拟频率。b、改变采样频率，f s=300Hz，观察 ，并做记录（打印曲线） ；进一步)(jeX降低采样频率 fs=200Hz，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的 曲线。)(jeX2、时域离散信号、系统和系统响应分析。a、观察信号 xb(n)和系统 hb(n)的时域和频域特性；利用线性卷积求信号 xb(n)通过系统 hb(n)的响应 y(n)，比较所求响应 y(n)和 hb(n)的时域和频域特性，注意它们之间有无差别，绘图说明，并用所学理论解释所得结果。b、观察系统 ha(n)对信号 xc

15、(n)的响应特性。利用线性卷积求系统响应 y(n)，并判断 y(n)的图形及其非零值序列长度是否与理论结果一致，对 xc(n)= ha(n)=R10(n)，说出一种定性判断 y(n)图形正确与否的方法，调用序列傅里叶变换数值计算子程序，求得 ，观察 特性曲线，定性判断结果的正确性。改变 xc(n)的长(kjeYkjeY度，取 N=5，重复改实验。注意参数变化的影响，说明变化前后的差异，并解释所得结果。3、卷积定理的验证。0 10 20 30 40 50 1 0.8 0.6 0.4 0.2 0 图 1.2 |xa(jf)|10将实验 2 中的信号换成 xa(n)，使 ，重复10374.24. T

16、A，实验 2a，打印 曲；对主程序做简单修改，按式 1.11 计算)(kjeY 1,)()( MkeHXkk jbjaj ， 所得结果之间有无差异？为什么？七、思考题1、在分析理想采样序列特性的实验中，采样频率不同时，相应的理想采样序列的傅里叶变换频谱的数字频率度量是否都相同？它们所对应的模拟频率是否相同？为什么？2、在卷积定理验证的实验中，如果选用不同的频域采样点数 M 值，例如，选 M=10 和 M=20，分别做序列的傅里叶变换，求得 1,0)()()( keHXeYkkk jbjaj ， 所得结果之间有无差异？为什么？八、实验报告1、简述实验目的及实验原理。2、按照实验步骤附上实验过程中

17、的信号序列、系统单位脉冲响应及系统响应序列的时域和幅频特性曲线，并对所得结果进行分析和解释。3、总结实验中的主要结论。4、简要回答思考题。九、其它说明11实验二：用 FFT 作谱分析 实验学时：2实验类型：验证实验要求：必修一、实验目的1、进一步加深 DFT 算法原理和基本性质的理解（因为 FFT 只是 DFT 的一种快速算法，所以 FFT 的运算结果必然满足 DFT 的基本性质）2、熟悉 FFT 算法原理和 FFT 子程序的应用3、学习用 FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用 FFT二、实验内容1、对以下信号逐个进行谱分析tt

18、ttxnxnnnnxR20cos168cos)(i4cs)(,0733,4)(,80)(654324其 他其 他（2.1）下面给出针对各信号的 FFT 变换区间 N 以及对连续信号 x6(t)的采样频率 fs，12供实验时参考。x1(n)，x 2(n)，x 3(n)，x 4(n)，x 5(n)：N=8，16x6(t)：f s=64(Hz)，N=16，32，642、令 x(n)=x4(n)+x5(n)，用 FFT 计算 8 点和 16 点离散傅里叶变换，X(k)=DFTx(n)并根据 DFT 的对称性，由 X(k)求出 X4(k)=DFTx4(n)和 X5(k)=DFTx5(n)，并与 1 中所

19、得的结果进行比较。提示：取 N=16 时，x 4(n) = x4(N-n)，x 5(n) = -x5(N-n)3、令 x(n)=x4(n)+jx5(n)，重复 2。三、实验原理、方法和手段参考教材第三章、第四章。四、实验组织运行要求根据本实验的特点、要求和具体条件，采用集中授课为主的方式进行。五、实验条件1、计算机2、matlab六、实验步骤1、复习 DFT 的定义、性质和用 DFT 作谱分析的有关内容。2、复习 FFT 算法原理与编程思想，并对照 DFT-FFT 运算流图和程序框图，读懂本实验提供的 FFT 子程序。3、编制信号产生子程序，产生式 2.1 中的典型信号供谱分析用。应当注意，如

20、果给出的是连续信号 xa(t)，则首先要根据其最高频率确定采样速率 fs 以及由频率分辨率选择采样点数 N，然后对其进行软件采样（即计算 x(n)= xa(nT)，0nN -1） ，产生对应序列 x(n)。对信号 x6 (t)，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。请实验者根据 DFT 的隐含周期性思考这个问题。4、编写主程序，图 2.1 给出了主程序框图，供参考。5、按实验内容要求，上机实验，并写出实验报告。本实验的主程序比较简单，直接根据图 2.1 所给的框图编写程序即可。编程13的难点是 FF

21、T 子程序。不过，各种语言的 FFT 子程序都可在有关信号处理的程序库中找到。由于 C 语言当前最普及，所以为了为实验提供方便，下面给出 C 语言FFT 函数，供参考。开始读入长度 N调用信号产生子程序产生实验信号调用绘图子程序 （ 函数 ） 绘制时间序列波形图调用 F F T 子程序 （ 函数 ） 计算信号的 D F T调用绘图子程序 （ 函数 ） 绘制 | X ( k ) | 的曲线结束图 2.1/*DITFFT 函数（C 语言）*/fft基 2DITFFT 函数要求：指向复数数组指针 X，FFT 长度为 2m，m 为正整数FFT 输出结果放在输入复数数组中。/*计算 N 点 FFT 子程

22、序*/*xr：=信号序列实部，xi：= 信号序列虚部，N=FFT 变换区间长度 N=2M*/*如果信号长度小于 N，应该给 xr，xi 后面补 0*/*计算结果 X(k)的实部和虚部分别存储在数组 xr 和 xi 中*/void Fft(double xr,double xi,int N,int M)int L,B,J,P,k,I;double rPartKB,iPartKB;/*分别代表 X（K+B ）的实部和虚部*/double rCf128,iCf128;/*rCf 存储旋转因子实部，iCf 存储旋转因子虚部*/*旋转因子数组长度应根据自己需要调整 */14/*计算旋转因子*/doubl

23、e PI2=8.0*aran(1.0)for(i=0;i=K)J=J-K;K=(int)(K/2+0.5);J=J+K;七、思考题1、在 N=8 时，x 2(n)和 x3(n)的幅频特性会相同吗？为什么？N=16 呢？2、如果周期信号的周期预先不知道，如何用 FFT 进行谱分析呢？八、实验报告1、简述实验原理及目的。2、结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用 FFT 作谱分析时有关参数的选择方法。3、总结实验所得主要结论。4、简要回答思考题。九、其它说明16实验三：用双线性变换法设计 IIR 数字滤波器 实验学时：2实验类型：验证实验要求：必修一

24、、实验目的1、熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法。2、掌握数字滤波器的计算机仿真方法。3、通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二、实验内容1、用双线性变换法设计一个巴特沃斯低通 IIR 数字滤波器。设计指标参数为在通频带内频率低于 0.2 时，最大衰减小于 1dB；在阻带内0.3, 频率区间上，最小衰减大于 15dB。2、以 0.02 为采样间隔，打印出数字滤波器在频率区间0，/2上的幅频响应特性曲线。3、用所设计的滤波器对实际心电图信号采样序列（在本实验后面给出）进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。三、

25、实验原理、方法和手段参考教材第六章四、实验组织运行要求根据本实验的特点、要求和具体条件，采用集中授课为主的方式进行。五、实验条件1、计算机172、matlab六、实验步骤1、复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计 IIR 数字滤波器的内容，用双线性变换法设计数字滤波器，使其系统函数 )215.094.1)(358.016.)(705.268.1( (7) 221 6 zzzzzzzH(3.1)式中 (3.2),2,1)()21kzCBAzHkk205.,904.38167.,28.331CBA由 3.1 和 3.2 式可见，滤波器 H(z)由三个二阶滤波器 H1(z)、 H2(z)

26、和 H3(z)级联组成，如下图所示H ( z )H1( z ) H3( z )H2( z )x ( n )y1( n ) y3( n ) = y ( n )y2( n )图 3.1 滤波器 H(z)的组成2、编写滤波器仿真程序，计算 H(z)对心电图信号采样序列 x(n)的响应序列y(n)。设 为第 k 级二阶滤波器 Hk(z)的输出序列， 为输入序列，如图)(k 1yk3.1 所示，由 3.2 式可得到差分方程)2()()2()1(2)()( 11 nyCnBnAynyAy kkkkkk(3.3)当 k=1 时， 。所以 H(z)对 x(n)的总响应序列 y(n)可以用顺序迭)(1xk代算法

27、得到。即依次对 k=1，2，3，求解差分方程 3.3，最后得到 。)(3y仿真程序就是实现上述求解差分方程和顺序迭代算法的通用程序。也可以直接调用 MATLAB filter 函数实现仿真183、在通用滤波器上运行仿真滤波程序，利用 MATLAB 绘图，完成实验内容2、和 3七、思考题1、用双线性变换法设计数字滤波器过程中，变换公式 12zTs中 T 的数值，对设计结果有无影响？为什么？八、实验报告1、简述实验目的及原理。2、由所打印的 特性曲线及设计过程简述双线性变换法的特点。)(jeH3、对比滤波前后的心电图信号波形，说明数字滤波器的滤波过程与滤波作用。4、简要回答思考题。九、其它说明心电

28、图信号采样序列 x(n)人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后才能作为判断心脏功能的有用信息。下面给出一实际心电图信号采样序列样本 x(n)，其中存在高频干扰。在实验中，以 x(n)作为输入序列，滤除其中的干扰成分。x(n)=-4, -2, 0, -4, -6, -4, -2, -4, -6, -6,-4, -4, -6, -6, -2, 6, 12, 8, 0, -16-38, -60, -84, -90, -66, -32, -4, -2, -4, -812, 12, 10, 6, 6, 6, 4, 0, 0, 00, 0, -2, -4, 0, 0,

29、0, -2, -2, 00, -2, -2, -2, -2, 0,19实验四： 用窗函数法设计 FIR 数字滤波器 实验学时：2实验类型：验证实验要求：必修一、实验目的1、掌握用窗函数法设计 FIR 数字滤波器的原理和方法。2、熟悉线性相位 FIR 数字滤波器特性。3、了解各种窗函数对滤波器特性的影响。二、实验内容1、用升余弦窗设计一线性相位低通 FIR 数字滤波器，截止频率 。radc4窗口长度 N=15，33。要求在两种窗口长度情况下，分别求出 h(n)，打印出相应的幅频特性和相频特性曲线，观察 3dB 带宽和 20dB 带宽。总结窗口长度 N 对滤波特性的影响。设计低通 FIR 数字滤波

30、器时，一般以理想低通滤波特性为逼近函数 ，)(jdeH即 ，其中ccje,021N)(sin21)(dedeHhc njnjdd c2、 ，用四种窗函数设计线性相位低通滤波器。绘制相应的幅4/,3c频特性曲线，观察 3dB 和 20dB 带宽以及阻带最小衰减，比较四种窗函数对滤波特性的影响。三、实验原理、方法和手段20如果所希望的滤波器理想频率响应函数为 ，则其对应的单位脉冲响)(jdeH应为 (4.1)nhnjdd)(21)(窗函数设计法的基本原理是用有限长单位脉冲响应序列 h(n)逼近 hd(n)。由于 hd(n)往往是无限长序列，且是非因果的，所以用窗函数 (n)将 hd(n)截断，并进

31、行加权处理，得到：h(n)= hd(n)(n) (4.2)h(n)就作为实际设计的 FIR 数字滤波器的单位脉冲响应序列，其频率响应函数 为： jeH 10)()(NnnjjeheH(4.3)式中，N 为所选窗函数 (n)的长度。用窗函数法设计的滤波器性能取决于窗函数 (n)的类型及窗口长度 N 的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度 N。这样选定窗函数类型和长度 N 后，求出单位脉冲响应 h(n)= hd(n)(n)，并按式 4.3 的要求求出 。 是否满足要求，要进行验算。一般在 h(n)尾)(jeH)(j部加零使长度满足 2 的整数次幂，

32、以便用 FFT 计算 。如果要观察细节，补)(jeH零点数增多即可。如果 不满足要求，则要重新选择窗函数类型和长度 N，)(j再次进行验算，直至满足要求。如果要求线性相位特性，则 h(n)还必须满足：h(n)=h(N-1-n)根据上式中的正、负号和长度 N 的奇、偶性又将线性相位 FIR 滤波器分成四类。要根据所设计的滤波器特性正确选择其中一类。例如，要设计线性相位低通特性，可选择 h(n)=h(N-1-n)一类，而不能选 h(n)=-h(N-1-n)一类。四、实验组织运行要求根据本实验的特点、要求和具体条件，采用集中授课为主的方式进行。五、实验条件1、计算机212、matlab六、实验步骤1

33、、复习用窗函数法设计 FIR 数字滤波器一节内容，阅读本实验原理，掌握设计步骤。2、编写程序。(1)编写能产生矩形窗、升余弦窗、改进升余弦窗和二阶升余弦窗的窗函数子程序。(2)编写主程序，主程序框图如图 4.1 所示，仅供参考，其中幅度特性要求用dB 表示。开始读入窗口长度 N计算 hd( n )调用窗函数子程序求 ( n )计算 h ( n ) = h d ( n ) * ( n )调用子程序 ( 函数 ) 计算 H ( k ) = D F T h ( n ) 调用绘图子程序 ( 函数 ) 绘制 H ( k ) 幅度相位曲线结束图 4.1 用窗函数法设计滤波器主程序框图设 )()()(22k

34、jHkHnhDFTIR画图时，用 20lg|H(k)|打印幅度特性。用 k 点对应的频率 。为使曲线包络更接Nk近 的幅度特性曲线， DFT 变换区间要选大些。例如窗口长度 N=33 时，可)(je通过在 h(n)末尾补零的方法，使长度变为 64，再进行 64 点 DFT，则可得到更精确的幅度衰减特性曲线。FFT 子程序参考实验一，也可以调用 MATLABA 函数 fft 实现。22七、思考题1、如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。2、如果要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为 1和 2，实求理想带通的单位脉冲响应 hd(n)。八、实验报告1、简述实验目的及原理。2、按照实验步骤及要求，比较各种情况下的滤波性能，说明窗口长度 N 和窗函数类型对滤波特性的影响。3、总结用窗函数法设计 FIR 滤波器的主要特点。4、简要回答思考题。九、其它说明