数字信号处理教程 程佩青 课后题答案
- 格式:pdf
- 大小:5.31 MB
- 文档页数:28
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
7.2 课后习题详解7-1 用冲激响应不变法将以下Ha (s )变换为H (z ),抽样周期为T 。
(1)H a (s )=(s +a )/[(s +a )2+b 2];(2)H a (s )=A/(s -s 0)n0,n 0为任意正整数。
解:(1)冲激响应不变法满足h (n )=h a (t )|t =nT =h a (nT ),T 为抽样间隔。
这种变换法必须让H a (s )先用部分分式展开。
由推出由冲激响应不变法可得(2)先引用拉氏变换的结论,可得按且可得可以递推求得7-2 设计一个模拟低通滤波器,要求其通带截止频率f p=20Hz,其通带最大衰减为R p=2dB,阻带截止频率为f st=40Hz,阻带最小衰减为A s=20dB,采用巴特沃思滤波器,画出滤波器的幅度响应。
解:巴特沃思模拟低通滤波器设计流程为:①利用教程(7-5-24)式求解滤波器阶次N;②利用教程(7-5-27a)式求解3dB截止频率Ωc;③查教程表7-2或表7-4获得归一化巴特沃思低通滤波器的系统函数H an(s);④将H an(s)根据Ωc的值去归一化求得所需的系统函数H a(s)。
已知Ωp=2π×20rad/s,Ωst=2π×40rad/s,R p=2dB,A s=20dB。
(1)按给定的参数由教程(7-5-24)式可求得取N=4。
(2)按教程(7-5-27a)式可求得巴特沃思滤波器3dB处的通带截止频率Ωc为(3)查教程表7-2可得N=4时归一化巴特沃思低通滤波器H an(s)(4)去归一化,求得所需的H a(s)为滤波器的幅度响应如图7-1所示。
图7-17-3 设计一个模拟高通滤波器,要求其阻带截止频率f st=30Hz,阻带最小衰减为A s=25dB,通带截止频率为f p=50Hz,通带最大衰减为R p=1dB。
(1)采用巴特沃思滤波器;(2)采用切比雪夫滤波器;(3)利用MATLAB工具箱函数设计椭圆函数滤波器。
目 录第1章 离散时间信号与系统1.1 复习笔记1.2 课后习题详解1.3 名校考研真题详解第2章 Z变换与离散时间傅里叶变换(DTFT)2.1 复习笔记2.2 课后习题详解2.3 名校考研真题详解第3章 离散傅里叶变换(DFT)3.1 复习笔记3.2 课后习题详解3.3 名校考研真题详解第4章 快速傅里叶变换(FFT)4.1 复习笔记4.2 课后习题详解4.3 名校考研真题详解第5章 数字滤波器的基本结构5.1 复习笔记5.2 课后习题详解5.3 名校考研真题详解第6章 几种特殊滤波器及简单一、二阶数字滤波器设计6.1 复习笔记6.2 课后习题详解6.3 名校考研真题详解第7章 无限长单位冲激响应(IIR)7.1 复习笔记7.2 课后习题详解7.3 名校考研真题详解第8章 有限长单位冲激响应(FIR)数字滤波器设计方法8.1复习笔记8.2 课后习题详解8.3 名校考研真题详解第9章 序列的抽取与插值——多抽样率数字信号处理基础9.1 复习笔记9.2 课后习题详解9.3 名校考研真题详解第10章 数字信号处理中的有限字长效应10.1 复习笔记10.2 课后习题详解10.3 名校考研真题详解第1章 离散时间信号与系统1.1 复习笔记一、离散时间信号——序列1.序列序列可以有三种表示法。
(1)函数表示法。
例如x(n)=a n u(n)。
(2)数列的表示法。
例如x(n)={...,-5,-3,-l,0,2,7,9,…)本书中,凡用数列表示序列时,都将n=0时x(o)的值用下划线(_)标注,这个例子中有z(-1)=-3,x(0)=-l,x(1)=0,…(3)用图形表示,如图l-1所示。
图1-1 离散时间信号的图形表示2.序列的运算(1)基于对序列幅度x(n)的运算序列的简单运算有①加法;②乘法;③累加;④序列绝对和;⑤序列的能量;⑥平均功率。
(2)基于对n的运算①移位,某序列为x(n)则x(n-m)就是x(n)的移位序列,当m=正数时,表示序列x(n)逐项依次右移(延时)m位;当m=负数时,表示序列 x(n)逐项依次左移(超前)m位;②翻褶,若序列为x(n),则x(-n)是以n=0为对称轴将x(n)序列加以翻褶;③时间尺度变换。
程佩青第四版答案【篇一:数字信号答案(第三版)程佩青 -需要的看看啊啊】数字信号处理教程课后习题及答案目录离散时间信号与系统 z变换离散傅立叶变换快速傅立叶变换数字滤波器的基本结构无限长单位冲激响应(iir)数字滤波器的设计方法有限长单位冲激响应(fir)数字滤波器的设计方法数字信号处理中有限字长效应第一章离散时间信号与系统1 .直接计算下面两个序列的卷积和y(n)?x(n)*h(n)h(n)???an , 0?n?n?1?0, 其他nn? x(n)????? n0,n0?n??0,n?n0请用公式表示。
分析:①注意卷积和公式中求和式中是哑变量m(n 看作参量),结果y(n)中变量是 n,??y(n)?x(m)h(n?m)??h(m)x(n?m) ; m????m???②分为四步(1)翻褶( -m ),(2)移位( n ),(3)相乘,(4)相加,求得一个 n 的 y(n) 值,如此可求得所有 n 值的y(n) ;③一定要注意某些题中在n 的不同时间段上求和范围的不同(3)?n0 ?n?1当n?n0?n?1n时 ,???n0?????n1?y(n)??x(m)h(n?m)m?n-n?1?n?1?n0??n?1?n0?,???nn?????m????m?n0?n?m?n0??m?n?n?1?m?n?n?1解:y(n)?x(n)*h(n)?m????x(m)h(n?m)y(n)?0?(1)(2)当n?n0时n当n0?n?n0?n?1时 ,部分重叠y(n)?nm?n0?x(m)h(n?m)m?n0?m?n0???n?m?n?n?m?n0???nmy(n)??n?n0?n?1?n0?,(???)(1)x(n)? ? (n),(2)x(n)? r3(n),如此题所示,因而要分段求解。
(3)x(n)? ? (n?2),(4)x(n)? 2nu(?n?1),h(n)?r5(n)h(n)?r4(n) h(n)?0.5nr3(n)h(n)?0.5nu(n)2 .已知线性移不变系统的输入为x(n),系统的单位抽样响应 ???n?n0??n?n?1????n?11???n?1?n?n0?n??n???,?????y(n)?n?n?n0,?????为h(n),试求系统的输出y(n),并画图。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
数字信号处理(程佩青)课后习题解答(1)题目1题目描述:给定一个长度为N的实数序列x(n),编写一段Python代码,实现求该序列的傅里叶变换,并将结果绘制成幅度谱和相位谱。
解答:import numpy as npimport matplotlib.pyplot as pltdef DFT(x):N = len(x)X = np.zeros(N, dtype=complex)for k in range(N):for n in range(N):X[k] += x[n] * np.exp(-2j * np.pi * k * n / N)return X# 输入实数序列x(n)x = [1, 2, 3, 4, 5]# 计算傅里叶变换X = DFT(x)# 计算幅度谱和相位谱X_mag = np.abs(X)X_phase = np.angle(X)# 绘制幅度谱plt.figure()plt.stem(range(len(X_mag)), X_mag)plt.xlabel('Frequency')plt.ylabel('Magnitude')plt.title('Magnitude Spectrum')# 绘制相位谱plt.figure()plt.stem(range(len(X_phase)), X_phase)plt.xlabel('Frequency')plt.ylabel('Phase')plt.title('Phase Spectrum')plt.show()题目2题目描述:根据以下信号的采样频率和采样点数,计算并绘制它的离散傅里叶变换(DFT)的幅度谱和相位谱。
采样频率:8kHz采样点数:1024信号:sin(2πf0n) + 0.5sin(2πf1n)其中,f0 = 500Hz,f1 = 1000Hz。
第四章 快速傅立叶变换运算需要多少时间。
计算需要多少时间,用,问直拉点的,用它来计算每次复加速度为平均每次复乘需如果一台通用计算机的FFT DFT[x(n)]512s 5 s50.1μμ解: 解: ⑴ 直接计算: 复乘所需时间:复加所需时间:⑵用FFT 计算: 复乘所需时间:复加所需时间:运算一次完成。
点试用一个为了提高运算效率值求今需要从值的点实序列是两个已知IFFTN n y n x k Y k X DFT n y n x N k Y k X ,,)(),()(),(,)(),()(),(.2sN T N01152.0512log105log105 2251262261=⨯⨯⨯=⨯⨯=--sT T T sN N T 013824.0 002304.0512log512105.0log105.0 2126262=+=∴=⨯⨯⨯=⨯⨯⨯=--sT T T sN N T 441536.1 130816.0)1512(512105.0)1(105.0 21662=+=∴=-⨯⨯⨯=-⨯⨯⨯=--sN T 31072.1 512105 105 26261=⨯⨯=⨯⨯=--值的过程。
)(),(完成计算点)可用一次()()(综上所述,构造序列)()()()(可得:)()()(再根据都是实序列,)(),(由原题可知:)()()()(()()(性质:又根据可得序列点作对取序列依据题意解 ]Im[ ]Re[ ][][ ][ ).()( )()()( )()();()( ::n y n x IFFT N k jY k X k Z n z n y n z n x n jy n x n z n y n x n jy n x k Y jIDFT k X IDFT k jY k X IDFT DFT n z IFFT N k Z k jY k X k Z k Y n y k X n x +===+=+=+=++=⇔⇔。
输出倒位序顺序频率抽取采用输入自然输出自然数顺序序时间抽取采用输入倒位流图抽取法的按时间抽取法及按频率画出基时), ,,( 2,16.3FFT N -=。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2x(m)()h n m -n1 1 1 0 0 0 0 y(n) 0 11 1 1 12 2 1 1 13 3 1 1 1 1 34 0 1 1 1 1 2 511111(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列, nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
第三章 离散傅立叶变换1.如下图,序列x(n)是周期为6的周期性序列,试求其傅立叶级数的系数。
∑∑=-===56265)(~)(~)(X ~:n nkj nkn e n x W n x k π解kj k j k j kj kj e e e e e 562462362262621068101214πππππ-----+++++=计算求得:。
339)5(~; 33)4(~ ; 0)3(~; 33)2(~;339)1(~;60)0(~j X j X X j X j X X +=-==+=-==。
并作图表示试求设)(~),(~)(~ .))(()(~),()(.264k X n x k X n x n x n R n x == ∑∑=-===56265)(~)(~)(~:n nk jnk n en x W n x k X π解k j kj k j e e e πππ---+++=3231。
计算求得: 3)5(~; 1)4(~ ; 0)3(~ ;1)2(~; 3)1(~ ; 4)0(~j X X X X j X X ====-==。
的周期卷积并作图与试求令其它,设 )(~)(~,))(()(~,))(()(~,)2()(,040,1)(.3464n h n x n h n h n x n x n R n h nn n n x ==-=⎩⎨⎧≤≤+=解:在一个周期内的计算值)(~)(~*)(~)(~m n h n h n x n y -==)(~)(~*)(~)(~m n h n h n x n y -==等各序列。
试画出所示如图已知)())((),())3((,))(()())((),())((,))((,13)(.47755633665n R n x n R n x n x n R n x n R n x n x P n x ----)()()5()(x(n)(4)N n 0 ),n -(n )()3()()()2()()(cos )()1()(52000n R n n x n nR n x n R a n x n R n a n x DFT N N N N n N ==<<===δω闭合形式表达式点试求以下有限长序列的])21sin()2sin()21sin()2sin([21])()()()([21)(]1111[)(][)(])([)()(cos )()()(cos )(:0)2(21020)2(2102)2(21)2(21)2(21222)2(21)2(21)2(21222)()(211)(10)(2110211000000000000000000002002002022002ϖπϖϖπϖωωϖπϖϖπϖϖπϖπϖπϖϖϖϖπϖπϖπϖϖϖωωωωωωωωππππππ-++⋅=--+--=--+--=+=+===---+---------+-++-----+---=---=+--=---=-∑∑∑∑k N e N e k N e N e a e ee e e e eeeee e a k R ee ee a k R eea k R e e e a k R en a k X n R n a n x k N j N j k Nj Njk Nj k Nj k Nj NjNjN jk Nj k N j k Nj NjNjNjN k j N j k j N j N N n nj N n nk j N N n nkj n j n j N n N nkj N N N N N N N 解)(111121)(21)()(21)()(cos )( )()(cos )( ) 1 (:)2()2(10)2(10)2(1020010200000k R e ee e a k R e e a k R e e e a k Re n a k X n R n a n x N k N j N j k Nj Nj N N n nN j N n nk N j N N n nk N j n j n j N n N nk N j N ⎥⎥⎦⎤⎢⎢⎣⎡--+--=⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡+===--+---=---=+--=---=-∑∑∑∑ωπωωπωωπωππωωπωω解⎥⎥⎥⎦⎤---⎢⎢⎢⎣⎡--=------+-++---)()()()(21)2(21)2(21)2(21222)2(21)2(21)2(212220000000ωπωπωπωωωωπωπωπωωωk Nj k N j k N j N jN jNjk N j k N j k N j Nj N j N j e e e eeee e e e e e a⎥⎥⎥⎥⎦⎤--⎢⎢⎢⎢⎣⎡+⋅=--+--)21sin()2sin()21sin()2sin(210)2(21020)2(21020000ωπωωπωωπωωπωk N e Nek N e Ne a k Nj Njk N j N jk Nj N N n nk NjnN n aea eak X n R a n x ππ210211)()()((2)--=---===∑)( )()( )()()( 0,)()( (3)02102010200k Re k Re n n k Re n x k X N n n n n x Nk n N j NN n nk N j N n NnkN j πππδδ--=--=-=-==<<-=∑∑)(1)( 11)1()())1(()(])1)2( 2[)1( 32()1)(()()()()( )()(411)1(32)1(321)1(110)1(1k R W Nk X N W W N k R W N k R N W N W W W N W W W nW nWW k X k R nW k X W k R nW k X n nR n x N kNkNkN N N n nk N N k N N k N k N k N N kN k N k N N n kn N N n nk Nk NN n N k n N k NN n N nkN N --=∴-=--+--=+--=-+-+++--++++=-=-==∴=∑∑∑∑∑-=---=+-=-=+-=∙∙∙∙∙∙)(kNN N n nkNN W Nk X n nR n x W n k X n R n n x --===∴=∑-=1)()()()4()( )()(5111022,则小题的结论根据第)(221111122)1(232)1(23210)1(2121)1(2)1()2()(12)2()(2)2(2)2()12()1(]1()2(4[)1(94)1)(()(k N kN kNN n nk NN n nkNk N N kN k N k N N k N k N k N N n kn NN n nk NkNN n kn N k NW N W N N k X W NN N k X N N nWN N W n N N W N W W W N W W W W n W nW k X W n k X W ---=∴----=+--=+--=-+--=-+-+++--++++=-=-=∑∑∑∑∑-=-=---=+-=-=+∙∙∙∙∙∙)∙∙∙±±±===∑-=,6,4,2,0)(~)3(?])0([)()2(?)()1(:;)(~1)(~).(~.61)/2(k k x X k X k X ek X Nn x n x N k nkN j 哪些序列列能做到成虚数外除时间原点使所有的哪些序列能够通过选择成为实数时间原点使所有的哪些序列能够通过选择问傅里叶级数这些序列可以表示成列如图画出了几个周期序π条件。
第七章 有限长单位冲激响应(FIR )数字滤波器的设计方法1. 用矩形窗设计一个FIR 线性相位低通数字滤波器。
已知 21,5.0==N c πω。
求出)(n h 并画出)(log 20ωj e H 曲线。
分析:此题给定的是理想线性相位低通滤波器,故⎪⎩⎪⎨⎧<<<<≤≤=-。
-- , , 0- , )(c c c c ωωππωωωωωωαωj j d e eH解:ωπππωωd eeH n h nj j d d ⎰-=)(21)()()](sin[21αωαωπωωπωωωωα--==⎰--n n d eec c c nj j cc⎪⎪⎩⎪⎪⎨⎧≤≤--====-=为其他故:其中n n n n n w n h n h N d c ,0200,)10(]2sin[)()()(5.0 102/)1( πππωα h( 0)= 9.7654073033E-4h( 1)= 3.5358760506E-2 h( 2)= -9.7657600418E-4 h( 3)= -4.5465879142E-2 h( 4)= 9.7651791293E-4 h( 5)= 6.3656955957E-2 h( 6)= -9.7658322193E-4 h( 7)= -1.0610036552E-1 h( 8)= 9.7643269692E-4 h( 9)= 3.1830877066E-1 h( 10)= 4.9902343750E-1 h( 11)= 3.1830900908E-1 h( 12)= 9.7669276875E-4 h( 13)= -1.0610023141E-1 h( 14)= -9.7654142883E-4 h( 15)= 6.3657015562E-2 h( 16)= 9.7660662141E-4 h( 17)= -4.5465819538E-2 h( 18)= -9.7654841375E-4 h( 19)= 3.5358794034E-2 h( 20)= 9.7658403683E-42.用三角形窗设计一个FIR 线性相位低通数字滤波器。