填空题(每空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 ) .
数 字信号处理试卷答案 完整版 一、填空题：（每空1分，共18分） 1、 数字频率ω是模拟频率Ω对采样频率s f 的归一化，其值是 连续 （连续还是离散？）。 2、 双边序列z 变换的收敛域形状为 圆环或空集 。 3、 某序列的DFT 表达式为∑-==1 0)()(N n kn M W n x k X ，由此可以看出，该序列时域的长度为 )1()(n N h n h --= ，此时对应系统的频率响应，则其对应的相位函数 为ωω?21)(--=N 。 8、请写出三种常用低通原型模拟滤波器 巴特沃什滤波器 、 切比雪夫滤波器 、 椭圆滤波器 。 二、判断题（每题2分，共10分） 1、 模拟信号也可以与数字信号一样在计算机上进行数字信号处理，只要加一道采样的工序就可 以了。 （╳） 2、 已知某离散时间系统为)35()]([)(+==n x n x T n y ，则该系统为线性时不变系统。(╳)
3、 一个信号序列，如果能做序列的傅里叶变换（DTFT ），也就能对其做DFT 变换。（╳） 4、 用双线性变换法进行设计IIR 数字滤波器时，预畸并不能消除变换中产生的所有频率点的非 线性畸变。 （√） 5、 阻带最小衰耗取决于窗谱主瓣幅度峰值与第一旁瓣幅度峰值之比。 （╳） 三、（15分）、已知某离散时间系统的差分方程为 系统初始状态为1)1(=-y ，2)2(=-y ，系统激励为)()3()(n u n x n =， 试求：（1）系统函数)(z H ，系统频率响应)(ωj e H 。 ??? ???=+=-241()2(2(2121c c y zi zi 解之得 31=c ，42-=c ， 故系统零输入响应为： k zi k y )2(43)(-= 0≥k 系统零状态响应为 即 3 21528123)(-+--+-=z z z z z z z Y zs 对上式取z 反变换，得零状态响应为 )(])3(2 15)2(823[)(k k y k k zs ε+-=
数字信号处理作业提交日期：2016年7月15日
实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ，当输入以随机信号()x n ，且 ()() () x n s n v n =+，其中()s n 表示原始信号，即期望信号。()v n 表示噪声，则输出()y n 为()=()()m y n h m x n m -∑，我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ，因此称()y n 为估计值，用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即 它的平方的统计期望最小：222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小，则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为：()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为：xs xx R R h =，可知：1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的
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、)125.0cos()(n n x π=的基本周期是 16 。 2、一个序列 )(n x 的离散傅里叶变换的变换定义为 ∑-=-= 1 0/2)()(N n N nk j e n x k X π 3、对于M 点的有限长序列，频域采样不失真恢复时域序列的条件是频域采样点数N 不小于M 4、有界输入一有界输出的系统称之为 稳定系统 三、填空题（本大题10分，每小题2分） 1、在对连续信号进行频谱分析时，频谱分析范围受 采样 速率的限制。 2、 ? ∞ ∞ -=ωωδd ( 1 。 3、对于一个系统而言，如果对于任意时刻0n ，系统在该时刻的响应仅取决于在时刻及其以前的输入，则称该系统为 因果 系统。。 4、对一个LSI 系统而言，系统的输出等于输入信号与系统单位采样响应的线性 卷积 。 5、假设时域采样频率为32kHz ，现对输入序列的32个点进行DFT 运算。此时，DFT 输出的各点频率间隔为 1000 Hz 。 四、计算题（本大题20分） 某两个序列的线性卷积为 ) 5(3)3(2)2(2)1()() ()()(-+-+-+-+=*=n n n n n n x n h n y l δδδδδ计算这两个序列的4点圆周卷积。 解：将序列)(rL n y l +的值列在表中，求n =0，1，2，3时这些值的和。只有序列)(n y l 和)4(+n y l 在30≤≤n 区间内有非零值，所以只需列 将30≤≤n 各列内的值相加，有)()(n h n y =④ )3(2)2(2)1(4)()(-+-+-+=n n n n n x δδδδ 五、分析推导题（本大题12分） 如果)(n x 是一个周期为N 的周期序列，则它也是周期为2N 的周期序列， 把 )(n x 看作周期为N 的周期序列，其DFT 为)(1k X ，再把)(n x 看作周 期为2N 的周期序列，其DFT 为)(2 k X ，试利用)(1k X 确定)(2k X 。 解： ∑ -== 1 1)()(N n nk N W n x k X ∑-== 1 20 22)()(N n nk N W n x k X 令n m 2=，则N M 2= ∑ -== 1 22)2 ( )(N m k m M W m x k X = )2/(1k X 六、证明题（本大题18分） 一个有限冲击响应滤波器，它的单位采样相应 )(n h 的长度为 )12(+N 。如果)(n h 为实偶序列，证明系统函数的零点对于单位圆成镜 像对出现。 证： )(n h 是偶序列，所以)()(n h n h -= ?? ? ??=z H z H 1)( ??? ? ??=-θ θ ρ ρj j e H e H 1)( 又因为)(n h 是实序列，故有)()(* * z H z H = ??? ? ??=-θ θ ρ ρ j j e H e H 1)1 ( * 所以 ??? ? ??=θ θ ρρj j e H e H 1)(* 当θ ρj e z =时 ()0)(== θ ρj e H z H 当θ ρ j e z 1 = 时 ()00 )1 ( )(* * ====θ θ ρρ j j e H e H z H 七、综合题（本大题20分） 已知连续时间信号)16000cos()(t t x a π=，用6000/1=T 对 其采样。 （1）求最小采样频率； （2）图示其频谱特性； （3）分析其频谱是否有混叠。 解：（1）信号的最高频率 π 160000=Ω， ππ12000/2==ΩT s （2） （3）πππ32000212000/20=Ω<==ΩT s 没有满足奈奎斯特定理，频谱有混叠。
实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形
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的围，频谱出现了镜像对称的部分。
第一章 离散时间信号与系统 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 313 si n()( )()8 73cos( )( )(πππ π-==-=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 时，不一定是周期序列， n m m m n n y n - - -∞ = - ? = = ≥ ∑ 2 3 1 2 5 . 0 ) ( 0 1 当 3 4 n m n m m n n y n 2 2 5 . 0 ) ( 1 ? = = - ≤ ∑ -∞ = - 当 a a a n y n a a a n y n n h n x n y a n u a n h n u n x m m n n m m n -= = ->-= = -≤=<<--==∑∑--∞ =---∞=--1)(11)(1) (*)()(1 0,)1()()()(:1 时当时当解
①当=0/2ωπ整数，则周期为0/2ωπ； ②； 为为互素的整数）则周期、（有理数当 , 2 0Q Q P Q P =ωπ ③当=0/2ωπ无理数 ，则)(n x 不是周期序列。 解：（1）014 2/3 πω=，周期为14 （2）06 2/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 括号内相等，所以是因果的。（x 括号内表达式满足小于等于y 括号内表达式，系统是因果的） │y(n)│=│g(n)x(n)│<=│g(n)││x(n)│x(n)有界，只有在g(n)有界时，y(n)有界，系统才稳定，否则系统不稳定 （3）T[x(n)]=x(n-n0) 线性，移不变，n-n0<=n 即n0>=0时系统是因果的，稳定 （5）线性，移变，因果，非稳定 （7）线性，移不变，非因果，稳定 （8）线性，移变，非因果，稳定 8.
武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02
一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分
实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');
实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');
正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');
随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');
第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）’）；
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
==============================绪论============================== 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:翻转)
信息科学与技术学院本科三年级 数字信号处理实验报告 2011 年12 月21日
实验一 序列的傅立叶变换 实验目的 进一步加深理解DFS,DFT 算法的原理;研究补零问题;快速傅立叶变换 （FFT ）的应用。 实验步骤 1. 复习DFS 和DFT 的定义，性质和应用； 2. 熟悉MATLAB 语言的命令窗口、编程窗口和图形窗口的使用；利用提供的 程序例子编写实验用程序；按实验内容上机实验，并进行实验结果分析；写出完整的实验报告，并将程序附在后面。 实验内容 1. 周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析补零后，对信号频谱的影响。 实验结果： 60 ,7)4(;60,5)3(; 40,5)2(;20,5)1()] (~[)(~,2,1,01 )1(,01,1)(~=========±±=???-+≤≤+-+≤≤=N L N L N L N L n x DFS k X m N m n L m N L m N n m N n x ) 52.0cos()48.0cos()(n n n x ππ+=
2. 有限长序列x(n)的DFT （1） 取x(n)(n=0:10)时，画出x(n)的频谱X(k) 的幅度； （2） 将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出 x(n)的频谱X(k) 的幅度； （3） 取x(n)(n:0~100)时，画出x(n)的频谱X(k) 的幅度。利用FFT 进行谱分析 已知：模拟信号 以t=0.01n(n=0:N-1)进行采样，求N 点DFT 的幅值谱。 请分别画出N=45; N=50;N=55;N=60时的幅值曲线。 实验结果： ) 8cos(5)4sin(2)(t t t x ππ+=
数字信号处理期末试卷(含答案) 填空题（每题2分，共10题） 1、 1、 对模拟信号（一维信号，是时间的函数）进行采样后，就是 信号，再 进行幅度量化后就是 信号。 2、 2、 )()]([ωj e X n x FT =，用)(n x 求出)](Re[ωj e X 对应的序列 为 。 ３、序列)(n x 的N 点DFT 是)(n x 的Z 变换在 的N 点等间隔采样。 ４、)()(5241n R x n R x ==，只有当循环卷积长度L 时，二者的循环卷积等于线性卷积。 ５、用来计算N ＝16点DFT ，直接计算需要_________ 次复乘法，采用基2FFT 算法，需要________ 次复乘法，运算效率为__ _ 。 ６、FFT 利用 来减少运算量。 ７、数字信号处理的三种基本运算是： 。 ８、FIR 滤波器的单位取样响应)(n h 是圆周偶对称的，N=6， 3)3()2(2 )4()1(5 .1)5()0(======h h h h h h ，其幅 度特性有什么特性？ ，相位有何特性？ 。 ９、数字滤波网络系统函数为 ∑=--= N K k k z a z H 111)(，该网络中共有 条反馈支路。 １０、用脉冲响应不变法将)(s H a 转换为)(Z H ，若)(s H a 只有单极点k s ，则系统)(Z H 稳定的条件是 （取s T 1.0=）。 一、 选择题（每题3分，共6题） 1、 1、 )6 3()(π-=n j e n x ，该序列是 。 A.非周期序列 B.周期 6π = N C.周期π6=N D. 周期π2=N 2、 2、 序列)1()(---=n u a n x n ，则)(Z X 的收敛域为 。 A.a Z < B.a Z ≤ C.a Z > D.a Z ≥ 3、 3、 对)70() (≤≤n n x 和)190()(≤≤n n y 分别作20点DFT ，得)(k X 和)(k Y ， 19,1,0),()()( =?=k k Y k X k F ，19,1,0)],([)( ==n k F IDFT n f ， n 在 范围内时，)(n f 是)(n x 和)(n y 的线性卷积。 A.70≤≤n B.197≤≤n C.1912≤≤n D.190≤≤n 4、 4、 )()(101n R n x =，)()(72n R n x =，用DFT 计算二者的线性卷积，为使计算量尽可 能的少，应使DFT 的长度N 满足 。 A.16>N B.16=N C.16 数字信号处理复习题 一、选择题 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(完整word版)数字信号处理题库(附答案)
相关主题
文本预览