应用_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≥，或频域间隔mt f 2121≤=πω（其中112T πω=）。

采样信号 的频谱是原信号频谱 的周期性重复，它每隔 重复出现一次。

当s ω＞2时，不会出现混叠现象，原信号的频谱的形状不会发生变化，从而能从采样信号中恢复原信号。

＞2的含义是：采样频率大于等于信号最高频率的2倍；这里的“不（注：s混叠”意味着信号频谱没有被破坏，也就为后面恢复原信号提供了可能！）(a)(b)(c)图* 抽样定理a)等抽样频率时的抽样信号及频谱（不混叠）b)高抽样频率时的抽样信号及频谱（不混叠）c) 低抽样频率时的抽样信号及频谱（混叠）2.1.2信号采样如图1所示，给出了信号采样原理图信号采样原理图（a ）由图1可见，)()()(t t f t f s T s δ⋅=，其中，冲激采样信号)(t s T δ的表达式为：∑∞-∞=-=n sT nT t t s)()(δδ其傅立叶变换为∑∞-∞=-n s s n )(ωωδω，其中ss T πω2=。

设)(ωj F ，)(ωj F s 分别为)(t f ，)(t f s 的傅立叶变换，由傅立叶变换的频域卷积定理，可得∑∑∞-∞=∞-∞=-=-=n ss n s s s n j F T n j F j F )]([1)(*)(21)(ωωωωδωωπω若设)(t f 是带限信号，带宽为m ω， )(t f 经过采样后的频谱)(ωj F s 就是将)(ωj F 在频率轴上搬移至 ,,,,,02ns 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 δδ其中)(nT t -δ是出现在时刻nT t =，强度为1的单位脉冲。

由于)(t e 的 数值仅在采样瞬时才有意义，同时，假设00)(<∀=t t e所以)(*t e 又可表示为*()()()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 。

设采样频率msωω2≥，则由式（9）知)(ωj F s是以sω为周期的谱线。

现选取一个频率特性⎪⎩⎪⎨⎧><=ccsT j H ωωωωω0)(（其中截止频率c ω满足2sc m ωωω≤≤）的理想低通滤波器与)(ωj F s 相乘，得到的频谱即为原信号的频谱)(ωj F 。

显然，)()()(ωωωj H j F j F s =，与之对应的时域表达式为)(*)()(t f t h t f s = （10）而∑∑∞-∞=∞-∞=-=-=n s s n s s nT t nT f nT t t f t f )()()()()(δδ)()]([)(1t Sa T j H F t h ccsωπωω==- 将)(t h 及)(t f s 代入式（10）得∑∞-∞=-==n scscsccssnT t Sa nT f T t Sa T t f t f )]([)()(*)()(ωπωωπω （11）式（11）即为用)(s nT f 求解)(t f 的表达式，是利用MATLAB 实现信号重构的基本关系式，抽样函数)(t Sa c ω在此起着内插函数的作用。

例：设ttt Sa t f sin )()(==，其)(ωj F 为： ⎪⎩⎪⎨⎧><=11)(ωωπωj F即)(t f 的带宽为1=m ω，为了由)(t f 的采样信号)(t f s 不失真地重构)(t f ，由时域采样定理知采样间隔πωπ=<ms T ，取π7.0=sT （过采样）。

利用MATLAB 的抽样函数tt 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 f t f )]([)()(*)()(πωπωωπω通过以上分析，得到如下的时域采样定理：一个带宽为w m 的带限信号f(t)，可唯一地由它的均匀取样信号fs(nTs)确定，其中，取样间隔Ts<π/w m, 该取样间隔又称为奈奎斯特间隔。

根据时域卷积定理，求出信号重构的数学表达式为：式中的抽样函数Sa(wct)起着内插函数的作用，信号的恢复可以视为将抽样函数进行不同时刻移位后加权求和的结果，其加权的权值为采样信号在相应时刻的定义值。

利用MATLAB 中的抽样函数来表示Sa(t)，有，，于是，信号重构的内插公式也可表示为：2.2设计的思路连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。

严格来说，MATLAB 并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。

当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。

时域对连续时间信号进行采样，是给它乘以一个采样脉冲序列，就可以得到采样点上的样本值，信号被采样前后在频域的变化，可以通过时域频域的对应关系分别求得了采样信号的频谱。

在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。

这样，抽样定理为连续时间信号与离散时间信号的相互转换提供了理论依据。

通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。

2.3详细设计过程2.3.1)(t Sa 的临界采样及重构1实现程序代码当采样频率小于一个连续的同信号最大频率的2倍，即m s ωω2=时，称为临界采样. 修改门信号宽度、采样周期等参数，重新运行程序，观察得到的采样信号时域和频域特性，以及重构信号与误差信号的变化。

Sa(t)的临界采样及重构程序代码；wm=1;%升余弦脉冲信号带宽wc=wm; %频率 Ts=pi/wm; %周期ws=2.4*pi/Ts; %理想低通截止频率 n=-100:100; %定义序列的长度是201nTs=n*Ts %采样点f=sinc(nTs/pi); %抽样信号Dt=0.005;t=-20:Dt:20;fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,le ngth(t)))); %信号重建t1=-20:0.5:20;f1=sinc(t1/pi);subplot(211);stem(t1,f1);xlabel('kTs');ylabel('f(kTs)');title('sa(t)=sinc(t/pi)的临界采样信号');subplot(212);plot(t,fa)xlabel('t');ylabel('fa(t)');title('由sa(t)=sinc(t/pi)的临界采样信号重构sa(t)');grid;2程序运行运行结果图与分析图3 )(tSa的临界采样及重构图运行结果分析：为了比较由采样信号恢复后的信号与原信号的误差，可以计算出两信号的绝对误差。

当t选取的数据越大，起止的宽度越大。

如有侵权请联系告知删除，感谢你们的配合！。