112-BF7-0103796-郝彩虹4.doc

本科学生综合性实验报告 学 院： 软件与通信工程 课程名称： 信号与系统 专业班级： 通信101 姓 名： 郝彩虹 学 号： 0103796 学生实验报告 学生姓名 郝彩虹 学号 0103796 同组人 实验项目 连续时间信号的采样与重构 □必修 □选修 □演示性实验 □验证性实验 □操作性实验 □综合性实验 实验地点 H123 实验仪器台号 8 指导教师 蒋娜 实验日期及节次 周四3、4 一、实验目的及要求： 1、目的 通过MATLAB编程实现对时域抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 2、内容及要求 二、仪器用具： 仪器名称 规格/型号 数量 备注 MATLAB 计算机 三、实验方法与步骤: 实验任务 连续信号f(t)=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t)，(1)计算该信号的奈奎斯特频率；（2）以不同的采样频率对该信号进行采样，画出采样前和采样后信号的频谱，对比观察和分析信号临界采样、过采样和欠采样情况下，信号频谱有何变化；（3）经过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号f(t)。 分析：该信号在w0=±40π(rad/s)的地方有近似于冲激信号的谱线。根据奈奎斯特采样定理，该信号的临界采样角频率为ws=2*w0=80π(rad/s)，临界采样周期Ts=0.025s。采样周期大于0.025s，则欠采样，频谱发生混叠；采样周期小于0.025s，则过采样，频谱不发生混叠，用低通滤波器可以将原信号的频谱滤出。 display(奈奎斯特周期0.025s， Ts0.025s过采样频谱不混叠； Ts0.025s欠采样频谱混叠); display(please input the value of sample period); Ts=input(Ts=); t=-0.4:0.01:0.4; y=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t); subplot(221); plot(t,y); xlabel(t 单位:s,Fontsize,8); title(f(t)); N=200; W=2*pi*5; k=-N:N; w=k*W/N; Y=0.01*y*exp(-j*t*w); Y=abs(Y); subplot(222) plot(w/pi,Y); title(F(j\omega)); xlabel(\omega 单位:pi); subplot(223); plot(t,y,b:); hold on; t2=-0.4:Ts:0.4; y2=cos(8*pi*t2)+2*sin(40*pi*t2)+cos(24*pi*t2); stem(t2,y2); xlabel(t 单位:s,Fontsize,8); title(fs(t)); hold off; Y2=Ts*y2*exp(-j*t2*w); Y2=abs(Y2); subplot(224); plot(w/pi,Y,b) xlabel(\omega 单位:pi); title(Fs(j\omega)); hold on; plot(w/pi,Y2,r); hold off; 当取Ts=0.01s，过采样 当Ts=0.03s，欠采样 Ts=0.025s，临界采样； 经过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号为f(t)， 程序如下： wm=40*pi; wc=1.2*wm; Ts=[0.03]; N=length(Ts); for k=1:N; n=-200:200; nTs=n*Ts(k); fs=(cos(8*pi*nTs)+2*sin(40*pi*nTs)+cos(24*pi*nTs)).*(Heaviside(nTs+pi)-Heaviside(nTs-pi)); t=-0.25:0.001:0.25; ft=f