数字信号处理a(双语)实验报告

数字信号处理a(双语)实验报告

ID:17683876

大小:819.35 KB

页数:23页

时间:2018-09-04

数字信号处理a(双语)实验报告_第1页
数字信号处理a(双语)实验报告_第2页
数字信号处理a(双语)实验报告_第3页
数字信号处理a(双语)实验报告_第4页
数字信号处理a(双语)实验报告_第5页
资源描述:

《数字信号处理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

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。