《数字信号处理实验指导书.docVIP

《数字信号处理》实验总结 《数字信号处理》实验指导书 温浩 吕秋霞 李继容编 五邑大学 第1章 MATLAB软件简介 MATLAB还精心设计了演示程序系统（Demo），内容包括MATLAB的内部主要函数和各个工具箱（Toolbox）的使用。一是在命令窗口中输入demo或demos命令；二是选择“help”→“Demos”菜单命令。 MATLAB软件的符号运算 各个变量名须用空格隔开。例如，定义x、y、z三个符号变量的语句格式为： syms x y z 另一种定义符号变量的语句格式为：sym(变量名)。例如，x、y、z三个符号变量定义的语句格式为：x=sym(x); y=sym(y); z=sym(z); sym语句还可以用来定义符号表达式，如，定义表达式x+1为符号表达式对象，语句为：sym(x+1);。 在符号运算中，可以用“simple”或者“simplify”函数来化简运算结果，例如： syms x; f1=sin(x)^2; f2=cos(x)^2; y=f1+f2 结果：y=sin(x)^2+cos(x)^2 y=simplify(y) 结果：y=1 ezplot：绘制符号表达式的曲线。ezplot(y, [a,b])，[a，b]表自变量取值范围，默认值为。 用syms x; y=-16*x^2+64*x+96;表示函数。 MATLAB为实现交互控制程序流程还提供了continue、break、pause、input、error、disp等命令。 stem函数：默认是空心的。 “fill”、“filled”，或者参数“.”。 第2章 离散时间信号的表示。 用function y=uDT(n); y=(n=0);表示单位阶跃序列。 矩形序列定义为：。 单边指数序列定义为：。 正弦序列定义为：。只有当为有理数时，正弦序列才具有周期性。 复指数序列定义为：。 欧拉公式：。 画出复指数序列的实部、虚部、模及相角随时间变化的曲线，并观察其时域特性。 n=0:30; A=2;a=-1/10;b=pi/6; x=A*exp((a+i*b)*n); subplot(2,2,1) stem(n,real(x),fill),grid on; title(实部),axis([0,30,-2,2]),xlabel(n) subplot(2,2,2) stem(n,imag(x),fill),grid on; title(虚部),axis([0,30,-2,2]) ,xlabel(n) subplot(2,2,3) stem(n,abs(x),fill),grid on; title(模),axis([0,30,0,2]) ,xlabel(n) subplot(2,2,4) stem(n,angle(x),fill),grid on; title(相角),axis([0,30,-4,4]) ,xlabel(n) 第3章 离散时间LTI系统的时域分析 离散时间LTI系统可用线性常系数差分方程来描述，即 可用y=filter(b,a,x)表示上式。其中，x为输入的离散序列；y为输出的离散序列；y的长度与x的长度一样。 求经过的零状态响应。 a=[3 -4 2]; b=[1 2]; n=0:30; x=(1/2).^n; y=filter(b,a,x); stem(n,y,fill),grid on; xlabel(n),title(系统响应y(n)) y=conv(x,h)。x与h表示离散时间信号值的向量。 第4章 z变换及离散时间LTI系统的z域 z变换： 单边z变换[单边：Z=ztrans(x), x=iztrans(z)]： 上式中的x和Z分别为时域表达式和z域表达式的符号表示，可通过sym函数来定义。 如：求的Z变换 x=sym(a^n*cos(pi*n)); Z=ztrans(x); simplify(Z) 结果：z/(z+a) 如果信号的z域表示式是有理函数，进行z反变换的另一个方法是对进行部分分式展开，然后求各简单分式的z反变换。设的有理分式表示为 [R,P,K]=residuez(B,A)：对进行部分分式展开，B，A分别表示X(z)的分子与分母多项式的系数向量；R为部分分式的系数向量；P为极点向量；K为多项式的系数。若X(z)为有理真分式，则K为零。 对函数进行部分分式展开，并求出其z反变换。 B = [18]; A = [18,3,-4,-1]; [R,P,K] = residuez(B,A) 结果： R= 0.3600 0.2400 0.4000 P= 0.5000 -0.3333 -0.3333 K= [] ，表示系统有一个二重极点。所以，X(z)的部分分式