《数字信号处理A双语实验报告.docVIP

实 验 报 告 实验名称 采样、系统性质及滤波 数字系统频率响应和样本处理算法实现 加窗和离散傅氏变换 数字滤波器设计 课程名称 数字信号处理A（双语） 班级学号 姓 名 周士亮 开课时间 2011/2012学年，第二学期 实验一：采样、系统性质及滤波 实验目的和任务 （1）熟悉MATLAB的主要操作命令。 （2）学会简单的矩阵输入和数据读写。 （3）掌握简单的绘图命令。 （4）用MATLAB编程并学会创建函数。 （5）观察采样引起的混叠。 （6）判别离散时间系统的时不变性。 （7）卷积计算 实验内容 A、观察采样引起的混叠 设模拟信号为，t的单位为毫秒(ms)。 设采样频率为3kHz，确定与混叠的采样重建信号。 画出和在范围内的连续波形。（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。） 分别用和在两信号波形上标记出3kHz采样点。 两信号波形是否相同？采样后的两序列是否相同？ 答：通过实验结果可以看出，两信号和的波形明显不同，但采样值是相同的。采样后的两序列和是相同。 B、判别离散时间系统的时不变性。 设输入序列为，系统实现对的抽取。 设。取延迟量D（例如D＝30）。记，画出、的序列波形。 编程求出系统对的响应以及对的响应 画出、的波形。 该系统是否为时不变的？ 答：根据系统时不变的性质，当时，有=。 从实验结果的图中可以看出，与的波形并不一致，即不重合。所以该系统为时变系统。 C、利用卷积计算出输入信号通过FIR滤波器的输出，并观察输出信号的input-on暂态、input-off暂态和稳态部分。（计算卷积可用conv命令） 考虑下面两个滤波器，第一个的单位脉冲响应为，另一个的单位脉冲响应为；输入为周期方波，在一个周期内。 1. 分别画出两个滤波器的输出的波形，并与书上p144例4.1.8的两幅图比较是否一致。 与课本上的图是一致的。 2. 计算出图中稳态部分的值。 答：第一个滤波器input-on的稳态值为 1, input-off的稳态值为0。 第二个滤波器input-on和input-off的稳态值均为 0。 三 、实验心得 本次实验首先熟悉了Matlab的运行环境,掌握了matlab的基本使用语句并学会了编写程序段来画图。观察到了采样引起的混叠，加深了对于相关知识的理解。并且使用matlab语句实现了卷积运算。 MATLAB实验程序 % ============= % problem 1 % ============= clear % estimate x(t) and xa(t) with a much higher sampling freq. fs1 time_period=6; % unit: ms fs1=50; % unit: kHz T1=1/fs1; % unit: ms n1=0:fix(time_period/T1); x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*sin(3*pi*n1*T1); xa=cos(pi*n1*T1); % obtain x(nT) and xa(nT) with given sampling freq. fs fs=3; T=1/fs; n=0:fix(time_period/T); x_sample=cos(5*pi*n*T)+4*sin(2*pi*n*T).*sin(3*pi*n*T); xa_sample=cos(pi*n*T); figure,plot(n1*T1,x,r,n1*T1,xa,b,n*T,x_sample,ro), hold on, stem(n*T,xa_sample,b:x) legend(x(t),xa(t),x(nT),xa(nT)),xlabel(t(ms)) % ============= % problem 2 % ============= clear % plot x(n) and x(n-D) D=30; N=500; n=1:N; x=sin(2*pi/100*n); for n=1:N+D, if (n-D)=0, xD(n)=0; else xD(n)=x(n-D); end