数字信号处理实验六时域采样与信号重建.doc

实验目的：
1。了解用MATLＡＢ语言进行时域抽样与信号重建的方法

3。掌握采样频率的确定方法和内插公式的编程方法.
实验内容
1认真阅读并输入实验原理与方法中介绍的例子，观察输出波形曲线，理解每一条语句的含义。．
２．已知一个连续时间信号f(t）=ｓｉnc(t)。取最高有限带宽频率fm=1Hz。
（１)分别显示原连续时间信号波形和Fｍ=ｆm、Fm＝2fm、Fｍ=3ｆm三种情况下抽样信号的波形。
实验程序：
dt＝0．１； f0＝1； Ｔ0=1/f0;
fm=f0； Tm=1/fｍ；
ｔ=—２:dｔ:2；
f=ｓinc(t）;
sｕbｐlot(４,１，1）,pｌot(ｔ，f,'k’）;
ａxｉs(［min（ｔ） max(ｔ）　1．1*miｎ（f）　1。1*max(f）］)；
titｌe（＇原连续信号和抽样信号＇）；
for i=1:3;
fs=ｉ＊ｆm; Ts=1/fｓ；
　 n=-2：Ｔs：２;
f=sinc(n）；
　 ｓｕbｐloｔ（４,1,i＋1），stem（n，ｆ,'fiｌｌed’，'k')；
　 　 axｉs（[mｉｎ(ｎ) max（n）　＊min（f)　*ｍax(f)］）;
eｎd
实验截图：
求解原连续信号波形和抽样信号所对应的幅度谱.
实验程序：
ｄｔ＝0。１;t=-4:ｄt：4;
N=lｅngtｈ（t）;f=sinc(ｔ)；Ｔm=１；fm=１/Ｔm;
ｗｍ=2＊pｉ*fｍ；k=1：N；
w1=k*wm/Ｎ；
F1=f＊ｅxp（－j*ｔ＇＊w1)*dt;
sｕbplｏt（4,1，1),ploｔ（ｗ1/（2＊ｐi),abs（F1))；griｄ
axiｓ([0　ｍax(4*fｍ) 1。1*miｎ(F１） ＊max（F1)］)；
foｒ ｉ=1：3；
if i〈= ２　ｃ=0 ，elsｅ ｃ＝，eｎｄ
　 fs＝（４－i+ｃ）*fm；　Ts=1／fs；
　 n＝-4：Ts：4;
f=sinc(ｎ）;
　 N=leｎgtｈ(n)； ｗm=2＊ｐi＊fs；
　 ｋ=１:N；
　 w＝ｋ*wｍ/N;
　　　 F＝f＊exｐ（—j＊n'*ｗ）*Ts;
sｕbpｌoｔ（4，１,5－i），plot(w/（2＊pi）,abｓ(Ｆ）,’k＇)；grid
　 axiｓ（［０ ｍax(4*fｍ） ＊min(F) 1．１＊max(F）］）;
ｅｎd