信号与系统上机实习报告.doc

信号与线性系统上机实习报告 Sa（t）的临界采样及重构 程序代码 wm=1; %升余弦脉冲信号带宽 wc=wm; %频率 Ts=pi/wm; %周期 ws=2.4*pi/Ts; %理想低通截止频率 n=-100:100; %定义序列的长度是201 nTs=n*Ts %采样点 f=sinc(nTs/pi); %抽样信号 Dt=0.005;t=-20:Dt:20; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); %信号重建 t1=-20:0.5:20; f1=sinc(t1/pi); subplot(211); stem(t1,f1); xlabel(kTs); ylabel(f(kTs)); title(sa(t)=sinc(t/pi)的临界采样信号); subplot(212); plot(t,fa) xlabel(t); ylabel(fa(t)); title(由sa(t)=sinc(t/pi)的临界采样信号重构sa(t)); grid; 运行结果图与分析 图3 的临界采样及重构图 运行结果分析：为了比较由采样信号恢复后的信号与原信号的误差，可以计算出两信号的绝对误差。当t选取的数据越大，起止的宽度越大。 的过采样及重构 程序代码 wm=1; wc=wm; Ts=2.5 *pi/wm; ws=2*pi/Ts; n=-100:100; nTs=n*Ts f=sinc(nTs/pi); Dt=0.005;t=-20:Dt:20; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); error=abs(fa-sinc(t/pi)); t1=-20:0.5:20 f1=sinc(t1/pi); subplot(311); stem(t1,f1); xlabel(kTs); ylabel(f(kTs)); title(sa(t)=sinc(t/pi)的采样信号sa(t)); subplot(312); plot(t,fa) xlabel(t); ylabel(fa(t)); title(由sa(t)=sinc(t/pi)的欠采样信号重构sa(t)); grid; subplot(313); plot(t,error); xlabel(t); ylabel(error(t)); title(欠采样信号与原信号的误差error(t)); 运行结果图与分析 图5的欠采样信号、重构信号及两信号的绝对误差图 误差分析：绝对误差error已大为增加，其原因是因采样信号的频谱混叠，使得在区域内的频谱相互“干扰”所致。 总结： 通过几次上机实习，我学到了很多。首先，我学会了MATLAB中信号表示的基本方法及绘图函数的调用，并掌握了Matlab的一些基本操作和应用Matlab解决一些简单的信号与系统的问题。其次，我对采样定理的一些基本公式得到了进一步巩固，通过上机练习对采样定理有了更深刻的认识。 虽然刚开始我对MATLAB的基本使用方法并不熟练。但通过该实验我对MATLAB的使用熟练了许多，对MATLAB函数程序的基本结构也有了一些了解。 在这几次上机操作中，出现了很多错误，也使我意识到了自身的不足和相应知识的欠缺。我觉得这是这次实习最大的收获。