取样定理的matlab实现,取样定理的MATLAB实现

第28卷第6期 200 8年1 2月 成宁学院学报 Journal of Xianning University Vol。28，No．秽Dec．2∞8 文章编号：1006—5342(2008)06—0039—02 取样定理的MATLAB实现’ 邢国泉 (成字学院生物藏学王程学院，湖托 成宁437100) 接 要：敬祥定理是把摸耘弦号突威数字信琴驭祥菝率选取爨一条譬要寝麓。本文逮取一模拟信号，誊凝艇西乎 台上用不同的取样频率对其采样，并利用三次样条函数对取样信号进行重建，将熏建出的遗续信号与原连续信号相比较。从而在Matlab平台上毵觋了取样定理的仿真说明．关键词：取棒定理；Matlab；重建 中图分类号：TP271+．8l 文献标识码：A 0 弓l言 从模拟信号到数字信号，邸扶连续信号勤离散信号的转换都是通过离散采样完成的，采样频率就是每秒种采样的个数．根据香农(Shannon)采样定联，要保证信号不失真，采样频率要大于信号最高频率的两倍． 1 采样定理简介 假设离散时间信号麓(￡)是通过对连续时阕信号髫。(t)取样获得⋯， “ +∞ +一 蕺(；)=∑牧(￡Ⅻ(t一哟=∑魏(摊D馥￡一狂∞= 靠(￡)∑8(t—nr) (1) 由卷积定理，则髫；(￡)的傅熏叶变换为： _ ^ +锥 ％(jto)=参。(如)堆等∑艿(∞一㈨。)= ' +∞ 寺∑菇。(jw-j黝o) (2) ‘^=一∞ ，’一 其中∞o=等 』 致样定理的叙述力：如果故(汹)在l∞l>l峨l 时，石。(如)=0，则连续时间信号扎(I)完全可以由离散时间信号瓤(￡)重建恢复，其充要条件是：取样羧率甜e≥2to。；鳓=2ta。被称为连续融阕信号弗。(f)的奈黛斯特(Nyquist)频糕 当∞。《2to。时，则会在进行茁(t)频谱复制时产生重叠，邋常稳为混叠(Aliasing)，则连续时闻信 号不能由它的取样信号重建恢复。 2 运用Matlabl编程举例分析旧J 以下选定一个连续信号故(￡)=芎。㈣“，选定 ·收稿日期：2008一07一09 不同的取样频率对其取样重建，以越来说骥取样定理． clear； close all； ． ％Analog Signal Dt=0．000 5；t=一O．005：Dt：0．005；xa=exp(一1 000簟abs(1))i ％Continuous—time Fourier Transfornl Wmax=2南pi宰2 000： K=5∞；l【=0：1：K；W=k事Wmax／K； Xa=xa宰exp(一j枣t’臻W)幸Dt； Xa=real(Xa)； W=[一fliplr(w)，w(2：501)]； Xa=[fliplr(Xa)，Xa(2：501)】； figure(1)，subplot(2，l，1)；plot(t木1 000，飙)；xlabel(t in msec．，；ylabel(K a(￡)，； title(primary continuous signal，；subplot(2，1力；plotCW／(2,pi木1仰),Xa宰l㈣； xlabel(Frequency in kHz)；ylabel(X_p(j ka№ga)事l 圈∞；title(Continuous—time Fourier Transform)； ％Discrete—time Fourier Transfo瑚 Ts=0．000 2；n=一25：i：25；x=exp(一l 000宰abs(n宰Ts))； K一500；k=0：1：K； W。pi木k／K；X=x木exp(一j事n’宰w