资源描述：
《数字信号处理a(双语)实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验报告实验名称采样、系统性质及滤波数字系统频率响应和样本处理算法实现加窗和离散傅氏变换数字滤波器设计课程名称数字信号处理A（双语）班级学号B09010622姓名周士亮开课时间2011/2012学年，第二学期实验一：采样、系统性质及滤波一、实验目的和任务（1）熟悉MATLAB的主要操作命令。（2）学会简单的矩阵输入和数据读写。（3）掌握简单的绘图命令。（4）用MATLAB编程并学会创建函数。（5）观察采样引起的混叠。（6）判别离散时间系统的时不变性。（7）卷积计算二、实验内容A、观察采样引起的混叠设模拟信号为，t的单位为毫秒(ms)。1.设采样频率为3k
2、Hz，确定与混叠的采样重建信号。2.画出和在范围内的连续波形。（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。）3.分别用和在两信号波形上标记出3kHz采样点。两信号波形是否相同？采样后的两序列是否相同？答：通过实验结果可以看出，两信号和的波形明显不同，但采样值是相同的。采样后的两序列和是相同。B、判别离散时间系统的时不变性。设输入序列为，系统实现对的抽取。1.设。取延迟量D（例如D＝30）。记，画出、的序列波形。2.编程求出系统对的响应以及对的响应3.画出、的波形。该系统是否为时不变的？答：根据系统时不变的性质，当时，有=。从实验结果
3、的图中可以看出，与的波形并不一致，即不重合。所以该系统为时变系统。C、利用卷积计算出输入信号通过FIR滤波器的输出，并观察输出信号的input-on暂态、input-off暂态和稳态部分。（计算卷积可用conv命令）考虑下面两个滤波器，第一个的单位脉冲响应为，另一个的单位脉冲响应为；输入为周期方波，在一个周期内。1.分别画出两个滤波器的输出的波形，并与书上p144例4.1.8的两幅图比较是否一致。与课本上的图是一致的。2.计算出图中稳态部分的值。答：第一个滤波器input-on的稳态值为1,input-off的稳态值为0。第二个滤波器input-on和i
4、nput-off的稳态值均为0。三、实验心得本次实验首先熟悉了Matlab的运行环境,掌握了matlab的基本使用语句并学会了编写程序段来画图。观察到了采样引起的混叠，加深了对于相关知识的理解。并且使用matlab语句实现了卷积运算。MATLAB实验程序%=============%problem1%=============clear%estimatex(t)andxa(t)withamuchhighersamplingfreq.'fs1'time_period=6;%unit:msfs1=50;%unit:kHzT1=1/fs1;%unit:msn1
5、=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);%obtainx(nT)andxa(nT)withgivensamplingfreq.'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
6、,xa,'b',n*T,x_sample,'ro'),holdon,stem(n*T,xa_sample,'b:x')legend('x(t)','xa(t)','x(nT)','xa(nT)'),xlabel('t(ms)')%=============%problem2%=============clear%plotx(n)andx(n-D)D=30;N=500;n=1:N;x=sin(2*pi/100*n);forn=1:N+D,if(n-D)<=0,xD(n)=0;elsexD(n)=x(n-D);endendfigure,subplot(2,1
7、,1),plot(1:N,x,'r:',1:length(xD),xD,'b'),legend('x(n)','xD(n)'),,xlabel('n')%ploty(n)andyD(n)forn=1:fix(N/2)y(n)=x(2*n);endforn=1:length(y)+D,if(n-D)<=0,y_delay(n)=0;elsey_delay(n)=y(n-D);endendforn=1:fix(length(xD)/2)yD(n)=xD(2*n);endsubplot(2,1,2),plot(1:length(y),y,'r:',1:leng
8、th(y_delay),y_delay,'r.:',1:length(yD),y
