数字信号处理实验教材正文.doc

实验一 离散时间信号与系统的实现与表示 由于MATLAB数值计算的特点，用它来分析离散时间信号与系统是很方便的。 一、离散时间信号 在MATLAB中，可以用一个列向量来表示一个有限长度的序列，例如： 在MATLAB中可以表示为： >>x = [2 1 0 2 3 -1 2 3]; 但是这种表示方法没有包含抽样的位置信息，要完全表示，要用和两个向量表示，例如： 在MATLAB中可以表示为： >> n=[-4 -3 -2 -1 0 1 2 3]; >> x = [2 1 0 2 3 -1 2 3]; 由于内存有限，MATLAB无法表示无限序列。 下面介绍几种典型离散信号的表示方法。 1、单位脉冲序列 （1）在MATLAB中zeros(1,N)函数来实现有限区间的，例如： >>x= zeros(1,N); >>x(1)=1; （2）我们还可以通过编写impseq函数来实现，该函数代码如下： function [x,n]= impseq(n0,n1,n2) %产生单位脉冲序列x(n)=δ(n-n0) %[x,n]= impseq(n0,n1,n2) % 1,n=n0 % x = % 1,n≠n0 %n=时间序列 n=n1:n2; x=[(n-n0)==0] 该函数产生如下信号，取值区间为，，得到序列 例1：输入命令： >>[x,n]= impseq(0,-5,5) 得到如下结果： x = 0 0 0 0 0 1 0 0 0 0 0 n = -5 -4 -3 -2 -1 0 1 2 3 4 5 2、单位阶跃序列 （1）在MATLAB中ones(1,N)函数来实现有限区间的，例如： >>x= ones(1,N); （2）我们还可以通过编写impseq函数来实现，该函数代码如下： function [x,n]= impseq(n0,n1,n2) %产生单位阶跃序列x(n)=u(n-n0);n1≤n0≤n2 %[x,n]= impseq(n0,n1,n2) % 1,n≥n0 % x = % 1,n<n0 %n=时间序列 n=n1:n2; x=[(n-n0)>=0] 该函数产生如下信号，取值区间为，，得到序列 例2：输入命令： >>[x,n]= impseq(0,-5,5) 得到如下结果： x = 0 0 0 0 0 1 1 1 1 1 1 n = -5 -4 -3 -2 -1 0 1 2 3 4 5 3、正弦序列 sin函数就可以产生正弦波，例如： >>n= 0:N-1; >>x=A*sin(n) 4、实指数序列 MATLAB实现： >>n= 0:N-1; >>x=a..^n 5、复指数序列 MATLAB实现： >>n= 0:N-1; >>x=exp((lu+j*w)*n); 6、随机序列 rand(1,N)：产生[0,1]上均匀分布的随机序列； randn(1,N)：产生均值为0，方差为1的高斯随机序列，即白噪声序列。 二、离散时间系统 离散时间系统的阶差分方程为： 在Matlab中，提供函数filter，用来求解给定输入和差分方程系数时，差分方程的数值解。 调用格式为： Y=filter(B,A,X) 其中：调用参数，分别为和各项差分项系数； 为输入信号序列，注意必须保证系数不为零。 输出序列与输入序列具有相同的长度。 返回参数为差分方程的数值解 例3：已知某线性时不变系统的差分方程如下： 当输入序列为时的输出结果。 MATLAB程序如下： N=41; a=[0.8 -0.44 0.36 0.22]; %输入序列X各项差分项系数 b=[1 0.7 -0.45 -0.6]; %输出序列Y各项差分项系数 x=[1 zeros(1,N-1)]; %单位脉冲序列 n=0:1:N-1; %时间序列 y=filter(a,b,x); %输入单位脉冲序列得到单位脉冲响应 stem(n,y);xlabel('n');ylabel('幅度'); 运行结果如下： 三、实验内容 1、编写程序用来产生下列基本序列，并绘制它们的波形图。 （1）单位脉冲序列，起点为0，终点为10