资源描述：
《数字信号处理报告—实习报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、数字信号处理报告中国地质大学（武汉）学校:姓名:班级:学号:指导老师:目录1-1:离散卷积的运算1-2：离散傅里叶变换及其应用1-3：IIR滤波器的设计1-4：FIR低通滤波器的设计1-5：用MATLAB设计一信号处理系统1-6:频率抽选1-7：实习心得体会KL离算卷积的应用1：基本要求设线性时不变(LTI)系统的冲激响应为h(n),输入序列为x(n)1>h(n)=(0.8)n,0WnW4；x(n)=u(n)-u(n-4)2、h(n)=(0.8)nu(n),x(n)=u(n)-u(n-4)3、h(n)=(0.8)nu(n),x(n)=u(n)求以上三种情况下系统的输出y(n),显示输入
2、和输出波形。2：实验目的及算法原理实验目的：(1)：掌握离散卷积计算机实现。(2)：进一步对离散信号卷积算法的理解。算法原理：把冲激响应h(n)与输入序列x(n)分别输入到程序川，然后调用离散卷积函数y=conv(x.,h)即可得到所要求的结果。8离散卷积定义为y{n)=工x(k)h(n_k)k=_8当序列为有限长时则y(n)=工x伙)/?(/?-k)k=03：程序代码xl=[l111];nxl=0:3;hl=[l0.80.640.8^30.8^4];nhl=0:4;yl二conv(xl,hl);subplot(3,3,1);stem(nxl,xl);title('序列xl');xla
3、bcl('n');ylabcl('xl(n)');subplot(3,3,2);stem(nhl,hl);title('序列hl');xlabel('n');ylabel('hl(n)，);subplot(3,3,3);stem(yl);title(，序列yl');xlabel('n');ylabel('yl(n)');x2=[l11l];nx2=0:3;nh2=0:l:20;h2=(0.8).八nh2;y2=conv(x2,h2);subplot(3,3,4);stem(nx2,x2);title(，序列x2');xlabcl('n');ylabcl('x2(n)');subplot
4、(3,3,5);stem(nh2,h2);title('序列h2');xlabel('n');ylabel('h2(n)，);subplot(3,3,6);stem(y2);title('序列y2');xlabel('n');ylabel('y2(n)，);nx3=0:1:20;x3=l.nx3;nh3=0:l:20;h3=(0.8).nh3;y3二conv(x3,h3);subplot(3,3,7);stem(nx3,x3);title(，序列x3');xlabcl('n');ylabcl('x3(n)');subplot(3,3,8);stem(nh3,h3);title('序列h
5、3,);xlabel('n');ylabel('h3(n)，);subplot(3,3,9);stem(y3);title('序列y3');xlabel('n);ylabel('y3(n)');4：结果分析序列x1Tn序列hinT序列y14iolTlIllk0510序列x2024序列h202040n序列y202040nn序列x3序列h3nn序列y30501-2离散傅里叶变换及应用1：基本要求设有离散序列x(n)=cos(0.48兀n)+cos(0.52nn)分析下列三种情况下的幅频特性。(1)采集数据长度N=16,分析16点的频谱，并画出幅频特性。采集数据长度N=16,并补零到64点，
6、分析其频谱，并画出幅频特性。(2)采集数据长度N=64,分析64点的频谱，并画出幅频特性。观察三幅不同的幅频特性图，分析和比较它们的待点及形成原因。2：实习目的及相关算法实习目的：(1)了解DFT及FFT的性质和特点。(2)利用FFT算法计算信号的频谱。相关算法：读入离散序列x(n)—A采集长度为N=16的数据，调用matlab中的函数fft(x,16)与fft(x,64)对其作离散傅里叶变换得到16点、64点的频谱一＞采集数据长度为N二64,,调用matlab中的函数fft(x,46)对其作离散傅里叶变换得到64点的频谱。3：相关实现代码n=0:l:15;nl二n;xl=cos(0.
7、48*pi*n)+cos(0.52*pi*n);gl=abs(fft(xl,16));subplot(3,2,l);stem(xl);title('xl');subplot(3,2,2);stem(gl);title(,gl,);n2二n;x2=cos(0.48*pi*n2)+cos(0.52*pi*n2);x2=[x2zeros(l,48)];g2=abs(fft(x2,64));subplot(323);stem(x2);title('x2
