基于MATLA数字信号处理实例分析.docVIP

湖北文理学院理工学院 学生结业论文 课程名称：MATLAB教程 结业论文名称：基于MATLAB的数字信号处理实例分析 专业名称：通信工程 班级：1011 学号： 学生姓名： 赵彦彦 教师姓名：李敏 2013年1月6日 基于MATLAB的数字信号处理实例分析 摘要 随着信息科学和计算技术的迅速发展，在人们的日常生活中，对信号的处理显得尤为重要，而计算机不能直接对模拟信号进行处理，使得人们对数字信号处理理论的认知与了解要求更为深入。由于计算机解决复杂的数字信号系统有一定的困难，而MATLAB的出现，解决了这一难题。MATLAB提供了用于数值运算和信号处理的数学计算软件包，同时可以实现系统级的通信系统设计与仿真。随着版本的不断升级，不同应用领域的专用库函数和模块汇集起来作为工具箱添加到软件包中，其功能越来越强大。本文是基于MATLAB的数字信号处理实例分析，主要介绍了用MATLAB对系统函数零点、极点分布图以及模拟周期信号的频谱分析（模拟信号x(t)等间隔T采样后x(nT)的N点DFT）。 关键字：MATLAB 数字信号 系统函数 频谱 系统函数零点、极点分布图 通过学习信号与系统、数字信号处理，掌握了传输函数和系统函数等，本文仅对系统函数X(z)零点和极点分布进行分析。 利用下面的程序段，观察系统函数X(z)零点和极点分布的特点 程序段如下：n=[1 0 0];m=[1 0.6 -0.16]; zplane(n,m); 执行结果如图： 改变系统函数X(z)，观察与上图的差异 程序段如下： n=[0 1 0 0];m=[1 -1 2 3]; zplane(n,m); 执行结果如下图： 程序及系统函数X(z)零点极点分布分析：通过以上两个程序，了解到该程序主要是对zplane(n,m)　绘制出系统函数H(z)的零极点图zplane(n,m)指令，绘制出系统函数H(z)的零极点图 N=16;n=0:N-1; %FFT的变换区间N=16 xnT=sin(2*pi*n*T)+cos(2*pi*n*T)+tan(2*pi*n*T); %对x(t)16点采样 Xk16=fft(xnT); %计算xnT的16点DFT Xk16=fftshift(Xk16); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心） subplot(3,1,1);stem(fk,abs(Xk16),.);box on %绘制8点DFT的幅频特性图 title((a) 16点|DFT[x(nT)]|);xlabel(f(Hz));ylabel(幅度); axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(Xk16))]) N=32;n=0:N-1; %FFT的变换区间N=16 xnT=sin(2*pi*n*T)+cos(2*pi*n*T)+tan(2*pi*n*T); %对x(t)32点采样 Xk32=fft(xnT); %计算xnT的32点DFT Xk32=fftshift(Xk32); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心） subplot(3,1,2);stem(fk,abs(Xk32),.);box on %绘制8点DFT的幅频特性图 title((b) 32点|DFT[x(nT)]|);xlabel(f(Hz));ylabel(幅度); axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(Xk32))]) N=64;n=0:N-1; %FFT的变换区间N=16 xnT=sin(2*pi*n*T)+cos(2*pi*n*T)+tan(2*pi*n*T); %对x(t)8点采样 Xk64=fft(xnT); %计算xnT的8点DFT Xk64=fftshift(Xk64); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心） subplot(3,1,3);stem(fk,abs(Xk64),.); box on%绘