数字信号模拟题6

数字信号模拟题6
数字信号模拟题6

模拟试题(六)

一、(20分)设()

0()()n

k n l k N l y n x l W --==∑是由长度为N 的序列()x n 激励单位取样响应为

, 0

()0, 0

kn

N k W n h n n -?≥=?

的因果系统的输出。证明:

(1) ()(),0,1,,1k y N X k k N ==- ,其中()X k 是序列()x n 的离散傅立叶变换(5分)。 (2) ()k h n 系统函数为1

()1()()1k k k N Y z H z X z W z

--=

=-和()(1)()k

k N k y n W y n x n -=-+(5分)。 (3) 如将1

1k N

W z --同时乘以()k H z 的分子和分母,则可得1

12

1()212cos k N k W z H z k z z

N π----=????-+ ? ????

?(5分)。

(4) 如果定义

12

()212cos ()k X z k z z V z N π--????=-+ ? ????

?,则有 2()2cos (1)(2)()()()(1)k

k k k k k N k k v n v n v n x n y n v n W v n N π

????=---+=--

? ????

?

, 再根据上面两式画出()x n 和()k y n 间的直接Ⅱ型结构图。由该图解释离散傅立叶变换(5分)。

二、(20分)一个理想Hilbert 变换器, 0(), 0j d j H e j ω

ωπ

πω-≤≤?=?-≤

。取一实序列()x n 作为

()j d H e ω的输入,()x n 的频谱表示为()()()j j j R I X e X e jX e ωωω=+。

(1) (5分)求()()()j j j d X e H e X e ωωω=。

(2) (5分)如果构造一复序列()()()y n x n jx n =+,求()j Y e ω,当()j X e ω的带宽为

[],B B ω∈-且0B π<<时,画出()j X e ω和()j Y e ω的示意图,比较它们的带宽变化。 (3) 求当

1212, 0(), 0N j j d N j je H e je ω

ω

ωωππω-??

- ???

-??- ????-≤≤?=?

?-≤

时的傅立叶反变换()d h n ,并验证10,()(1)2d d d N h h n h N n -??

==--- ???

(10分)

三、(20分)已知一理想模拟带限微分器,其频率响应为, ()0, a j H j T π?

ΩΩ≤?

Ω=??Ω?其他,

T 是取样周期。

(1) 用12()j a r H e H j j r T T T ωω

π∞

=-∞??=- ???∑乘模拟低通滤波器, ()0, LP T H j T π?Ω≤?Ω=??Ω

?其他,

得到数字带限微分器()()()j j d LP H e H e H j ωω=Ω。求()j d H e ω傅立叶反变换()d h n ,其中T ω=Ω(10分)

。 (2) 画出()j d H e ω的示意图。(提示:2(1)ax ax xe dx a ax e c -=-+?)(10分)。 四、(20分)已知一线性移不变因果系统的差分方程

()(1)()x n ax n u n =--+

其中()u n 是均值为零、方差为2δ的白噪声过程,()x n 是当()u n 为输入时的系统输出,它也是一随机过程,a 是实常数。已知输入和输出间的互相关序列

2, 0

()0, 0xu m R m m δ?==?>?

(1) 求随机信号()x n 的自相关序列()xx R m (5分)。

(2) 先求出由()xx R m 表示的a 和2δ;再由随机过程{}n x 的N 个取样(0),(1),,(

x x x N - 构成的()xx R m 估计

1*0

1

()()(), 11N m

N n R m x n x n m N m N N

--==

+-+≤≤-∑

求出a 和2δ的估计?a 和2?δ,说明a 是否为正值、其绝对值是否大于1,为什么(5分)?

(3) 求出系统函数()H z 、功率谱()j xx S e ω和由?a 和2?δ表示的谱估计?()j xx

S e ω(10分)。 五、(20分)(1) 证明周期图21

()()j j per N S e X e N

ωω=在z 平面单位圆j z e ω=上的等间隔取样为

2

21()

(), 0,1,,1j per k N

S e X k k N N

ωπω=

=

=- 其中10

()()N j j n N n X e x n e ω

ω--==∑,而()X k 是随机过程{}n x 的N 个取样(0),(1),,(1)

x x x N - 构成序列的离散傅立叶变换(5分)。

(2) 因为{}()X k 是一随机过程。当已知[]E ()x n m =时,其中m 是常数,求[]E ()X k (5分)。

(3) 设某随机过程{}()x n 的功率谱

1, ()0,

j xx c S e ω

ωωωωπ

?≥?=?

≤?? 其中()()j j m

xx xx

m S e R

m e

ω

ω∞

-=-∞

=

∑。若取傅立叶变换()()()j j m N R

xx

m I e w m R

m e ω

ω∞

-=-∞

=

∑,其中

()R w m 是长度为N 的矩形窗。利用1

()()()2N xx

R I S

W d π

πωθωθθπ-

=

-?和窗谱幅度函数

sin(/2)()sin(/2)R N W ωωω=

,用20,,c c N πωωωωω===-

和2c N

π

ωω=+四个点做出()j N I e ω示意图,然后分析()j N I e ω和()j xx S e ω的区别,(其中2/c N ωπ )(10分)。

数字信号处理实验一

一、实验目的 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;

数字信号处理期末试卷!

数字信号处理模拟试题一 一、单项选择题(本大题共10小题,每小题2分,共20分) 在每小题列出的四个备选项中只有一个是符合题目要求的,请将其代码填写在题后的括号内。错选、多选或未选均无分。 1.在对连续信号均匀采样时,要从离散采样值不失真恢复原信号,则采样角频率Ωs与信号最高截止频率Ωc应满足关系(A ) A.Ωs>2Ωc B.Ωs>Ωc C.Ωs<Ωc D.Ωs<2Ωc 2.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?(D) A.y(n)=y(n-1)x(n) B.y(n)=x(n)/x(n+1) C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1) 3.已知某序列Z变换的收敛域为5>|z|>3,则该序列为(D ) A.有限长序列 B.右边序列 C.左边序列 D.双边序列 4.实偶序列傅里叶变换是(A ) A.实偶序列 B.实奇序列 C.虚偶序列 D.虚奇序列 5.已知x(n)=δ(n),其N点的DFT[x(n)]=X(k),则X(N-1)=(B) A.N-1 B.1 C.0 D.-N+1 6.设两有限长序列的长度分别是M与N,欲通过计算两者的圆周卷积来得到两者的线性卷积,则圆周卷积的点数至少应取(B ) A.M+N B.M+N-1 C.M+N+1 D.2(M+N) 7.下面说法中正确的是(C) A.连续非周期信号的频谱为周期连续函数 B.连续周期信号的频谱为周期连续函数 C.离散非周期信号的频谱为周期连续函数 D.离散周期信号的频谱为周期连续函数 8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构?(C ) A.直接型 B.级联型 C.频率抽样型 D.并联型 9.下列关于FIR滤波器的说法中正确的是(C) A.FIR滤波器容易设计成线性相位特性

数字信号处理试题

一、 单 项选择题 1. 序列x(n)=Re(e jn π/12 )+I m (e jn π/18 ),周期为( )。 A. 18π B. 72 C. 18π D. 36 2. 设C 为Z 变换X(z)收敛域内的一条包围原点的闭曲线,F(z)=X(z)z n-1 ,用留数法求X(z)的反变换时( )。 A. 只能用F(z)在C 内的全部极点 B. 只能用F(z)在C 外的全部极点 C. 必须用收敛域内的全部极点 D. 用F(z)在C 内的全部极点或C 外的全部极点 3. 有限长序列h(n)(0≤n ≤N-1)关于τ= 2 1 -N 偶对称的条件是( )。 A. h(n)=h(N-n) B. h(n)=h(N-n-1) C. h(n)=h(-n) D. h(n)=h(N+n-1) 4. 对于x(n)= n )21(u(n)的Z 变换,( )。 A. 零点为z=21,极点为z=0 B. 零点为z=0,极点为z=21 C. 零点为z=21,极点为z=1 D. 零点为z=2 1 ,极点为z=2 5、)()(101n R n x =,)()(72n R n x =,用DFT 计算二者的线性卷积,为使计算量尽可能的少,应使DFT 的长度N 满足 。 A.16>N B.16=N C.160,Z 变换的收敛域为( )。 A. 0<|z|<∞ B. |z|>0 C. |z|<∞ D. |z|≤∞ 9.在对连续信号均匀采样时,要从离散采样值不失真恢复原信号,则采样角频率Ωs 与信号最高截止频率Ωc 应满足关系( ) A. Ωs>2Ωc B. Ωs>Ωc C. Ωs<Ωc D. |Ωs<2Ωc 10.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?( ) A.y(n)=y(n-1)x(n) B.y(n)=x(n)/x(n+1) C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1)

数字信号处理实验1,2,3,4

实验一 连续时间系统的时域和频域分析相关MATLAB 函数1.设描述连续时间系统的微分方程为:)()()()()()()()(01)1(1)(01)1(1)(t f b t f b t f b t f b t y a t y a t y a t y a m m m m n n n n +'+++=+'+++---- 则可用向量和表示该系统,即 a b ] ,,,,[011a a a a a n n -=],,,,[011b b b b b m m -=注意,向量和的元素一定要以微分方程时间求导的降幂次序排列,且缺项要用0补齐。a b 如微分方程)()()(2)(3)(t f t f t y t y t y +''=+'+''表示该系统的向量为 ]2 3 1[=a ]1 0 1[=b (1)求解冲激响应:impulse()函数impulse()函数有以下四种调用格式: ① impulse(b,a) 该调用格式以默认方式绘制由向量和定义的连续时间系统的冲激响应的时域波形。a b ② impulse(b,a,t)该调用格式绘制由向量和定义的连续时间系统在时间范围内的冲激响应的时a b t ~0域波形。③ impulse(b,a, t1:p:t2)该调用格式绘制由向量和定义的连续时间系统在时间范围内,且以时间间a b 21~t t 隔均匀抽样的冲激响应的时域波形。p ④ y=impulse(b,a,t1:p:t2)该调用格式并不绘制系统冲激响应的波形,而是求出由向量和定义的连续时间系a b 统在时间范围内以时间间隔均匀抽样的系统冲激响应的数值解。21~t t p (2)求解阶跃响应:step()函数 step()函数也有四种调用格式:① step(b,a) ② step(b,a,t) ③ step(b,a, t1:p:t2) ④ y=step(b,a,t1:p:t2) 上述调用格式的功能与impulse()函数完全相同。

数字信号处理期末试卷(含答案)全..

数字信号处理期末试卷(含答案) 一、单项选择题(在每小题的四个备选答案中,选出一个正确答案,并将正确答案的序号填在括号内。 1.若一模拟信号为带限,且对其抽样满足奈奎斯特采样定理,则只要将抽样信号通过( )即可完全不失真恢复原信号。 A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2.下列系统(其中y(n)为输出序列,x(n)为输入序列)中哪个属于线性系统?( ) A.y(n)=x 3(n) B.y(n)=x(n)x(n+2) C.y(n)=x(n)+2 D.y(n)=x(n 2) 3..设两有限长序列的长度分别是M 与N ,欲用圆周卷积计算两者的线性卷积,则圆周卷积的长度至少应取( )。 A .M+N B.M+N-1 C.M+N+1 D.2(M+N) 4.若序列的长度为M ,要能够由频域抽样信号X(k)恢复原序列,而不发生时域混 叠现象,则频域抽样点数N 需满足的条件是( )。 A.N ≥M B.N ≤M C.N ≤2M D.N ≥2M 5.直接计算N 点DFT 所需的复数乘法次数与( )成正比。 A.N B.N 2 C.N 3 D.Nlog 2N 6.下列各种滤波器的结构中哪种不是FIR 滤波器的基本结构( )。 A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR 滤波器的幅度响应H(w)特点( ): A 关于0=w 、π、π2偶对称 B 关于0=w 、π、π2奇对称 C 关于0=w 、π2偶对称 关于=w π奇对称 D 关于0=w 、π2奇对称 关于=w π偶对称 8.适合带阻滤波器设计的是: ( ) A )n N (h )n (h ---=1 N 为偶数 B )n N (h )n (h ---=1 N 为奇数 C )n N (h )n (h --=1 N 为偶数

数字信号处理试题及答案

数字信号处理试题及答案 一、填空题:(每空1分,共18分) 1、 数字频率ω是模拟频率Ω对采样频率s f 的归一化,其值是 连续 (连续还是离散?)。 2、 双边序列z 变换的收敛域形状为 圆环或空集 。 3、 某序列的 DFT 表达式为∑-==1 0)()(N n kn M W n x k X ,由此可以看出,该序列时域的长度为 N ,变换后数字频域上相邻两个频率样点之间的间隔是 M π 2 。 4、 线性时不变系统离散时间因果系统的系统函数为2 52) 1(8)(2 2++--=z z z z z H ,则系统的极点为 2,2 1 21-=-=z z ;系统的稳定性为 不稳定 。系统单位冲激响应)(n h 的初值 4)0(=h ;终值)(∞h 不存在 。 5、 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点 的有限长序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为 64+128-1=191点 点的序列,如果采用基FFT 2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为 256 点。 6、 用冲激响应不变法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的 映射变换关系为T ω = Ω。用双线性变换法将一模拟滤波器映射为数字滤波器时,模拟频率Ω 与数字频率ω之间的映射变换关系为)2 tan(2ω T =Ω或)2arctan(2T Ω=ω。 7、当线性相位 FIR 数字滤波器满足偶对称条件时,其单位冲激响应)(n h 满足的条件为 )1()(n N h n h --= ,此时对应系统的频率响应)()()(ω?ω ωj j e H e H =,则其对应的相位函数 为ωω?2 1 )(-- =N 。 8、请写出三种常用低通原型模拟滤波器 巴特沃什滤波器 、 切比雪夫滤波器 、 椭圆滤波器 。 二、判断题(每题2分,共10分) 1、 模拟信号也可以与数字信号一样在计算机上进行数字信号处理,只要加一道采样的工序就可 以了。 (╳) 2、 已知某离散时间系统为)35()]([)(+==n x n x T n y ,则该系统为线性时不变系统。(╳)

数字信号处理实验4-6

实验4 离散系统的变换域分析 一、实验目的 1、熟悉对离散系统的频率响应分析方法; 2、加深对零、极点分布的概念理解。 二、实验原理 离散系统的时域方程为 其变换域分析方法如下: 频域: 系统的频率响应为: Z域: 系统的转移函数为:

分解因式: , 其中和称为零、极点。 三、预习要求 1.在MATLAB中,熟悉函数tf2zp、zplane、freqz、residuez、zp2sos的使用,其中:[z, p,K]=tf2zp(num,den)求得有理分式形式的系统转移函数的零、极点;zplane(z,p)绘制零、极点分布图;h=freqz(num,den,w)求系统的单位频率响应;[r,p,k]=residuez (num,den)完成部分分式展开计算;sos=zp2sos(z,p,K)完成将高阶系统分解为2阶系统的串联。 2.阅读扩展练习中的实例,学习频率分析法在MATLAB中的实现; 3.编程实现系统参数输入,绘出幅度频率响应和相位响应曲线和零、极点分布图。 四、实验内容 求系统 的零、极点和幅度频率响应和相位响应。 解析: 【代码】 num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336]; [z,p,k]=tf2zp(num,den); disp('零点');disp(z); disp('极点');disp(p); disp('增益系数');disp(k); figure(1) zplane(num,den)

figure(2) freqz(num,den,128) 【图形】 -2 -1.5 -1 -0.500.5 1 1.5 -1.5 -1 -0.5 0.51 1.5 Real Part I m a g i n a r y P a r t 0.1 0.2 0.30.40.50.60.70.80.9 1 -800 -600-400-2000 Normalized Frequency (?π rad/sample) P h a s e (d e g r e e s ) 0.1 0.2 0.30.40.50.60.70.80.9 1 -40-2002040Normalized Frequency (?π rad/sample) M a g n i t u d e (d B ) 【结果】 零点 -1.5870 + 1.4470i

数字信号实验报告

北京科技大学 《信号系统与信号处理综合实验》实验 报告 学号:__________ 姓名:_____________________ 专业:____________ 年月日

目录: 1实验一CCS使用实验 2实验二、SEED-DTK6446 Linux开发环境搭建3实验三、Linux平台实验 4二、音频采集回放实验 5三、视频采集回放实验 6OSD图像叠加实验 7图像边缘检测实验

课程实验目的 1.数字信号处理是一门理论与实践并重的课程,在学习理论知识的同时再配合经典DSP实验,可以加深对数字信号处理软、硬件的理解与掌握。 2.接触并了解SEED-DTK6446实验箱,学会通过Linux操作平台,利用SEED-DTK6446实验箱完成一些经典的实验历程,加深对数字信号处理的了解。 3. 学习并掌握SEED-DTK6446 CCS开发环境的搭建,建立好所有编译测试环境,为下面的实验做好准备工作。 实验一 CCS使用实验 一、实验目的 1.熟悉CCS3.3集成开发环境,掌握工程的生成方法; 2.熟悉SEED-DTK6446实验环境; 3. 学习用标准C 语言编制程序; 4.掌握CCS3.3集成开发环境的调试方法; 二、实验内容 1.DSP源文件的建立; 2.DSP程序工程文件的建立; 3. 学习使用CCS3.3集成开发工具的调试工具。 三、实验步骤 1.创建源文件:选择File →New →Source File 命令;打开配套光盘\03. Examples of program\01.SEEE-DTK6446 CCS Examples\examples\3.1.1 math。 2.创建工程文件:点击Project-->New,创建新工程;点击Project选择add files to project,添加源程序math.c。 3. 设置编译与连接选项:点击Project选择Build Opitions; 4. 工程编译与调试:点击Project →Build all,对工程进行编译;点击File →load program,在弹出的对话框中载入debug 文件夹下的.out可执行文件;点击debug →Go Main回到C程序的入口;运行程序并观察输出结果。 四.实验要求:

数字信号处理期末试题及答案汇总

数字信号处理期末试题及答案汇总

数字信号处理卷一 一、填空题(每空1分, 共10分) 1.序列()sin(3/5)x n n π=的周期为 。 2.线性时不变系统的性质有 律、 律、 律。 3.对4 ()()x n R n =的Z 变换为 ,其收敛域为 。 4.抽样序列的Z 变换与离散傅里叶变换DFT 的关系为 。 5.序列x(n)=(1,-2,0,3;n=0,1,2,3), 圆周左移2位得到的序列为 。 6.设LTI 系统输入为x(n) ,系统单位序列响应为h(n),则系统零状态输出y(n)= 。 7.因果序列x(n),在Z →∞时,X(Z)= 。 二、单项选择题(每题2分, 共20分) 1.δ(n)的Z 变换是 ( )A.1 B.δ(ω) C.2πδ(ω) D.2π 2.序列x 1(n )的长度为4,序列x 2(n )的长度为3,则它们线性卷积的长度是 ( ) A. 3 B. 4 C. 6 D. 7 3.LTI 系统,输入x (n )时,输出y (n );输入 为3x (n-2),输出为 ( ) A. y (n-2) B.3y (n-2) C.3y (n ) D.y (n ) 4.下面描述中最适合离散傅立叶变换DFT 的是 ( ) A.时域为离散序列,频域为连续信号 B.时域为离散周期序列,频域也为离散周期序列

C.时域为离散无限长序列,频域为连续周期信号 D.时域为离散有限长序列,频域也为离散有限长序列 5.若一模拟信号为带限,且对其抽样满足奈奎斯特条件,理想条件下将抽样信号通过即可完全不失真恢复原信号()A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 6.下列哪一个系统是因果系统()A.y(n)=x (n+2) B. y(n)= cos(n+1)x (n) C. y(n)=x (2n) D.y(n)=x (- n) 7.一个线性时不变离散系统稳定的充要条件是其系统函数的收敛域包括()A. 实轴 B.原点C.单位圆 D.虚轴 8.已知序列Z变换的收敛域为|z|>2,则该序列为()A.有限长序列 B.无限长序列C.反因果序列 D.因果序列 9.若序列的长度为M,要能够由频域抽样信号X(k)恢复原序列,而不发生时域混叠现象,则频域抽样点数N需满足的条件是( ) A.N≥M B.N≤M C.N≤2M D.N≥2M 10.设因果稳定的LTI系统的单位抽样响应h(n),在n<0时,h(n)= ( )

数字信号处理实验报告(全)

实验一、离散时间系统及离散卷积 1、单位脉冲响应 源程序: function pr1() %定义函数pr1 a=[1,-1,0.9]; %定义差分方程y(n)-y(n-1)+0.9y(n-2)=x(n) b=1; x=impseq(0,-40,140); %调用impseq函数 n=-40:140; %定义n从-40 到140 h=filter(b,a,x); %调用函数给纵座标赋值 figure(1) %绘图figure 1 (冲激响应) stem(n,h); %在图中绘出冲激 title('冲激响应'); %定义标题为:'冲激响应' xlabel('n'); %绘图横座标为n ylabel('h(n)'); %绘图纵座标为h(n) figure(2) %绘图figure 2 [z,p,g]=tf2zp(b,a); %绘出零极点图 zplane(z,p) function [x,n]=impseq(n0,n1,n2)%声明impseq函数 n=[n1:n2]; x=[(n-n0)==0]; 结果: Figure 1: Figure 2:

2、离散系统的幅频、相频的分析 源程序: function pr2() b=[0.0181,0.0543,0.0543,0.0181]; a=[1.000,-1.76,1.1829,-0.2781]; m=0:length(b)-1; %m从0 到3 l=0:length(a)-1; %l从0 到3 K=5000; k=1:K; w=pi*k/K; %角频率w H=(b*exp(-j*m'*w))./(a*exp(-j*l'*w));%对系统函数的定义 magH=abs(H); %magH为幅度 angH=angle(H); %angH为相位 figure(1) subplot(2,1,1); %在同一窗口的上半部分绘图 plot(w/pi,magH); %绘制w(pi)-magH的图形 grid; axis([0,1,0,1]); %限制横纵座标从0到1 xlabel('w(pi)'); %x座标为 w(pi) ylabel('|H|'); %y座标为 angle(H) title('幅度,相位响应'); %图的标题为:'幅度,相位响应' subplot(2,1,2); %在同一窗口的下半部分绘图 plot(w/pi,angH); %绘制w(pi)-angH的图形 grid; %为座标添加名称

数字信号处理试卷及答案

A 一、 选择题(每题3分,共5题) 1、)6 3()(π-=n j e n x ,该序列是 。 A.非周期序列 B.周期6 π = N C.周期π6=N D. 周期π2=N 2、序列)1()(---=n u a n x n ,则)(Z X 的收敛域为 。 A.a Z < B.a Z ≤ C.a Z > D.a Z ≥ 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、)()(101n R n x =,)()(72n R n x =,用DFT 计算二者的线性卷积,为使计算量尽可能的少,应使DFT 的长度N 满足 。 A.16>N B.16=N C.16

2017数字信号处理模拟题a答案

1. 两个有限长序列x1(n),0≤n ≤33和x2(n),0≤n ≤36,做线性卷积后结果的长度是 70 , 若对这两个序列做64点循环卷积,则圆周卷积结果中n= 6 至 64 为线性卷积结果。 2. 一线性时不变系统,输入为 x (n )时,输出为y (n ) ;则输入为2x (n )时,输出为 ; 输入为x (n-3)时,输出为 3. 若正弦序列x(n)=sin(30n π/120)是周期的,则周期是N= 8 4. 如果一台计算机的速度为平均每次复乘5μS ,每次复加0.5μS ,用它来计算512 点的DFT[x(n)],问直接计算需要多少时间,用FFT 运算需要多少时间。 1、 直接计算 复乘所需时间 62621510510512 1.31072T N s --=??=??= 复加所需时间()6610.51010.5105125110.130816T N N s --=???-=???= 所以12 1.441536T T T s =+= 2、用FFT 计算 复乘所需时间 66122512510log 510log 5120.0115222 N T N s --=?? =??= 复加所需时间662220.510log 0.510512log 5120.002304T N N s --=??=??= 所以120.013824T T T s =+=

6.设系统差分方程 y(n)=ay(n-1)+x(n) 其中x(n)为输入,y(n)为输出。当边界条件选为y(-1)=0时,是判断系统是否线性的、移不变的

7.用级联型结构实现以下系统函数,试问一共能构成几种级联型网络,并画出其中一种的信号流图。 ()() ()() 22 41 1.41()0.50.90.8Z Z Z H z Z Z Z +-+= -++

数字信号处理实验六-时域采样与信号的重建

实验目的: 1.了解用MATLAB语言进行时域抽样与信号重建的方法 2.进一步加深对时域信号抽样与恢复的基本原理的理解 3.掌握采样频率的确定方法和内插公式的编程方法。 二.实验内容 1认真阅读并输入实验原理与方法中介绍的例子,观察输出波形曲线,理解每一条语句的含义。. 2.已知一个连续时间信号f(t)=sinc(t)。取最高有限带宽频率fm=1Hz。(1)分别显示原连续时间信号波形和Fm=fm、Fm=2fm、Fm=3fm三种情况下抽样信号的波形。 实验程序: dt=0.1; f0=1; T0=1/f0; fm=f0; Tm=1/fm; t=-2:dt:2; f=sinc(t); subplot(4,1,1),plot(t,f,'k'); 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=sinc(n); subplot(4,1,i+1),stem(n,f,'filled','k'); axis([min(n) max(n) 1.1*min(f) 1.1*max(f)]); end 实验截图: (2)求解原连续信号波形和抽样信号所对应的幅度谱。实验程序: dt=0.1;t=-4:dt:4;

N=length(t);f=sinc(t);Tm=1;fm=1/Tm; wm=2*pi*fm;k=1:N; w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt; subplot(4,1,1),plot(w1/(2*pi),abs(F1));grid axis([0 max(4*fm) 1.1*min(F1) 1.1*max(F1)]); for i=1:3; if i<= 2 c=0 ,else c=0.2,end fs=(4-i+c)*fm; Ts=1/fs; n=-4:Ts:4; f=sinc(n); N=length(n); wm=2*pi*fs; k=1:N; w=k*wm/N; F=f*exp(-j*n'*w)*Ts; subplot(4,1,5-i),plot(w/(2*pi),abs(F),'k');grid axis([0 max(4*fm) 1.1*min(F) 1.1*max(F)]); end 实验截图:

数字信号 实验四

一、 实验目的和要求: (1)进一步掌握线性卷积的计算机编程方法,利用卷积的方法观察系统响应的时域特性。 (2)掌握循环卷积的计算机编程方法,并比较与线性卷积的差别,验证二者之间的关系。利用循环卷积的方法观察、分析系统响应的时域特性。 二、 实验内容与原理: 1.实验原理: (1)线性卷积: 线性时不变系统(Linear Time-Invariant System, or L. T. I 系统)输入、输出间的关系为:当系统输入序列为)(n x ,系统的单位脉冲响应为)(n h ,输出序列为)(n y ,则系统输出为:∑∞ -∞ =-=*=m m n h m x n h n x n y )()()()()(;上式称为线性卷积。 (2)循环卷积 设两个有限长序列)(1n x 和)(2n x ,长度分别为1N 和2N ,)()(11k X n x DFT N ?? ?→←点 )()(22k X n x DFT N ???→←点。如果)()()(21k X k X k X ?=,则 ∑---==1 21)())(()()]([)(N m N N n R m n x m x k X IDFT n x 上式称为)(1n x 和)(2n x 的循环卷积。 (3)两个有限长序列的线性卷积 序列)(1n x 和)(2n x ,长度分别为L 点和P 点,)(3n x 为这两个序列的线性卷积,则 )(3n x 为∑∞ -∞ =-= *=m m n x m x n x n x n x )()()()()(2 1 213且线性卷积)(3n x 的非零值长度为 1-+P L 点。 (4)循环卷积与线性卷积的关系 序列)(1n x 为L 点长,序列)(2n x 为P 点长,若序列)(1n x 和)(2n x 进行N 点的循环卷积)(n x c ,其结果是否等于该两序列的线性卷积)(n x l ,完全取决于循环卷积的长度。 由教材相关推导,得∑∞ -∞ =+= q N l c n R qN n x n x )()()(,也就是说,循环卷积是线性卷积 的周期延拓序列再取主值区间。 当1-+≥P L N 时循环卷积等于线性卷积,即)()(n x n x l c =;

(完整版)数字信号处理试卷及答案

江 苏 大 学 试 题 课程名称 数字信号处理 开课学院 使用班级 考试日期

江苏大学试题第2A页

江苏大学试题第3A 页

江苏大学试题第页

一、填空题:(每空1分,共18分) 8、 数字频率ω是模拟频率Ω对采样频率s f 的归一化,其值是 连续 (连续还是离散?)。 9、 双边序列z 变换的收敛域形状为 圆环或空集 。 10、 某序列的DFT 表达式为∑-== 10 )()(N n kn M W n x k X ,由此可以看出,该序列时域的长度为 N , 变换后数字频域上相邻两个频率样点之间的间隔是 M π 2 。 11、 线性时不变系统离散时间因果系统的系统函数为2 52) 1(8)(22++--=z z z z z H ,则系统的极点为 2,2 1 21-=-=z z ;系统的稳定性为 不稳定 。系统单位冲激响应)(n h 的初值4)0(=h ; 终值)(∞h 不存在 。 12、 如果序列)(n x 是一长度为64点的有限长序列)630(≤≤n ,序列)(n h 是一长度为128点的有限长 序列)1270(≤≤n ,记)()()(n h n x n y *=(线性卷积),则)(n y 为 64+128-1=191点 点的序列,如果采用基FFT 2算法以快速卷积的方式实现线性卷积,则FFT 的点数至少为 256 点。 13、 用冲激响应不变法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之间的映射变换 关系为T ω = Ω。用双线性变换法将一模拟滤波器映射为数字滤波器时,模拟频率Ω与数字频率ω之 间的映射变换关系为)2tan(2ωT = Ω或)2 arctan(2T Ω=ω。 当线性相位FIR 数字滤波器满足偶对称条件时,其单位冲激响应)(n h 满足的条件为)1()(n N h n h --= ,

数字信号处理实验报告 六

程序P6.1 % 程序 P6_1 % 将一个有理数传输函数 % 转换为因式形式 num = input('分子系数向量 = '); den = input('分母系数向量 = '); [z,p,k] = tf2zp(num,den); sos = zp2sos(z,p,k) Q6.1 使用程序p6.1,生成如下有限冲激响应传输函数的一个级联实现:H1(Z)=2+10Z^-1+23Z^-2+34Z^-3+31Z^-4+16Z^-5+4Z^-6 画出级联实现的框图。H1(Z)是一个线性相位传输函数吗? 分子系数向量 = [2,10,23,34,31,16,4] 分母系数向量 = [1,0,0,0,0,0,0] sos = 2.0000 6.0000 4.0000 1.0000 0 0 1.0000 1.0000 2.0000 1.0000 0 0 1.0000 1.0000 0.5000 1.0000 0 0 Y[k] 2 11 X[k] Q6.2 使用程序p6.1,生成如下有限冲激响应传输函数的一个级联实现:H2(Z)=6+31Z^-1+74Z^-2+102Z^-3+74Z^-4+31Z^-5+6Z^-6 画出级联实现的框图。H2(Z )是一个线性相位传输函数吗?只用4个乘法器生成H2(Z)的一个级联实现。显示新的级联结构的框图。 分子系数向量 = [6,31,74,102,74,31,6] 分母系数向量 = [1,0,0,0,0,0,0] sos =

6.0000 15.0000 6.0000 1.0000 0 0 1.0000 2.0000 3.0000 1.0000 0 0 1.0000 0.6667 0.3333 1.0000 0 0 Y[k] 6 11 X[k] Q6.3 使用程序 6.1生成如下因果无限冲激响应传输函数的级联实现: H1(Z)=(3+8Z^-1+12Z^-2+7Z^-3+2Z^-4-2Z^-5)/(16+24Z^-1+24Z^-2+14Z^-3+5Z^-4+Z^-5),画出级联实现的框图。 分子系数向量 = [3,8,12,7,2,-2] 分母系数向量 = [16,24,24,14,5,1] sos = 0.1875 -0.0625 0 1.0000 0.5000 0 1.0000 2.0000 2.0000 1.0000 0.5000 0.2500 1.0000 1.0000 1.0000 1.0000 0.5000 0.5000 Y [k ] 0.1875 11 11X [k ] Q6.4使用程序6.1生成如下因果无限冲激响应传输函数的级联实现:

数字信号实验处理实验报告

中国石油大学(北京) China University of Petroleum 基于MATLAB的数字 信号处理 实验名称:DFT和DCT的应用 姓名: 学号: 日期:2017.4.10

1. 用FFT 计算线性卷积(预习实验指导书p9内容) 两个序列:x1=[1 2 3 4 5 6]; x2=[3 2 1 1 2 3 4],用DFT 求两者的循环卷积, (1) 取FFT 长度L=10,问哪些点上的循环卷积=线性卷积? (2) 取FFT 长度L=12,结果又如何? 程序代码: x1=[1 2 3 4 5 6];x2=[3 2 1 1 2 3 4];Y1=conv(x1,x2); N=10;m=length(x1)+length(x2)-1;n=0:m-1; subplot(311); stem(n,Y1); title('线性卷积结果'); x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x1))]; xk1=fft(x1,N); xk2=fft(x2,N); Yk=xk1.*xk2; y1=ifft(Yk,N); subplot(312); stem(0:N-1,y1); title('长度取10的循环卷积结果'); N=12; x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x1))]; xk1=fft(x1,N); xk2=fft(x2,N); Yk1=xk1.*xk2; y2=ifft(Yk1,N); subplot(313); stem(0:N-1,y2); title('长度取12的循环卷积结果'); 结果: 2. 探地雷达信号谱分析 线性卷积结果 长度取10的循环卷积结 果 长度取12的循环卷积结果

数字信号处理试卷大全..

北京信息科技大学 2010 ~2011 学年第一学期 《数字信号处理》课程期末考试试卷(A) 一、填空题(本题满分30分,共含4道小题,每空2分) 1.两个有限长序列x1(n),0≤n≤33和x2(n),0≤n≤36,做线性卷积 后结果的长度是,若对这两个序列做64点圆周卷积,则圆周卷积结果中n= 至为线性卷积结果。 W的、和三个固有特性来实现2.DFT是利用nk N FFT快速运算的。 3.IIR数字滤波器设计指标一般由、、和等 四项组成。 4.FIR数字滤波器有和两种设计方法,其结构 有、和等多种结构。 二、判断题(本题满分16分,共含8道小题,每小题2分,正 确打√,错误打×) 1.相同的Z变换表达式一定对应相同的时间序列。() 2.Chirp-Z变换的频率采样点数M可以不等于时域采样点数N。() 3.按频率抽取基2 FFT首先将序列x(n)分成奇数序列和偶数序列。() 4.冲激响应不变法不适于设计数字带阻滤波器。() 5.双线性变换法的模拟角频率Ω与数字角频率ω成线性关系。() 6.巴特沃思滤波器的幅度特性必在一个频带中(通带或阻带)具有等

波纹特性。( ) 7. 只有FIR 滤波器才能做到线性相位,对于IIR 滤波器做不到线性相 位。( ) 8. 在只要求相同的幅频特性时,用IIR 滤波器实现其阶数一定低于 FIR 阶数。( ) 三、 综合题(本题满分18分,每小问6分) 若x (n)= {3,2,1,2,1,2 },0≤n≤5, 1) 求序列x(n)的6点DFT ,X (k)=? 2) 若)()]([)(26k X W n g DFT k G k ==,试确定6点序列g(n)=? 3) 若y(n) =x(n)⑨x(n),求y(n)=? 四、 IIR 滤波器设计(本题满分20分,每小问5分) 设计一个数字低通滤波器,要求3dB 的截止频率f c =1/π Hz ,抽样频率f s =2 Hz 。 1. 导出归一化的二阶巴特沃思低通滤波器的系统函数H an (s)。 2. 试用上述指标设计一个二阶巴特沃思模拟低通滤波器,求其系 统函数H a (s),并画出其零极点图。 3. 用双线性变换法将H a (s)转换为数字系统的系统函数H(z)。 4. 画出此数字滤波器的典范型结构流图。 五、 FIR 滤波器设计(本题满分16分,每小问4分)

数字信号处理完整试题库

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分

相关文档
最新文档