MATLAB实现抽样定理探讨及仿真

1、应用MATLAB实现抽样定理探讨及仿真课程设计的目的利用MATLAB,仿模信号抽样与恢复系统的实际实现，探讨过抽样和欠抽样的信号以及抽样与恢复系统的性能。课程设计的原理抽样信号模拟信号经过(A/D)变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs,重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。时域采样定理从采样信号丸(£)恢复原信号处)必需满足两个条件：(1)处)必须是带限信号，其频谱函数在|如吐各处为零；对信号的要求，即只有带限信号才能适用采样定理。(2)取样频率不能过低，必须尖2气

2、或£。对取样频率的要求，即取样频率要足够大,采得的样值要足够多，才能恢复原信号。如果采样频率殴工5大于或等于，即叫丑2%2为连续信号氾的有限频谱，那么采样离散信号E能无失真地恢复到原来的连续信号冷。一个频谱在区间-战，瑶以外为零的频带有限信号了们,可唯一地由其在均匀间隔雾需上的样点值)所确定。根据时域与频域的对称性，可以由时域采样定理直接推出频域采样定理。四小5(a)图2.1抽样定理等抽样频率时的抽样信号及频谱不混叠高抽样频率时的抽样信号及频谱不混叠低抽样频率时的抽样信号及频谱混叠2.1信号采样如图1所示，给出了信号采样原理图信号采样原理图a由图1可见，fs(t)f(t)Ts(t)，

3、其中，冲激采样信号启(t)的表达式为:Ts(t)(tnTs)其傅立叶变换为s(nnns)，其中s。设F(j),Fs(j)分别为f(t),fs(t)的T傅立叶变换，由傅立叶变换的频域卷积定理，可得1Fs(j)土F(j)*s(2n1ns)yFj(ns)Tsn假设设f(t)是带限信号，带宽为m,f(t)经过采样后的频谱Fs(j)就是将F(j)在频率轴上搬移至0,s,2s,ns,处幅度为原频谱的1Ts倍。因此，当s2皿时，频谱不发生混叠；而当S2m时，频谱发生混叠。2.1.3信号重构设信号f(t)被采样后形成的采样信号为fs(t)，信号的重构是指由fs(t)经过插处理后，恢复出原来信号f(t)的过程。

4、又称为信号恢复。假设设f(t)是带限信号，带宽为经采样后的频谱为Fs(j)。设采样频率由式9知Fs(j)是以s为周期的谱线。Ts现选取一个频率特性H(j)0其中截止频率c满足成的理想低通滤波器与FJj)相乘，得到的频谱即为原信号的频谱F(j)。显然，F(jFs(j)H(j),与之对应的时域表达式为f(t)h(t)*fs(t)10fs(t)f(t)(tnTs)f(nTs)(tnTs)h(t)_1FH(j)Ts-Sa(ct)将h(t)及fs(t)代入式10得f(t)fs(t)*Ts二Sa(ct)Tsf(nTs)Sac(tnTs)11MATLAB实现信号重构的根本关系式，抽样函式11即为用f(nTs

5、)求解f(t)的表达式，是利用数Sa(ct)在此起着插函数的作用。三、抽样定理的仿真和探讨3.1.1Sa(t)的临界采样及重构图当采样频率小于一个连续的同信号最大频率的2倍，即s2m时，称为临界采样.修改门信号宽度、采样周期等参数，重新运行程序，观察得到的采样信号时域和频域特性，以及重构信号与误差信号的变化。程序运行结果:由sa(t)=sinc(V同的1临界采祥信号重构sap)I.!JI.i.I.J,'-20-15F书051015203.1.2Sa(t)的过采样及重构当采样频率大于一个连续的同信号最大频率的2倍，即s2m时，称为过采样的变化。在不同采样频率的条件下，观察对应采样信号的时

6、域和频域特性，以及重构信号与误差信程序运行结果：WLLsim.啊的原信号节1IIr泛1.1I'.1甲平TfT"?甲ftII二】E已C.3J=3-T=J，、】-10-3-64-2024kTg631004r0.2-0L由sa(t)=sinci：t/pi)的过采祥信号重构“0)IIT-0-6-202469103.1.3Sa(t)的欠采样及重构当采样频率小于一个连续的同信号最大频率的2倍，即s2m时，称为过采样。利用频域滤波的方法修改实验中的局部程序，完成对采样信号的重构。程序运行结果：sa(t>=sjnc(t/pi)的原信号由sa(t)=sinc(Vpi)l欠采样信号重构飒)

7、误差分析：绝对误差error已大为增加，其原因是因采样信号的频谱混叠，使得在。区域的频谱相互“干扰"所致。四、课题研讨的小结该课程设计使我们对采样定理的一些根本公式得到了进一步稳固。在整个实验过程中，我们查阅了很多相关知识，从这些书籍中我们受益良多。虽然学习过采样过程和恢复过程，但是认识不深，实践能力也有所欠缺，通过这次实验对采样过程和恢复过程有了进一步掌握。通过实验的设计使我们对采样定理和信号的重构有了深一步的掌握，也让我们在实践的过程中了解到团队合作的重要性。虽然在实验过程中出现很多错误，但是在教师的帮助和团队成员的齐心协力下，不断的修正错误，同时也学会了MATLAB号表示的根本

8、方法及绘图函数的调用。虽然刚开场我们对MATLAB的根本使用方法没有太深刻的认识，但是该实验使我们对MATLAB函数程序的根本构造有所了解，也提高了我们独立完成实验的能力和理论联系实际的应用能力。通过这次课程设计，我们不仅学到了学科知识，锻炼了实践能力，更重要的是学到了学习的方法和团队合作的重要性。我们团队分工有序，每个人都能按时完成各自的任务。在遇到问题时，大家都能够互相理解，互相帮助，最后圆满完成课题！附录:一、Sa(t)的临界采样及重构1.Sa(t*勺临界采样及重构程序代码；wm=1;wc=wm;Ts=pi/wm;ws=2.4*pi/Ts;n=-100:100;nTs=n*Ts;f=si

9、nc(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,length(t);subplot(311);plot(t,fa)xlabel('t');ylabel('fa(t)');title('sa(t)=sinc(t/pi)的原信号');grid;t1=-20:0.5:20;f1=sinc(t1/pi);subplot(312);stem(t1,f1);xlabel('kTs');ylabel(

10、'f(kTs)');title('sa(t)=sinc(t/pi)的临界采样信号');grid;subplot(313);plot(t,fa)xlabel('t');ylabel('fa(t)');title('由sa(t)=sinc(t/pi)的临界采样信号重构sa(t)');grid;2程序运行运行结果图与分析图3.1.1Sa(t)的临界采样及重构图运行结果分析：为了比拟由采样信号恢复后的信号与原信号的误差，可以计算出两信号的绝对误差。当t选取的数据越大，起止的宽度越大。二、Sa(t)的过采样及重构1.Sa(t

11、*勺过采样及重构程序代码；wm=1;wc=1.1*wm;Ts=1.1*pi/wm;ws=2*pi/Ts;n=-100:100;nTs=n*Ts;f=sinc(nTs/pi);Dt=0.005;t=-10:Dt:10;fa=f*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t);subplot(411);plot(t,fa)xlabel('t');ylabel('fa(t)');title('sa(t)=sinc(t/pi)的原信号');grid;error=a

12、bs(fa-sinc(t/pi);t1=-10:0.5:10;f1=sinc(t1/pi);subplot(412);stem(t1,f1);xlabel('kTs');ylabel('f(kTs)');title('sa(t)=sinc(t/pi)的采样信号');grid;subplot(413);plot(t,fa)xlabel('t');ylabel('fa(t)');title('由sa(t)=sinc(t/pi)的过采样信号重构sa(t)');grid;subplot(414);plot(

13、t,error);xlabel('t');ylabel('error(t)');title('过采样信号与原信号的误差error(t)');grid;2.程序运行运行结果图与分析。24&s?(t)=sinc(t/pi|的采，号.£V.：v20-2G一小EmeA2L-J罕#5686。，卜20年1IIIQ-S102Q£4kTs由ss(tsinctVpi)的过采样信号重构”tJJIxIL-4过采杵信号亍原信号的误箜5顷-108-6-A202d图3.1.2Sa(t)的过采样信号、重构信号及两信号的绝对误差图运行分析：将原始信号

14、分别修改为抽样函数Sa(t)正弦信号sin(20*pi*t)+cos(20*pi*t)、指数信号e-2tu(t)时，在不同采样频率的条件下，可以观察到对应采样信号的时域和频域特性，以及重构信号与误差信号的变化。三、Sa(t；的欠采样及重构1.Sa(t*勺欠采样及重构程序代码；wm=1;wc=wm;Ts=2.5*pi/wm;ws=2*pi/Ts;n=-100:100;nTs=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,length

15、(t);error=abs(fa-sinc(t/pi);subplot(411);plot(t,fa)xlabel('t');ylabel('fa(t)');title('sa(t)=sinc(t/pi)的原信号');grid;t1=-20:0.5:20;f1=sinc(t1/pi);subplot(412);stem(t1,f1);xlabel('kTs');ylabel('f(kTs)');title('sa(t)=sinc(t/pi)的采样信号sa(t)');grid;subplot(413);plot(t,fa)xlabel('t');ylabel('fa(t)');title('由sa(t)=sinc(t/pi)的欠采样信号重构sa(t)');grid;