原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数字信号处理实验六时域采样与信号的重建
实验目的：
1.了解用MATLAB语言进行时域抽样与信号重建的方法
2.进一步加深对时域信号抽样与恢复的基本原理的理解
3.掌握采样频率的确定方法和内插公式的编程方法。
实验内容
1认真阅读并输入实验原理与方法中介绍的例子，观察输出波形曲线，理解每一条语句的含义。.
2．已知一个连续时间信号f(t)=sinc(t)。取最高有限带宽频率fm=1Hz。
（1）分别显示原连续时间信号波形和Fm=fm、Fm=2fm、Fm=3fm三种情况下抽样信号的波形。
实验程序：
dt=0.1; f0=1; T0=1/f0;
fm=f0; Tm=1/fm;
t=-2:dt:2;
f=sinc(t);
subplot(4,1,1),plot(t,f,k);
axis([min(t) max(t) 1.1*min(f) 1.1*max(f)]);
title(原连续信号和抽样信号);
for i=1:3;
fs=i*fm; Ts=1/fs;
n=-2:Ts:2;
f=sinc(n);
subplot(4,1,i+1),stem(n,f,filled,k);
axis([min(n) max(n) 1.1*min(f) 1.1*max(f)]);
end
实验截图：
求解原连续信号波形和抽样信号所对应的幅度谱。
实验程序：
dt=0.1;t=-4:dt:4;
N=length(t);f=sinc(t);Tm=1;fm=1/Tm;
wm=2*pi*fm;k=1:N;
w1=k*wm/N;
F1=f*exp(-j*t*w1)*dt;
subplot(4,1,1),plot(w1/(2*pi),abs(F1));grid
axis([0 max(4*fm) 1.1*min(F1) 1.1*max(F1)]);
for i=1:3;
if i= 2 c=0 ,else c=0.2,end
fs=(4-i+c)*fm; Ts=1/fs;
n=-4:Ts:4;
f=sinc(n);
N=length(n); wm=2*pi*fs;
k=1:N;
w=k*wm/N;
F=f*exp(-j*n*w)*Ts;
subplot(4,1,5-i),plot(w/(2*pi),abs(F),k);grid
axis([0 max(4*fm) 1.1*min(F) 1.1*max(F)]);
end
实验截图：
用时域卷积方法（内插公式）重建信号。
实验程序：
f0=1; T0=1/f0; dt=0.01;
fm=f0; Tm=1/fm;
t=-4:dt:4*f0;
x=sinc(t);
subplot(4,1,1),plot(t,x);
axis([min(t),max(t),1.1*min(x),1.1*max(x)]);
title(用时域卷积重建抽样信号)
for i=1:3;
fs=i*fm; Ts=1/fs;
n=-4:(4*f0)/Ts;
t1=-4:Ts:4*f0;
x1=sinc(n/fs);
T_N=ones(length(n),1)*t1-n*Ts*ones(1,length(t1))
xa=x1*sinc(fs*pi*T_N);
subplot(4,1,i+1),plot(t1,xa);
axis([min(t1),max(t1),1.1*min(xa),1.1*max(xa)]);
end
实验截图：
文档评论（0）