DSP论文第一章.docxVIP

基于MATLAB的时域采样定理杨凤凤严亭关珊珊谌瑶谢叶青 (湖北文理学院物理与电子工程学院湖北襄阳 441053)摘要：本次应用MATLAB验证时域采样定理，加深理解时域采样定理的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采样、频谱分析和采样信号恢复的方法。计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。关键字：时域采样定理；连续信号；MATLAB；采样信号从模拟信号到数字信号，即从连续信号到离散信号的转换都是通过离散采样完成的，采样频率即每秒钟采样的个数。要保证信号不失真，则采样频率要大于信号最高频率的两倍。对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。1 时域采样定理的概念假设离散时间信号 是通过对连续时间信号 取样获得(1)由 卷积定理 , 则 xi(t) 的傅里叶变换为 :(2)其中，取样定理的叙述为: 如果 在 时 , = 0, 则连续时间信号 完全可以由离散时间信号重建恢复 , 其充要条件是:取样频率 ≥2 ; = 2被称为连续时间信 号的奈奎斯特频率。当 2时 , 则会在进行 频谱复制时产生重叠,通常称为混叠,则连续时间信号不能由它的取样信号重建恢复。2 时域采样定理的验证2.1验证方法（步骤）（1）画出连续时间信号的时域波形及其幅频特性曲线，信号为　fy= sin(2*pi*60*t)+cos(2*pi*30*t)对信号进行采样，得到采样序列,画出采样频率分别为80Hz,120Hz，150Hz时的采样序列波形；对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。（4）对信号进行谱分析，观察与3中结果有无差别。（5）由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。2.2 MATLAB实现编程 %实现采样频谱分析绘图函数 function fz=caiyang(fy,fs) %第一个输入变量是原信号函数，信号函数fy以字符串的格式输入 %第二个输入变量是采样频率 fs0=10000; tp=0.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(fx1)]) % 画原信号幅度频谱 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),stem(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))%信号的恢复及频谱函数 function fh=huifu