数字信号处理第2章习题解答
2.1 今对三个正弦信号1()cos(2)a x t t π=，2()cos(6)a x t t π=-，3()cos(10)a x t t π=进行理想采样，采样频率为8s πΩ=，求这三个序列输出序列，比较其结果。画出
1()a x t 、2()a x t 、3()a x t 的波形及采样点位置并解释频谱混淆现象。
解：采样周期为2184
T ππ=
= 三个正弦信号采样得到的离散信号分别表示如下：
1()cos(2)cos()42a n x n n π
π=?=
2()cos(6)cos()42a n x n n π
π=-?=-
3()cos(10)cos()42
a n x n n π
π=?=
输出序列只有一个角频率2π
，其中1()a x n 和3()a x n 采样序列完全相同，2()a x n 和
1()a x n 、3()a x n 采样序列正好反相。
三个正弦信号波形及采样点位置图示如下：
t
x a 1(t )
t
x a 2(t )
t
x a 3(t )
三个正弦信号的频率分别为1Hz 、3Hz 和5Hz ，而采样频率为4Hz ，采样频率大于第一个正弦信号频率的两倍，但是小于后两个正弦信号频率的两倍，因而由第一个信号的采样能够正确恢复模拟信号，而后两个信号的采样不能准确原始的模拟信号，产生频谱混叠现象。
2.3 给定一连续带限信号()a x t 其频谱当f B >时，()a X f 。求以下信号的最低采样频率。
（1）2()a x t （2）(2)a x t （3）()cos(7)a x t Bt π
解：设()a x t 的傅里叶变换为()a X j Ω
（1）2
()a x t 的傅里叶变换为
22()[()]B
a a B
X j X j d ππ
ωωω-?Ω-?
因为22,22B B B B πωππωπ-≤≤-≤Ω-≤ 所以44B B ππ-≤Ω≤
即2()a x t 带限于2B ，最低采样频率为4B 。 （2）(2)a x t 的傅里叶变换为
1
(/2)2
a X j Ω 2/22B B ππ-≤Ω≤，即44B B ππ-≤Ω≤
即(2)a x t 带限于2B ，最低采样频率为4B 。 （3）()771
()cos(7)()2
j Bt j Bt a a x t Bt x t e e πππ-=
+ 根据傅里叶变换的频移性质，()cos(7)a x t Bt π的傅里叶变换为
[]1
((7)((7)2
a a X j B X j B ππΩ-+Ω+ 它为一个带宽为2B 的带通信号，其通带范围为59
22
B f B ≤≤。
根据带通模拟信+号的采样定理，最小采样频率为1/4
4(1) 4.52
B B ?+=。
补充知识：带通模拟信号的采样定理
设带通模拟信号的频带限制在L f 和H f 之间，其频谱最低频率大于L f ，最高频率小于H f ，信号带宽H L B f f =-。此带通模拟信号所需最小抽样频率s f 等于
21s k f B n ??
=+ ???
式中，B 为信号带宽；n 为商(
H f B )的整数部分，1,2,n =L ；为商(H f
B
)的小数部分，01k <<。
2.5 一带通模拟信号如图所示，现用以下采样频率对其采样。 （1）25 Hz （2）50 Hz （3）100 Hz 求采样后的频谱。
解：采样后的频谱分别如下图所示：
25 Hz 采样的频谱（注意：每一个三角形频谱都产生了混叠，以幅度的增加表示）
50 Hz 采样的频谱（没有混叠）
100 Hz 采样的频谱（没有混叠）
2.6 一带通模拟信号如图所示，求不产生混叠的最低采样频率。
解：由图知，60H f =，40L f =，20B =，
3H
f B
=的整数部分为3，小数部分为0，根据带通模拟信号的采样定理，最小采样频率为
240B Hz =。
0 40 60 -40 -60 f /Hz
0 25 50 -25 -50 f /Hz
75 100 -75 -100 0 25 50 -25 -50 f /Hz
75 100 -75 -100 0 25 50 -25 -50 f /Hz
填空题(每空2分，共20分) 信号与系统的时域分析与处理 1.序列x(n)的能量定义为__________。 2.线性移不变系统是因果系统的充分必要条件是__________。 3.设两个有限长序列的长度分别为N 和M ，则它们线性卷积的结果序列长度为__________。 4.线性系统同时满足_____和_____两个性质。 5.某线性移不变系统当输入x(n) =δ(n-1)时输出y(n) =δ(n -2) + δ(n -3)，则该系统的单位冲激响应h(n) =__________。 6.序列x(n) = cos (3πn)的周期等于__________。 7.线性移不变系统的性质有______、______和分配律。 8. 已知系统的单位抽样响应为h(n)，则系统稳定的充要条件是__________。 9.线性移不变系统是因果系统的充分必要条件是________。 10．序列x(n) = nR 4(n -1)，则其能量等于 _______ 。 11.两序列间的卷积运算满足_______，_______与分配率。 12信号处理有两种形式;其中一种是（ASP 模拟信号处理）;另一种是（DSP ：数字信号处理）。 13数字信号处理可以分为两类：信号(分析)和信号 (过滤) . 14数字信号是指 (时间) 和 (幅度)都离散的信号. 15．一个离散LTI 系统稳定的充要条件是系统的脉冲响应 h(n)满足关系式: ( ()h n ∞-∞<∞∑).LTI 离散系 统因果的充要条件是当且仅当 (h(n)=0,n<0). 16．互相关 ryx(l) 可以用卷积运算表示为（ryx(l)=y(l)*x(-l)）, 自相关 rxx(l)可写为 (rxx(l)=x(l)*x(-l) ) 17．若 LTI 系统的脉冲响应是有限长的,则该系统可称为(FIR:有限长脉冲响应) 滤波器, 否则称为 (IIR ：无 限长脉冲响应) 滤波器. 18．2n u(n)*δ(n-1)=( ). 0.8 n u(n)* 0.8 n u(n)=( ) 离散时间傅里叶变换（DTFT ） 1. 输入x(n)=cos(ω0n)中仅包含频率为ω0的信号，输出y(n)=x(n)cos(4 πn)中包含的频率为__________。 2．输入x(n)=cos(ω0n)中仅包含频率为ω0的信号，输出y(n)=x 2(n)中包含的频率为__________。 3．系统差分方程为y(n)=x(n)-x(n-1) 的系统被称为 (数字微分器). 4．实序列的DTFT 有两个重要属性:(周期性)和 (对称性), 根据这两个性质,我们只需要考虑[0,π]频率范围上的X(ejw) . 5．若DTFT[x(n)]= X(ejw), 则 DTFT[x*(n)]=(X*(e-jw)), DTFT[x(-n)]=( X(e-jw)); DTFT[x(n-k)]=( X(ejw) e-jwk). 6．DTFT[ (0.5)n u(n)]=（1 10.5jw e --）; 7．x(n)={ 1,2,3,4},DTFT[x(n)]=(1+2 e-jw+3 e-j2w+4 e-j3w ) .
第二章 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. 有一个线性移不变的系统，其系统函数为： 2z 2 1 )21)(2 11(2 3)(11 1<<-- - = ---z z z z H 1）用直接型结构实现该系统 2)讨论系统稳定性，并求出相应的单位脉冲响应)(n h 4．试用冲激响应不变法与双线性变换法将以下模拟滤波器系统函数变换为数字滤波器系统函数： H(s)= 3) 1)(s (s 2 ++其中抽样周期T=1s 。 三、有一个线性移不变的因果系统，其系统函数为： ) 21)(2 1 1(2 3)(111------= z z z z H 1用直接型结构实现该系统 2)讨论系统稳定性，并求出相应的单位脉冲响应)(n h 七、用双线性变换设计一个三阶巴特沃思数字低通虑波器，采样频率为kHz f s 4=（即采样周期为s T μ250=）,其3dB 截止频率为kHz f c 1=。三阶模拟巴特沃思滤波器为： 3 2 ) ()(2)(211)(c c c a s s s s H Ω+Ω+Ω+= 解1)2 111112 5 12 3) 21)(2 1 1(2 3)(------+-- = --- = z z z z z z z H …………………………….. 2分 当2 1 2> >z 时： 收敛域包括单位圆……………………………6分 系统稳定系统。……………………………….10分 1111 1211 2 111)21)(2 11(2 3)(------- -= -- - = z z z z z z H ………………………………..12分 )1(2)()2 1 ()(--+=n u n u n h n n ………………………………….15分 4.（10分）解： 3 1 11)3)(1(1)(+- +=++= s s s s s H ………………1分 1 311)(------ -= Z e s T Z e T z H T T ……………………3分
实验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、一线性时不变系统，输入为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 n /120)是周期的，则周期是N二8 。 10、用窗函数法设计FIR数字滤波器时，过渡带的宽度不但与窗的类型有关，还与窗的采样点数有关 11、DFT与DFS有密切关系，因为有限长序列可以看成周期序列的主值区间截断，而周期序列可以看成有限长序列的周期延拓。 12、对长度为N的序列x(n)圆周移位m位得到的序列用Xn(n)表示，其数学表达式为x m(n)= x((n-m)) N R(n)。 13、对按时间抽取的基2-FFT流图进行转置，并将输入变输出，输出变输入即可得到按频率抽取的基 2-FFT流图。 14、线性移不变系统的性质有交换率、结合率和分配律。
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
广西大学计电学院《数字信号处理》课程实验 适用专业：电信和通信工程专业；实验学时：9 学时 一、实验的性质、任务和基本要求 （一）本实验课的性质、任务 数字信号处理课程实验是数字信号处理课程的有效的补充部分，通过实验，使学生巩固和加深数字信号处理的理论知识的理解和掌握，在实验过程中了解简单但是完整的数字信号处理的工程实现方法和流程。通过实践进一步加强学生独立分析问题和解决问题的能力、实际动手能力、综合设计及创新能力的培养。 （二）基本要求 掌握数字信号处理基本理论知识和滤波器设计及应用。 （三）实验选项
二、实验教学内容 实验1 数字信号处理在音乐信号延时和混响处理中的应用（时域处理）－p278 1、实验目的和内容 掌握数字信号处理方法在音乐信号混响处理中的应用（在时域处理）。 2、实验内容 按照p278的内容要求，采用MATLAB设计相应滤波器及小型应用软件APP，实现对一段音乐信号进行延时和混响处理和输出，能实现混响程度的调节。 3、实验要求 1）提供MATLAB程序； 2）写实验报告，对设计的思路和步骤结果进行分析和说明。 实验2 数字信号处理在音乐信号均衡处理中的应用（频域处理）－p281 1、实验目的 掌握数字信号处理方法在音乐信号均衡处理中的应用（在频域处理）。 2、实验内容 按照p281的内容要求，采用MATLAB设计相应滤波器，及小型应用软件APP,实现对一段音乐信号进行均衡处理和输出，更高要求是能实现均衡程度的调节。 3、实验要求 1）提供MATLAB程序； 2）写实验报告，对设计的思路和步骤结果进行分析和说明。 实验3 双音多频（DTMF）通信设计的MATLAB仿真 1、实验目的
数字信号处理作业实验题报告 第一章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)')
第1章选择题 1．信号通常是时间的函数，数字信号的主要特征是：信号幅度取 ；时间取 B 。 A.离散值；连续值 B.离散值；离散值 C.连续值；离散值 D.连续值；连续值 2．数字信号的特征是（ B ） A ．时间离散、幅值连续 B ．时间离散、幅值量化 C ．时间连续、幅值量化 D ．时间连续、幅值连续 3．下列序列中属周期序列的为（ D ） A ．x(n) = δ(n) B ．x(n) = u(n) C ．x(n) = R 4(n) D ．x(n) = 1 4．序列x(n)=sin ??? ??n 311的周期为（ D ） A ．3 B ．6 C ．11 D ．∞ 5. 离散时间序列x (n )=cos(n 73π-8π )的周期是 ( C ) A. 7 B. 14/3 C. 14 D. 非周期 6．以下序列中（ D ）的周期为5。 A ．)853cos( )(ππ+=n n x B. )853sin()(ππ+=n n x C. )852()(π+=n j e n x D. )852()(ππ+=n j e n x 7．下列四个离散信号中，是周期信号的是（ C ）。 A ．sin100n B. n j e 2 C. n n ππ30sin cos + D. n j n j e e 5431 π - 8．以下序列中 D 的周期为5。 A.)853cos( )(π+=n n x B.)853sin()(π+=n n x C.)852 ()(π +=n j e n x D.)852 ()(ππ+ =n j e n x 9．离散时间序列x (n )=cos ??? ??+353ππ n 的周期是( C ) A.5 B.10/3 C.10 D.非周期 10.离散时间序列x(n)=sin （ 5n 31π+）的周期是( D ) A.3 B.6 C.6π D.非周期 11.序列x (n )=cos ? ?? ??n 5π3的周期为( C ) A.3 B.5 C.10 D.∞ 12．下列关系正确的为（ C ） A ．u(n)=∑=n k 0 δ (n) B ．u(n)=∑∞=0k δ (n) C ．u(n)=∑-∞=n k δ (n) D ．u(n)=∞-∞=k δ (n)
实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）；
成绩： 《数字信号处理》作业与上机实验 （第三四章） 班级： 13-电信 学号： 姓名： 任课老师：李宏民 完成时间： 15.11.08 信息与通信工程学院 2015—2016学年第1学期
第3-4章 离散傅里叶变换及快速算法与应用 一、实验内容及要求： 1、对信号在x(n)={1,2,3,4,5,6,7,8}，n=0，1，2....7,求其N=8点的DFT 。要求采用基于时间抽取算法编写FFT 实现程序，画出DFT 幅度谱与相位谱,并将计算结果与用MATLAB 自带的FFT 函数计算结果进行比较。 2、一个由40Hz 和100 Hz 正弦信号构成的信号，受零均值随机噪声的干扰（噪声服从标准正态分布，由randn 函数产生，n(t)= randn(m,n) 返回一个m*n 的随机矩阵），即()sin(2100)sin(240)()s t t t n t ππ=++g g 。数据 采样率为500Hz ，试用FFT 函数来分析其信号频率成分。 ① 求其幅度频谱，从频谱图（横坐标以HZ 为单位，用plot 函数画图）中能否观察出信号的2个频率分量？ ② 提高采样点数，再求该信号的幅度频谱图，此时幅度频谱发生了什么变化？信号的2个模拟频率和数字频率各为多少？FFT 频谱分析结果与理论上是否一致？ 3、研究高密度频谱与高分辨率频谱。频率分辨率是指所用的算法能将信号中两个靠得很近的谱峰分开的能力。信号末尾补零由于没有对原信号增加任何新的信息，因此不能提高频率分辨率，但可以减小栅栏效应，所得到的频谱称为高密度频谱。在维持采样频率不变的情况下，为提高分辨率只能增加采样点数N ，此时所得到的频谱称为高分辨率频谱。设有连续信号 ) 1092cos()1072cos()105.62cos()(333t t t t x a ??π+??π+??π=
==============================绪论============================== 1. A/D 8bit 5V 00000000 0V 00000001 20mV 00000010 40mV 00011101 29mV ==================第一章 时域离散时间信号与系统================== 1. ①写出图示序列的表达式 答：3)1.5δ(n 2)2δ(n 1)δ(n 2δ(n)1)δ(n x(n)-+---+++= ②用δ(n) 表示y (n )={2,7,19,28,29,15} 2. ①求下列周期 ②判断下面的序列是否是周期的; 若是周期的， 确定其周期。 （1）A是常数 8ππn 73Acos x(n)??? ? ??-= （2）)8 1 (j e )(π-=n n x 解： （1） 因为ω= 73π, 所以314 π2=ω, 这是有理数， 因此是周期序列， 周期T =14。 （2） 因为ω= 81, 所以ω π2=16π, 这是无理数， 因此是非周期序列。 ③序列)Acos(nw x(n)0?+=是周期序列的条件是是有理数2π/w 0。 3.加法 乘法 序列{2，3，2，1}与序列{2，3，5，2，1}相加为__{4，6，7，3，1}__，相乘为___{4，9，10，2} 。 移位 翻转：①已知x(n)波形，画出x(-n)的波形图。 ② 尺度变换：已知x(n)波形，画出x(2n)及x(n/2)波形图。 卷积和：①h(n)*求x(n)，其他0 2 n 0n 3，h(n)其他03n 0n/2设x(n) 例、???≤≤-=???≤≤= ②已知x (n )={1,2,4,3},h (n )={2,3,5}, 求y (n )=x (n )*h (n ) x (m )={1,2,4,3},h (m )={2,3,5}，则h (-m )={5,3,2}(Step1:翻转)
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名：
学号： 组员： 摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯
目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (3) 2 课程设计原理及流程 (5) 3.课程设计原理过程 (6) 3.1 语音信号的采集 (6) 3.2 语音信号的时频分析 (6) 3.3合成后语音加噪声处理 (9) 3.3.1 噪声信号的时频分析 (9) 3.3.2 混合信号的时频分析 (10) 3.4滤波器设计及消噪处理 (13) 3.4.1 设计IIR和FIR数字滤波器 (13) 3.4.2 合成后语音信号的消噪处理 (16) 3.4.3 比较滤波前后语音信号的波形及频谱 (17) 3.4.4回放语音信号 (18) 3.5结果分析 (18) 4 结束语 (19)
【最新整理，下载后即可编辑】 第二章 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)
数字信号处理复习题 一、选择题 1、某系统)(),()()(n g n x n g n y =有界，则该系统（ A ）。 A.因果稳定 B.非因果稳定 C.因果不稳定 D. 非因果不稳定 2、一个离散系统（ D ）。 A.若因果必稳定 B.若稳定必因果 C.因果与稳定有关 D.因果与稳定无关 3、某系统),()(n nx n y =则该系统（ A ）。 A.线性时变 B. 线性非时变 C. 非线性非时变 D. 非线性时变 4.因果稳定系统的系统函数)(z H 的收敛域是（ D ）。 A.9.0 计电学院《数字信号处理》课程实验 适用专业：电子通信工程专业；实验学时：9 学时 一、实验的性质、任务和基本要求 (一)本实验课的性质、任务 数字信号处理课程实验是数字信号处理课程的有效的补充部分，通过实验，使学生巩固和加深数字信号处理的理论知识的理解和掌握，在实验过程中了解简单但是完整的数字信号处理的工程实现方法和流程。通过实践进一步加强学生独立分析问题和解决问题的能力、实际动手能力、综合设计及创新能力的培养。 （二）基本要求 掌握数字信号处理基本理论知识和滤波器设计及应用。 （三）实验选项 二、实验教学内容 实验一 1、实验目的和要求 1）加深理解时域采样定理、体会使用MATLAB的离散FT函数fft( )来解决涉及模拟信号的问题； 2）加深理解对带通信号的采样特性，学会采用MATLAB解决该问题； 3）加深理解在频率采样法中，过渡点对所设计滤波器特性的影响。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第一组：张毅雷凌峰白法聪覃昱滔刘强何新文 第二组：邓志强林盛勇李日胜黎少锋梁聪杨晨 实验二 1、实验目的和要求 （1）加深理解采用数字信号处理方法对模拟信号处理的过程、掌握使用MATLAB处理的方法；对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线； （2）加深对截断效应的理解； （3）掌握使用MATLAB设计滤波器，并对语音信号处理的方法。对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第九组:汪涛张汉毅巫金敏张经中柳泽举 第六组：罗涛梁乐杰黄乃生 实验三 1、实验目的和要求 掌握采用MATLAB数字滤波器设计软件编制方法。软件要求在界面内有不同类型（高通低通带通带阻）滤波器的选择、或者只对低通滤波器采用不同方法设 数字信号处理期末复习题 一、单项选择题（在每个小题的四个备选答案中选出一个正确答案，并将正确答案的号码写在题干后面的括号内，每小题1分，共20分） 1.要从抽样信号不失真恢复原连续信号，应满足下列条件的哪几条( ① )。 (Ⅰ)原信号为带限 (Ⅱ)抽样频率大于两倍信号谱的最高频率 (Ⅲ)抽样信号通过理想低通滤波器 ①.Ⅰ、Ⅱ②.Ⅱ、Ⅲ ③.Ⅰ、Ⅲ④.Ⅰ、Ⅱ、Ⅲ 2.在对连续信号均匀采样时，若采样角频率为Ωs，信号最高截止频率为Ωc，则折叠频率为( ④ )。 ①Ωs ②.Ωc ③.Ωc/2 ④.Ωs/2 3.若一线性移不变系统当输入为x(n)=δ(n)时输出为y(n)=R3(n)，则当输入为u(n)-u(n-2)时输出为( ② )。 ①.R3(n) ②.R2(n) ③.R3(n)+R3(n-1) ④.R2(n)-R2(n-1) 4.已知序列Z变换的收敛域为｜z｜>1，则该序列为( ② )。 ①.有限长序列②.右边序列 ③.左边序列④.双边序列 5.离散系统的差分方程为y(n)=x(n)+ay(n-1)，则系统的频率响应( ③ )。 ①当｜a｜<1时，系统呈低通特性 ②.当｜a｜>1时，系统呈低通特性 ③.当0 6.序列x(n)=R5(n)，其8点DFT记为X(k)，k=0,1,…,7，则X(0)为( ④ )。 ①.2 ②.3 ③.4 ④.5 7.下列关于FFT的说法中错误的是( ① )。 ①.FFT是一种新的变换 ②.FFT是DFT的快速算法 ③.FFT基本上可以分成时间抽取法和频率抽取法两类 ④.基2 FFT要求序列的点数为2L(其中L为整数) 8.下列结构中不属于FIR滤波器基本结构的是( ③ )。 ①.横截型②.级联型 ③.并联型④.频率抽样型 9.已知某FIR滤波器单位抽样响应h(n)的长度为(M+1)，则在下列不同特性的单位抽样响应中可以用来设计线性相位滤波器的是( ④ )。 ①.h［n］=-h［M-n］ ②.h［n］=h［M+n］ ③.h［n］=-h［M-n+1］ ④.h［n］=h［M-n+1］ 10.下列关于用冲激响应不变法设计IIR滤波器的说法中错误的是( ④ )。 ①.数字频率与模拟频率之间呈线性关系 ②.能将线性相位的模拟滤波器映射为一个线性相位的数字滤波器 ③.容易出现频率混叠效应 ④.可以用于设计高通和带阻滤波器 11.利用矩形窗函数法设计FIR滤波器时，在理想特性的不连续点附近形成的过滤带的宽度近似等于( ① )。 ①.窗函数幅度函数的主瓣宽度 ②.窗函数幅度函数的主瓣宽度的一半 数字信号处理实验17 例题1： b=[8,-4,11,-2]; a=[1,-1.25,0.75,-0.125]; [sos,g]=tf2sos(b,a) [r,p,k]=residuez(b,a) 运行结果： sos = 1.0000 -0.1900 0 1.0000 -0.2500 0 1.0000 -0.3100 1.3161 1.0000 -1.0000 0.5000 g = 8 r =-8.0000 -12.0000i -8.0000 +12.0000i 8.0000 p = 0.5000 + 0.5000i 0.5000 - 0.5000i 0.2500 k = 16 例题2： sos=[1 0.9 0 1 -0.25 0;1 -3 2 1 1 0.5]; g=0.5; [b,a]=sos2tf(sos,g) [C,B,A]=dir2par(b,a) 子函数：dir2par（b，a）； function[C,B,A]=dir2par(num,den) M=length(num); N=length(den); [r1,p1,C]=residuez(num,den); p=cplxpair(p1,10000000*eps); I=cplxcomp(p1,p); r=r1(I); K=floor(N/2); B=zeros(K,2); A=zeros(K,3); if K*2==N; for i=1:2:N-2 Brow=r(i:1:i+1,:); Arow=p(i:1:i+1,:); [Brow,Arow]=residuez(Brow,Arow,[]); B(fix((i+1)/2),:)=real(Brow); A(fix((i+1)/2),:)=real(Arow); end; [Brow,Arow]=residuez(r(N-1),p(N-1),[]); B(K,:)=[real(Brow),0]; A(K,:)=[real(Arow),0]; else for i=1:2:N-1 Brow=r(i:1:i+1,:); Arow=p(i:1:i+1,:); [Brow,Arow]=residuez(Brow,Arow,[]); B(fix(i+1)/2,:)=real(Brow); A(fix(i+1)/2,:)=real(Arow); end end 子函数：cplxcomp(p1,p2) function I=cplxcomp(p1,p2) I=[]; for j=1:length(p2) for i=1:length(p1) if(abs(p1(i)-p2(j))<0.0001) I=[I,i]; end; end end; I=I'; 运行结果： b = 0.5000 -1.0500 -0.3500 0.9000 a = 1.0000 0.7500 0.2500 -0.1250 C = -7.2000 B = 3.9846 1.6308 3.7154 0 A = 1.0000 1.0000 0.5000 1.0000 -0.2500 0 例题3： b=[8,-4,11,-2]; a=[1,-1.25,0.75,-0.125]; [K,C]=tf2latc(b,a) [b,a]=latc2tf(K,C) 运行结果： K = -0.7327 0.6032 -0.1250《数字信号处理》课程实验题目
数字信号处理试题及参考答案
数字信号处理实验十七
相关主题
文本预览