实验一 序列的产生、运算和Z反变换.docVIP

实验一 序列的产生、运算和z反变换 （2课时） 实验目的 学习用MATLAB进行数字信号处理的基本方法； 掌握用MATLAB生成常用序列的方法； 掌握用MATLAB进行序列的运算和z反变换的方法。 实验设备 计算机一台，安装了MATLAB软件 实验原理和基本操作 常用序列的生成 （1）、单位脉冲序列δ(n)及迟延的单位脉冲序列δ(n-np) 用MATLAB编写的用来生成单位脉冲序列及其移位序列的函数为impseq.m，其输入参数为序号起始值ns，序号终止值nf及脉冲位置np。函数主要语句如下： function [x,n] = impseq (np,ns,nf) if ((np ns) | (np nf) | (ns nf)) % 检查输入参数正确性 error(参数必须满足 ns = np = nf) end n = [ns:nf]; % 生成位置向量 x = [(n-np) == 0];’r’) 例1 》[x,n]=impseq(2,-1,5) function [x,n] = stepseq (np,ns,nf) if ((np ns) | (ns nf) | (np nf)) % 检查输入参数正确性 error(参数必须满足 ns = np = nf) end n = [ns:nf]; % 生成位置向量 x = [(n-np) = 0]; % 生成阶跃脉冲序列 stem(n,x,’r’) 例2 》[x,n]=stepseq(0,-2,4) 运行结果： （3）、矩形序列 用MATLAB编写的生成矩形序列的程序如下：其中N是矩形序列的长度. function [x,n] = rn(N) %n = [0:N]; % 生成位置向量 %x = [n= 0]; n=[-N:2*N-1]; x=[n=0 n=N-1]; stem(n,x,’r’) 例3：[x,n]=rn(4); 结果： 序列的运算 （1）、序列的移位 y(n)=x(n-k). 完全地表示序列x(n)要用x和n两个向量，n是位置向量（序号向量）。 原始的序列x(n)用x和nx（nx是位置向量）表示，移位后的序列用y和ny（ny为位置向量）表示。移位不影响向量x的值，因此y=x。移位体现为位置向量的改变，nx的每个元素都加了一个k后成为ny，即ny=nx+k。把移位函数命名为seqshift.m，可以写出： function [y,ny] = seqshift(x,nx,k) ny = nx + k; % 位置向量移位 y = x; % 序列的值不变 例4 》x=[1,2,3,4,5] 》nx=[-1,0,1,2,3] 》[y,ny] =seqshift(x,nx,2) 运行结果： （2）、序列的翻褶 y(n)=x(-n). 在MATLAB中，这一运算由两次调用fliplr函数来实现：y=fliplr(x)的作用是把行向量x中元素排列的次序左右翻转；ny=-fliplr(nx)既将序号向量左右翻转，又改变其各元素的正负号，这样就得到了翻褶后的向量y和ny，写成函数seqfold.m，其主要语句如下： function [y,n] = seqfold(x,nx) y = fliplr(x); % 将序列数值左右翻转 ny = -fliplr(nx); % 将序列位置对零位置左右翻转，故同时改变正负号 例5 》x=[1,2,3,4,5] 》nx=[-1,0,1,2,3] 》[y,ny] =seqfold(x,nx) 运行结果： (3)、序列的和与积 序列的加法： y(n)=x1(n)+x2(n)；序列的乘法： y(n)=x1(n) .*x2(n)； 序列的加法可用算术运算符“+”实现，序列相乘要用点乘运算符“.*”实现。加法和乘法要求x1和x2有同样的位置向量n。先要把x1和x2的位置向量n1和n2扩展成包括两者最左和最右端点的公共位置向量n，显然min(n)应是min(n1)和min(n2)中的小者，而max(n)应是max(n1)和max(n2)中的大者。然后把x1和x2放到扩展后的位置向量n上，因为n可能比n1和n2长，必须把x1和x2的长度加大，把加长后的x1和x2分别命名为y1和y2。这需要把x1和x2放到位置向量n的对应位置上，并把多余的元素置零。程序中先把y1和y2全置零，然后用find函数来寻找对应的下标区间。因为find函数的输出是满足条件的输