数字信号处理实验指导书(15-16-2).docVIP

数字信号处理 实验讲义 二O一六年三月  目 录 实验一 离散时间信号的时域分析 3 实验二 离散时间系统的时域分析 6 实验三 离散时间信号的频域分析 9 实验四 线性时不变离散时间系统的频域分析 13 实验五 IIR数字滤波器的设计 17 实验六 FIR数字滤波器的设计 24 附录A MATLAB系统的常用概念 28 附录B 信号处理工具箱函数 32 前言 数字信号处理研究数字序列信号的表示方法，并对信号进行运算，以提取包含在其中的特殊信息。近几十年来，由于在研究及应用两方面均取得了进展，数字信号处理领域已日趋成熟。由于计算机的大量使用，从而很容易向学生提供实际环境，以验证所学的概念和算法。 本指导书编程语言是MATLAB，它广泛应用于高性能数值计算和可视化。本书假定读者已具备MATLAB基础知识。前面的一些实验帮助学生理解信号处理的重要概念，后面以设计性实验项目为主，有利于加强对原理的理解并且加强对技术的应用。附录中给出了本书中用到的MATALB函数及简要解释。 实验一 离散时间信号的时域分析 一．实验目的 熟悉MATLAB中产生信号和绘制信号的基本命令；熟悉序列的简单运算，如：加法、标量乘法、时间反转、延时、乘法等。 二．实验相关知识准备 1．用到的MATLAB命令 运算符号： ： . + - * / ; % sawtooth square 离散时间信号用数字序列x[n]来表示，常用的信号有单位冲激信号，单位阶跃信号，正弦信号，指数信号等 三．实验内容 1．离散时间序列的产生 (1) 单位抽样序列（单位冲激信号）的产生和绘制 %program 1 clf; % 清除图形窗口 n=-10:20; % 产生向量n，取值-10 - 20，间隔为1 u=[zeros(1,10) 1 zeros(1,20)]; % 生成单位冲激信号，前面10个0，后面20个0 stem(n,u); % 以n为横坐标，u为纵坐标画杆状图 xlabel(Time index n);ylabel(Amplitude);%定义坐标轴名称 title(Unit Sample Sequence); %标注图名 axis([-10 20 0 1.2]); %定义坐标轴范围 S=ones(1, N)(表示产生长度为N的一维行向量), 请采用与程序1相似的过程产生单位阶跃信号并画图。 指数信号的产生可以用命令 .^ 和 exp。 %program 2 %生成一个复指数信号 clf; % c=-(1/12)+(pi/6)*i; %生成一个复数 k=2;%复数 n=0:40; %横坐标（41个点） x=k*exp(c*n); %生成一个复指数信号； subplot(2,1,1); %子图1 stem(n,real(x)); %复数的实部图 xlabel(‘时间序n’)；%坐标轴 ylabel(‘振幅’);%坐标轴 title(‘实部’)；%标注图名 subplot(2,1,2); %子图2 stem(n,imag(x)); %虚部图 xlabel(‘时间序n’)； ylabel(‘振幅’); title(‘虚部’)； %program 3 %生成一个实指数信号 clf; n=0:35; a=1.2; %指数函数的底 k=0.2; %指数函数的系数 x=k*a.^n; %指数函数表达式 stem(n,x); %画出实指数函数的杆状图 xlabel(‘时间序n’)； ylabel(‘振幅’); (4) 正余弦信号的产生和绘制，可使用函数sin和cos %program 4 % 产生余弦信号 n=0:40; %横坐标（41个点） f=0.1; %余弦信号的频率 phase=0; %余弦信号的初相 A=1.5; %余弦信号的振幅 arg=2*pi*f*n-phase; %余弦信号的相位 x=A*cos(arg); %余弦信号 clf; %Clear old