信号与系统实验报告信号采样与重构声音的延时与混响.doc

《信号与系统》 ——课程设计 实验一 信号的采样与重构 一、实验内容： 1.应用MATLAB实现连续信号的采样与重构仿真，了解MATLAB软件，学习应用MATLAB软件的仿真技术。 2.加深理解采样与重构的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采用与重构的方法。计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差。 3. 加深对采样定理的理解和掌握，以及对信号恢复的必要性；掌握对连续信号在时域的采样与重构的方法。 二、实验原理 （1）连续时间信号 连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。 (2)采样定理 模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率 fs，重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。 时域采样定理从采样信号恢复原信号必需满足两个条件： a、必须是带限信号，其频谱函数在 ＞ 各处为零；（对信号的要求，即只有带限信号才能适用采样定理。） b、 取样频率不能过低，必须 ＞2 （或 ＞2）。 一个理想采样器可以看成是一个载波为理想单位脉冲序列的幅值调制器。 图2 信号的采样 （4） 信号重构 设信号被采样后形成的采样信号为，信号的重构是指由经过内插处理后，恢复出原来信号的过程,又称为信号恢复。 三、实验步骤及代码 (一).%%%%%%%%%%% 产生一个连续sin（）信号 %%%%%%%%%%%%%%%%%% f=100; t=(1:50)/2000; %时间轴步距 x=sin(2*pi*t*f); figure subplot(211); plot(x); %绘制x(t)的图形图片号加底框 xlabel(t);ylabel(x(t)); title(连续时间信号sin（）的波形); %图片命名 grid; n=0:255; %长度 N=256; %设采样点的N值 Xk=abs(fft(x,N)); subplot(212); %频域波形 plot(n,Xk); axis([0 N 1.2*min(Xk) 1.2*max(Xk)]); %可用axis函数来调整图轴的范围 xlabel(时域频谱波形图); ylabel(|Xk|); title(信号sin（）的频谱波形); (二)%%%%%%%%%%%% 对原始信号进行采样并滤波重构 %%%%%%%%%%%% t1=3*t; f1=sin(2*pi*t1*f); figure subplot(211); stem(t1,f1); xlabel(kTs); ylabel(f(kTs)); title(欠采样的信号波形); [B,A]=butter(2,450/500); %设置低通滤波器参数 [H,w]=freqz(B,A,512,2000); fa=filter(B,A,f1); subplot(212); plot(fa) xlabel(t); ylabel(fa(t)); title(欠采样信号重构后的波形); t2=0.5*t;f2=sin(2*pi*t2*f); Figure,subplot(211); stem(t2,f2);xlabel(kTs); ylabel(f(kTs));title(临界采样的信号波形); [B,A]=butter(2,450/500); %设置低通滤波器参数 [H,w]=freqz(B,A,512,2000); fb=filter(B,A,f2);subplot(212); plot(fb),xlabel(t),ylabel(fb(t));title(临界采样信号重构后的波形); t3=0.2*t; f3=sin(2*pi*t3*f); figure subplot(2