matlab实现信号的混叠现象,应用matlab实现连续信号的采样与重构

应用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 2 1  m t f 2 1 2 1     1 1 2 T   是原信号频谱的周期性重复，它每隔重复出现一次。当 ＞2s  时，不会出现混叠现象，原信号的频谱的形状不会发生变化，从而能从采样 信号中恢复原信号。 (注： ＞2的含义是：采样频率大于等于 s  信号最高频率的2倍；这里的“不混叠”意味着信号频谱没有被破坏，也就为 后面恢复原信号提供了可能！) (a) (b)(c) 图* 抽样定理 a) 等抽样频率时的抽样信号及频谱(不混叠) b) 高抽样频率时的抽样信号及频谱(不混叠) c) 低抽样频率时的抽样信号及频谱(混叠) 2.1.2信号采样 如图1所示，给出了信号采样原理图信号采样原理图(a) 由图1可见， ，其中，冲激采样信号 的表达式 ) ( ) ( ) ( t t f t f s T s    ) (t s T  为：       n s T nT t t s ) ( ) (   其傅立叶变换为 ，其中 。设 ， 分别为       n s s n ) (     s s T   2  ) (  j F ) (  j F s ， 的傅立叶变换，由傅立叶变换的频域卷积定理，可得 ) (t f ) (t f s              n s s n s s s n j F T n j F j F )] ( [ 1 ) ( * ) ( 2 1 ) (          若设 是带限信号，带宽为 ， 经过采样后的频谱 就是将 ) (t f m  ) (t f ) (  j F s 在频率轴上搬移至 处(幅度为原频谱的 倍) 。因 ) (  j F   , , , , , 0 2 ns s s       s T 1 此，当 时，频谱不发生混叠；而当 时，频谱发生混叠。 m s   2  m s   2  一个理想采样器可以看成是一个载波为理想单位脉冲序列 的幅值调制 ) (t T  器，即理想采样器的输出信号 ，是连续输入信号 调制在载波 上的 ) ( * t e ) (t e ) (t T  结果，如图2所示。 图2 信号的采样 用数学表达式描述上述调制过程，则有 ) ( ) ( ) ( * t t e t e T  理想单位脉冲序列 可以表示为 ) (t T      0 ) ( ) ( n T nT t t   其中 是出现在时刻 ，强度为1的单位脉冲。由于 的 ) ( nT t   nT t  ) (t e 数值仅在采样瞬时才有意义，同时，假设 0 0 ) (    t t e 所以 又可表示为 ) ( * t e * 0 ( ) ( ) ( ) n e t e nT t nT       2.1.3信号重构 设信号 被采样后形成的采样信号为 ，信号的重构是指由 经过 ) (t f ) (t f s ) (t f s 内插处理后，恢复出原来信号 的过程。又称为信号恢复。 ) (t f 若设 是带限信号，带宽为 ，经采样后的频谱为 。设采样频率 ) (t f m  ) (  j F s ，则由式(9)知 是以 为周期的谱线。现选取一个频率特性 m s   2  ) (  j F s s  (其中截止频率 满足 )的理想低通滤波         c c s T j H      0 ) ( c  2 s c m      器与 相乘，得到的频谱即为原信号的频谱 。 ) (  j F s ) (  j F 显然， ，与之对应的时域表达式为 ) ( ) ( ) (    j H j F j F s (10) ) ( * ) ( ) ( t f t h t f s  而               n s s n s s nT t nT f nT t t f t f ) ( ) ( ) ( ) ( ) (   ) ( )] ( [ ) ( 1 t Sa T j H F t h c c s        将 及 代入式(10)得 ) (t h ) (t f s(11)         n s c s c s c c s s nT t Sa nT f T t Sa T t f t f )] ( [ ) ( ) ( * ) ( ) (      式(11)即为用 求解 的表达式，是利用 MATLAB 实现信号重构的 ) ( s nT f ) (t f 基本关系式，抽样函数 在此起着内插函数的作用。 ) ( t Sa c  例：设 ，其 为： t t t Sa t f sin ) ( ) (   ) (  j F         1 0 1 ) (     j F 即 的带宽为 ，为了由 的采样信号 不失真地重构 ，由时 ) (t f 1  m  ) (t f ) (t f s ) (t f 域采样定理知采样间隔 ，取 (过采样) 。利用 MATLAB 的      m s T  7 . 0  s T 抽样函数 来表示 ，有 。据此可知： t t t Sinc   ) sin( ) (  ) (t Sa ) / ( ) (  t Sinc t Sa          n s c s c s c c s s nT t Sinc nT f T t Sa T t