matlab验证时域采样定理.docx

1、matlab验证时域采样定理第1章摘要 1第2章基本原理 2第3章实验步骤 5第4章MATLAB实现编程 5第5章实验结果与分析 85.1程序分析 85.2信号的波形及幅度频谱 85.3结果分析 9第6章总结 12参考文献 13第1章摘要一、 数字信号处理数字信号处理是将信号以数字方式表示并处理的理论和技术。 数字信号处理与模拟信号处理是信号处理的子集。数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。 因此 在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数 转换器实现。而数字信号处理的输出经常也要变换到模拟域， 这是通过数模转 换器实现的。数字信号处理的算法需要

2、利用计算机或专用处理设备。 数字信号处理技术 及设备具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等突出优点， 这些都是模拟信号处理技术与设备所无法比拟的。数字信号处理的核心算法是离散傅立叶变换 (DFT)，是DFT使信号在数字 域和频域都实现了离散化，从而可以用通用计算机处理离散信号。 而使数字信 号处理从理论走向实用的是快速傅立叶变换 (FFT) ,FFT的出现大大减少了 DFT 的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。随着大规模集成电路以及数字计算机的飞速发展，加之从 60年代末以来 数字信号处理理论和技术的成熟和完善， 用数字方法来处理信号，即数字信号 处

3、理，已逐渐取代模拟信号处理。随着信息时代、数字世界的到来，数字信号处理已成为一门极其重要的学 科和技术领域。二、 实验目的本次课程设计应用 MATLA验证时域采样定理。了解MATLA软件，学习应 用MATLA软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一 些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立 工作能力。加深理解时域采样定理的概念，掌握利用 MATLA分析系统频率响应的方法和掌握利用MATLA实现连续信号采样、频谱分析和采样信号恢复的方法。 计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总 结采样频率对信号恢复产生误差的影响，从而验

4、证时域采样定理。三、MATLA啲介绍MATLAE是 一套功能十分强大的工程计算及数据分析软件，广泛应用于各 行各业。MATLAB是矩阵实验室之意。除具备卓越的数值计算能力外，它还提 供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB勺基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形 式十分相似,故用MATLA来解算问题要比用C,FORTRA等语言完全相同的事情 简捷得多在新的版本中也加入了对C,FORTRAN,c+JAVA的支持.可以直接调 用,用户也可以将自己编写的实用程序导入到 MATLAB数库中方便自己以后 调用。第2章基本原理(1)时域采样定理1、

5、对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续 信号的频谱以采样频率为周期进行周期性的延拓形成的。2、 设连续信号的的最高频率为 Fmax，如果采样频率Fs 2Fmax，那么采 样信号可以唯一的恢复出原连续信号，否则 Fs 2Fmax会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号(2)设计原理图fa(t) 一连续信号取样脉冲信号(3)信号的时域采样与频谱分析对一个连续信号fa (t)进行理想采样的过程可以用下式表示fa(t) fa(t)S(t) ( 1)A 其中fa(t)为fa(t)的理想采样，S(t)为周期脉冲信号，即s(t) (t nT) (2)nA Afa(t

6、)的傅里叶变换Fa(j )为A 1Fa(j ) Faj( m s) (3)I mA上式表明，Fa(j )为Fa(j )的周期延拓，其延拓周期为采样角频率(s=2n /T )。只有满足采样定理时，才不会发生频率混叠失真。AA在计算机上用高级语言编程，直接按照(3)式计算fa(t)的频谱Fa(j )很不方便，下面导出用序列的傅里叶变换来计算 Fa(j )的公式。将(2)式代入(1)式，并进行傅里叶变换,AFa(j )fa(t) (t n T)ejtdtnnfa(t) (t nT)e j dt (4)nfa(nT)e j nT式中的fa (nT)就是采样后的序列f (n)，即：f (n)= fa (

7、nT)， f (n)的傅里叶变换为 F(ej ) f(n)e j n (5)nA比较和可知 Fa(j ) F(ej ) T (6)说明两者之间只在频率度量上差一个常数因子 To实验过程中应注意这一差别离散信号和系统在时域均可以用序列来表示， 序列图形给人以形象直观的 印象，它可以加深我们对信号和系统的时域特征的理解。 本实验还将观察分析几种信号及系统的时域特性。为了观察分析各种序列的频域特性，通常对 F(ej )在0 , 2n 上进行M点采样来观察分析。对长度为 N的有限长序列f (n)，有k k, k 0,1,.,M 1其中 M通常m应取的大一些，以便观察谱的细节变化。取模 F(ej k)可

8、绘出幅频特性曲线。 采样信号的恢复(内插函数法)设信号f(t)被采样后形成的采样信号为fs(t)，信号的重构是指由fs(t)经过内插处理后，恢复出原来信号f (t)的过程。又称为信号恢复。信号恢复的时域表达式f(t)h(t)* fs(t)(8)而fs(t)f(t) (tnnTs) f(nTs) (tnnTs)h(t)1F H(j )TsSa( ct)将h(t)及fs(t)代入式(8)得：c T c cf(t) fs(t)*Ts Sa( ct) 一 f(nTs)Sa c(t nTs) (9)n式(9)即为用f(nTs)求解f(t)的表达式，是利用MATLA实现信号恢复的 基本关系式，抽样函数Sa

9、( ct)在此起着内插函数的作用。内插公式表明模拟信号f(t)等于各采样点函数f (nTs)乘以对应内插函数 的总和，即只要采样频率高于信号频率的两倍，模拟信号就可用它的采样信号 代表，而不会丢失任何信息。这种理想低通滤波器的模拟信号完全等于模拟信 号f(t)，是一种无失真的恢复。第3章实验步骤1、 画出连续时间信号的时域波形及其幅频特性曲线，信号为f (x)= sin(2*pi*60*t)+cos(2*pi*25*t)+cos(2*pi*30*t) ；2、 对信号进行采样，得到采样序列，画出采样频率分别为80Hz, 120 Hz, 150 Hz时的采样序列波形；3、 对不同采样频率下的采样序

10、列进行频谱分析，绘制其幅频曲线，对比 各频率下采样序列和的幅频曲线有无差别。4、 对信号进行谱分析，观察与3中结果有无差别。5、 由采样序列恢复出连续时间信号 ，画出其时域波形，对比与原连续时 间信号的时域波形。第4章 MATLA实现编程澈现采样频谱分析绘图函数fun cti on fz=caiya ng(fy,fs)%第一个输入变量是原信号函数，信号函数 fy以字符串的格式输入%第二个输入变量是采样频率fs0=10000; tp=0.1;t=-tp:1/fs0:tp;k仁0:999; k2=-999:-1;m1=length(k1); m2=length(k2);f=fs0*k2/m2,fs

11、0*k1/m1; % 设置原信号的频率数组w=-2*pi*k2/m2,2*pi*k1/m1;fx仁 eval(fy);FX1=fx1*exp(-j*1:le ngth(fx1)*w);%求原信号的离散时间傅里叶变换figure%画原信号波形subplot(2,1,1),plot(t,fx1,r)title( 原信号),xlabel( 时间 t (s)axis(min(t),max(t),min(fx1),max(fx1) %FZ=fz*exp(-j*1:le ngth(fz)*w);%米样信号的离散时间傅里叶变换figure%画米样序列波形subplot(2,1,1),stem(t,fz,.)

12、.title( 取样信号),xlabel( 时间 t (s)lin e(mi n(t),max(t),0,0)%画采样信号幅度频谱subplot(2,1,2),plot(f1,abs(FZ),m)title( 取样信号幅度频谱),xlabel( 频率f (Hz)%信号的恢复及频谱函数fun cti on fh=huifu(fz,fs)%第一个输入变量是采样序列 .专业.整理.%第二个输入变量是得到采样序列所用的采样频率T=1/fs; dt=T/10; tp=0.1;t=-tp:dt:tp; n=-tp/T:tp/T;TMN=o nes(le ngth( n),1)*t -n *T*o nes(

13、1,le ngth(t);fh=fz*si nc(fs*TMN); % 由采样信号恢复原信号k1=0:999; k2=-999:-1;m1=length(k1); m2=length(k2);w=-2*pi*k2/m2,2*pi*k1/m1;FH=fh*exp(-j*1:le ngth(fh)*w);%恢复后的信号的离散时间傅里叶变换figure%画恢复后的信号的波形subplot(2,1,1),plot(t,fh,g),st仁sprintf( 由取样频率 fs=%d,fs);时间t (s)画重构信号的幅度频谱设置频率数组频率 f (Hz)st2=恢复后的信号；st=st1,st2; titl

14、e(st) , xlabel( axis(mi n( t),max(t),mi n( fh),max(fh) lin e(mi n( t),max(t),0,0) %f=10*fs*k2/m2,10*fs*k1/m1; %subplot(2,1,2),plot(f,abs(FH),g) title( 恢复后信号的频谱),xlabel( axis(-100,100,0,max(abs(FH)+2);%t函数输入一个信号f1=si n( 2*pi*60*t)+cos(2*pi*25*t)+cos(2*pi*30*t);% fs0=caiyang(f1,80); % 频率 fs 2fmax，即 欠采样fr0=huifu(fs0,80);fs1= caiyang(f1,120);% 频率 2也，临 界采样fr仁 huifu(fs1,120);fs2=caiyang(f1,150);% 频率 fs 2仏，即 过采样fr2=huifu(fs2,150);第5章实验结果与分析1、程序分析TMN=o nes(le ngth( n),1)*t- n*T*o nes(1,le ngth(t);fh=fz*s in c(fs*TMN);%由米样信号恢复原信号plot(t,f)%绘制fx的波形stem(t,f)%绘制一个二维杆图（画离散波形）subpolt(,)%在一个窗口画多个波形图f