资源描述：
《matlab时域采样定理》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、一、实验名称：matlab验证时域采样定理二、实验目的 本次课程设计应用MATLAB验证时域采样定理。了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立工作能力。 加深理解时域采样定理的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采样、频谱分析和采样信号恢复的方法。计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。三、实验原理 MATLAB是一套功能十分强大
2、的工程计算及数据分析软件，广泛应用于各行各业。MATLAB是矩阵实验室之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。 1、对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。 2、设连续信号的的最高频率为fc，如果采样频率fs，那么采样信号可以唯一的恢复出原连续信号，否则会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。四、实验步骤 1、画出连续时间信号的时域波形及其幅频特性曲线，信号为 (x)=sin(2*pi*60*t)+cos(2*pi*25*t)+
3、cos(2*pi*30*t)； 2、对信号进行采样，得到采样序列，画出采样频率分别为80Hz，120Hz，150Hz时的采样序列波形； 3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。 4、对信号进行谱分析，观察与3中结果有无差别。 5、由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。 五、MATLAB实现编程%实现采样频谱分析绘图函数functionfz=caiyang(fy,fs)%第一个输入变量是原信号函数，信号函数fy以字符串的格式输入%第二个输入变量是采样频率fs0=10000;tp=0.
4、1;t=[-tp:1/fs0:tp];k1=0:999;k2=-999:-1;m1=length(k1);m2=length(k2);f=[fs0*k2/m2,fs0*k1/m1];%设置原信号的频率数组w=[-2*pi*k2/m2,2*pi*k1/m1];fx1=eval(fy);FX1=fx1*exp(-j*[1:length(fx1)]'*w);%求原信号的离散时间傅里叶变换figure%画原信号波形subplot(2,1,1),plot(t,fx1,'r')title('原信号'),xlabel('时间t(s)')axis([min(t),max(t),min(fx1),max(fx
5、1)])%画原信号幅度频谱subplot(2,1,2),plot(f,abs(FX1),'r')title('原信号幅度频谱'),xlabel('频率f(Hz)')axis([-100,100,0,max(abs(FX1))+5])%对信号进行采样Ts=1/fs;%采样周期t1=-tp:Ts:tp;%采样时间序列f1=[fs*k2/m2,fs*k1/m1];%设置采样信号的频率数组t=t1;%变量替换fz=eval(fy);%获取采样序列FZ=fz*exp(-j*[1:length(fz)]'*w);%采样信号的离散时间傅里叶变换figure%画采样序列波形subplot(2,1,1),st
6、em(t,fz,'.'),title('取样信号'),xlabel('时间t(s)')line([min(t),max(t)],[0,0])%画采样信号幅度频谱subplot(2,1,2),plot(f1,abs(FZ),'m')title('取样信号幅度频谱'),xlabel('频率f(Hz)')%信号的恢复及频谱函数functionfh=huifu(fz,fs)%第一个输入变量是采样序列%第二个输入变量是得到采样序列所用的采样频率T=1/fs;dt=T/10;tp=0.1;t=-tp:dt:tp;n=-tp/T:tp/T;TMN=ones(length(n),1)*t-n'*T*ones
7、(1,length(t));fh=fz*sinc(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:length(fh)]'*w);%恢复后的信号的离散时间傅里叶变换figure%画恢复后的信号的波形subplot(2,1,1),plot(t
