数字信号处理2

一、填空题

1.若初始信号的频带宽度是有限的，经过抽样后能够不失真的还原出初始信号，那

么抽样频率应该满足什么样的数量关系？

2.序列，则的收敛域为 。

)1()(---=n u a n x n )(Z X 3.一个线性移不变因果系统稳定的充分必要条件是 。

4.序列的傅里叶变换的定义即正变换的表达式为 ，离散傅里叶变换的

x(n)定义即其正变换的表达式为 ，反变换的表达式为 。

5.已知序列：x(n),0≤n≤15; g(n),0≤n≤19. X(k)、G(k)分别是它们的32 点DFT.

令y(n)=IDFT[X(k)G(k)],0≤n≤31,则y(n)中相等于x(n)与g(n)线性卷积中的有 点，其序号从 到 。

6.利用DFT 计算连续时间信号（模拟信号）时可能出现的问题 、 、

、 。

7.快速傅里叶变换算法可以分为两大类，即为 和 。

8.无限长单位冲激响应滤波器有 、 、 、

四种基本网络结构。

9.数字滤波器按照频率特性分可以划分为 、 、 、 、 等类型。

10.抽样序列在单位圆上的z 变换，就等于其理想抽样信号的 。

二、判断题

1、模拟信号也可以与数字信号一样在计算机上进行数字信号处理，只要加一道采样的工序就

可以了。

2、已知某离散时间系统为，则该系统为线性时不变系统。 )35()]([)(+==n x n x T n y

3、一个信号序列，如果能做序列的傅里叶变换（），也就能对其做变换。DTFT DFT

4、用双线性变换法进行设计数字滤波器时，预畸并不能消除变换中产生的所有频率点的IIR 非线性畸变。

5、在IIR 数字滤波器的设计中，用脉冲响应不变法设计时，从模拟角频率向数字角频率转换时，转换关系是线性的。6． 在时域对连续信号进行抽样，在频域中，所得频谱是原信号频谱的周期延拓。7、x(n)=cos （w 0n)所代表的序列一定是周期的。8、y(n)=x 2(n)+3所代表的系统是时不变系统。

9、 用窗函数法设计FIR 数字滤波器时，改变窗函数的类型可以改变过渡带的宽度。10、阻带最小衰耗取决于窗谱主瓣幅度峰值与第一旁瓣幅度峰值之比。

三、简答题

1. 试写出利用FFT 计算IFFT 的步骤

行检查术资中资料

2.画出按时域抽取N=4点基2FFT 的信号流图3、用DFT 对连续信号进行谱分析的误差问题有哪些？ 4、简要分析FFT 的计算量和算法特点。四、证明题

设表示长度为N 的有限长序列的DFT 。)(k X )(n x （1）

证明如果满足关系式)(n x )1()(n N x n x ---=则

0)0(=X （2）证明当N 为偶数时，如果 )1()(n N x n x --=则0)2(

=N X 五、实验题已知A=[1 2 3 4],B=[3 4 5 6].C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B.求出C 、D 、E 、F 、G ，并用绘图函数stem 绘出A 、B 、C 、D 、E 、F 、G 的序列图。请写出其matlab 程序六、设计题系统函数为

，

H (z )=1s +1?1s +3试用冲激响应不变法，设计IIR 数字滤波器.

参考答案

一、填空题

1.若初始信号的频带宽度是有限的，经过抽样后能够不失真的还原出初始信号，那么抽样频率应该满足什么样的数量关系？ ,即抽样频率大于或等于两倍信号谱的f s ≥2f ?最高频率

2序列，则的收敛域为 。

)1()(---=n u a n x n )(Z X a Z ≤3.一个线性移不变因果系统稳定的充分必要条件是 系统函数的全部极点必须在单位圆内 。4.序列的傅里叶变换的定义即正变换的表达式为

x(n)X （e jω）=∑∞n =?∞x(n)e ?jωn ，离散傅里叶变换的定义即其正变换的表达式为

，反变换X (k )=∑N ?1n =0x(n)e ?j 2πN nk 的表达式为 。x (n )=1N ∑N ?1k =0X(k)e j 2πN nk 5.已知序列：x(n),0≤n≤15; g(n),0≤n≤19. X(k)、G(k)分别是它们的32 点DFT.令y(n)=IDFT[X(k)G(k)],0≤n≤31,则y(n)中相等于x(n)与g(n)线性卷积中的有 29点，其序号从 3 到 31 。6.利用DFT 计算连续时间信号（模拟信号）时可能出现的问题 频率相应的混叠失真 、频谱泄露 、 栅栏效应 、频率分辨力 。7.快速傅里叶变换算法可以分为两大类，即为 按时间抽取法 和 按频率抽取法 。8.无限长单位冲激响应滤波器有 直接Ⅰ型 、 直接Ⅱ型 、 级联型 、 并联型 四种基本网络结构。9.数字滤波器按照频率特性分可以划分为低通 、高通 、带通 、带阻 、全通 等类型。10.抽样序列在单位圆上的z 变换，就等于其理想抽样信号的傅里叶变换。二、判断题5、模拟信号也可以与数字信号一样在计算机上进行数字信号处理，只要加一道采样的工序就可以了。 （╳）6、已知某离散时间系统为，则该系统为线性时不变系统。(╳)

)35()]([)(+==n x n x T n y 7、一个信号序列，如果能做序列的傅里叶变换（），也就能对其做变换。（╳）

DTFT DFT 8、用双线性变换法进行设计数字滤波器时，预畸并不能消除变换中产生的所有频率点的

IIR 非线性畸变。 （√）

5、在IIR 数字滤波器的设计中，用脉冲响应不变法设计时，从模拟角频率向数字角频率转换时，转换关系是线性的。（ √ ）

6． 在时域对连续信号进行抽样，在频域中，所得频谱是原信号频谱的周期延拓。（ √ ）接地线课件断习题装置高中资具高中处理，因此电源高

7、x(n)=cos （w 0n)所代表的序列一定是周期的。（　×　）8、y(n)=x 2(n)+3所代表的系统是时不变系统。 （　√ ）9、 用窗函数法设计FIR 数字滤波器时，改变窗函数的类型可以改变过渡带的宽度。（ √ ）

10、阻带最小衰耗取决于窗谱主瓣幅度峰值与第一旁瓣幅度峰值之比。 （╳）三、简答题1. 试写出利用FFT 计算IFFT 的步骤1）对X(k)取共轭，得 (k);2）对 (k)做N 点FFT ；3）对2）的结果取共轭并除以N 2.画出按时域抽取N=4点基2FFT

的信号流图3、用DFT 对连续信号进行谱分析的误差问题有哪些？ 答：混叠失真；频谱泄漏；栅栏效应。4、简要分析FFT 的计算量和算法特点。答：计算量：复数乘法 次，复数加法 次。2log 22N N N N 2log 算法特点：（1）以碟形运算为基础进行组合计算，计算因子W K 的指数K 与运算所在的级数和组内位置有关。（2）中间数据的存储，可采用原位存储法。即每次碟形运算的结果可以存储在原数据的同一个存储单元。这样在高速硬件实现时，可节省存储器。（3）输入序列的混序。因为DFT 输入序列是顺序采样的，所以在计算FFT 之前需要进行序列按混序要求排序。排序算法很多，较常用的计算混序号的方法有二进制序号反转算法。四、证明题设表示长度为N 的有限长序列的DFT 。)(k X )(n x （3）证明如果满足关系式)(n x )

1()(n N x n x ---=则0)0(=X

（4）证明当N 为偶数时，如果 )1()(n N x n x --=则0)2

(=N X 解 （1）∑∑∑∑∑-=-=-=-=-=---====121201010010)1()()()()0()()(N N n N n N n N n N N n nk N n N x n x n x W n x X W n x k X 令m n N =--1∑∑-=-=-=012120)()()0(N n N n m x n x X 显然可得 0)0(=X （2） （将n 分为奇数和偶数两部分表示）

∑∑-=-=-==1010)1)(()()2(N n n N n jk n x e n x N X π ∑∑-=+-=-++-=

120121202)1)(12()1)(2(N r r N r r r x r x ∑∑-=-=+-=120120)12()2(N r N r r x r x ()1221)12()21(120120+=--+---=∑∑-=-=k r N r x r N x N r N r 令 ∑∑-==+-+=1200

2)12()12(N r N k r x r x 显然可得 02(=N X 五、实验题

已知A=[1 2 3 4],B=[3 4 5 6].C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B.

求出C 、D 、E 、F 、G ，并用绘图函数stem 绘出A 、B 、C 、D 、E 、F 、G 的序列图。请写出其

matlab 程序

Close all;Clear all;clc;A=[1 2 3 4];B=[3 4 5 6];C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B;subplot(2,4,1);stem(A,'.');subplot(2,4,2);stem(B,'.');subplot(2,4,3);stem(C,'.');subplot(2,4,4);stem(D,'.');subplot(2,4,5);stem(E,'.');subplot(2,4,6);stem(F,'.');subplot(2,4,7);stem(G,'.');

六、设计题

系统函数为 ，

H (z )=1s +1?1s +3试用冲激响应不变法，设计IIR 数字滤波器.

H (z )=T 1?z ?1e ?T ?T 1?z ?1e ?3T

设T=1，则有H (z )=

0.318z ?1

1?0.4177z ?1+0.01831z ?2

壁薄、开处理卷电气要在事力高中压器组

数字信号处理实验作业

实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。

数字信号处理答案第二章

第二章 2.1 判断下列序列是否是周期序列。若是，请确定它的最小周期。 （1）x(n)=Acos(6 85ππ+n ) （2）x(n)=)8( π-n e j （3）x(n)=Asin(3 43π π+n ) 解 (1)对照正弦型序列的一般公式x(n)=Acos(?ω+n )，得出=ω85π。因此5 16 2=ωπ是有理数，所以是周期序列。最小周期等于N= )5(165 16 取k k =。 （2）对照复指数序列的一般公式x(n)=exp[ωσj +]n,得出81=ω。因此πωπ162=是无理数，所以不是周期序列。 （3）对照正弦型序列的一般公式x(n)=Acos(?ω+n )，又x(n)=Asin(343ππ+n )＝Acos(-2π3 43π π-n )＝Acos(6143-n π)，得出=ω43π。因此3 8 2=ωπ是有理数，所以是周期序列。最小周期等于 N= )3(83 8 取k k = 2.2在图2.2中，x(n)和h(n)分别是线性非移变系统的输入和单位取样响应。计算并列的x(n)和h(n)的线性卷积以得到系统的输出y(n)，并画出y(n)的图形。 (a) 1 11 1 (b) (c) 11 1 11 0 0 -1-1 -1 -1 -1 -1 -1 -1 2 2 2 22 2 3 3 3 3 34 44 … … …n n n n n n x(n)x(n) x(n) h(n)h(n) h(n)2 1 u(n) u(n) u(n)a n ===2 2

解利用线性卷积公式 y(n)=∑∞ -∞ =- k k n h k x) ( ) ( 按照折叠、移位、相乘、相加、的作图方法，计算y(n)的每一个取样值。 (a) y(0)=x(O)h(0)=1 y(l)=x(O)h(1)+x(1)h(O)=3 y(n)=x(O)h(n)+x(1)h(n-1)+x(2)h(n-2)=4,n≥2 (b) x(n)=2δ(n)-δ(n-1) h(n)=-δ(n)+2δ(n-1)+ δ(n-2) y(n)=-2δ(n)+5δ(n-1)= δ(n-3) (c) y(n)= ∑∞ -∞ =-- k k n k n u k u a) ( ) (=∑∞ -∞ = - k k n a= a a n - -+ 1 11 u(n) 2.3 计算线性线性卷积 (1) y(n)=u(n)*u(n) (2) y(n)=λn u(n)*u(n) 解：(1) y(n)= ∑ ∞ -∞ = - k k n u k u) ( ) ( =∑ ∞ = - ) ( ) ( k k n u k u=(n+1),n≥0 即y(n)=(n+1)u(n) (2) y(n)=∑∞ -∞ = - k k k n u k u) ( ) ( λ

数字信号处理实验一

一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.');

xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid;

数字信号处理作业答案

数字信号处理作业

DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~1k X 是周期性的，周期为N ，而)(~2k X 也是周期性的，周期为N 2。试利用)(~1k X 确定)(~2k X 。（76-4）

2. 研究两个周期序列)(~n x 和)(~n y 。)(~n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~k W 。（76-5）

3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000) ()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）

2020年数字信号处理大作业新版修订

2019~2020年度《数字信号处理》大作业题目与要求 大作业要求： 本学期大作业总分40分，学生可选择任意数量的题目完成，只要所选题目总分达到40分即可，所选题目总分如果超过40分，超过的部分不计入大作业总分。大作业以电子版的形式提交，内容应包括详细的程序设计思路与题目分析（题目分析指的是对该题目中所用到的知识点的说明，不要照搬书上或网上的内容，写出你自己对该知识点的理解。），程序截图，程序源码，其中设计思路和程序截图可写在同一个文档中，程序源码可以是.txt或.m 文件，并在源码中标注代码注释。另：题目中有GUI设计要求的部分占该题目分值的20%，功能实现部分占该题目分值的80%。 注：以下题目均用MATLAB完成。 大作业题目： 1、实现有限长序列的基本运算（包括：加法、乘法、累加、移位、翻褶、抽取、插值、卷积和），并以GUI的形式将这些运算整合起来，使用者可通过向GUI输入任意有限长序列得到对应的运算结果。（5分） 2、设计一个GUI，实现奈奎斯特采样定理，要求：1、在GUI中输入任意一个模拟信号，显示该模拟信号的时域和频域谱图；2、在GUI中设置任意采样频率，对输入的模拟信号进行采样处理，显示采样信号的时域和频域谱图； 3、在GUI中实现采样信号向模拟信号的恢复功能，要求显示恢复后的模拟信号的时域和频域谱图。（10分） 3、通过GUI动态展示z变换与s变换之间的所有关系。（5分） 4、设计一个GUI，通过向GUI输入任意系统函数，得到其对应系统的相关信息（包括：系统频率响应中的幅度响应和相位响应、系统零极点的分布、系统的稳定性判定）。（10分） 5、设计一个GUI，实现利用DFT（或FFT）完成任意时域信号的频谱分析，要求：1、可在GUI中输入时域数字或模拟信号；2、可设置DFT点数；3、在GUI中显示输入信号经DFT（或FFT）处理后的频谱图；3、若输入信号为模拟信号，需完成对该模拟信号的采样，采样频率可在GUI中设置。（10分） 6、在GUI中，实现IIR滤波器的直接型、级联型和并联型三种结构之间的任意转换，要求：在GUI中输入任意一型的系统函数后可在该GUI中显示出对应的另外两型的系统函数。（10分） 7、实现巴特沃斯样本模拟低通滤波器及其对应的数字低通滤波器的设计，以GUI的形式给出。要求：输入所需的模拟低通滤波器参数指标后，程序能将该指标转化为数字低通滤波器指标（在GUI中应能选择转化方式：冲激响应不变法、双线性变换法），并在GUI中显示出所给参数下巴特沃斯样本模拟低通滤波器及其对应的数字低通滤波器的频率响应中幅度响应的频谱图。（15分） 8、已知某组数字信号（见大作业数据压缩包中HWDATA.mat文件），该信号中除了目标信号之外还掺杂有强噪声，但噪声与目标信号的频率不重叠，要求采用本学期已学的知识对该信

数字信号处理实验二

实验报告（本科） 学号 2015141443002 姓名柏冲 专业通信工程 日期 2017/12/4 实验题目时域采样和频域采样 一、实验目的

时域采样理论与频域采样理论是数字信号处理中重要的理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使得采样后的信号不丢失信息；要求掌握频率采样会引起时域周期化的概念，以及频域采样定理及其对频域采样点数选择的指导作用。 二、实验过程 附：源程序 (1)时域采样 Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=1000; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); %M点FFT[(xnt)] subplot(3,2,1); stem(xnt,'.'); %调用编绘图函数stem绘制序列图 box on;title('(a) Fs=1000Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,2);stem(fk,abs(Xk),'.');title('(a) T*FT[xa(nT)],Fs=1000Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]); % Fs=300Hz和 Fs=200Hz的程序与上面Fs=1000Hz完全相同。 Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=300; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M1=fix(M); Xk=T*fft(xnt,M1); %M点FFT[(xnt)] subplot(3,2,3); stem(xnt,'.'); %调用自编绘图函数stem绘制序列图 box on;title('(b) Fs=300Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,4);stem(fk,abs(Xk),'.');title('(b) T*FT[xa(nT)],Fs=300Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]); Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=200; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M2=fix(M);

数字信号处理第二章上机题作业

数字信号处理作业实验题报告 第一章16.(1) 实验目的： 求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。 实验要求： 运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。 源程序： B1=1;A1=[1, -0.6, 0.08]; ys=2; %设差分方程 xn=[1, zeros(1, 20)]; %xn=单位脉冲序列，长度N=31 xi=filtic(B1, A1, ys); hn1=filter(B1, A1, xn, xi); %求系统输出信号hn1 n=0:length(hn1)-1; subplot(2, 1, 1);stem(n, hn1, '.') title('单位脉冲响应'); xlabel('n');ylabel('h(n)') xn=ones(1, 20); sn1=filter(B1, A1, xn, xi); %求系统输出信号sn1 n=0:length(sn1)-1; Subplot(2, 1, 2); stem(n, sn1, '.') title('单位阶跃响应'); xlabel('n'); ylabel('s(n)')

运行结果： 实验分析： 单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。

第二章31题 实验目的： 用matlab判断系统是否稳定。 实验要求： 用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。 源程序： A=[2, -2.98, 0.17, 2.3418, -1.5147]; B=[0, 0, 1, 5, -50]; subplot(2,1,1); zplane(B,A); %求H(z)的极点 p=roots(A); %求H(z)的模 pm=abs(p); if max(pm)<1 disp('系统因果稳定'), else,disp('系统因果不稳定'),end un=ones(1,800); sn=filter(B, A, un); n=0:length(sn)-1; subplot(2, 1, 2);plot(n, sn) xlabel('n');ylabel('s(n)')

西电数字信号处理大作业

第二章 2.25 已知线性时不变系统的差分方程为 若系统的输入序列x(x)={1,2,3,4,2,1}编写利用递推法计算系统零状态响应的MATLAB程序，并计算出结果。 代码及运行结果： >> A=[1,-0.5]; >> B=[1,0,2]; >> n=0:5; >> xn=[1,2,3,4,2,1]; >> zx=[0,0,0];zy=0; >> zi=filtic(B,A,zy,zx); >> yn=filter(B,A,xn,zi); >> figure(1) >> stem(n,yn,'.'); >> grid on;

2.28图所示系统是由四个子系统T1、T2、T3和T4组成的，分别用单位脉冲响应或差分方程描述为 T1： 其他 T2: 其他 T3： T4: 编写计算整个系统的单位脉冲响应h(n)，0≤n≤99的MATLAB程序，并计算结果。 代码及结果如下： >> a=0.25;b=0.5;c=0.25; >> ys=0; >> xn=[1,zeros(1,99)]; >> B=[a,b,c]; >> A=1; >> xi=filtic(B,A,ys); >> yn1=filter(B,A,xn,xi); >> h1=[1,1/2,1/4,1/8,1/16,1/32]; >> h2=[1,1,1,1,1,1]; >> h3=conv(h1,h2); >> h31=[h3,zeros(1,89)]; >> yn2=yn1+h31; >> D=[1,1];C=[1,-0.9,0.81]; >> xi2=filtic(D,C,yn2,xi); >> xi2=filtic(D,C,ys); >> yn=filter(D,C,yn2,xi); >> n=0:99; >> figure(1) >> stem(n,yn,'.'); >> title('单位脉冲响应'); >> xlabel('n');ylabel('yn');

数字信号处理上机作业

数字信号处理上机作业 学院：电子工程学院 班级：021215 组员：

实验一：信号、系统及系统响应 1、实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 (2) 熟悉时域离散系统的时域特性。 (3) 利用卷积方法观察分析系统的时域特性。 (4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 2、实验原理与方法 (1) 时域采样。 (2) LTI系统的输入输出关系。 3、实验内容及步骤 (1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。 (2) 编制实验用主程序及相应子程序。 ①信号产生子程序，用于产生实验中要用到的下列信号序列： a. xa(t)=A*e^-at *sin(Ω0t)u(t) b. 单位脉冲序列：xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③有限长序列线性卷积子程序 用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数conv。 conv 用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。调用格式如下： y=conv (x, h) 4、实验结果分析 ①分析采样序列的特性。 a. 取采样频率fs=1 kHz,，即T=1 ms。 b. 改变采样频率，fs=300 Hz，观察|X(e^jω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的|X(e^j ω)|曲线。 程序代码如下： close all;clear all;clc; A=50; a=50*sqrt(2)*pi; m=50*sqrt(2)*pi; fs1=1000; fs2=300; fs3=200; T1=1/fs1; T2=1/fs2; T3=1/fs3; N=100;

数字信号处理_吴镇扬_第二版_第三章习题答案

∞ ∑ 3.3 x %1(n) = x (n + 8r) 说明x%1(n)的周期是8 r=?∞ 2π ? j π kn N ?1 7 7 ? j kn = e ∴X% 1(k) = x %1(n)W Nkn = e ∑ ∑ ∑ 8 4 n=0 n=0 n=0 ? j π 0×n 7 X 1(0) = e ∑ 4 = 3 n=0 ? j π n 7 X 1(1) = e ∑ 4 = (1+ 2 / 2) (1? j ) n=0 ? j π 2n ? j π 3n 7 X 1(2) = e 7 = ? j ;X 1(3) = e ∑ ∑ 4 4 = (1? 2 / 2) (1+ j ) n=0 n=0 ? j π 4n ? j π 5n 7 X 1(4) = e 7 =1;X 1(5) = e ∑ ∑ 4 4 =(1? 2 / 2) (1? j ) = (1+ 2 / 2) (1+ j ) n=0 n=0 ? j π 6n ? j π 7n 7 X 1(6) = e 7 = j ;X 1(7) = e ∑ ∑ 4 4 n=0 n=0

3.6 解：(1)x(n) =δ(n ) N ?1 根据定义有：X(k) = x (n)W = δ(n )W Nkn =1 N ?1 ∑ ∑ kn N n=0 n=0 (2)x (n) =δ(n ?n 0),0 < n 0 < N N ?1 X (k ) = x (n)W = δ (n ? n 0)W N kn =W N kn N ?1 ∑ ∑ kn 0< n 0 < N 0 N n=0 n=0 (3)x (n) = a n 0 < n < N ?1 1?(a W N k )N 1?a N 1?(a W N k ) 1?aW N k N ?1 X (k ) = x (n)W N kn = ∑ = n=0

数字信号处理作业+答案讲解

数字信号处理作业 哈尔滨工业大学 2006.10

DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~ 2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~ 1k X 是周期性的，周期为N ，而)(~ 2k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2k X 。（76-4）

2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列 )(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~ k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~ n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~ k X 和)(~ k Y 求)(~ k W 。（76-5）

3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）

数字信号处理作业-答案

数字信号处理作业-答案

数字信号处理作业

DFT 习题 1. 如果)(~ n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1 k X 表示)(~ n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2 k X 表示)(~ n x 的离散傅里叶级数之系数。当然，)(~ 1 k X 是周期性的，周期为N ，而)(~ 2 k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2 k X 。（76-4）

2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~ n w 定义为)()()(~~ ~ n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~ n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地，由于)(~ n y 的周期为M ，其离散傅里叶级数之系数)(~ k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~ k W 。（76-5）

3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=0 0)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）

数字信号处理 答案 第二章(精编文档).doc

【最新整理，下载后即可编辑】 第二章 2.1 判断下列序列是否是周期序列。若是，请确定它的最小周期。 （1）x(n)=Acos(6 85ππ+n ) （2）x(n)=)8( π-n e j （3）x(n)=Asin(3 43π π+ n ) 解 (1)对照正弦型序列的一般公式x(n)=Acos(?ω+n )，得出= ω8 5π 。因此 5 16 2= ωπ 是有理数，所以是周期序列。最小周期等于N= )5(165 16 取k k =。 （2）对照复指数序列的一般公式 x(n)=exp[ωσj +]n,得出8 1 =ω。因此 πω π 162=是无理数，所以不是周期序列。 （3）对照正弦型序列的一般公式x(n)=Acos(?ω+n )，又x(n)=Asin(3 43ππ+n ) ＝Acos( -2π343ππ-n )＝Acos(6143-n π)，得出= ω43π。因此3 8 2=ωπ是有理数，所以是周期序列。最小周期等于N=)3(83 8 取k k = 2.2在图2.2中，x(n)和h(n)分别是线性非移变系统的输入和单位取样响应。计算并列的x(n)和h(n)的线性卷积以得到系统的输出y(n)，并画出y(n)的图形。

(a) 1 11 1 (b) (c) 11 1 11 0 0 -1 -1 -1 -1 -1 -1 -1 2 2 2 22 2 3 3 3 34 44 … … …n n n n n n x(n)x(n) x(n) h(n)h(n) h(n)2 1 u(n) u(n) u(n)a n ===2 2 解 利用线性卷积公式 y(n)=∑∞ -∞=-k k n h k x )()( 按照折叠、移位、相乘、相加、的作图方法，计算y(n)的每一个取样值。 (a) y(0)=x(O)h(0)=1 y(l)=x(O)h(1)+x(1)h(O)=3 y(n)=x(O)h(n)+x(1)h(n-1)+x(2)h(n-2)=4,n ≥2 (b) x(n)=2δ(n)-δ(n-1) h(n)=-δ(n)+2δ(n-1)+ δ(n-2) y(n)=-2δ(n)+5δ(n-1)= δ(n-3) (c) y(n)= ∑∞ -∞=--k k n k n u k u a )()(= ∑∞ -∞ =-k k n a = a a n --+111u(n)

数字信号处理作业-2012

《数字信号处理Ⅰ》作业 姓名： 学号： 学院： 2012 年春季学期

第一章 时域离散信号和时域离散系统 月 日 一 、判断： 1、数字信号处理和模拟信号处理在方法上是一样的。（ ） 2、如果信号的取值和自变量都离散，则称其为模拟信号。（ ） 3、如果信号的取值和自变量都离散，则称其为数字信号。（ ） 4、时域离散信号就是数字信号。（ ） 5、正弦序列都是周期的。（ ） 6、序列)n (h )n (x 和的长度分别为N 和M 时，则)n (h )n (x *的长度为N+M 。（ ） 7、如果离散系统的单位取样响应绝对可和，则该系统稳定。（ ） 8、若满足采样定理，则理想采样信号的频谱是原模拟信号频谱以s Ω（采样频率）为周期进行周期延拓的结果。（ ） 9、序列)n (h )n (x 和的元素个数分别为21n n 和，则)n (h )n (x *有（1n n 21-+）个元素。（ ） 二、选择 1、R N (n)和u(n)的关系为（ ）： A. R N (n)=u(n)-u(n-N) B. R N (n)=u(n)+u(n-N) C. R N (n)=u(n)-u(n-N-1) D. R N (n)=u(n)-u(n-N+1) 2、若f(n)和h(n)的长度为别为N 、M ，则f(n)*h(n)的长度为 （ ）： A.N+M B.N+M-1 C.N-M D.N-M+1 3、若模拟信号的频率范围为[0，1kHz],对其采样，则奈奎斯特速率为（ ）： A.4kHz B. 3kHz C.2kHz D.1kHz 4、LTIS 的零状态响应等于激励信号和单位序列响应的（ ）： A.相乘 B. 相加 C.相减 D.卷积 5、线性系统需满足的条件是（ ）： A.因果性 B.稳定性 C.齐次性和叠加性 D.时不变性 6、系统y(n)=f(n)+2f(n-1)(初始状态为0)是（ ）： A. 线性时不变系统 B. 非线性时不变系统 C. 线性时变系统 D. 非线性时变系统

长沙理工数字信号处理大作业数字滤波器设计

IIR及FIR数字滤波器 一题干 对模拟信号进行低通滤波处理，要求通带0≤f≤4kHz，通带衰减小于0.5dB，阻带4.5k Hz≤f<∞，阻带衰减大于50dB，设采样频率Fs=20kHz。 (1)设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。 (2)分别用脉冲响应不变法和双线性变换法设计IIR低通数字滤波器，求出Ha(z) 的分子、分母多项式系数Bz和Az，并画出幅频响应损耗函数曲线 (3)采用窗函数法（分别用汉宁窗、哈明窗、布莱克曼窗函数）设计满足要求的FIR 低通滤波器，求出h(n)，并画出幅频响应损耗函数曲线. (4)用频率采样法设计满足要求的FIR低通滤波器，求出h(n),并画出幅频响应损耗函数曲线。

二求解过程 具体内容如下: （1）设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。 程序： wp=2*pi*4000; ws=2*pi*5800; Rp=0.5; As=50; [N,wc]=buttord(wp,ws,Rp,As,'s'); [B,A]=butter(N,wc,'s'); k=0:511; fk=0:20000/512:20000; wk=2*pi*fk; Hk=freqs(B,A,wk); plot(fk/1000,20*log10(abs(Hk))); grid on xlabel('频率/kHz'); ylabel('幅度/dB'); axis([0,6,-65,5]); 波形图：

A = 1.0e+207 * 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0020 2.1576 B = 1.0e+207 * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.1576 N = 46

数字信号处理第三章作业.pdf

数字信号处理第三章作业 1．(第三章习题3)在图P3-2中表示了两个周期都为6的周期性序列，确定这个两个序列的周期卷积的结果3()x n ，并画出草图。 2．(第三章习题5)如果()x n 是一个具有周期为N 的周期性序列，它也是具有周期为2N 的周期性序列。令~1()X k 表示当()x n 看做是具有周期为N 的周期性序列的DFS 系数。而~2()X k 表示当()x n 看作是具有周期为2N 的周期性序列的DFS 系数。当然~1()X k 是具有周期为N 的周期性序列，而~2()X k 是具有周期为2N 的周期性序列，试根据~1()X k 确定~2()X k 。 3．(第三章习题6) （a ）试证明下面列出的周期性序列离散傅里叶级数的对称特性。在证明中，可以利用离散傅里叶级数的定义及任何前面的性质，例如在证明性质③时可以利用性质①和②。 序列 离散傅里叶级数 ① *()x n ~*()X k - ②*()x n - ~*()X k ③Re ()x n ???? ~ e ()X k ④Im ()j x n ???? ~()o X k

（b ）根据已在（a ）部分证明的性质，证明对于实数周期序列()x n ，离散傅里叶级数的下列对称性质成立。 ①~~Re ()Re ()X k X k ????=-???????? ②~~Im ()Im ()X k X k ????=--???????? ③~~()()X k X k =- ④~~arg ()arg ()X k X k ????=--???????? 4．(第三章习题7)求下列序列的DFT (a) {}11 1-，，，-1 (b) {}1 j 1j -，，，- (c) ()cn 0n 1x n N =≤≤-， (d) 2n ()sin 0n 1x n N N π??=≤≤- ??? ， 5．(第三章习题8)计算下列各有限长序列的离散傅立叶变换（假设长度为N ） 1 0)()(0) ()()() ()()(00-≤≤=<<-==N n a n x c N n n n n x b n n x a n δδ 6．(第三章习题9)在图P3-4中表示了一有限长序列)(n x ，画出序列)(1n x 和)(2n x 的草图。（注意：)(1n x 是)(n x 圆周移位两个点） )())(()() ())2(()(442441n R n x n x n R n x n x -=-=

DSP大作业(哈工程)

DSP原理与应用 学号： 姓名： 日期：2017年5月23日星期二

1.DSP的生产厂商主要有哪些？分别有什么系列？ 答： ①德州仪器公司（最有名的DSP芯片厂商）。TI公司在市场上主要的三个系 列产品： （1）面向数字控制、运动控制的TMS320C2000系列，主要包括TMS320C24x/F24x、TMS320LC240x/LF240x、TMS320C24xA/LF240xA、TMS320C28xx等； （2）面向低功耗、手持设备、无线终端应用的TMS320C5000系列，主要包括TMS320C54x、TMS320C54xx、TMS320C55x等； （3）面向高性能、多功能、复杂应用领域的TMS320C6000系列，主要包括TMS320C62xx、TMS320C64xx、TMS320C67xx等。 ②美国模拟器件公司。其主要的系列： （1）定点DSP芯片有ADSP2101/2103/2105、ADSP2111/2115、ADSP2126/2162/2164、ADSP2127/2181、ADSP-BF532以及Blackfin系列； （2）浮点DSP芯片有ADSP21000/21020、ADSP21060/21062，以及虎鲨TS101、TS201S。 ③Motorola公司（发布较晚）。其主要的系列包括： （1）定点DSP 处理器MC56001； （2）与IEEE浮点格式兼容的的浮点DSP芯片MC96002； （3）DSP53611、16位DSP56800、24位的DSP563XX和MSC8101等产品。 ④杰尔公司。主要系列有： 嵌入式DSP内核的SC1000和SC2000系列，主要面向电信基础设施、移动通信、多媒体服务器及其它新兴应用。 2.浮点DSP和定点DSP各自有什么特点？ 答： 浮点DSP和定点DSP在宏观上有很大的特点区别，包括动态范围、速度、价格等等。 （1）动态范围：定点DSP的字长每增加1bit,动态范围扩大6dB。16bit字长的动态范围为96dB。程序员必须时刻关注溢出的发生。例如，在作图像处理时，图像作旋转、移动等，就很容易产生溢出。这时，要么不断地移位定标，要么作截尾。前者要耗费大量的程序空间和执行时间，后者则很快带来图像质量的劣化。总之，是使整个系统的性能下降。在处理低信噪比信号的场合，例如进行语音识别、雷达和声纳信号处理时，也会发生类似的问题。 32bit浮点运算DSP的动态范围可以作到1536dB，这不仅大大扩大了动态范围，提高了运算精度，还大大节省了运算时间和存储空间，因为大大减少了定标，移位和溢出检查。 由于浮点DSP的浮点运算用硬件来实现，可以在单周期内完成，因而其处理速度大大高于定点DSP。这一优点在实现高精度复杂算法时尤为突出，为复杂算法的实时处理提供了保证。 32bit浮点DSP的总线宽度较定点DSP宽得多，因而寻址空间也要大得多。这一方面为大型复杂算法提供了可能、因为省的DSP目标子程序已使用到几十MB存储器或更多；另一方面也为高级语言编译器、DSP操作系统等高级工具软件的应用提供了条件。DSP的进一步发展，必然是多处理器的应用。新型的浮点DSP已开始在通信口的设置和强化、资源共享等方面有所响应。

DSP与数字信号处理作业

1、什么是DSP？简述DSPs的特点？简述DSPs与MCU、FPGA、ARM的区别？学习DSP开发需要哪些知识？学习DSP开发需要构建什么开发环境？(15分) 答：（1）DSP是Digital Signal Processing(数字信号处理的理论和方法）的缩写，同时也是Digital Signal Processor（数字信号处理的可编程微处理器）的缩写。通常流过器件的电压、电流信号都是时间上连续的模拟信号，可以通过A/D器件对连续的模拟信号进行采样，转换成时间上离散的脉冲信号，然后对这些脉冲信号量化、编码，转化成由0和1构成的二进制编码，也就是常说的数字信号。DSP能够对这些数字信号进行变换、滤波等处理，还可以进行各种各样复杂的运算，来实现预期的目标。 （2）DSP既然是特别适合于数学信号处理运算的微处理器，那么根据数字信号处理的要求，DSP芯片一般具有下面所述的主要特点：1）程序空间和数据空间分开，CPU可以同时访问指令和数据； 2）在一个指令周期内可以完成一次乘法和一次加法运算； 3）片内具有快速RAM，通常可以通过独立的数据总线在程序空间和数据空间同时访问； 4）具有低开销和无开销循环及跳转的硬件支持； 5）具有快速的中断处理和硬件I/O支持； 6）可以并行执行多个操作； 7）支持流水线操作，使得取址、译码和执行等操作可以重复执行。（3）DSP采用的是哈佛结构，数据空间和存储空间是分开的，通过

独立的数据总线在数据空间和程序空间同时访问。而MCU采用的是冯·诺依曼结构，数据空间和存储空间共用一个存储器空间，通过一组总线（地址总线和数据总线）连接到CPU）。很显然，在运算处理能力上，MCU不如DSP；但是MCU价格便宜，在对性能要求不是很高的情况下，还是很具有优势的。 ARM是Advanced RISC（精简指令集）Machines的缩写是面向低运算市场的RISC微处理器。ARM具有比较强的事务管理功能，适合用来跑跑界面、操作系统等，其优势主要体现在控制方面，像手持设备90%左右的市场份额均被其占有。而DSP的优势是其强大的数据处理能力和较高的运算速度，例如加密/解密、调制/解调等。 FPGA是Field Programmable Gate Array（现场可编程门阵列）的缩写，它是在PAL、GAL、PLD等可编程器件的基础上进一步发展的产物，是专用集成电路中集成度最高的一种。FPGA采用了逻辑单元阵列LCA（Logical Cell Array）的概念，内部包括了可配置逻辑模块CLB、输入/输出模块IOB、内部连线三个部分。用户可以对FPGA内部的逻辑模块和I/O模块进行重置配置，已实现用户自己的逻辑。它还具有静态可重复编程和动态在系统重构的特性，使得硬件的功能可以像软件一样通过编程来修改。使用FPGA来开发数字电路，可以大大缩短设计时间，减少PCB面积，提高系统的可靠性；同时FPGA可以用VHDL或Verilog HDL来编程，灵活性强。由于FPGA能够进行编程、除错、再编程和重复操作，因此可以充分地进行设计开发和验证。当电路有少量改动时，更能显示出FPGA的优势，其现场编程能力可