数字信号处理实验报告95152new

《数字信号处理实验报告95152new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字信号处理课程设计报告通信11级16班学号：52111604姓名：薛焕福实验一离散时间系统及离散卷积一、实验目的1、熟悉MATLAB软件的使用方法。2、熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。3、利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲响应。4、熟悉离散卷积的概念，并利用MATLAB计算离散卷积。二、实验原理本实验利用MATLAB的库函数，实验冲击响应零极点、幅度相位谱的绘制。熟悉MATLAB的使用。并利用离散卷积的概念编写一个计算离散卷积的程序。三、实验内容及要求1.离散时间系统的单位冲击响应（1）选择一个离散时间系

2、统；（2）用笔进行差分方程的递推计算；（3）编制差分方程的递推计算程序；（4）在计算机上实现递推运算；（5）将程序计算结果与笔算的计算结果进行比较，验证程序运行的正确性；2.离散系统的幅频、相频的分析方法（1）离散系统的幅频、相频的分析方法（2）用笔计算几个特殊的幅频、相频的值，画出示意曲线图；（3）编制离散系统的幅频、相频的分析程序；（4）在计算机上进行离散系统的幅频、相频特性计算，并画出曲线；（5）通过比较，验证程序的正确性；3.离散卷积的计算（1）选择两个有限长序列，用笔计算其线性卷积；（2）编制有限长序列线性卷积程序；（3）利用计算程序对（1）选择的有限长序

3、列进行卷积运算；（4）比较结果验证程序的正确性。四、实验过程1、离散时间系统的单位脉冲响应程序:a=[1,-1,0.91];%系统的系数b=[1,2,3];n0=0;n=-20:120;x=[(n-n0)==0];%利用filter函数进行滤波，从而求得系统的单位冲击响应h=filter(b,a,x);figure(1)stem(n,h);%利用stem绘制曲线图title('冲击响应');xlabel('n');ylabel('h(n)');hfigure(2)[z,p,g]=tf2zp(b,a);%求传递函数zplane(z,p)%求零极点（2）实验结果：2、离

4、散系统的幅频、相频的分析程序：b=[0.6,1.2,0.4,0.9];%分子的系数a=[1.0,-1.7,1.1,-0.5];%分母的系数m=0:length(b)-1;l=0:length(a)-1;K=500;k=1:K;w=pi*k/K;H=(b*exp(-j*m'*w))./(a*exp(-j*l'*w));%求传递函数magH=abs(H);%求传递函数的幅值angH=angle(H);%求传递函数的相位figure(1)subplot(2,1,1);plot(w/pi,magH);%定义坐标轴的区间grid;axis([0,1,0,1]);xlabel(

5、'w(pi)');ylabel('

6、H

7、');title('·幅度，相位响应');subplot(2,1,2);plot(w/pi,angH);grid;xlabel('w(pi)');ylabel('angle(H)');（2）实验结果：3、离散卷积的计算程序：functionpr13()n=-5:50;u1=stepseq(0,-5,50);u2=stepseq(10,-5,50);x=u1-u2;h=((0.9).^n).*u1;figure(1)subplot(3,1,1);stem(n,x);axis([-5,50,0,2]);title('输入序列');

8、xlabel('n');ylabel('x(n)');subplot(3,1,2);stem(n,h);axis([-5,50,0,2]);title('冲击响应序列');xlabel('n');ylabel('h(n)');[y,ny]=conv_m(x,n,h,n);subplot(3,1,3);stem(ny,y);axis([-5,50,0,8]);title('输出响应');xlabel('n');ylabel('y(n)');function[x,n]=stepseq(n0,n1,n2)n=n1:n2;x=[(n-n0)>=0];function[y,n

9、y]=conv_m(x,nx,h,nh)nyb=nx(1)+nh(1);nye=nx(length(x))+nh(length(h));ny=[nyb:nye];y=conv(x,h);（2）实验结果：实验二离散傅立叶变换与快速傅立叶变换一、实验目的：1、加深理解离散傅立叶变换及快速傅立叶变换概念；2、学会应用FFT对典型信号进行频谱分析的方法；3、研究如何利用FFT程序分析确定性时间连续信号；4、熟悉应用FFT实现两个序列的线性卷积的方法。二、实验原理在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以使用离散Fouier变换(DFT)。