第数字信号处理讲义--3章_连续时间信号的采样剖析

第数字信号处理讲义--3章_连续时间信号的采样剖析
第数字信号处理讲义--3章_连续时间信号的采样剖析

第3章连续时间信号的采样

[教学目的]

1.理解周期采样的原理,掌握采样的频域表示法,奈奎斯特采样定理;

2.掌握样本重构带限信号的原理与条件;

3.掌握连续信号转换成离散信号的方法,理想低通滤波器特点,冲激响应的概念;

4.掌握离散时间信号的连续时间处理方法;

5.了解量化误差产生的原因和影响。

[教学重点与难点]

重点:

1.采样的频域表示法,奈奎斯特采样定理;

2.样本重构带限信号;

难点:

1.采样的频域表示法;

2.样本重构带限信号;

3.1周期采样

在某些合理条件限制下,一个连续时间信号能用其采样序列来完全给予表示,连续时间信号的处理往往是通过对其采样得到的离散时间序列的处理来完成的。本节将详细讨论采样过程,包括信号采样后,信号的频谱将发生怎样的变换,信号内容会不会丢失,以及由离散信号恢复成连续信号应该具备哪些条件等。采样的这些性质对离散信号和系统的分析都是十分重要的。要了解这些性质,让我们首先从采样过程的分析开始。

采样器可以看成是一个电子开关,它的工作原理可由图3-1(a)来说明。设开关每隔T秒短暂地闭合一次,将连续信号接通,实现一次采样。如果开关每次闭合的时间为τ秒,那么采样器的输出将是一串周期为T,宽度为τ的脉冲。而脉冲的幅度却是重复着在这段τ时间内信号的幅度。如果以xa(t)代表输入的连续信号,如图3-1(b)所示,以xp(t)表示采样输出信号,它的结构如图3-1(d)所示。显然,这个过程可以把它看作是一个脉冲调幅过程。被调制的脉冲载波是一串周期为T、宽度为τ的矩形脉冲信号,如图3-1(c)所示,并以p(t)表示,而调制信号就是输入的连续信号。因而有

t

x

)(

x

)(

t

)(t

p

p

a

一般开关闭合时间都是很短的,而且τ越小,采样输出脉冲的幅度就越准确地反映输入信号在离散时间点上的瞬时值。当τ<

图 3-1 连续时间信号的采样过程

3.2 采样的频域表示 1.理想采样

理想采样就是假设采样开关闭合时间无限短,即τ→0的极限情况。此时,采样脉冲序列p(t)变成冲激函数序列s(t),如图 3-1(e)所示。这些冲激函数准确地出现在采样瞬间,面积为1。采样后,输出理想采样信号的面积(即积分幅度)则准确地等于输入信号xa(t)在采样瞬间的幅度。理想采样过程如图3-1(f )所示。 冲激函数序列s(t)为 3-1式

以 表示理想采样的输出,以后我们都以下标a 表示连续信

号(或称模拟信号),如x a (t);而以它的顶部符号(∧)表示它的理

想采样,如

。这样我们就可将理想采样表示为 3-2式

代入3-1式得: 3-3式

由于δ(t-nT)只在t=nT 时不为零,故

(a )(b )

x a

)

(c )(e )(d )( f )

)

()(nT t t s n -=∑∞

-∞=δ)(?t x a )(?t x

a )()()(?t s t x t x a a =)()()(?nT t t x t x

a n a -=∑∞

-∞=δ)

()()(?nT t nT x t x a

n a

-=∑

-∞

2 理想采样信号的频谱

我们首先看看通过理想采样后信号频谱发生了什么变化。 由于在连续时间信号与系统中已学过,式(1-17)表示时域相乘, 则频域(傅里叶变换域)为卷积运算。 所以由式(1-17)可知, 若各个信号的傅里叶变换分别表示为:

则应满足

3-5式

现在来求S(j Ω)=F [s(t)]。 由于s(t)是以采样频率重复的冲激脉冲,因此是一个周期函数,可表示为傅里叶级数,即

此级数的基频为采样频率,即: 一般称fs 为频率,单位为赫兹(Hz),Ωs 为角频率,单位为弧度/秒; 习惯上都统称为“频率”。 它们的区别由符号f 及Ω来识别。 根据傅氏级数的知识,系数ak 可以通过以下运算求得

以上结果的得出是考虑到在|t|≤T/2的积分区间内,只有一个冲激脉冲δ(t),其他冲激δ(t-nT),n ≠0 都在积分区间之外,且利用了以下关系: 因而

由此得出

dt e t x j X dt

e t s j S dt

e t x j X t j a a t j t j a a Ω-∞∞

-Ω-∞

∞-Ω-∞

-???=Ω=Ω=Ω)(?)(?)()()()()

()(21)(?Ω*Ω=Ωj S j X j X a a π∑∞

-∞=Ω=k t

jk k s

e a t s )(T

f s 1=s

s f T ππ

22==ΩT dt e t T dt e nT t T dt e t s T a T T t jk n t jk T T t

jk T T k s

s

s

1

)(1)(1)(12/2/2/2/2/2/?∑??-Ω-∞-∞

=Ω--Ω--==-==δδ?∞

∞-=dt

t t f f )()()0(δ∑∞-∞

=Ω=k t jk s

e T t s 1)([]

∞-∞

=Ω∞-∞

=Ω=??

???

?==Ωk t jk k t jk s

s

e F T

e T F t s F j S 11)]([)(

所以

3-6式

将式(3-6)代入式(3-5)可得

根据冲激函数的性质,可得 3-7式 或者

由此看出,一个连续时间信号经过理想采样后,其频谱将沿着频率轴以采样频率Ωs=2π/T 为间隔而重复,这就是说频谱产生了周期性延拓,如图3-2所示。也就是说,理想采样信号的频谱,是Xa(j Ω)的周期延拓函数,其周期为Ωs ,而频谱的幅度则受1/T 加权,由于T 是常数,所以除了一个常数因子外,每一个延拓的谱分量都和原频谱分量相同。因此只要各延拓分量与原频谱分量不发生频率混叠,则有可能恢复出原信号。也就是说,如果xa(t)是限带信号,其频谱如图3-2(a )所示,且最高频谱分量Ωh 不超过Ωs/2,即

那么原信号的频谱和各次延拓分量的谱彼此不重叠,如图3-2(c )所示。 这时采用一个截止频率为Ωs/2的理想低通滤波器, 就可得到不失真的原信号频谱。也就是说,可以不失真地还原出原来的连续信号。

)(2][s jk k e F s

Ω-Ω=Ωπδ∑∑∞

-∞=∞-∞=Ω-ΩΩ=Ω-Ω=Ωk s

s s k k k T j S )()(2)(δδπθθδθθ

θδθδππd k j X T d k j X T j X k T j X s a k k s a k a s a )()(1)()(1)()(221)(?-Ω-Ω=-Ω-Ω=??

????Ω*Ω-Ω=Ω∑?∑?∑∞-∞=∞∞-∞-∞=∞

∞-∞

-∞=∑∞

-∞=Ω-Ω=Ωk s a a jk j X T j X )

(1)(?∑∞-∞=??

? ??-Ω=Ωk a a T jk j X T j X π21)(???????

?

Ω≥ΩΩ<ΩΩ=Ω2||02

||)()(s s a a j X j X

图 3-2 时域采样后, 频谱的周期延拓

(a )原始限带信号频谱; (b ) 采样函数频谱

;

(c )已采样信号频谱(Ωs>2Ωh );(d )已采样信号频谱(Ωs<2Ωh ) 如果信号的最高频谱Ωh 超过Ωs/2,则各周期延拓分量产生频谱的交叠,称为混叠现象,如图1-10(d )所示。由于Xa(j Ω)一般是复数,所以混叠也是复数相加。为了简明起见,在图1-10中我们将Xa(j Ω)作为标量来处理。

我们将采样频率之半(Ωs/2)称为折叠频率,即 它如同一面镜子,当信号频谱超过它时,就会被折叠回来,造成频谱的混叠。

图3-3说明了在简单余弦信号情况下频谱混叠的情况。在图3-3(a )中,给出该余弦信号

的傅里叶变换Xa(j Ω)。

图(b )是在Ω0<Ωs/2时, 的傅里叶变换。图(c )是在Ω0>

Ωs/2时, 的傅里叶变换。(d )和(e )则分别对应于Ω0<Ωs/2=π/T 和Ω0>π/T 时低通滤波器输出的傅里叶变换,在没有混叠时((b)和(d)),恢复出的输出ya(t)为 在有混叠时,则是

s

s

s

s s s

(a )(b )

(c )

(d )s s s

…T

s π

=

Ω2t t x a 0cos )(Ω=)(?t x a )(?t x a t t y a 0cos )(Ω=t

t y s a )cos()(0

Ω-Ω=

这就是说,作为采样和恢复的结果,高频信号cos Ω0t 已经被当作和低频信号cos(Ωs-Ω0)t 是一样的东西被冒名顶替了。这个讨论就是奈奎斯特采样定理的基础。

图 3-3 一个余弦信号采样中的混叠效果

由此得出结论:要想采样后能够不失真地还原出原信号, 则采样频率必须大于两倍信号谱的最高频率(Ωs>2Ωh ),这就是奈奎斯特采样定理。 即

fs>2fh

频率Ωh 一般称为奈奎斯特频率,而频率2Ωh 称为奈奎斯特率。 采样频率必须大于奈奎斯特率。

在实际工作中,为了避免频谱混淆现象发生,采样频率总是选得比奈奎斯特频率更大些,例如选到(3~4)Ωh 。同时为了避免高于折叠频率的杂散频谱进入采样器造成频谱混淆,一般在采样器前加入一个保护性的前置低通滤波器,称为防混叠滤波器,其截止频率为Ωs/2,以便滤除掉高于Ωs/2 的频率分量。

00

s 2s

s

(c )(b )

(a )00

s 0(e )

(d )s 0

3.3 由样本重构带限信号

1 原理:如果理想采样满足奈奎斯特定理,即模拟信号谱的最高频率小于折叠频率

则采样后不会产生频谱混叠,由式(3-7)知

故将 通过一个理想低通滤波器,这个理想低通滤波器应该只让基带频谱通过,因而其带宽应该等于折叠频率,它的特性如图3-4所示。

图3-4

采样信号通过这个滤波器后,就可滤出原模拟信号的频谱

因此,在输出端可以得到原模拟信号 理想低通滤波器虽不可实现,但是在一定精度范围内,可用一个可实现的滤波器来逼近它。

2.由采样信号序列重构带限信号 理想低通滤波器的冲激响应为 ???

????Ω≥ΩΩ<ΩΩ=Ω2||02

||)()(s

s a a j X j X 2

||)(1)(?s

a

a j X T j X Ω<ΩΩ=Ω)(

?Ωj X a s

??????

?

Ω≥ΩΩ<Ω=Ω2||02

||)(s s T j H )()()(?)(Ω=ΩΩ=Ωj X j H j X j Y a a a )

()(t x t y a a =T

t T t t d e T d e j H t h s s t j t

j s

s

/)/sin(2

/)2/sin(2)(21)(2/2

/ππππ=

ΩΩ=Ω=ΩΩ=ΩΩΩ-Ω∞∞-??

由 与h(t)的卷积积分,即得理想低通滤波器的输出为

这里h(t-nT)称为内插函数:

它的波形如图3-5所示,其特点为:在采样点nT 上,函数值为1; 其余采样点上,函数值都为零。

图 3-5 内插函数

由于ya(t)=xa(t),因此以上卷积结果也可以表示为 3-8式

式(3-8)称为采样内插公式,即信号的采样值xa(nT)经此公式而得到连续信号xa(t)。 也就是说,xa(t)等于各xa(nT)乘上对应的内插函数的总和。在每一采样点上,只有该点所对应的内插函数不为零,这使得各采样点上信号值不变,而采样点之间的信号则由加权内插函数波形的延伸叠加而成,如图3-6所示。这个公式说明了,只要采样频率高于两倍信号最高频率,则整个连续信号就可以完全用它的采样值来代表,而不会丢掉任何信息。这就是奈奎斯特采样定理的意义。由上面讨论可看出采样内插公式只限于使用到限带(频带有限)信号上。 )(?t x a )()()()()()()()()()(?)(nT t h nT x d nT t h x d t h nT t x d t h x t y a n a n n a a a -=--=-??????-=-=∑?∑∑??∞

-∞

=∞

-∞-∞

=∞-∞=∞∞

-∞

∞-τ

τδττττδττττT nT t T nT t nT t h /)(]

/)(sin[)(--=

-ππ

T nT t T

nT t nT x t x n a a /)(]

/)(sin[)

()(--=∑∞

-∞=ππ

图 3-6 采样内插恢复

3.4 连续时间信号的离散时间处理 随着信号传输和处理手段的数字化发展,越来越有必要将连续信号转化为离散信号处理。

一、

C/D 转换 时域分析

∑+∞-∞

=-=n nT t t p )()(δ)(][nT x n x

c =)(t x c -∞=-=n nT

j c P e nT x j X ω)()(∑+∞-∞=-=k s c P k j X T j X ))((1)(ωωω

)

()(ωωj X T X p d =)

()()(T H j X T Y d p p ωωω=)

()()(ΩΩΩ=j d j d j d e H e X e Y )()()(T H T X T Y d d p ωωω=)(?)()(ωωωj H

j X j Y p c =

)

(][nT x n x c d =)(][nT y n y c d = +∞)(x t x c p = ??

?<>=2

/||)

()(2/||0s d P s T H j TX ωωωωωω??

?<>=2

/||)

()(2

/||0s d c s T H j X ωωωωωω

即: 例1:数字微分器 带限微分

3.6 利用离散时间信号处理改变采样频率 3.6.1 脉冲串采样

从时域角度看: )(ωj H c ??

?<>=2/||)(2

/||0s d s T H ωωωωω?][n x p π<ΩΩ=Ω||)()(T H H c

??

?<<=ωωωωωω||0||)(c c j j H π<ΩΩ=Ω||)()(T j e H j d

2

c 2||0||)(s

c c c j c e j H ωωωωωωωω=???><=?-θθωθω

d

e X e P e X j j j )()(1

)()(-?=-)s k ωω-)()

(k j s

e ωω∑+∞δT nT t T nT t n x t x n c /)(]

/)(sin[)

()(--=∑∞

-∞

=ππT nT t T nT t n y t y n c /)(]

/)(sin[)

()(--=

∑∞

-∞

=ππ

从频域角度看:

3.6.2 离散信号抽取与内插 抽取——从序列中提取每第N 个点上样本的过程。

2.内插

抽取又称为减抽样, 内插又称为增抽样。

减抽样使信号的频带扩展,但提高了数据的传输率。 增抽样虽降低了信息的传输率,但节省了传输频带。

对 以T 抽样再以N 抽取,则相当于对

以NT 为间隔来抽取。

例:某一离散时间序列 ,其傅立叶变换如图a 所示。现采用

抽取、内插等手段对信号进行处理,欲使处理后信号 的频谱

如图b 所示,请给出信号处理过程的系统框图,并画出各处相应的频谱图。

π

ωωω

<=)/(1)(T j Y T e Y j T

e TX j X T j c /)

()(π<Ω=ΩΩT j X j H j Y c c c /)()()(π<ΩΩΩ=Ω][][][][][nN x nN x n x n x n x p b p ==→→∑+∞-∞=-=n n

j b j b e n x e X ωω][)(∑+∞

-∞

=-=n n j p e nN x ω][)

(][)(N j

p m N m j p j b e X e m x e X ωωω==∑+∞-∞

=-nN

m =?????=其他的整倍数为0][][N n N n x n x d p )()(N

j d j p e X e X ωω

=)(t x c )(t x c )(n x )

(n

y

3.7 实际应用考虑的问题

3.7.1 消除混叠的预滤波

在采样前加一低通滤波器,以滤除高于2倍采样频率成分,以避免高频成分的干扰。

3.7.2 A/D转换中的量化误差

数字信号不仅在时间上是离散的,而且在取值上也不连续,即数字信号的取值必须为某个规定的最小数量单位的整数倍。

因此,为了将模拟信号转换成数字信号,还必须将采样/保持电路输出的采样值按照某种近似方式归并到相应的离散电平上,也就是将模拟信号在取值上离散化,我们把这个过程称为量化。将量化后的结果(离散电平)用数字代码来表示,称为编码。于单极性模拟信号,一般采用自然二进制编码;对于双极性模拟信号,则通常采用二进

制补码。经过编码后得到的代码就是A/D转换器输出的数字量。

在量化过程中,量化结果(离散电平)都是其中一个最小离散电平的整数倍,我们将这个最小离散电平值称为量化单位,记作Δ。它也

就是数字量的最低有效位LSB的1所代表的模拟电压值。

由于采样/保持电路输出的信号在取值上是连续的,不可能所有的采样值都恰好为量化单位的整数倍,所以量化前后不可避免地存在着误差,这种误差称为量化误差,用ε表示。

量化可以按两种近似方式进行:只舍不入量化方式和有舍有入(四舍五入)量化方式。以3位A/D转换器为例来说明这两种方式,设采样值的变化范围是0~8 V。只舍不入量化方式如图3-8所示。当采样/保持电路输出的电压uS介于两个量化电平之间时,采用取整的方法将其归并为较低的量化电平。

例如,无论uS = 5.9 V还是uS = 5.1 V,都将其归并为5 V的量化电平,输出的编码都为101。可见,采用只舍不入量化方式,最大误差εmax等于量化单位Δ。

四舍五入量化方式如图3-9所示。当采样/保持电路输出的电压uS介于两个量化电平之间时,采用四舍五入的方式将其归并为最相近那个量化电平。例如,若uS = 5.49 V,就将其归并为5 V的量化电平,输出的编码为101;若uS = 5.50 V,就将其归并为6 V的量化电平,

输出的编码为110。可见,采用四舍五入量化方式,最大量化误差εmax只有量化单位的一半(Δ/2),比只舍不入量化方式的最大

转换器都采用这种量化方式。

图3-9 四舍五入量化方式

量化误差是A/D转换的固有误差,只能减小,不能完全消除。减小量化误差的主要措施就是减小量化单位。但是当输入模拟电压的变化范围一定时,量化单位越小就意味着量化电平的个数越多,编码的位数越大,电路也就越复杂。

另外:

D/A转换的误差以及与A/D转换和D/A转换中抽取与内插的应用等也要引起注意.这里从略

本章小结:

本章建立并研究了连续信号现经周期采样而得到的离散时间序列之间的关系.

奈奎斯特定理表明对一个限带信号而言,只要采样的频率相对于连续信号中的最高频率的来说足够高的话,则周期采样信号可以无失真的恢复原信号.

相应于采样频率的变化将使信号在频域频谱重复并重新归一化.

最后研究了在实际处理中要注意的几个实际问题.

000

001

010

011

100

101

110

111

编码

S

000

001

010

011

100

101

110

111

编码

S

利用MATLAB实现连续信号的采样与重构仿真课程设计

目录 1、摘要 (1) 2、正文 (2) 2.1、设计目的 (2) 2.2、设计原理 (2) (1)、MTLAB简介 (2) (2)、连续时间信号 (2) (3)、采样定理 (3) (4)、信号重构 (5) 2.3、信号采样和恢复的程序 (5) (1)设计连续信号 (6) (2)设计连续信号的频谱 (7) (3)设计采样信号 ........................................错误!未定义书签。 (4)设计采样信号的频谱图 (9) (5)设计低通滤波器 (10) (6)恢复原信号 (12) 3、总结和致谢........................... 错误!未定义书签。

4、参考文献 (15) 1.摘要 本次课程设计使用MATLAB实现连续信号的采样和重构仿真,了解MATLAB软件,学习使用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用,以及一些关键命令的掌握,理解,分析等。初步掌握线性系统的设计方法,培养独立工作能力。 加深理解采样和重构的概念,掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采用和重构的方法。计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差,并由此总结采样频率对信号重构误差的影响。 要做到以下基本要求: 1. 掌握利用MATLAB分析系统频率响应的方法,增加对仿真软件MATLAB的感性认识,学会该软件的操作和使用方法。 2. 掌握利用MATLAB实现连续信号采用和重构的方法,加深理解采样和重构的概念。 3 . 初步掌握线性系统的设计方法,培养独立工作能力。 4. 学习MATLAB中信号表示的基本方法及绘图函数的调用,实现对常用连续时间信号的可视化表示,加深对各种电信号的理解。 5. 加深理解采样对信号的时域和频域特性的影响;验证信号和系统的基本概念、基本理论,掌握信号和系统的分析方法。 6. 加深对采样定理的理解和掌握,以及对信号恢复的必要性;掌握对连续信号在时域的采样和重构的方法。

实验连续时间信号与系统的基本表示与分析

实验1 连续时间信号与系统的基本表示与 分析 实验目的: 1. 了解并掌握用Matlab 分析信号与系统时用到的主要基本知识; 2. 掌握用matlab 表示常见的连续信号,并进行信号的基本运算 3. 熟悉使用Matlab 分析系统的方法。 实验内容: 1. 画出下列连续信号的波形 (1)()() 2()t x t e u t -=- (2)[](1) ()(1)(2)t x t e u t u t --=--- (3)()cos(10)()t x t e t u t π-= 2. 判断系统是否为线性系统,是否具有时不变特性。 (1)()(2)y t x t = (2)()2 ()y t x t = 指导资料: 1.1MATLAB 的基础知识 1.1.1数据的表示和运算 MATLAB 中,变量可以通过变大时直接赋值,例如直接输入 >>a=2+2 得到的结果为 a=4 如果输入的表达式后面加上分号“;”,那么结果就不会显示出来。由于MA TLAB 的变量名对字母大小写敏感,因此“a ”和“A ”是两个不同的变量名。 MATLAB 主要用到以下数值运算符: + 加 - 减 * 乘 / 除 ^ 乘方(幂)

' (矩阵)转置 这些符号可以对数值或已经定义过的变量进行运算,并给变量直接赋值。例如,假设变量“a”在上面已经定义过,则 >>b=2*a 得到的结果为 b=8 MATLAB中有一些预定义的变量可以直接使用。信号与系统中常用的变量有: i或j 1 pi π(3.1415926……) 在信号与系统中,常用以下函数进行计算和对变量的赋值: abs 数值的大小(实数的绝对值) angel 复数的角度,以弧度表示 real 求复数的实部 imag 求复数的虚部 cos 余弦函数,假设角度是弧度值 sin 正弦函数,假设角度是弧度制 exp 指数海曙 sqrt 求平方根 例如: >>y=2*(1+4*i) y=2.000+8.000i >>c=abs(y) c=8.2462 >>d=angle(y) d=1.3258 1.1.2矩阵的表示和运算 MATLAB是基于矩阵和向量的代数运算,甚至标量也可以看做是1×1 的矩阵,因此,MATLAB中对矩阵和向量的操作比较简单。 向量可以用两种方法定义。第一种是指定元素建立向量: v = [1 3 5 7]; 这个命令创建了一个1×4 的行向量,元素为1,3,5和7.可以用逗号代替空格来分隔元素: v = [1,3,5,7]; 如果要增加向量的元素,可以表示为 v(5) = 8 得到的向量为v = [1 3 5 7 8]。前面定义过的向量还可以用来定义新的向量,例如,前面已经定义过的向量v ,再定义向量a和b: a = [9 10]; b = [ v a]; 得到向量b为

连续时间信号傅里叶级数分析及MAtlAB实现

课程设计任务书 学生姓名:专业班级: 指导教师:工作单位: 题目: 连续时间信号傅里叶级数分析及MATLAB实现 初始条件: MATLAB 6.5 要求完成的主要任务: 深入研究连续时间信号傅里叶级数分析的理论知识,利用MA TLAB强大的图形处理功能,符号运算功能以及数值计算功能,实现连续时间周期信号频域分析的仿真波形。 1.用MATLAB实现周期信号的傅里叶级数分解与综合。 2.用MATLAB实现周期信号的单边频谱及双边频谱。 3.用MATLAB实现典型周期信号的频谱。 4.撰写《MATLAB应用实践》课程设计说明书。 时间安排: 学习MATLAB语言的概况第1天 学习MATLAB语言的基本知识第2、3天 学习MATLAB语言的应用环境,调试命令,绘图能力第4、5天 课程设计第6-9天 答辩第10天 指导教师签名:年月日 系主任(或责任教师)签名:年月日

目录 摘要................................................................................................................................................ I Abstract .......................................................................................................................................... II 绪论. (1) 1 MATLAB简介 (2) 1.1 MATLAB语言功能 (2) 1.2 MATLAB语言特点 (2) 2 傅里叶级数基本原理概要 (4) 2.1 周期信号的傅里叶分解 (4) 2.2 三角形式和指数形式傅里叶级数及各系数间的关系 (4) 2.3 周期信号的频谱 (5) 3 用MATLAB实现周期信号的傅立叶级数分解与综合 (6) 3.1 合成波形与原波形之间的关系 (6) 3.2 吉布斯现象 (6) 4 用MATLAB实现周期信号的单边频谱及双边频谱。 (8) 4.1 单边,双边(幅度,相位)频谱及其关系 (8) 4.1.1单边,双边(幅度,相位) (8) 4.1.2 单边,双边频谱关系 (9) 4.2以单边幅度频谱为例,研究脉冲宽度与频谱的关系 (10) 4.3以单边幅度频谱为例,研究脉冲周期与频谱的关系 (11) 5用MATLAB实现典型周期信号的频谱 (13) 5.1 周期方波脉冲频谱的MATLAB实现 (13) 5.2 周期三角波脉冲频谱的MATLAB 实现 (14) 6 小结及心得体会 (17) 参考文献 (18) 附录: (19)

信号与系统 连续时间信号卷积运算

连续时间信号的卷积运算的MATILAB实现 薛皓20091453 例1:已知两连续时间信号如图9-3所示,试用matlab求f(t)=f1(t)*f2(t),并绘出f(t)的时域波形图。 图1-1 连续时间信号波形图示例 实现上述过程的matlab命令如下: p=0.5; k1=0:p:2; f1=0.5*k1; k2=k1; f2=f1; [f,k]=sconv(f1,f2,k1,k2,p) 上述命令绘制的波形图也在图9-3中示出。图9-3中给出了抽样时间间隔p=0.5时的处理效果。而图9-4给出了抽样时间间隔p=0.01时的处理效果。

图1-2 例1的连续时间信号波形图 习题1:已知f1(t)=1(2t 1≤≤),f2(t)=1(3t 2≤≤),用matlab 实现其卷积并绘制出卷积曲线。 解:程序代码如下: >> p=0.01; k1=1:p:2; f1=ones(size(k1)).*(k1>1); k2=2:p:3; f2=ones(size(k2)).*(k2>2); f=conv(f1,f2); f=f*p; k0=k1(1)+k2(1); k3=k1(length(k1))+k2(length(k2)); subplot(2,2,1) plot(k1,f1) title('f1(t)') xlabel('t') ylabel('f1(t)') subplot(2,2,2) plot(k2,f2)

title('f2(t)') xlabel('t') ylabel('f2(t)') subplot(2,2,3) plot(k,f); h=get(gca,'position'); h(3)=2.5*h(3); 0 set(gca,'position',h) title('f(t)=f1(t)*f2(t)') xlabel('t') ylabel('f(t)') 绘制图形如图2-1所示。 图2-1 习题2:)1()2/1t ()t (2f ),1t ()t ()t (1f δ-+δ=-ε-ε=,求其卷积。 程序代码: p=0.01; t1=0:p:1; f1=ones(size(t1)).*(t1>0); t2=-0.5:p:1; f2=(t2==-0.5)-(t2==1); f=conv(f1,f2); f=f*p; t=-0.5:p:2;

(完整word版)连续时间信号分析答案

实验一 连续时间信号分析 一、实验目的 (一)掌握使用Matlab 表示连续时间信号 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉常用信号的波形和特性 (二)掌握使用Matlab 进行连续时间信号的相关运算 1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换 2、学会运用Matlab 进行连续时间信号微分、积分运算 3、学会运用Matlab 进行连续时间信号相加、相乘运算 4、学会运用Matlab 进行连续时间信号卷积运算 二、实验条件 一台电脑、winXP 系统、matlab7.0软件 三、实验内容 1、利用Matlab 命令画出下列连续信号的波形图。 (1))4/3t (2cos π+ 代码: clear all;close all;clc; K=2;a=3; t=0:0.01:3; ft=K*cos(a*t+pi/4); plot(t,ft),grid on axis([-5,5,-2.2,2.2]) title('2cos(3t+4π)')

-5 -4 -3 -2 -1 1 2 3 4 5 -2-1.5-1-0.500.511.5 22cos(3t+4π) (2) )t (u )e 2(t -- -3 -2-10123 -3 -2 -1 1 2 3 指数信号与阶跃信号的乘积

代码: 函数文件: function f=uCT(t) f=(t>=0); 命令文件: clear all;close all;clc; a=-1; t=-5:0.01:5; ft=(2-exp(a*t)).*uCT(t); %y=2-exp(a*t); %plot(t,y),grid on plot(t,ft),grid on axis([-3,3,-3,3]); title('指数信号与阶跃信号的乘积') (3))]2()(u )][t (cos 1[--+t u t π

信号与系统课设 常用连续时间信号的可视化及微积分运算

成绩评定表

课程设计任务书

目录 一、引言 (1) 二、Matlab入门 (2) 2.1 Matlab7.0介绍 (2) 2.2利用Matlab7.0编程完成习题设计 (3) 三、Matlab7.0实现连续时间信号微积分运算的设计 (4) 3.1常用连续时间信号的类别及原理 (4) 3.2编程设计及实现 (4) 3.3运行结果及其分析 (7) 四、结论 (16) 五、参考文献 (17)

一、引言 人们之间的交流是通过消息的传播来实现的,信号则是消息的表现形式,消息是信号的具体内容。 《信号与系统》课程是一门实用性较强、涉及面较广的专业基础课,该课程是将学生从电路分析的知识领域引入信号处理与传输领域的关键性课程,对后续专业课起着承上启下的作用. 该课的基本方法和理论大量应用于计算机信息处理的各个领域,特别是通信、数字语音处理、数字图像处理、数字信号分析等领域,应用更为广泛。 近年来,计算机多媒体教序手段的运用逐步普及,大量优秀的科学计算和系统仿真软件不断涌现,为我们实现计算机辅助教学和学生上机实验提供了很好的平台。通过对这些软件的分析和对比,我们选择MATLAB语言作为辅助教学工具,借助MATLAB强大的计算能力和图形表现能力,将《信号与系统》中的概念、方法和相应的结果,以图形的形式直观地展现给我们,大大的方便我们迅速掌握和理解老师上课教的有关信号与系统的知识。 MATLAB是当前最优秀的科学计算软件之一,也是许多科学领域中分析、应用和开发的基本工具。MATLAB全称是Matrix Laboratory,是由美国Mathworks公司于20世纪80年代推出的数学软件,最初它是一种专门用于矩阵运算的软件,经过多年的发展,MATLAB已经发展成为一种功能全面的软件,几乎可以解决科学计算中的所有问题。而且MATLAB编写简单、代码效率高等优点使得MATLAB在通信、信号处理、金融计算等领域都已经被广泛应用。它具有强大的矩阵计算能力和良好的图形可视化功能,为用户提供了非常直观和简洁的程序开发环境,因此被称为第四代计算机语言。MATLAB 强大的图形处理功能及符号运算功能,为我们实现信号的可视化及系统分析提供了强有力的工具。MATLAB 强大的工具箱函数可以分析连续信号、连续系统,同样也可以分析离散信号、离散系统,并可以对信号进行各种分析域计算,如相加、相乘、移位、反折、傅里叶变换、拉氏变换、Z 变换等等多种计算。 作为信号与系统的基本分析软件之一,利用MATLAB进行信号与系统的分析与设计是通信以及信息工程学科的学生所要掌握的必要技能之一。通过学习并使用MATLAB语言进行编程实现课题的要求,对学生能力的培养极为重要。尤其会提高综合运用所学理

连续时间信号的抽样及频谱分析-时域抽样信号的频谱--信号与系统课设

1 引言 随着科学技术的迅猛发展,电子设备和技术向集成化、数字化和高速化方向发展,而在学校特别是大学中,要想紧跟技术的发展,就要不断更新教学和实验设备。传统仪器下的高校实验教学,已严重滞后于信息时代和工程实际的需要。仪器设备很大部分陈 旧,而先进的数字仪器(如数字存储示波器)价格昂贵不可能大量采购,同时其功能较为单一,与此相对应的是大学学科分类越来越细,每一专业都需要专用的测量仪器,因此仪器设备不能实现资源共享,造成了浪费。虚拟仪器正是解决这一矛盾的最佳方案。基于PC 平台的虚拟仪器,可以充分利用学校的微机资源,完成多种仪器功能,可以组合成功能强大的专用测试系统,还可以通过软件进行升级。在通用计算机平台上,根据测试任务的需要来定义和设计仪器的测试功能,充分利用计算机来实现和扩展传统仪器功能,开发结构简单、操作方便、费用低的虚拟实验仪器,包括数字示波器、频谱分析仪、函数发生器等,既可以减少实验设备资金的投入,又为学生做创新性实验、掌握现代仪器技术提供了条件。 信号的时域分析主要是测量测试信号经滤波处理后的特征值,这些特征值以一个数值表示信号的某些时域特征,是对测试信号最简单直观的时域描述。将测试信号采集到计算机后,在测试VI 中进行信号特征值处理,并在测试VI 前面板上直观地表示出信号的特征值,可以给测试VI 的使用者提供一个了解测试信号变化的快速途径。信号的特征值分为幅值特征值、时间特征值和相位特征值。 尽管测量时采集到的信号是一个时域波形,但是由于时域分析工具较少,所以往往把问题转换到频域来处理。信号的频域分析就是根据信号的频域描述来估计和分析信号的组成和特征量。频域分析包括频谱分析、功率谱分析、相干函数分析以及频率响应函数分析。 信号在时域被抽样后,他的频谱X(j )是连续信号频谱X(j )的形状以抽样频率为间隔周期重复而得到,在重复过程中幅度被p(t)的傅里叶级数Pn加权。因为Pn只是n的函数,所以X(j )在重复的过程中不会使其形状发生变化。假定信号x(t)的频谱限制在- m~+ m的范围内, 若以间隔Ts对xa(t)进行抽样,可知抽样信号X^(t)的频谱X^(j )是以s为周期重复。显然,若在抽样的过程中s<2 m,则X^(j )将发生频谱混叠现象,只有在抽样的过程中满足s>=2 m条件,X^(j )才不会产生频谱的混叠,接收端完全可以由x^(t)恢复原连续信号xa(t),这就是低通信号抽样定理的核心内容。

连续时间信号的采样实验

实验一 连续时间信号的采样 一、 实验目的 进一步加深对采样定理和连续信号傅立叶变换的理解。 二、实验步骤 1.复习采样定理和采样信号的频谱 采样定理 如果采样频率s F 大于有限带宽信号)(t x a 带宽0F 的两倍,即 02F F s > (1) 则该信号可以由它的采样值)()(s a nT x n x =重构。否则就会在)(n x 中产生混叠。该有限带宽模拟信号的02F 被称为乃魁斯特频率。 必须注意,在)(t x a 被采样以后,)(n x 表示的最高模拟频率为2/s F Hz (或πω=)。 2.熟悉如何用MATLAB 语言实现模拟信号表示 严格地说,除了用符号处理工具箱(Symbolics)外,不可能用MATLAB 来分析模拟信号。然而如果用时间增量足够小的很密的网格对)(t x a 采样,就可得到一根平滑的曲线和足够长的最大时间来显示所有的模态。这样就可以进行近似分析。令t ?是栅网的间隔且s T t <

连续时间信号和系统时域分析及MATLAB实现[精编版]

连续时间信号和系统时域分析及MATLAB实 现[精编版] MATLAB课程设计任务书 姓名:王** 学号:2010******010 题目:

连续时间信号和系统时域分析及MATLAB实现 初始条件: MATLAB 7.5.0 ,Windows XP系统 实验任务: 一、用MATLAB实现常用连续时间信号的时域波形(通过改变参数,分析其时域特性)。 1、单位阶跃信号, 2、单位冲激信号, 3、正弦信号, 4、实指数信号, 5、虚指数信号, 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加, 2、相乘, 3、数乘, 4、微分, 5、积分 三、用MATLAB实现信号的时域变换(参数变化,分析波形变化) 1、反转, 2、使移(超时,延时), 3、展缩, 4、倒相, 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解, 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子,对每个例子,要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子,四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 目录

1.1 MATLAB设计目的 (1) 1.2 MATLAB语言特点 (1) 2常用连续时间信号的时域波形 (1) 2.1单位阶跃信号 (1) 2.2单位冲激信号 (2) 2.3正弦信号 (3) 2.4实指数信号 (4) 2.5虚指数信号 (5) 2.6复指数信号 (6) 3 连续时间信号的时域运算 (7) 3.1相加 (7) 3.2相乘 (8) 3.3数乘 (9) 3.4微分 (10) 3.5积分 (11) 4.1反转 (12) 4.2时移 (13) 4.3展缩 (14) 4.4倒相 (15) 4.5综合变化 (16) 5连续时间信号简单的时域分解 (17) 5.1信号的交直流分解 (17)

实验二_连续时间信号的频域分析

实验二 连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs 现象”,了解其特点以及产生的原因; 3、掌握连续时间傅里叶变换的分析方法及其物理意义; 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质; 5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序,并能利用这些程序对一些典型信号进行频谱分析,验证CTFT 、DTFT 的若干重要性质。 基本要求:掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用Matlab 编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号,只要满足狄利克利条件,就可以展开成傅里叶级数。 三角傅里叶级数为: ∑∞ =++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 2.1 或: ∑∞ =++ =1 00)cos()(k k k t k c a t x ?ω 2.2 其中1 02T π ω= ,称为信号的基本频率(Fundamental frequency ),k k b a a ,和,0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度,k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位,它们都是频率0ωk 的函数,绘制出它们与0ωk 之间的图像,称为信号的频谱图(简称“频谱”),k c -0ωk 图像为幅度谱,k ?-0ωk 图像为相位谱。 三角形式傅里叶级数表明,如果一个周期信号x(t),满足狄里克利条件,就可以被看作是由很多不同频率的互为谐波关系(harmonically related )的正弦信号所组成,其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component),其幅度(amplitude )为k c 。

连续时间信号的频域分析(信号与系统课设).

福建农林大学计算机与信息学院 信息工程类 课程设计报告 课程名称:信号与系统 课程设计题目:连续时间信号的频域分析 姓名: 系:电子信息工程 专业:电子信息工程 年级:2008 学号: 指导教师: 职称: 2011 年 1 月10 日

福建农林大学计算机与信息学院信息工程类 课程设计结果评定

目录 1课程设计的目的 (1) 2课程设计的要求 (1) 3课程设计报告内容.....................................................................1-13 3.1连续信号的设计..................................................................1-11 3.2验证傅里叶变换的调制定理 (11) 3.3周期信号及其频谱 (12) 4总结 (13) 参考文献 (14)

连续时间信号的频域分析 1.课程设计的目的 (1)熟悉MATLAB语言的编程方法及MATLAB指令; (2)掌握连续时间信号的基本概念; (3)掌握门函数、指数信号和抽样信号的表达式和波形; (4)掌握连续时间信号的傅里叶变换及其性质; (5)掌握连续时间信号频谱的概念以及幅度谱、相位谱的表示; (6)掌握利用MATLAB进行信号的傅里叶变换以及时域波形和频谱的表示;(7)通过连续时间信号的频域分析,更深刻地理解了连续时间信号的时域和频域间的关系,加深了对连续时间信号的理解。 2.课程设计的要求 (1)自行设计以下连续信号:门函数、指数信号和抽样信号。要求:(a)画出以上信号的时域波形图; (b)实现以上信号的傅里叶变换,画出以上信号的幅度谱及相位谱,并对相关结果予以理论分析; (c)对其中一个信号进行时移和尺度变换,分别求变换后信号的傅里叶变换,验证傅里叶变换的时移和尺度变换性质。 (2)自行设计信号,验证傅里叶变换的调制定理。 (3)自行设计一个周期信号,绘出该信号的频谱,并观察周期信号频谱的特点。 3.课程设计报告内容 3.1(a)①门函数(矩形脉冲): MATLAB中矩形脉冲信号用rectpuls函数表示: y=rectpuls (t,width) %width缺省值为1 >> t=-2:0.001:2; T=2; yt=rectpuls (t,T); plot(t,yt); axis([-2,2,0,1.5]); grid on; %显示格线

离散系统与连续时间系统的根本差别是:离散系统(图)有采样开

离散系统与连续时间系统的根本差别是:离散系统(图3)有采样开关存在,而连续系统则无。连续信号经过采样开关变成离散信号(图4),采样开关起这理想脉冲发生器的作用,通过它将连续信号调制成脉冲序列。 图3 离散系统方块图 图4 离散型时间函数 调制之后的信号中,包含与脉冲频率相关的高频频谱(图5),相邻两频谱不相重叠的条件是: max 2f f s 其中: s f ---采样开关的采样频率 m ax f ---连续信号频谱中的最高频率 这就是采样定理,通常选择采样频率时取四倍连续信号的最大频率。实验中,信号源产生频率可调的周期性信号,计算机通过A/D 板将信号采集入内存,通过软件示波器显示出来,调整采样频率,可以得到不同的采样结果,以波形图直观显示出来。由此,可考察波形失真程度。 三、实验使用的仪器设备及实验装置 1. 装有LabVIEW 软件和PCI-1200数据采集卡的计算机一台 2. 频率计或信号发生器一台 3. 外接端子板、数据采集板、计算机、组态软件 基于LabVIEW 的信号测试系统主要包括信号发生器、DAQ 数据采集卡和计算机软件三部分组成。A/D 数据采集采用NI 公司PCMCIA 接口的PCI-1200型多功能数据采集卡;L abVIEW 7.1软件。 将PCI-1200数据采集卡插到计算机主板上的一个空闲的PCI 插槽中,接好各种附件,其驱动程序就是NI-DAQ 。附件包括一条50芯的数据线,一个型号为CB-50LP 的转接板,转接板直接与外部信号连接。 图5 信号频谱图

四、具体实验步骤 (一)通过LabVIEW进行模拟信号的数据采集 1. 安装数据采集卡,根据数据采集卡接线指示(图6)连接线路,并检查测试。 2. 熟悉LabVIEW软件中与数据采集相关的控件与设置项。 3. 编制DAQ程序,并调试数据采集组态。 4. 应用该组态软件进行波形数据采集并存储,信号种类设置为正弦波,分别设置 信号发生器频率为50,100Hz,观察并记录波形变化。 5. 设置信号种类为方波或锯齿波,重复上述实验。 (二)采样定理验证实验 1. 按图8连接线路,并检查测试。 2. 熟悉GeniDAQ软件中与数据采集相关的控件与设置项。 3. 编制、调试数据采集组态。 4. 应用该组态软件进行波形数据采集并存储,信号种类设置为正弦波,分别设置 信号发生器频率为50,100Hz,采集频率设置为50、100、150、200、300、500Hz,观察并记录波形变化,体验采样定理的正确性。 五、实验准备及预习要求 1.认真阅读实验指导书,在老师答疑和同学讨论的基础上,完成实验准备任务: 1).了解数据采集及其硬件(A/D变换器和数据采集卡)选择的基本知识; 2).熟悉G语言编程环境和虚拟仪器的含义; 1.理解采样定理的意义;

连续时间信号的抽样及频谱分析-时域抽样信号的频谱__信号与系统课设

连续时间信号的抽样及频谱分析-时域抽样信号的频谱__信号与系统课设

1 引言 随着科学技术的迅猛发展,电子设备和技术向集成化、数字化和高速化方向发展,而在学校特别是大学中,要想紧跟技术的发展,就要不断更新教学和实验设备。传统仪器下的高校实验教学,已严重滞后于信息时代和工程实际的需要。仪器设备很大部分陈旧,而先进的数字仪器(如数字存储示波器)价格昂贵不可能大量采购,同时其功能较为单一,与此相对应的是大学学科分类越来越细,每一专业都需要专用的测量仪器,因此仪器设备不能实现资源共享,造成了浪费。虚拟仪器正是解决这一矛盾的最佳方案。基于PC 平台的虚拟仪器,可以充分利用学校的微机资源,完成多种仪器功能,可以组合成功能强大的专用测试系统,还可以通过软件进行升级。在通用计算机平台上,根据测试任务的需要来定义和设计仪器的测试功能,充分利用计算机来实现和扩展传统仪器功能,开发结构简单、操作方便、费用低的虚拟实验仪器,包括数字示波器、频谱分析仪、函数发生器等,既可以减少实验设备资金的投入,又为学生做创新性实验、掌握现代仪器技术提供了条件。

信号的时域分析主要是测量测试信号经滤波处理后的特征值,这些特征值以一个数值表示信号的某些时域特征,是对测试信号最简单直观的时域描述。将测试信号采集到计算机后,在测试VI 中进行信号特征值处理,并在测试VI 前面板上直观地表示出信号的特征值,可以给测试VI 的使用者提供一个了解测试信号变化的快速途径。信号的特征值分为幅值特征值、时间特征值和相位特征值。 尽管测量时采集到的信号是一个时域波形,但是由于时域分析工具较少,所以往往把问题转换到频域来处理。信号的频域分析就是根据信号的频域描述来估计和分析信号的组成和特征量。频域分析包括频谱分析、功率谱分析、相干函数分析以及频率响应函数分析。 信号在时域被抽样后,他的频谱X(j )是连续信号频谱X(j )的形状以抽样频率为间隔周期重复而得到,在重复过程中幅度被p(t)的傅里叶级数Pn加权。因为Pn只是n的函数,所以X(j )在重复的过程中不会使其形状发生变化。假定信号x(t)的频谱限制在- m~+ m

应用_MATLAB实现连续信号的采样与重构

抽样定理及应用 2.1课程设计的原理 2.1.1连续信号的采样定理 模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样,信号采样后其频谱产生了周期延拓,每隔一个采样频率 fs ,重复出现一次。为保证采样后信号的频谱形状不失真,采样频率必须大于信号中最高频率成分的两倍,这称之为采样定理。时域采样定理从采样信号 恢复原信号 必需满足两个条件: (1) 必须是带限信号,其频谱函数在 > 各处为零;(对信号的要求, 即只有带限信号才能适用采样定理。) (2) 取样频率不能过低,必须 >2 (或 >2)。(对取样频率的要 求,即取样频率要足够大,采得的样值要足够多,才能恢复原信号。)如果采样频率 大于或等于 ,即 ( 为连续信号 的有限频谱),则采样离散信号能无失真地恢复到原来的连续信号 。一个频 谱在区间(- , )以外为零的频带有限信号,可唯一地由其在均匀 间隔 ( < )上的样点值 所确定。根据时域与频域的对称性, 可以由时域采样定理直接推出频域采样定理。一个时间受限信号()t f ,它集中在(m m ωω+-,)的时间范围内,则该信号的频谱()ωj F 在频域中以间隔为1ω的冲激序列进行采样,采样后的频谱)(1ωj F 可以惟一表示原信号的条件为重复周期 m t T 21≥,或频域间隔m t f 21 21≤ = πω(其中112T πω=)。采样信号 的频谱是原 信号频谱 的周期性重复,它每隔 重复出现一次。当s ω>2 时, 不会出现混叠现象,

原信号的频谱的形状不会发生变化,从而能从采样信号中恢复原信号。 >2的含义是:采样频率大于等于信号最高频率的2倍;这里的“不(注: s 混叠”意味着信号频谱没有被破坏,也就为后面恢复原信号提供了可能!) (a) (b) (c) 图* 抽样定理 a)等抽样频率时的抽样信号及频谱(不混叠) b)高抽样频率时的抽样信号及频谱(不混叠) c) 低抽样频率时的抽样信号及频谱(混叠) 2.1.2信号采样 如图1所示,给出了信号采样原理图

5.连续时间信号的抽样与量化资料

第5章 连续时间信号的抽样与量化 5.1 学习要求 1. 掌握时域抽样过程及时域抽样定理,会求已知信号的奈奎斯特频率;了解抽样信号的频 谱及其求解方法。 2. 掌握抗混叠滤波处理 3. 深刻理解连续时间信号的内插恢复过程; 4. 理解频域抽样定理; 5. 了解连续时间信号的离散处理过程。 5.2 学习重点 1. 时域抽样定理。 5.3知识结构 5.4 内容摘要 5.4.1 时域抽样定理 1. 时域抽样 就是利用抽样脉冲序列)(t p 从时域连续信号)(t f 中抽取一系列的离散样值,这种离散信号通常称为抽样信号,以)(t f s 表示 ,抽样信号傅里叶变换为: ()s n n FT s n F P t p t f t f ωω-? =∑∞ -∞ =)()()( ()()dt e t p T n P t jn T T s s s s ω--? = 22 1 ,称为)(t p 的傅里叶级数的系数。 n P 取决于抽样脉冲序列的形状,可以是,也可以是矩形脉冲抽样。 (1) 冲激抽样 设单位冲激序列)(t T δ为: ∑∞ -∞ =-= n s T nT t t )()(δδ ,

()()dt e t T n P t jn T T T s s s s ωδ--? = 221 = s T 1 抽样信号为:()()()()()s T s s n f t f t t f nT t nT δδ∞ =-∞ =?=?-∑ 则抽样信号)(t f s 的频谱为:∑∞ -∞ =-= n s s s n F T F )(1 )(ωωω (2) 矩形脉冲序列的抽样 如果抽样脉冲序列是周期为s T ,幅度为1,宽度为τ的矩形脉冲序列)(t p , 则抽样信号)(t f s 的频谱为: )()2 ( )](*)([21)(s s n s s n F n Sa T p F F ωωτ ωτ ωωπω-= =∑∞ -∞ = 2. 时域抽样定理 时域抽样定理是指一个频谱受限的信号)(t f ,如果频谱只占据m ω-到m ω的范围,则信号)(t f 可以用等间隔的抽样值唯一地表示,而抽样间隔m s f T 21 ≤(其中m m f πω2=),或者说,最低抽样频率为m f 2。 抽样的最大间隔m s f T 21=称为奈奎斯特间隔,最低抽样速率m s f f 2=称为奈奎斯特速率,最低抽样频率m s ωω2=称为奈奎斯特频率。 5.4.2 频率混叠效应和信号抽样频率的选择 由于信号在时域上的抽样而造成信号在频域上的频谱混叠称作频率混叠效应。减小频率混叠效应有两种途径:一是提高信号的抽样频率s ω,即缩小抽样周期T 。二是对被抽样的信号)(t f 预先进行抗混叠滤波处理,将非带限信号变成带限信号,然后按抽样定理抽样。后一种方法虽然使信号丢失了部分高频分量,但可以有效地保护信号()x t 中低频分量将不因抽样而受到干扰。 当信号有效带宽a f 已知时,若取抗混叠滤波器截止频率c a f f ≈,当滤波器具有-50~-60dB/倍频程衰减率,那么滤波后的信号以3s a f f =抽样即可。

实验二 连续时间信号的频域分析

实验二连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因; 3、掌握连续时间傅里叶变换的分析方法及其物理意义; 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质; 5、学习掌握利用Matlab语言编写计算CTFS、CTFT和DTFT的仿真程序,并能利用这些程序对一些典型信号进行频谱分析,验证CTFT、DTFT的若干重要性质。 基本要求:掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用Matlab编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS分析 任何一个周期为T 1 的正弦周期信号,只要满足狄利克利条件,就可以展开成傅里叶级数。 三角傅里叶级数为: ∑∞=+ + = 1 0 0 )] sin( ) cos( [ )( k k k t k b t k a a t xω ω 2.1

或: ∑∞ =++=100)cos()(k k k t k c a t x ?ω 2.2 其中1 02T πω=,称为信号的基本频率(Fundamental frequency ),k k b a a ,和,0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度,k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位,它们都是频率0ωk 的函数,绘制出它们与0ωk 之间的图像,称为信号的频谱图(简称“频谱”),k c -0ωk 图像为幅度谱,k ?-0ωk 图像为相位谱。 三角形式傅里叶级数表明,如果一个周期信号x(t),满足狄里克利条件,就可以被看作是由很多不同频率的互为谐波关系(harmonically related )的正弦信号所组成,其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component),其幅度(amplitude )为k c 。也可以反过来理解三角傅里叶级数:用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。 指数形式的傅里叶级数为: ∑∞-∞== k t jk k e a t x 0)(ω 2.3 其中,k a 为指数形式的傅里叶级数的系数,按如下公式计算: ? --=2/2/1110)(1T T t jk k dt e t x T a ω 2.4 指数形式的傅里叶级数告诉我们,如果一个周期信号x(t),满足狄里克利条件,那么,它就可以被看作是由很多不同频率的互为谐波关系(harmonically related )的周期复指数信号所组成,其中每一个不同频率的周期复指数信号称为基本频率分量,其复幅度

信号与系统 连续时间信号的采样 实验

实验一连续时间信号的采样 一.简述实验原理的目的; 实验原理: 采样定理:如果采样频率大于有限带宽信号带宽的两倍,即,则该信号可以由它的采样值重构。否则就会在中产生混叠。该有限带宽模拟信号的被称为奈奎斯特频率。 实验目的: 用不同的采样频率采样并进行比较,验证采样定理。 对不同信号采样验证时域与频域的性质。 二.结合实验中得到的实验结果曲线与理论结果比较,并分析说明误差产生的原因; matlab是用网栅间隔来表示模拟信号的。其实现积分是离散的点的求和,与理论值存在误差。 题目采用数值方法估计,是把用一个栅格序列来近似,由于,题中用之间的有限长度信号来近似表示。 三.总结实验所得主要结论 采样频率必须大于2即奈奎斯特采样频率(为有限带宽信号的带宽),该信号才可以由它的采样值重构。 若小于2则会在中产生混叠,为欠采样。产生的傅立叶变换图形是把互相交叠的的复制品叠加的结果,所以产生了失真。 四.简要回答思考题。 1.通过实验说明信号的时域与频域成反比的关系。 从实验一做出的图形看出: 的时域是时域的1/2倍,而的频域却是的频域的2倍。 2.分别求出奈奎斯特采样间隔,并与例一的信号的奈奎斯特采样间隔比较。 的奈奎斯特采样频率为8000样本/秒;则奈奎斯特采样间隔为1/8000;

是例一的信号的1/2;的奈奎斯特采样频率为2000样本/秒;则奈奎斯特采样间隔为1/2000,是例一的信号的2倍。 五.实验内容 用MATLAB语言实现对连续信号 进行采样;并验证采样定理。 1.对于 a.以样本/秒采样得到。求并画出。 % 模拟信号 Dt=0.00005; t=-0.0025:Dt:0.0025; %t在-0.0025到0.0025之间; xa=exp(-1000*abs(2*t)); %离散时间信号 Ts=0.0001;n=-25:1:25;x=exp(-1000*abs(2*n*Ts)); %离散时间傅立叶变换 K=500; k=0:1:K; w=pi*k/K; X=x*exp(-j*n'*w); X=real(X); w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(1,1,1) subplot(2,1,1);plot(t*1000,xa); xlabel('t 毫秒'); ylabel('x1(n)'); title('离散信号');hold on stem(n*Ts*1000,x);gtext('Ts=0.2毫秒');hold off subplot(2,1,2); plot(w/pi,X); xlabel('以pi为单位的频率'); ylabel('X1(w)'); title('离散时间傅立叶变换');

连续时间信号和系统时域分析及MATLAB实现

MATLAB课程设计任务书 姓名:王** 学号:2010******010 题目: 连续时间信号和系统时域分析及MATLAB实现 初始条件: MATLAB 7.5.0 ,Windows XP系统 实验任务: 一、用MATLAB实现常用连续时间信号的时域波形(通过改变参数,分析其时域特性). 1、单位阶跃信号, 2、单位冲激信号, 3、正弦信号, 4、实指数信号, 5、虚指数信号, 6、复指数信号. 二、用MATLAB实现信号的时域运算 1、相加, 2、相乘, 3、数乘, 4、微分, 5、积分 三、用MATLAB实现信号的时域变换(参数变化,分析波形变化) 1、反转, 2、使移(超时,延时), 3、展缩, 4、倒相, 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解, 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子,对每个例子,要求画出对应波形. 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形. 给出几个典型例子,四种调用格式. 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形.

目录1 MATLAB简介1 1.1 MATLAB设计目的1 1.2 MATLAB语言特点1 2常用连续时间信号的时域波形1 2.1单位阶跃信号1 2.2单位冲激信号2 2.3正弦信号3 2.4实指数信号4 2.5虚指数信号5 2.6复指数信号6 3 连续时间信号的时域运算7 3.1相加7 3.2相乘8 3.3数乘9 3.4微分10 3.5积分11 4.1反转12 4.2时移13 4.3展缩14

相关文档
最新文档