最新《数字信号处理》期末试题库(有答案)
一. 填空题
1、一线性时不变系统，输入为 x（n）时，输出为y（n）；则输入为2x（n）时，输出为 2y(n) ；输入为x（n-3）时，输出为 y(n-3) 。
2、从奈奎斯特采样定理得出，要使实信号采样后能够不失真还原，采样频率fs与信号最高频率f max关系为： fs>=2f max。
3、已知一个长度为N的序列x(n)，它的离散时间傅立叶变换为X（e jw），它的N点离散傅立叶变换X（K）是关于X（e jw）的 N 点等间隔采样。
4、有限长序列x(n)的8点DFT为X（K），则X（K）= 。
5、用脉冲响应不变法进行IIR数字滤波器的设计，它的主要缺点是频谱的交叠所产生的现象。
6．若数字滤波器的单位脉冲响应h（n）是奇对称的，长度为N，则它的对称中心是 (N-1)/2 。
7、用窗函数法设计FIR数字滤波器时，加矩形窗比加三角窗时，所设计出的滤波器的过渡带比较窄，阻带衰减比较小。8、无限长单位冲激响应（IIR）滤波器的结构上有反馈环路，因此是递归型结构。
9、若正弦序列x(n)=sin(30nπ/120)是周期的,则周期是N= 8 。
10、用窗函数法设计FIR数字滤波器时，过渡带的宽度不但与窗的类型有关，还与窗的采样点数有关
11．DFT与DFS有密切关系，因为有限长序列可以看成周期序列的主值区间截断，而周期序列可以看成有限长序列的周期延拓。12．对长度为N的序列x(n)圆周移位m位得到的序列用xm(n)表示，其数学表达式为xm(n)= x((n-m))NRN(n)。
13．对按时间抽取的基2-FFT 流图进行转置，并 将输入变输出，输出变输入 即可得到按频率抽取的基2-FFT 流图。
14.线性移不变系统的性质有 交换率 、 结合率 和分配律。 15.用DFT 近似分析模拟信号的频谱时，可能出现的问题有混叠失真、 泄漏 、 栅栏效应 和频率分辨率。
16.无限长单位冲激响应滤波器的基本结构有直接Ⅰ型，直接Ⅱ型， 串联型 和 并联型 四种。
17.如果通用计算机的速度为平均每次复数乘需要5μs ，每次复数加需要1μs ，则在此计算机上计算210点的基2 FFT 需要 10 级蝶形运算，总的运算时间是______μs 。 二．选择填空题
1、δ(n)的z 变换是 A 。
A. 1
B.δ(w)
C. 2πδ(w)
D. 2π
2、从奈奎斯特采样定理得出，要使实信号采样后能够不失真还原，采样频率f s 与信号最高频率f max 关系为： A 。
A. f s ≥ 2f max
B. f s ≤2 f max
C. f s ≥ f max
D. f s ≤f max
3、用双线性变法进行IIR 数字滤波器的设计，从s 平面向z 平面转换的关系为s= C 。
A. 1111z z z --+=- B . 1111z z z ---=+s C. 11211z z T z ---=+ D. 11
211z z T z
--+=- 4、序列x 1（n ）的长度为4，序列x 2（n ）的长度为3，则它们线性卷积的长度是 B ，5点圆周卷积的长度是 。
A. 5, 5 B . 6, 5 C. 6, 6 D. 7, 5
5、无限长单位冲激响应（IIR ）滤波器的结构是 C 型的。
A. 非递归
B. 反馈
C.递归
D. 不确定
6、若数字滤波器的单位脉冲响应h（n）是对称的，长度为N，则它的
对称中心是 B 。
A. N/2
B.（N-1）/2
C. （N/2）-1
D. 不确
定
7、若正弦序列x(n)=sin(30nπ/120)是周期的,则周期是N= D 。
A. 2π
B. 4π
C. 2
D. 8
8、一LTI系统，输入为 x（n）时，输出为y（n）；则输入为2x（n）时，输出为 A ；输入为x（n-3）时，输出为。A. 2y（n），y（n-3） B. 2y（n），y（n+3） C. y（n），y（n-3）
D. y（n），y（n+3）
9、用窗函数法设计FIR数字滤波器时，加矩形窗时所设计出的滤波器，
其过渡带比加三角窗时 A ，阻带衰减比加三角窗时。
A.窄，小
B. 宽，小
C. 宽，大
D. 窄，大
10、在N=32的基2时间抽取法FFT运算流图中，从x(n)到X(k)需 B 级蝶形运算
过程。A. 4 B. 5 C. 6 D. 3
11．X(n)=u(n)的偶对称部分为（ A ）。
A． 1/2+δ(n)/2 B. 1+δ(n) C. 2δ(n) D. u(n)- δ(n)
12. 下列关系正确的为（ B ）。
A．∑
=-
=
n
k
k n
n
u
) (
)
(δ B.∑∞
=-
=
) (
)
(
k
k n
n
uδ
C．∑
-∞
=-
=
n
k
k n
n
u)
(
)
(δ D. ∑∞
-∞
=-
=
k
k n
n
u)
(
)
(δ
13．下面描述中最适合离散傅立叶变换DFT的是（ B ）A．时域为离散序列，频域也为离散序列
B．时域为离散有限长序列，频域也为离散有限长序列
C．时域为离散无限长序列，频域为连续周期信号
D．时域为离散周期序列，频域也为离散周期序列
14．脉冲响应不变法（ B ）
A．无混频，线性频率关系B．有混频，线性频率关系C．无混频，非线性频率关系D．有混频，非线性频率关系
15．双线性变换法（ C ）
A．无混频，线性频率关系B．有混频，线性频率关系C．无混频，非线性频率关系D．有混频，非线性频率关系
16．对于序列的傅立叶变换而言,其信号的特点是（ D ）A．时域连续非周期，频域连续非周期B．时域离散周期，频域连续非周期
C．时域离散非周期，频域连续非周期D．时域离散非周期，频域连续周期
17．设系统的单位抽样响应为h(n)，则系统因果的充要条件为
（ C ）
A．当n>0时，h(n)=0 B．当n>0时，h(n)≠0
C．当n<0时，h(n)=0 D．当n<0时，h(n)≠0
18.若一模拟信号为带限，且对其抽样满足奈奎斯特条件，则只要将抽样信号通过( A )即可完全不失真恢复原信号。
A.理想低通滤波器
B.理想高通滤波器
C.理想带通滤波器
D.理想带阻滤波器
19.若一线性移不变系统当输入为x(n)=δ(n)时输出为y(n)=R3(n)，则当输入为u(n)-u(n-2)时输出为( C )。
A.R3(n)
B.R2(n)
C.R3(n)+R3(n-1)
D.R2(n)+R2(n-1)
20.下列哪一个单位抽样响应所表示的系统不是因果系统?( D )
A.h(n)=δ(n)
B.h(n)=u(n)
C.h(n)=u(n)-u(n-1)
D.h(n)=u(n)-u(n+1)
21.一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包括( A )。
A.单位圆
B.原点
C.实轴
D.虚轴
22.已知序列Z变换的收敛域为｜z｜<1，则该序列为( C )。
A.有限长序列
B. 无限长右边序列
C.无限长左边序列
D. 无限长双边序列
23.实序列的傅里叶变换必是( A )。
A.共轭对称函数
B.共轭反对称函数
C.奇函数
D.偶函数
24.若序列的长度为M，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N需满足的条件是( A )。
A.N≥M
B.N≤M
C.N≤2M
D.N≥2M
25.用按时间抽取FFT计算N点DFT所需的复数乘法次数与( D )成正比。
A.N
B.N2
C.N3
D.Nlog2N
26.以下对双线性变换的描述中不正确的是( D )。
A.双线性变换是一种非线性变换
B.双线性变换可以用来进行数字频率与模拟频率间的变换
C.双线性变换把s平面的左半平面单值映射到z平面的单位圆内
D.以上说法都不对
27.以下对FIR和IIR滤波器特性的论述中不正确的是( A )。
A.FIR滤波器主要采用递归结构
B.IIR滤波器不易做到线性相位
C.FIR滤波器总是稳定的
D.IIR滤波器主要用来设计规格化的频率特性为分段常数的标准滤波器
28、设系统的单位抽样响应为h(n)=δ(n-1)+δ(n+1)，其频率响应为（ A ）
A．H(e jω)=2cosω B. H(e jω)=2sinω C. H(e jω)=cosω D. H(e j ω)=sinω
29. 若x(n)为实序列，X(e jω)是其离散时间傅立叶变换，则（ C ）A．X(e jω)的幅度合幅角都是ω的偶函数
B．X(e jω)的幅度是ω的奇函数，幅角是ω的偶函数
C．X(e jω)的幅度是ω的偶函数，幅角是ω的奇函数
D．X(e jω)的幅度合幅角都是ω的奇函数
30. 计算两个Ｎ１点和Ｎ2点序列的线性卷积，其中Ｎ１>Ｎ２，至少要做( B )点的ＤＦＴ。
A. Ｎ１
B. Ｎ１+Ｎ２-１
C. Ｎ１+Ｎ２+１
D. N2
31. y(n)+0.3y(n-1) = x(n)与y(n) = -0.2x(n) + x(n-1)是( C )。
A. 均为IIR
B. 均为FIR
C. 前者IIR，后者FIR
D. 前者FIR, 后者IIR
三．判断题
1、在IIR数字滤波器的设计中，用脉冲响应不变法设计时，从模拟角频率向数字角频率转换时，转换关系是线性的。（√）
2．在时域对连续信号进行抽样，在频域中，所得频谱是原信号频谱
的周期延拓。（√）
3、x(n)=cos（w0n)所代表的序列一定是周期的。（×）
4、y(n)=x2(n)+3所代表的系统是时不变系统。（√）
5、用窗函数法设计FIR数字滤波器时，改变窗函数的类型可以改变
过渡带的宽度。（√）
6、有限长序列的N点DFT相当于该序列的z变换在单位圆上的N点等
间隔取样。（√）
7、一个线性时不变离散系统是因果系统的充分必要条件是：系统函数
H(Z)的极点在单位圆内。（×）
8、有限长序列的数字滤波器都具有严格的线性相位特性。（×）
9、x(n) ,y(n)的线性卷积的长度是x(n) ,y(n)的各自长度之和。（×）
10、用窗函数法进行FIR数字滤波器设计时，加窗会造成吉布斯效应。（√）
12、在IIR数字滤波器的设计中，用双线性变换法设计时，从模拟角
频率向数字角频率转换时，转换关系是线性的。（×）
13．在频域中对频谱进行抽样，在时域中，所得抽样频谱所对应的序
列是原序列的周期延拓。（√）
14、有限长序列h(n)满足奇、偶对称条件时，则滤波器具有严格的线
性相位特性。（√）
15、y(n)=cos[x(n)]所代表的系统是线性系统。（×）
16、x(n) ,y(n)的循环卷积的长度与x(n) ,y(n)的长度有关；x(n) ,y(n)
的线性卷积的长度与x(n) ,y(n)的长度无关。（×）
17、在N=8的时间抽取法FFT运算流图中，从x(n)到x(k)需3级蝶形运算过程。（√）
18、用频率抽样法设计FIR数字滤波器时，基本思想是对理想数字滤波器的频谱作抽样，以此获得实际设计出的滤波器频谱的离散值。（√）
19、用窗函数法设计FIR数字滤波器和用频率抽样法设计FIR数字滤波器的不同之处在于前者在时域中进行，后者在频域中进行。（√）20、用窗函数法设计FIR数字滤波器时，加大窗函数的长度可以减少过渡带的宽度，改变窗函数的种类可以改变阻带衰减。（√）21、一个线性时不变的离散系统，它是因果系统的充分必要条件是：系统函数H(Z)的极点在单位圆外。（×）
22、一个线性时不变的离散系统，它是稳定系统的充分必要条件是：系统函数H(Z)的极点在单位圆内。（√）
23.对正弦信号进行采样得到的正弦序列必定是周期序列。( × )
24.常系数差分方程表示的系统必为线性移不变系统。( × )
25.序列的傅里叶变换是周期函数。( √ )
26.因果稳定系统的系统函数的极点可能在单位圆外。( × )
27.FIR滤波器较之IIR滤波器的最大优点是可以方便地实现线性相位。(√ )
28. 用矩形窗设计FIR滤波器，增加长度N可改善通带波动和阻带衰减。（×）
29. 采样频率fs=5000Hz，DFT的长度为2000，其谱线间隔为2.5Hz。（√）
三、计算题
一、设序列x(n)={4，3，2，1} ，另一序列h(n) ={1，1，1，1}，n=0,1,2,3
（1）试求线性卷积 y(n)=x(n)*h(n)（2）试求6点循环卷积。（3）试
求8点循环卷积。
二．数字序列 x(n)如图所示. 画出下列每个序列时域序列： (1) x(n-2); (2)x(3-n); (3)x[((n-1))6],(0≤n ≤5);
(4)x[((-n-1))6],(0≤n ≤5);
x[((n-1))]
n
n
三．已知一稳
定的LTI 系统的H(z)为
)
21)(5.01()
1(2)(111------=z z z z H 试确定该系统H(z)的收敛域和脉冲响应h[n]。
解：系统有两个极点，其收敛域可能有三种形式，|z|<0.5, 0.5<|z|<2, |z|>2
因为稳定，收敛域应包含单位圆，则系统收敛域为：0.5<|z|<2
1
1
111213/25.013/4)21)(5.01()1(2)(--------=---=z z z z z
z H )1(23
2)
()5.0
(34)(--+=n u n u n h n n
n
四．设x(n)是一个10点的有限序列
x （n ）={ 2,3,1,4,-3,-1,1,1,0,6}，不计算DFT ，试确定下列表达式的值。
(1) X(0), (2) X(5), (3) ∑=9
0)(k k X ,（4）∑=-9
5/2)(k k j k X e π
解：（1）
（2）
（3） （4）
五． x(n)和h(n)是如下给定的有限序列 x(n)={5, 2, 4, -1, 2}， h(n)={-3, 2, -1 }
(1) 计算x(n)和h(n)的线性卷积y(n)= x(n)* h(n)；(2) 计算x(n)和h(n)的6 点循环卷积y 1(n)= x(n)⑥h(n)；(3) 计算x(n)和h(n)的8 点循环卷积y 2(n)= x(n)⑧h(n)；比较以上结果，有何结论？
14
][]0[19
0===∑=n N
n x X W 12
][][]5[1
19
180510
-=-=
==???-=∑∑====奇
偶
奇数
偶数n n n n n n x n x X n n W
20
]0[*10][][101]0[9
9
===∑∑==x k X k X x k k 0
]8[*10][]
[101]))210[((]
[]
))[((2
)10/2(9
2
)10/2(9
10)/2(===-?
--=-=-∑∑x k X e
k X e
x k X e m n x k j k k j k m N k j N πππ
解：（1）y(n)= x(n)* h(n)={-15,4,-3,13,-4,3,2}
5 2 4 -1 2
-3 2 1
5 2 4 -1 210 4 8 -2 4-15 -
6 -12 3 -6
-15 4 -3 13 -4 3 2
(2)y 1(n)= x(n)⑥h(n)= {-13,4,-3,13,-4,3} (3)因为8>(5+3-1),
所以y 3(n)= x(n)⑧h(n)＝{-15,4,-3,13,-4,3,2，0} y 3(n)与y(n)非零部分相同。
六．用窗函数设计FIR 滤波器时，滤波器频谱波动由什么决定 _____________，滤波器频谱过渡带由什么决定_______________。 解：窗函数旁瓣的波动大小，窗函数主瓣的宽度
七．一个因果线性时不变离散系统，其输入为x[n]、输出为y[n]，系统的差分方程如下：
y （n ）-0.16y(n-2)= 0.25x(n-2)＋x(n) (1)
求系统的系统函数 H(z)=Y(z)/X(z);系统稳定吗?画出系统直接
型II 的信号流图;
(2) 画出系统幅频特性。
解：(1)方程两边同求Z 变换： Y(z)-0.16z -2Y(z)= 0.25z -2X(z)＋X(z)
2
216.0125.01)()()(---+=
=z z z X z Y z H
5 2 4 -1 2
-3 2 15 2 4 -1 210 4 8 -2 4
-15 -6 -12 3 -6
-15 4 -3 13 -4 3 2
2
-13 4 -3 13 -4 3 2
(2)系统的极点为：0.4和－0.4,在单位圆内，故系统稳定。
(3)
(4)
2
2
八．如果需要设计FIR 低通数字滤波器，其性能要求如下： (1)阻带的衰减大于35dB, (2)过渡带宽度小于π/6.
请选择满足上述条件的窗函数，并确定滤波器h(n)最小长度N
解：根据上表，我们应该选择汉宁窗函数，
十．已知 FIR DF 的系统函数为H(z)=3-2z -1+0.5z -2-0.5z -4＋2z -5-3z -6,试分
48
6
8≥≤
N N π
π()
()
x n y n
别画出直接型、线性相位结构量化误差模型。
123456
十一．两个有限长的复序列x [n ]和h [n ]，其长度分别为N 和M ，设两序列的线性卷积为y [n ]=x [n ]*h [n ]，回答下列问题：. (1) 序列y [n ]的有效长度为多长？
(2) 如果我们直接利用卷积公式计算y [n ] ，那么计算全部有效y [n ]的需要多少次复数乘法？
(3) 现用FFT 来计算y [n ]，说明实现的原理，并给出实现时所需满足的条件，画出实现的方框图，计算该方法实现时所需要的复数乘法计算量。
解：(1) 序列y [n ]的有效长度为：N+M-1；
(2) 直接利用卷积公式计算y[n]， 需要MN 次复数乘法
(3) 需要L L 2log 3次复数乘法。
十二．用倒序输入顺序输出的基2 DIT-FFT 算法分析一长度为N 点的复
补零
补零
L点-DFT
L点-DFT
L点-IDFT
y(n)1e 2(n)
3
序列x[n] 的DFT，回答下列问题：
(1) 说明N所需满足的条件，并说明如果N不满足的话，如何处理？
(2)如果N=8, 那么在蝶形流图中，共有几级蝶形？每级有几个蝶形？
确定第2级中蝶形的蝶距(d m)和第2级中不同的权系数(W N r )。
(3)如果有两个长度为N点的实序列y1[n]和y2 [n]，能否只用一次N点的
上述FFT运算来计算出y1[n]和y2 [n]的DFT，如果可以的话，写出实现的原理及步骤，并计算实现时所需的复数乘法次数；如果不行，说明理由。
解(1)N应为2的幂，即N＝2m，（m为整数）；如果N不满足条件，可以补零。
(2)3级，4个，蝶距为2，W N0，W N2
(3) y[n]=y1[n]+jy2[n]
十三．考虑下面4个8点序列，其中 0≤n ≤7，判断哪些序列的8点DFT 是实数，那些序列的8点DFT 是虚数，说明理由。 (1) x 1[n ]={-1, -1, -1, 0, 0, 0, -1, -1}, (2) x 2[n ]={-1, -1, 0, 0, 0, 0, 1, 1}, (3) x 3[n ]={0, -1, -1, 0, 0, 0, 1, 1}, (4) x 4[n ]={0, -1, -1, 0, 0, 0, -1, -1},
解：
DFT[x e （n ）]=Re[X （k ）]
DFT[x 0（n ）]=jIm[X （k ）]
x 4[n ]的DFT 是实数 , 因为它们具有周期性共轭对称性；x 3[n ] 的DFT 是
虚数 , 因为它具有周期性共轭反对称性 十四. 已知系统函数2
11
3.025.0125.02)(---+-+=
z z z z H ，求其差分方程。
)
()()(*
n N X n N x n x e e e -=-=)
()()(*n N X n N x n x o o o --=--
=
解：
2
11
3.025.0125.02)(---+-+=
z z z z H
2
11
3.025.0125.02)()(---+-+=
z z z z X z Y
)25.02)(()3.025.01)((121---+=+-z z X z z z Y
)1(25.0)(2)2(3.0)1(25.0)(-+=-+--n x n x n y n y n y
十五.已知)1)(()8
14
31)((121---+=+-z z X z z z Y ，画系统结构图。
解：)1)(()8
14
31)((121---+=+-z z X z z z Y
1
1111
2
11
25.015
5.016)25.01)(5.01(1125.075.011)()()(----------
-=--+=
+-+=
=z z z z z z z z z X z Y z H
直接型I ： 直接型II ：
级联型：
并联型：
x [n
y [n ]
x [n ]
y [n ]
x [n y [n ]
n ]
数字信号处理实验作业
实验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所示）转换为级联型和并联型。
数字信号处理知识点总结
《数字信号处理》辅导 一、离散时间信号和系统的时域分析 （一） 离散时间信号 （1）基本概念 信号：信号传递信息的函数也是独立变量的函数，这个变量可以是时间、空间位置等。 连续信号：在某个时间区间，除有限间断点外所有瞬时均有确定值。 模拟信号：是连续信号的特例。时间和幅度均连续。 离散信号：时间上不连续，幅度连续。常见离散信号——序列。 数字信号：幅度量化，时间和幅度均不连续。 （2）基本序列（课本第7——10页） 1）单位脉冲序列 1,0()0,0n n n δ=?=?≠? 2）单位阶跃序列 1,0 ()0,0n u n n ≥?=?≤? 3）矩形序列 1,01 ()0,0,N n N R n n n N ≤≤-?=?<≥? 4）实指数序列 ()n a u n 5）正弦序列 0()sin()x n A n ωθ=+ 6）复指数序列 ()j n n x n e e ωσ= （3）周期序列 1）定义：对于序列()x n ，若存在正整数N 使()(),x n x n N n =+-∞<<∞ 则称()x n 为周期序列，记为()x n ，N 为其周期。 注意正弦周期序列周期性的判定（课本第10页） 2）周期序列的表示方法： a.主值区间表示法 b.模N 表示法 3）周期延拓 设()x n 为N 点非周期序列，以周期序列L 对作()x n 无限次移位相加，即可得到周期序列()x n ，即 ()()i x n x n iL ∞ =-∞ = -∑ 当L N ≥时，()()()N x n x n R n = 当L N <时，()()()N x n x n R n ≠ （4）序列的分解 序列共轭对称分解定理：对于任意给定的整数M ，任何序列()x n 都可以分解成关于/2c M =共轭对称的序列()e x n 和共轭反对称的序列()o x n 之和，即
数字信号处理实验一
实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结
数字信号处理作业答案
数字信号处理作业
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文件），该信号中除了目标信号之外还掺杂有强噪声，但噪声与目标信号的频率不重叠，要求采用本学期已学的知识对该信
数字信号处理实验
实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=
1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）
数字信号处理复习总结-最终版
绪论：本章介绍数字信号处理课程的基本概念。 0.1信号、系统与信号处理 1.信号及其分类 信号是信息的载体，以某种函数的形式传递信息。这个函数可以是时间域、频率域或其它域，但最基础的域是时域。 分类： 周期信号/非周期信号 确定信号/随机信号 能量信号/功率信号 连续时间信号/离散时间信号/数字信号 按自变量与函数值的取值形式不同分类： 2.系统 系统定义为处理（或变换）信号的物理设备，或者说，凡是能将信号加以变换以达到人们要求的各种设备都称为系统。 3.信号处理 信号处理即是用系统对信号进行某种加工。包括：滤波、分析、变换、综合、压缩、估计、识别等等。所谓“数字信号处理”，就是用数值计算的方法，完成对信号的处理。 0.2 数字信号处理系统的基本组成 数字信号处理就是用数值计算的方法对信号进行变换和处理。不仅应用于数字化信号的处理，而且
也可应用于模拟信号的处理。以下讨论模拟信号数字化处理系统框图。 （1）前置滤波器 将输入信号x a(t)中高于某一频率（称折叠频率，等于抽样频率的一半）的分量加以滤除。 （2）A/D变换器 在A/D变换器中每隔T秒（抽样周期）取出一次x a(t)的幅度，抽样后的信号称为离散信号。在A/D 变换器中的保持电路中进一步变换为若干位码。 （3）数字信号处理器（DSP） （4）D/A变换器 按照预定要求，在处理器中将信号序列x(n)进行加工处理得到输出信号y(n)。由一个二进制码流产生一个阶梯波形，是形成模拟信号的第一步。 （5）模拟滤波器 把阶梯波形平滑成预期的模拟信号；以滤除掉不需要的高频分量，生成所需的模拟信号y a(t)。 0.3 数字信号处理的特点 （1）灵活性。（2）高精度和高稳定性。（3）便于大规模集成。（4）对数字信号可以存储、运算、系统可以获得高性能指标。 0.4 数字信号处理基本学科分支 数字信号处理（DSP）一般有两层含义，一层是广义的理解，为数字信号处理技术——DigitalSignalProcessing，另一层是狭义的理解，为数字信号处理器——DigitalSignalProcessor。 0.5 课程内容 该课程在本科阶段主要介绍以傅里叶变换为基础的“经典”处理方法，包括：（1）离散傅里叶变换及其快速算法。（2）滤波理论（线性时不变离散时间系统，用于分离相加性组合的信号，要求信号频谱占据不同的频段）。 在研究生阶段相应课程为“现代信号处理”（AdvancedSignalProcessing）。信号对象主要是随机信号，主要内容是自适应滤波（用于分离相加性组合的信号，但频谱占据同一频段）和现代谱估计。 简答题： 1．按自变量与函数值的取值形式是否连续信号可以分成哪四种类型? 2．相对模拟信号处理，数字信号处理主要有哪些优点? 3．数字信号处理系统的基本组成有哪些？
西电数字信号处理大作业
第二章 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');
数字信号处理实验答案完整版
数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入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并用stem语句画出A、B、C、D、E、F、G。 clear all; 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; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。
数字信号处理上机作业
数字信号处理上机作业 学院：电子工程学院 班级：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;
数字信号处理实验三
实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期
x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response
数字信号处理作业+答案讲解
数字信号处理作业 哈尔滨工业大学 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）
什么是数字信号处理
什么是数字信号处理？有哪些应用？ 利用数字计算机或专用数字硬件、对数字信号所进行的一切变换或按预定规则所进行的一切加工处理运算。 例如：滤波、检测、参数提取、频谱分析等。 对于DSP:狭义理解可为Digital Signal Processor 数字信号处理器。广义理解可为Digital Signal Processing 译为数字信号处理技术。在此我们讨论的DSP的概念是指广义的理解。 数字信号处理是利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以得到符合人们需要的信号形式。 信号处理的实质是对信号进行变换。 信号处理的目的是获取信号中包含的有用信息，并用更直观的方式进行表达。 DSP的应用几乎遍及电子学每一个领域。 ▲通用数字信号处理器：自适应滤波，卷积，相关，数字滤波，FFT, 希尔伯特变换，波形生成，窗函数等等。 ▲语音信号处理：语音增强、识别、合成、编码、信箱等，文字/语音转换 ▲图形/图像处理：三维动画，图象鉴别/增强/压缩/传输，机器人视觉等等图 ▲特殊应用数字信号处理：振动和噪声分析与处理，声纳和雷达信号处理， 通信信号处理, 地震信号分析与处理，汽车安全及全球定位，生物医学工程等等。 在医疗、军事、汽车等行业，以及通信市场、消费类电子产品等中具有广阔的市场前景。 数字信号处理系统的基本组成：前置预滤波器（PrF）、a/d变换器(ADC)、数字信号处理器(DSP)、d/a变换器(DAC)、模拟滤波器(PoF) 数字信号处理特点: 1.大量的实时计算（FIR IIR FFT）， 2.数据具有高度重复（乘积和操作在滤波、卷积和FFT中等常见） 数字信号处理技术的意义、内容 数字信号处理技术是指数字信号处理理论的应用实现技术，它以数字信号处理理论、硬件技术、软件技术为基础和组成，研究数字信号处理算法及其实现方法。 意义： 在21世纪，数字信号处理是影响科学和工程最强大的技术之一 它是科研人员和工程师必须掌握的一门技巧 DSP芯片及其特点 ▲采用哈佛结构体系：独立的程序和数据总线，一个机器周期可同时进行程序读出和数据存取。对应的：冯·诺依曼结构。 ▲采用流水线技术： ▲硬件乘法器：具有硬件连线的高速“与或”运算器 ▲多处理单元：DSP内部包含多个处理单元。 ▲特殊的DSP指令：指令具有多功能，一条指令完成多个动作；如：倒位序指令等 ▲丰富的外设▲功耗低：一般DSP芯片功耗为0.5～4W。采用低功耗技术的DSP芯片只有0.1W/3.3V、1.6V （电池供电） DSP芯片的类别和使用选择 ▲按特性分：以工作时钟和指令类型为指标分类▲按用途分：通用型、专用型DSP芯片 ▲按数据格式分：定点、浮点各厂家还根据DSP芯片的CPU结构和性能将产品分成若干系列。 TI公司的TMS320系列DSP芯片是目前最有影响、最为成功的数字信号处理器，其产品销量一直处于领先地位，公认为世界DSP霸主。 ?目前市场上的DSP芯片有： ?美国德州仪器公司(TI):TMS320CX系列占有90%
数字信号处理实验4
数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线
（2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
数字信号处理作业-答案
数字信号处理作业-答案
数字信号处理作业
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）
数字信号处理作业-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 -=-=
数字信号处理
Matlab上机实验 报告 ； 学院：理学院 专业：10 电信 姓名：贺茂海 学号：2010142110 完成日期：2012.10.20
matlab上机实验 实验内容：1）阅读例子程序，观察输出波形，理解每条语句的含义。 （2）已知有限长序列x(n)=[7,6,5,4,3,2],求DFT和IDFT，要求：画出序列傅立叶变换对应的幅度谱和相位谱；画出原信号与傅立叶逆变换IDFT[X(k)]的图形进行比较。 （3）已知周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期重复次数为3次时的DFS和IDFS。要求：画出原信号序列的主值和周期序列的图形；画出离散傅立叶变换对应的幅度谱和相位谱。 （4）求x(n)=[7,6,5,4,3,2], 0=