当前位置:文档之家› 数字信号处理实验六-时域采样与信号的重建

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

实验目的:
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
实验截图:
(2)求解原连续信号波形和抽样信号所对应的幅度谱。

实验程序:
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
实验截图:
(3)用时域卷积方法(内插公式)重建信号。

实验程序:
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
实验截图:。

相关主题