信号与系统实验四 答案

合集下载

信号与系统实验(MATLAB 西电版)实验4 离散时间信号的时域基本运算_OK

信号与系统实验(MATLAB 西电版)实验4  离散时间信号的时域基本运算_OK

图 4.5 序列及其平移
2021/7/3
16
实验4 离散时间信号的时域基本运算
2.
已知序列f(k)={2,3,1,2,3,4,3,1},对应的k值为 -3≤k≤4 f1(k)=f(k-2),f2(k)=f(-k),f3(k)=f(k-1)ε(k), f4(k)=f(-k+2),f5(k)=f(k+1), f6(k)=f(k-2)ε(k),f7(k)=f(k+2)ε(k)
5) MATLAB x1=-2:2; %序列1 k1=-2:2; k0=2; k=k1+k0; f=x1; stem(k,f,′filled′); axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]); 序列及其平移如图4.5
2021/7/3
15
实验4 离散时间信号的时域基本运算
2021/7/3
17
实验4 离散时间信号的时域基本运算
(1) 在计算机中输入程序,验证并记录实验结果,经过 (2) 对于设计性实验,应自行编制完整的实验程序,重复 验证性实验的过程,并在实验报告中给出完整的自编程序。
2021/7/3
18
axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]);
序列及其翻转如图4.3所示。
2021/7/3
11
实验4 离散时间信号的时域基本运算
图 4.3 序列及其翻转
2021/7/3
12
实验4 离散时间信号的时域基本运算
4)
MATLAB
x1=-2:2;
%序列1
两个序列的乘法如图4.2
2021/7/3

信号与系统实验四实验报告

信号与系统实验四实验报告

实验四 时域抽样与频域抽样一、实验目的加深理解连续时间信号的离散化过程中的数学概念和物理概念,掌握时域抽样定理的基本内容。

掌握由抽样序列重建原连续信号的基本原理与实现方法,理解其工程概念。

加深理解频谱离散化过程中的数学概念和物理概念,掌握频域抽样定理的基本内容。

二、 实验原理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件:对于基带信号,信号抽样频率sam f 大于等于2倍的信号最高频率m f ,即m sam f f 2≥。

时域抽样是把连续信号x (t )变成适于数字系统处理的离散信号x [k ] ;信号重建是将离散信号x [k ]转换为连续时间信号x (t )。

非周期离散信号的频谱是连续的周期谱。

计算机在分析离散信号的频谱时,必须将其连续频谱离散化。

频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。

三.实验内容1. 为了观察连续信号时域抽样时抽样频率对抽样过程的影响,在[0,0.1]区间上以50Hz 的抽样频率对下列3个信号分别进行抽样,试画出抽样后序列的波形,并分析产生不同波形的原因,提出改进措施。

)102cos()(1t t x ⨯=π答: 函数代码为: t0 = 0:0.001:0.1;x0 =cos(2*pi*10*t0);plot(t0,x0,'r')hold onFs =50;t=0:1/Fs:0.1;x=cos(2*pi*10*t); stem(t,x); hold offtitle('连续信号及其抽样信号')函数图像为:)502cos()(2t t x ⨯=π同理,函数图像为:)0102cos()(3t t x ⨯=π同理,函数图像为:由以上的三图可知,第一个图的离散序列,基本可以显示出原来信号,可以通过低通滤波恢复,因为信号的频率为20HZ,而采样频率为50>2*20,故可以恢复,但是第二个和第三个信号的评论分别为50和100HZ,因此理论上是不能够恢复的,需要增大采样频率,解决的方案为,第二个信号的采样频率改为400HZ,而第三个的采样频率改为1000HZ,这样可以很好的采样,如下图所示:2. 产生幅度调制信号)200cos()2cos()(t t t x ππ=,推导其频率特性,确定抽样频率,并绘制波形。

信号与系统课后习题参考答案

信号与系统课后习题参考答案

1试分别指出以下波形是属于哪种信号?题图1-11-2 试写出题1-1 图中信号的函数表达式。

1-3 已知信号x1(t)与x2(t)波形如题图1-3 中所示,试作出下列各信号的波形图,并加以标注。

题图1-3⑴x1(t2)⑵ x1(1 t)⑶ x1(2t 2)⑷ x2(t 3)⑸ x2(t 2) ⑹x2(1 2t)2⑺x1(t) x2( t)⑻x1(1 t)x2(t 1)⑼x1(2 t) x2(t 4)21- 4 已知信号x1(n)与x2 (n)波形如题图1-4中所示,试作出下列各信号的波形图,并加以标注。

题图1-4⑴x1(2n 1) ⑵ x1(4 n)⑶ x1(n)2⑷ x2 (2 n)⑸ x2(n 2) ⑹ x2(n 2) x2( n 1)⑺x1(n 2) x2(1 2n)⑻x1(1 n) x2(n 4)⑼ x1(n 1) x2(n 3)1- 5 已知信号x(5 2t )的波形如题图1-5 所示,试作出信号x(t)的波形图,并加以标注。

题图1-51- 6 试画出下列信号的波形图:1⑴ x(t) sin( t) sin(8 t)⑵ x(t) [1 sin( t )] sin(8 t)21⑶x(t) [1 sin( t)] sin(8 t)⑷ x(t) sin( 2t )1-7 试画出下列信号的波形图:⑴ x(t)1 e t u(t) ⑵ x(t) e t cos10 t[u(t 1) u(t 2)]⑶ x(t)(2 e t)u(t)⑷ x(t) e (t 1)u(t)⑸ x(t)u(t22 9) ⑹ x(t)(t2 4)1-8 试求出以下复变函数的模与幅角,并画出模与幅角的波形图1j2 ⑴ X (j ) (1 e j2)⑵ X( j1 e j4⑶ X (j ) 11 ee j ⑷ X( j )试作出下列波形的奇分量、偶分量和非零区间上的平均分量与交流分量。

题图 1-10形图。

题图 1-141-15 已知系统的信号流图如下,试写出各自系统的输入输出方程。

(仅供参考)信号与系统第四章习题答案

(仅供参考)信号与系统第四章习题答案

e −sT
=
−sT
2 − 4e 2
+ 2e −sT
Ts 2
(f) x(t) = sin πt[ε (t)− ε (t − π )]
sin π tε (t ) ↔
π s2 + π 2
L[sin
πtε (t
−π
)]
=
L e jπt
− 2
e− jπt j
ε (t
−π
)
∫ ∫ =
1 2j
∞ π
e
jπt e−st dt
4.3 图 4.2 所示的每一个零极点图,确定满足下述情况的收敛域。
(1) f (t) 的傅里叶变换存在
(2) f (t )e 2t 的傅里叶变换存在
(3) f (t) = 0, t > 0
(4) f (t) = 0, t < 5
【知识点窍】主要考察拉普拉斯变换的零极点分布特性。 【逻辑推理】首先由零极点写出拉普拉斯变换式,再利用反变换求取其原信号,即可求取其收
= cosϕ eω0tj + e−ω0tj − sin ϕ eω0tj − e−ω0tj
2
2j
=
cos 2
ϕ

sin 2
ϕ j
e
ω0 t j
+
cosϕ 2
+
sin ϕ 2j
e −ω 0tj
F(s) =
L
cosϕ 2

sin ϕ 2j
eω0tj
+
cos 2
ϕ
+
sin ϕ 2j
e
−ω0
t
j
ε
(t
)
∫ ∫ =

信号与系统第四章习题参考答案13

信号与系统第四章习题参考答案13

《信号与系统》第四章习题参考答案4-1 解 (1)111()ataL es s a s s a -⎡⎤-=-=⎣⎦++ (2)[]2221221sin 2cos 111s s L t t s s s ++=+++++ (3)()2212tL te s -⎡⎤=⎣⎦+(4)[]21sin(2)4L t s =+,由S 域平移性质,得 ()21s i n (2)14tL e t s -⎡⎤=⎣⎦++ (5)因为1!nn n L t s +⎡⎤=⎣⎦,所以 []2211212s L t s s s++=+= 由S 域平移性质,得 ()()23121ts L t e s -+⎡⎤+=⎣⎦+(6)()2211cos sL at s s a -=-⎡⎤⎣⎦+,由S 域平移性质,得 (){}()2211cos ts L at e s s aβββ-⎡⎤-=-⎣⎦+++ (7)232222L t t s s ⎡⎤+=+⎣⎦ (8)732()327tL t es δ-⎡⎤-=-⎣⎦+ (9)[]22sinh()L t s βββ=-,由S 域平移性质,得()22sinh()atL e t s a βββ-⎡⎤=⎣⎦+-(10)由于()211cos ()cos 222t t Ω=+Ω 所以 222221111c o s ()22424ss L t s s s s ⎛⎫⎡⎤Ω=+∙=+ ⎪⎣⎦+Ω+Ω⎝⎭(11)()()()11111at t L e e a a s a s s a s βββββ--⎡⎤⎛⎫-=-= ⎪⎢⎥--++++⎣⎦⎝⎭ (12)由于()221cos()1ts L e t s ωω-+⎡⎤=⎣⎦++所以 ()()()221cos()1a t a s e L et s ωω--++⎡⎤=⎣⎦++(13)因为(2)(1)(1)(1)(1)(1)t t t te u t e t e e u t ------⎡⎤-=-+-⎣⎦且()(1)(1)2(1)(1)(1)11sst t e e L t eu t L eu t s s ------⎡⎤⎡⎤--=-=⎣⎦⎣⎦++所以 ()(1)(2)2211(2)(1)(1)11s t s s e L teu t e e s s s -----⎡⎤+⎡⎤-=+=⎢⎥⎣⎦+++⎣⎦(14)()(1)tL e f t F s -⎡⎤=+⎣⎦,由尺度变换性质,得(1)ta t L e f aF as a -⎡⎤⎛⎫=+⎢⎥ ⎪⎝⎭⎣⎦(15)()t L f aF as a ⎡⎤⎛⎫=⎪⎢⎥⎝⎭⎣⎦,再由s 域平移性质,得 []2()()at t L e f aF a s a aF as a a -⎡⎤⎛⎫=+=+ ⎪⎢⎥⎝⎭⎣⎦(16)31cos(6)cos (3)cos(3)2t t t -=∙13cos(9)cos(3)44t t =+32213cos (3)48149s s L t s s ⎡⎤=+⎣⎦++由s 域微分性质,得()()22322222213181327cos (3)481494819d s s s s L t t ds s s s s ⎡⎤--⎛⎫⎢⎥⎡⎤=-+=+ ⎪⎣⎦⎢⎥++⎝⎭++⎣⎦(17)[]2cos(2)4sL t s =+,连续两次应用s 域微分性质,有 []()2224cos(2)4s L t t s-=+,()3232224cos(2)4s sL t t s-⎡⎤=⎣⎦+(18)111atL es s a -⎡⎤-=-⎣⎦+,由s 域积分性质,得111111(1)at sL e ds t s s a ∞-⎛⎫⎡⎤-=- ⎪⎢⎥+⎣⎦⎝⎭⎰ln()ln ln s s a s s a ⎛⎫=+-=- ⎪+⎝⎭ (19)351135tt L ee s s --⎡⎤-=-⎣⎦++,由s 域积分性质,得 33111115ln 353t t s e e s L ds t s s s --∞⎛⎫⎡⎤-+⎛⎫=-= ⎪ ⎪⎢⎥+++⎝⎭⎣⎦⎝⎭⎰(20)()22sin aL at s a =⎡⎤⎣⎦+,由s 域积分性质,得()1122211sin 1arctan 21s s at s a s L ds d t s a a a s a π∞∞⎡⎤⎛⎫⎛⎫===-⎢⎥ ⎪ ⎪+⎝⎭⎝⎭⎛⎫⎣⎦+ ⎪⎝⎭⎰⎰ 4-2 解(1)因为()()sin ()2T f t t u t u t ω⎡⎤⎛⎫=--⎪⎢⎥⎝⎭⎣⎦()sin ()sin 22T T t u t t u t ωω⎡⎤⎛⎫⎛⎫=+-- ⎪ ⎪⎢⎥⎝⎭⎝⎭⎣⎦ 所以可借助延时定理,得()()sin ()sin 22T T L f t L t u t L t u t ωω⎧⎫⎡⎤⎛⎫⎛⎫=+--⎡⎤⎡⎤⎨⎬ ⎪ ⎪⎢⎥⎣⎦⎣⎦⎝⎭⎝⎭⎣⎦⎩⎭222222221sT T s ee S S S ωωωωωω--⎛⎫=+=+ ⎪+++⎝⎭(2)因为()()()sin sin cos cos sin t t t ωϕωϕωϕ+=+ 所以()222222cos sin cos sin sin s s L t s s s ωϕϕωϕϕωϕωωω++=+=⎡⎤⎣⎦+++ 4-3 解此题可巧妙运用延时性质。

信号与系统课后习题答案

信号与系统课后习题答案

习 题 一 第一章习题解答基本练习题1-1 解 (a) 基频 =0f GCD (15,6)=3 Hz 。

因此,公共周期3110==f T s 。

(b) )30cos 10(cos 5.0)20cos()10cos()(t t t t t f ππππ+==基频 =0f GCD (5, 15)=5 Hz 。

因此,公共周期5110==f T s 。

(c) 由于两个分量的频率1ω=10π rad/s 、1ω=20 rad/s 的比值是无理数,因此无法找出公共周期。

所以是非周期的。

(d) 两个分量是同频率的,基频 =0f 1/π Hz 。

因此,公共周期π==01f T s 。

1-2 解 (a) 波形如图1-2(a)所示。

显然是功率信号。

t d t f TP T TT ⎰-∞→=2)(21lim16163611lim 22110=⎥⎦⎤⎢⎣⎡++=⎰⎰⎰∞→t d t d t d T T T W(b) 波形如图1.2(b)所示。

显然是能量信号。

3716112=⨯+⨯=E J (c) 能量信号 1.0101)(lim101025=-===⎰⎰∞∞---∞→T t ttT e dt edt eE J(d) 功率信号,显然有 1=P W1-3 解 周期T=7 ,一个周期的能量为 5624316=⨯+⨯=E J 信号的功率为 8756===T E P W 1-5 解 (a) )(4)2()23(2t tt δδ=+; (b) )5.2(5.0)5.2(5.0)25(5.733-=-=----t e t e t et tδδδ(c) )2(23)2()3sin()2()32sin(πδπδπππδπ+-=++-=++t t t t 题解图1-2(a) 21题解图1-2(b) 21(d) )3()3()(1)2(-=----t e t t et δδε。

1-6 解 (a) 5)3()94()3()4(2-=+-=+-⎰⎰∞∞-∞∞-dt t dt t t δδ(b) 0)4()4(632=+-⎰-dt t t δ(c) 2)]2(2)4(10[)]42(2)4()[6(63632=+++-=+++-⎰⎰--dt t t dt t t t δδδδ(d)3)3(3)(3sin )(1010=⋅=⎰⎰∞-∞-dt t Sa t dt ttt δδ。

信号与系统(第四版)第四章课后答案

信号与系统(第四版)第四章课后答案

第5-10页

©西安电子科技大学电路与系统教研中心
信号与系统 电子教案
4.1 拉普拉斯变换
四、常见函数的单边拉普拉斯变换
1. (t ) 1, 2.( t) 或1 3. ( t ) s, 4. 指数信号e
1
s
, 0

1 s s0
s0t
(t 2)
f1(t) 1 0 1 f2(t) 1 t
例1:e (t 2) e
-t
2
e
(t 2)
e
2

1 s 1
e
2s
-1 0
第5-17页

1
t
©西安电子科技大学电路与系统教研中心
信号与系统 电子教案
4.2 拉普拉斯变换性质
1 1e sT
例2: 单边冲激 T(t ) 1 e sT e s 2T 例3: 单边周期信号 fT(t ) (t ) f1(t ) f1(t T ) f1(t 2T ) F1(s )(1 e sT e s 2T )
8 e 2 s
s
f(t ) 1 0 1 y(t ) 2 4 t
二、尺度变换
2s
2
(1 e 2 s 2s e 2 s )
2 e 2 s 2 (1 e 2 s 2s e 2 s ) s
第5-16页

0
2
4
t
©西安电子科技大学电路与系统教研中心
信号与系统 电子教案
拉氏逆变换的物理意义
f (t )
2 j 1

j
j
F (s)est ds

信号与系统实验四-信号的采样及恢复

信号与系统实验四-信号的采样及恢复

信号与系统实验四-信号的采样及恢复实验四信号的采样及恢复⼀、实验⽬的1、加深理解连续时间信号离散化过程中的数学概念和物理概念;2、掌握对连续时间信号进⾏抽样和恢复的基本⽅法;3、通过实验验证抽样定理。

⼆、实验内容1、为了观察连续信号时域抽样时,抽样频率对抽样过程的影响,在[0,0.1]区间上以50Hz 的抽样频率对下列3个信号分别进⾏抽样,试画出抽样后序列的波形,并分析产⽣不同波形的原因,提出改进措施。

(1))102cos()(1t t x ?=π(2))502cos()(2t t x ?=π(3))1002cos()(3t t x ?=π2、产⽣幅度调制信号)200cos()2cos()(t t t x ππ=,推导其频率特性,确定抽样频率,并绘出波形。

3、对连续信号)4cos()(t t x π=进⾏抽样以得到离散序列,并进⾏重建。

(1)⽣成信号)(t x ,时间t=0:0.001:4,画出)(t x 的波形。

(2)以10=sam f Hz 对信号进⾏抽样,画出在10≤≤t 范围内的抽样序列)(k x ;利⽤抽样内插函数)/1()(sam r f T T t Sa t h =??=π恢复连续信号,画出重建信号)(t x r 的波形。

)(t x 与)(t x r 是否相同,为什么?(3)将抽样频率改为3=sam f Hz ,重做(2)。

4、利⽤MATLAB 编程实现采样函数Sa 的采样与重构。

三、实验仪器及环境计算机1台,MATLAB7.0软件。

四、实验原理对连续时间信号进⾏抽样可获得离散时间信号,其原理如图8-1。

采样信号)()()(t s t f t f s ?=,)(t s 是周期为s T 的冲激函数序列,即)()()(∑∞-∞=-==n sT nT t t t s sδδ则该过程为理想冲激抽样。

其中s T 称为采样周期,ss T f 1=称为抽样频率, ss s T f π⼤于等于2倍的原信号频率m f 时,即m s f f 2≥(抽样时间间隔满⾜ms f T 21≤),抽样信号的频谱才不会发⽣混叠,可⽤理想低通滤波器将原信号从采样信号中⽆失真地恢复。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验四 基于窗函数的FIR DF 的设计
提示:
1. Matlab 中提供了很多常用的窗函数,其中一些窗函数的调用形式为:
矩形窗:w=boxcar(N)
三角形窗:w=bartlett(N) 汉宁窗:w=hanning(N) 哈明窗:w=hamming(N) 布莱克曼窗:w=blackman(N)
其中,输入参数N 表示窗口的长度,返回的变量w 是一个长度为N 的列向量,表示窗函数在这N 点的取值。

2. b=fir1(N,Wc,'ftype',Window)
fir1函数用来设计FIR 滤波器。

其中N 为滤波器的阶数;Wc 是截止频率,其取值在0~1之间,它是以π为基准频率的标称值,设计低通和高通滤波器时,Wc 是标量,设计带通和带阻滤波器时,Wc 是1×2的向量;设计低通和带通滤波器时,无需 'ftype',当ftype=high 时,设计高通滤波器,当ftype=stop 时,设计带阻滤波器;Window 表示设计滤波器所采用的窗函数类型,Window 的长度为N+1,若Window 缺省,则fir1默认使用哈明窗;b 对应设计好的滤波器的系数h(n),即单位冲激响应,h(n)的长度为N+1。

需注意)(n h 的长度与滤波器的阶数间的关系。

FIR 滤波器的系统函数可表示为:
∑-=-=1
)()(N n n z n h z H
)(n h 的长度为N ,而滤波器的阶数为1-N 阶。

3. 求数字滤波器的频率响应
h=freqz(b,a,w)
其中,b 和a 分别为系统函数)(z H 的分子多项式和分母多项式的系数。

对于FIR 滤波器,此处的b 即为h(n),a 可看作1。

实验题目:
1. 分别用矩形窗和哈明窗设计FIR 低通滤波器,设窗宽11=N ,截止频率rad c πω
2.0=,要求绘出两种窗函数设计的滤波器幅频曲线,并进行比较。

2. 设计一个线性相位FIR 低通滤波器,通带截止频率为rad p πω2.0=,阻带截止频率为
rad s πω4.0=,阻带最小衰减为dB s 50=α。

要求分别绘制理想冲激响应,窗函数的时域
波形及幅频特性,实际冲激响应,FIR 滤波器的幅频特性和相频特性。

0.51
0.2
0.4
0.6
0.8
1
1.2
w/pi
幅频响应
0.51
-120
-100-80-60-40
-20
020w/pi
d B
N=11;
h1=fir1(N-1,0.2,boxcar(N)); h2=fir1(N-1,0.2,hamming(N)); w=0:0.01:pi; H1=freqz(h1,1,w);
H1db=20*log10(abs(H1)/max(abs(H1))); H2=freqz(h2,1,w);
H2db=20*log10(abs(H2)/max(abs(H1))); subplot(1,2,1);
plot(w/pi,abs(H1),'-.',w/pi,abs(H2)); legend('矩形窗','哈明窗'); xlabel('w/pi');ylabel('幅频响应');
subplot(1,2,2);
plot(w/pi,H1db,'-.',w/pi,H2db); xlabel('w/pi');ylabel('dB'); legend('矩形窗','哈明窗');
2040
理想单位响应n h d (n )
2040哈明窗
n
w (n )
2040
实际单位响应
n
h (n )
0.51
-150
-100
-50
0幅频响应
w/pi
分贝数
0.5
1
-200
-150
-100
-500哈明窗幅频特性00.51
-1500
-1000
-500
滤波器相频特性
wp=0.2*pi; ws=0.4*pi; tr_width=ws-wp; N=ceil(8*pi/tr_width); n=0:N-1; wc=(ws+wp)/2; alpha=(N-1)/2; m=n-alpha;
hd=sin(wc*m)./(pi*m); w_ham=(hamming(N))';
h=hd.*w_ham;
w=0:0.01:pi;
H=freqz(h,1,w);
H1=freqz(hd,1,w);
dbH=20*log10(abs(H)/max(abs(H)));
dbH1=20*log10(abs(H1)/max(abs(H1)));
subplot(231);stem(n,hd,'.');title('理想单位响应');
xlabel('n');ylabel('hd(n)');
subplot(232);stem(n,w_ham,'.');title('哈明窗');
xlabel('n');ylabel('w(n)');
subplot(233);stem(n,h,'.');title('实际单位响应');
xlabel('n');ylabel('h(n)');
subplot(234);plot(w/pi,dbH);title('幅频响应');
xlabel('w/pi');ylabel('分贝数');
subplot(235);plot(w/pi,dbH1);title('哈明窗幅频特性');
subplot(236);plot(w/pi,180/pi*unwrap(angle(H)));title('滤波器相频特性');。

相关文档
最新文档