matlab数字信号处理实验报告.pdf

数字信号处理实验报告 基础实验篇 实验一 离散时间系统及离散卷积 一、实验原理 利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统 频率响应等的图像并于笔算结果进行比较，找出异同。编译合适程序能计 算取值范围不同的离散卷积。 二、实验目的 （1）熟悉MATLAB软件的使用方法。 （2）熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。 （3）利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲 响应。 三、实验步骤 （1）自编并调试实验程序，并且，给实验程序加注释； （2）按照实验内容完成笔算结果； （3）验证计算程序的正确性，记录实验结果。 （4）至少要求一个除参考实例以外的实验结果，在实验报告中，要描述清楚 实验结果对应的系统，并对实验结果进行解释说明。 四、实验源程序及实验结果 实验1-1 运行结果 xlabel ('n'); ylabel ('h (n)'); figure (2) [z,p,g] tf2zp (b,a); zplane (z,p) title ('零极点'); function [x,n] chongji (n1,n2) n [n1:n2]; x [n 0]; function shiyan1() a [1,-1,0.9]; b 1; x chongji (-20,120); n -20:120; h filter (b,a,x); figure (1) stem (n,h); title ('冲击响应'); 实验1-2 运行结果 b [0.0181,0.0543, 0.0543,0.0181]; a [1.000,-1.76, 1.1829,-0.2781]; w pi*freqspace (500); H freqz(b,a,w); MH abs(H); AH angle (H); subplot (2,1,1); plot (w/pi,MH); grid; axis([0,1,0,1]); xlabel ('w (pi)'); ylabel ('|H |'); title ('幅度、相位响应'); subplot (2,1,2); plot (w/pi,AH); grid; xlabel ('w (pi)'); ylabel ('angle (H)'); 实验1-3 运行结果 n 0:30; %输入x (n)和冲激响应h (n) x zeros(1,length (n)); h zeros(1,length (n)); x ([find ((n> 0)&(n< 4))]) 1; h ([find ((n> 0)&(n< 8))]) 0. 5; figure (1) subplot (3,1,1); stem (n,x); axis([0,30,0,2]); title ('输入序列'); xlabel ('n'); ylabel ('x (n)'); subplot (3,1,2); stem (n,h); axis([0,30,0,2]); title ('冲激响应序列'); xlabel ('n'); ylabel ('h (n)'); %输出响应 y conv (x,h); subplot (3,1,3); n 0:length (y)-1; stem (n,y); title ('输出响应'); xlabel ('n'); ylabel ('y (n)'); 实验二 离散傅立叶变换与快速傅立叶变换 一、 实验原理 对有限长序列使用离散Fouier变换(DFT)可以很好的反映序列的频谱特性， 而且易于用快速算法在计算机上实现，当序列x(n)的长度为N 时，它的DFT定 义为 N1