数字信号处理实验指导手册
数字信号处理实验指导手册【模板】
数字信号处理实验指导手册西安文理学院机械电子工程系目录实验一离散时间信号 (2)实验二时域采样定理 (7)实验三离散时间系统 (10)实验四线性卷积与圆周卷积 (13)实验五用FFT作谱分析 (16)实验六用双线性变换法设计IIR数字滤波器 (18)实验七 FIR滤波器设计 (20)实验一离散时间信号【实验目的】用MATLAB实现离散时间信号的表示和运算，掌握MATLAB的基本命令和编程方法，为后续实验打基础。
【实验原理】在数字信号处理中，所有的信号都是离散时间信号，因此应首先解决在MATLAB中如何表示离散信号。
设一模拟信号经A/D变换后，得到序列信号由于MATLAB对下标的约定为从1开始递增，因此要表示，一般应采用两个矢量，如：这表示了一个含9个采样点的矢量：【实验内容】熟悉下面序列（信号）的产生方法及相关运算1、单位采样序列2、单位阶跃序列3、信号翻转4、信号相加5、信号折叠6、信号移位【参考程序】单位采样序列1、impluse1.m （图1-1）n=10;x=zeros(1,n);x(1)=1;plot(x,'*');2、impluse2.m（图1-2）n=-5:5;x=[n==0];stem(x,'*');3、impluse3.m（图1-3）n=1:10;n0=3;x=[(n-n0)==1];plot(x,'*');单位阶跃序列1、steps1.m（图1-4）n=10;x=ones(1,n);plot(x,'*');2、steps2.m（图1-5）n=10;x=ones(1,n);x(1)=0;x(2)=0;plot(x,'*');3、steps3.m（图1-6）n=0:10;n0=4;x=[(n-n0)>=0];plot(x,'*');信号翻转（中心对称）flip.m（图1-7）x=0:6;figure(1);subplot(211);plot(x,'*');y=fliplr(x);subplot(212);plot(y,'*');信号相加Showsigadd.m（图1-8）clear allfigure(1);[x1 n1]=impseq(-4,-5,5); subplot(3,1,1);stem(n1,x1,'*');title('信号相加演示');ylabel('x1');[x2 n2]=stepseq(0,-3,7); subplot(3,1,2);stem(n2,x2,'*');ylabel('x2');[y n]=sigadd(x1,n1,x2,n2); subplot(3,1,3);stem(n,y,'*');ylabel('y=x1+x2)')信号折叠%showfold.m（图1-9）clear all[x n]=impseq(2,-5,3); subplot(2,1,1);stem(n,x);title('信号折叠演示'); [y,n]=sigfold(x,n); subplot(2,1,2);stem(n,y);信号移位%showsigshift.m（图1-10）figure(1)clear all[x,n]=stepseq(0,-5,5);subplot(211);stem(n,x,'*');title('信号移位演示’);[y,n]=sigshift(x,n,3);subplot(212)stem(n,y);【思考题】1、用MATLAB实现信号的相加、移位、反转运算。
数字信号处理实验指导书(M)
数字信号处理实验电子信息科学与技术实验室2007年7月目录实验一离散时间信号的时域表示 (3)实验二离散信号的卷积和 (6)实验三离散傅立叶变换及其特性验证 (8)实验四信号处理中FFT的应用 (11)实验五离散系统的Z域分析 (15)实验六无限冲激响应(IIR)数字滤波器的三种结构 (19)实验七冲激响应不变法IIR数字滤波器设计 (23)实验八双线性变换法IIR数字滤波器设计 (26)实验一 离散时间信号的时域表示一、实验目的1、熟悉Matlab 命令，掌握离散时间信号－序列的时域表示方法。
2、掌握用Matlab 描绘二维图像的方法。
3、掌握用Matlab 对序列进行基本的运算和时域变换的方法。
二、实验原理与计算方法（一）序列的表示方法 序列的表示方法有列举法、解析法和图形法，相应的用Matlab 也可以有这样几种表示方法，分别介绍如下：1、列举法 在Matlab 中，用一个列向量来表示一个有限长序列，由于一个列向量并不包含位置信息，因此需要用表示位置的n 和表示量值的x 两个向量来表示任意一个序列，如：例1.1：>>n=[-3,-2,-1,0,1,2,3,4]; >>x=[2,1,-1,0,1,4,3,7];如果不对向量的位置进行定义，则Matlab 默认该序列的起始位置为n=0。
由于内存有限，Matlab 不能表示一个无限序列。
2、解析法对于有解析表达式的确定信号，首先定义序列的范围即n 的值，然后直接写出该序列的表达式，如：例1.2：实现实指数序列nn x )9.0()(=，100≤≤n 的Matlab 程序为：>>n=[0:10]; >>x=(0.9).^n;例 1.3：实现正余弦序列)5.0sin(2)31.0cos(3)(n n n x πππ++=，155≤≤n 的Matlab 程序为：>>n=[5:15];>>x=3*cos(0.1*pi*n+pi/3)+2*sin(0.5*pi*n); 3、图形法在Matlab 中用图形法表示一个序列，是在前两种表示方法的基础上将序列的各个量值描绘出来，即首先对序列进行定义，然后用相应的画图语句画图，如：例1.4：绘制在例1.1中用列举法表示的序列的图形，则在向量定义之后加如下相应的绘图语句：>>stem(n,x);此时得到的图形的横坐标范围由向量n 的值决定，为-3到4，纵坐标的范围由向量x 的值决定，为-1到7。
《数字信号处理》上机实验指导书
《数字信号处理》上机实验指导书实验1 离散时间信号的产生1．实验目的数字信号处理系统中的信号都是以离散时间形态存在，所以对离散时间信号的研究是数字信号处理的基本所在。
而要研究离散时间信号，首先需要产生出各种离散时间信号。
MATLAB 是一套功能强大的工程计算及数据处理软件，广泛应用于工业，电子，医疗和建筑等众多领域。
使用MATLAB软件可以很方便地产生各种常见的离散时间信号，而且它还具有强大的绘图功能，便于用户直观地输出处理结果。
通过本实验，学生将学习如何用MATLAB产生一些常见的离散时间信号，并通过MATLAB中的绘图工具对产生的信号进行观察，加深对常用离散信号的理解。
2．实验要求本实验要求学生运用MATLAB编程产生一些基本的离散时间信号，并通过MATLAB的几种绘图指令画出这些图形，以加深对相关教学内容的理解，同时也通过这些简单的函数练习了MATLAB的使用。
3．实验原理（1）常见的离散时间信号1）单位抽样序列，或称为离散时间冲激，单位冲激：?(n)???1?0n?0 n?0如果?(n)在时间轴上延迟了k个单位，得到?(n?k)即：?1n?k ?(n?k)??0n?0?2）单位阶跃序列n?0?1 u(n)?n?0?0在MATLAB中可以利用ones( )函数实现。
x?ones(1,N);3）正弦序列x(n)?Acos(?0n??)这里，A,?0,和?都是实数，它们分别称为本正弦信号x(n)的振幅，角频率和初始相位。
f0??02?为频率。
x(n)?ej?n4）复正弦序列5）实指数序列x(n)?A?n（2）MATLAB编程介绍MATLAB是一套功能强大，但使用方便的工程计算及数据处理软件。
其编程风格很简洁，没有太多的语法限制，所以使用起来非常方便，尤其对初学者来说，可以避免去阅读大量的指令系统，以便很快上手编程。
值得注意得就是，MATLAB中把所有参与处理的数据都视为矩阵，并且其函数众多，希望同学注意查看帮助，经过一段时间的训练就会慢慢熟练使用本软件了。
数字信号处理实验指导书(需印刷)
实验一 常见离散信号的MATLAB 产生和图形显示1、实验目的：（１）熟悉MATLAB 应用环境，常用窗口的功能和使用方法。
（２）加深对常用离散时间信号的理解。
（３）掌握简单的绘图命令。
（４）掌握线性卷积的计算机编程方法。
2、实验原理：（１）单位抽样序列⎩⎨⎧=01)(n δ0≠=n n如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即：⎩⎨⎧=-01)(k n δ≠=n k n（２）单位阶跃序列⎩⎨⎧=01)(n u 00<≥n n（３）矩形序列 ⎩⎨⎧=01)(n R N 其他10-≤≤N n（４）正弦序列)sin()(ϕ+=wn A n x（５）复正弦序列jwnen x =)(（６）指数序列na n x =)(（７）线性时不变系统的响应为如下的卷积计算式：∑∞-∞=-=*=m m n h m x n h n x n y )()()()()(3、实验内容及步骤：（１）复习常用离散时间信号的有关内容。
（２）编制程序产生上述6种序列（长度可输入确定，对(4) (5) (6)中的参数可自行选择），并绘出其图形。
（３）已知系统的单位脉冲响应),(9.0)(n u n h n=输入信号)()(10n R n x =，试用卷积法求解系统的输出)(n y ，并绘出n n x ~)(、n n h ~)(及n n y ~)(图形。
4、实验用MATLAB 函数介绍（１）数字信号处理中常用到的绘图指令(只给出函数名，具体调用格式参看help)figure(); plot(); stem(); axis(); grid on; title(); xlabel(); ylabel(); text(); hold on; subplot()（２）离散时间信号产生可能涉及的函数zeros(); ones(); exp(); sin(); cos(); abs(); angle(); real(); imag(); （３）卷积计算可能涉及的函数conv(); length()注：实验过程中也可以使用自己编制的自定义函数，如impseq()、stepseq()等。
数字信号处理实验指导书--Matlab版
下：
由于 e jω 是频率的周期函数，所以系统的频率特性也是频率的周期函数，且周期为 2π ，因 此研究系统频率特性只要在[ −π ,π ]范围内就可以了。
三、实验内容及步骤
1、在实验上机前，认真复习离散卷积和差分方程的有关内容，仔细阅读本实验原理和 步骤。准备实验所需数据。
2、离散卷积计算。 1）在主界面下进入“实验五” 的“离散卷积”子系统，选定几组不同的 x(n)、h(n) 代入离散卷积程序计算 y(n)，记录实验结果。 3、差分方程迭代解法 1）¨在主界面下进入“实验五” 的“差分方程”子系统，本实验要求首先确定系统方 程系数个数和需要输出 y(n)的样点个数，然后依次输入 a[0]~a[N]、b[0]~b[N]这些系数的 数值和 y(n)的初始值 y(0)~y(N-1)，h(n)的初始值 h(0)~h(N-1)，运行后即得冲激响应 h(n) 和阶跃信号激励下的响应 y(n)的结果和图形显示。 2）自己选择一个离散时间系统，写出其差分方程，设好初始条件进行实验并记录实验 结果。 4、在主界面下进入“实验六”即“离散系统 Z 域分析”，本实验中给出了计算系统
明了系统结构、参数、特性三者之间的关系，即同一结构，参数不同其特性也不同。 例如，下图所示离散系统:
数学模型由下列差分方程描述： y(n)=ay(n-1)+x(n)
系统函数 系统频率特性
H (z) = z , | z |> a z−a
幅度特性为
H (e jω )
=
e jω e jω − a
《数字信号处理》实验指导书(完整)
《数字信号处理》实验指导书通信教研室安阳工学院二零零九年三月第1章 系统响应及系统稳定性1.1 实验目的● 学会运用MATLAB 求解离散时间系统的零状态响应； ● 学会运用MATLAB 求解离散时间系统的单位取样响应； ● 学会运用MATLAB 求解离散时间系统的卷积和。
1.2 实验原理及实例分析1.2.1 离散时间系统的响应离散时间LTI 系统可用线性常系数差分方程来描述，即∑∑==-=-Mj jN i i j n x b i n y a 00)()( （1-1） 其中，i a （0=i ，1，…，N ）和j b （0=j ，1，…，M ）为实常数。
MATLAB 中函数filter 可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter 的语句格式为y=filter(b,a,x)其中，x 为输入的离散序列；y 为输出的离散序列；y 的长度与x 的长度一样；b 与a 分别为差分方程右端与左端的系数向量。
【实例1-1】 已知某LTI 系统的差分方程为试用MATLAB 命令绘出当激励信号为)()2/1()(n u n x n=时，该系统的零状态响应。
解：MATLAB 源程序为>>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)')程序运行结果如图1-1所示。
1.2.2 离散时间系统的单位取样响应系统的单位取样响应定义为系统在)(n 激励下系统的零状态响应，用)(n h 表示。
MATLAB 求解单位取样响应可利用函数filter ，并将激励设为前面所定义的impDT 函数。
《数字信号处理》实验指导书(正文)
实验一 离散时间信号分析一、实验目的1．掌握各种常用的序列，理解其数学表达式和波形表示。
2．掌握在计算机中生成及绘制数字信号波形的方法。
3．掌握序列的相加、相乘、移位、反褶等基本运算及计算机实现与作用。
4．掌握线性卷积软件实现的方法。
5．掌握计算机的使用方法和常用系统软件及应用软件的使用。
6．通过编程，上机调试程序，进一步增强使用计算机解决问题的能力。
二、实验原理1．序列的基本概念离散时间信号在数学上可用时间序列来表示，其中代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为∞<<∞-n 的整数，n 取其它值)(n x 没有意义。
离散时间信号可以是由模拟信号通过采样得到，例如对)(t x a 模拟信号进行等间隔采样，采样间隔为T ，得到一个{})(nT x a 有序的数字序列就是离散时间信号，简称序列。
2．常用序列常用序列有：单位脉冲序列（单位采样）)(n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。
3．序列的基本运算序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。
4．序列的卷积运算∑∞∞-*=-=)()()()()(n h n x m n h m x n y上式的运算关系称为卷积运算，式中代表两个序列卷积运算。
两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。
其计算的过程包括以下4个步骤。
(1)反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
(2)移位：将)(m h -移位n ，得)(m n h -。
当n 为正数时，右移n 位；当n 为负数时，左移n 位。
(3)相乘：将)(m n h -和)(m x 的对应点值相乘。
(4)求和：将以上所有对应点的乘积累加起来，即得)(n y 。
《数字信号处理》实验指导书
R6(n)
0.5
0 0.2
0
5
10
15
20
25 n
30
35
40
45
ቤተ መጻሕፍቲ ባይዱ
50
Magnitude
0.1
0 4
-3
-2
-1
0 Frequency(rad)
1
2
3
Phase
2 0 -2 -4 -3 -2 -1 0 Frequency(rad) 1 2 3
k＝0 时的直流分量及其合成的波形：
0.5
the Kth harmonic
the Kth harmonic
0
the Kth harmonic
0 5 10 15 20 25 n 30 35 40 45 50
0
-0.5
-0.5
0
5
10
15
20
sum of the first K+1 harmonics
1
sum of the first K+1 harmonics
25 n
30
35
0
-0.5
-0.5
0
5
10
15
20
sum of the first K+1 harmonics
1
sum of the first K+1 harmonics
25 n
30
35
40
45
50
1
0.5
0.5
0 0 5 10 15 20 25 n 30 35 40 45 50
0 0 5 10 15 20 25 n 30 35 40 45 50
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理实验指导手册西安文理学院机械电子工程系目录实验一离散时间信号 (2)实验二时域采样定理 (7)实验三离散时间系统 (10)实验四线性卷积与圆周卷积 (13)实验五用FFT作谱分析 (16)实验六用双线性变换法设计IIR数字滤波器 (18)实验七 FIR滤波器设计 (20)实验一 离散时间信号【实验目的】用MATLAB 实现离散时间信号的表示和运算，掌握MATLAB 的基本命令和编程方法，为后续实验打基础。
【实验原理】在数字信号处理中，所有的信号都是离散时间信号，因此应首先解决在MATLAB 中如何表示离散信号。
设一模拟信号经A/D 变换后，得到序列信号}),1(),0(),1(,{)}({)( x x x n x n x -==由于MATLAB 对下标的约定为从1开始递增，因此要表示)(n x ，一般应采用两个矢量，如：]5,4,3,2,1,0,1,2,3[---=n ]1,2,5,4,0,2,3,1,1[-=y这表示了一个含9个采样点的矢量：)}5(,),1(),2(),3({)(x x x x n y ---=【实验内容】熟悉下面序列（信号）的产生方法及相关运算 1、 单位采样序列 2、 单位阶跃序列 3、 信号翻转 4、 信号相加 5、 信号折叠 6、 信号移位【参考程序】单位采样序列1、impluse1.m （图1-1） n=10; x=zeros(1,n); x(1)=1;plot(x,'*');2、 impluse2.m （图1-2）n=-5:5; x=[n==0]; stem(x,'*');3、impluse3.m （图1-3） n=1:10; n0=3;x=[(n-n0)==1]; plot(x,'*');单位阶跃序列1、steps1.m （图1-4） n=10; x=ones(1,n); plot(x,'*');2、steps2.m （图1-5） n=10; x=ones(1,n); x(1)=0; x(2)=0;图1-1 单位采样序列1图1-2 单位采样序列2图1-3 单位采样序列3plot(x,'*');3、steps3.m （图1-6） n=0:10; n0=4;x=[(n-n0)>=0]; plot(x,'*');信号翻转（中心对称）flip.m （图1-7）x=0:6; figure(1); subplot(211); plot(x,'*'); y=fliplr(x); subplot(212); plot(y,'*');信号相加Showsigadd.m （图1-8） clear all figure(1);[x1 n1]=impseq(-4,-5,5);图1-4 单位阶跃序列1图1-5 单位阶跃序列2图1-6 单位阶跃序列3图1-7 信号翻转subplot(3,1,1);stem(n1,x1,'*');title('信号相加演示');ylabel('x1');[x2 n2]=stepseq(0,-3,7);subplot(3,1,2);stem(n2,x2,'*');ylabel('x2');[y n]=sigadd(x1,n1,x2,n2);subplot(3,1,3);stem(n,y,'*');ylabel('y=x1+x2)')图1-8 信号相加信号折叠%showfold.m（图1-9）clear all[x n]=impseq(2,-5,3);subplot(2,1,1);stem(n,x);title('信号折叠演示');[y,n]=sigfold(x,n);subplot(2,1,2);stem(n,y);信号移位%showsigshift.m（图1-10）figure(1)clear all[x,n]=stepseq(0,-5,5);subplot(211); stem(n,x,'*');title('信号移位演示’); [y,n]=sigshift(x,n,3); subplot(212) stem(n,y);【思考题】1、用MATLAB 实现信号的相加、移位、反转运算。
2、用MATLAB 计算信号的能量。
图1-9 信号折叠图1-10 信号移位实验二 时域采样定理【实验目的】1、掌握奈奎斯特抽样定理的含义，以及在实际应用中需要注意的问题。
2、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。
【实验原理】1、理想抽样信号连续信号)(t x a 经过时间间隔为T 的均匀采样后，得到的理想抽样信号)(t x a可表示为)()()(t t x t x T a a δ⋅=其傅里叶变换为∑∞-∞=-=k a a Tjkj XTj X )2(1)(πΩΩ 可见，一个连续时间信号经过理想抽样后，其频谱将以抽样频率s πΩ2=为间隔而重复，这就是频谱产生周期延拓。
如果信号的最高频谱h Ω超过2s Ω，则各周期延拓分量产生频谱的交叠，称为混叠现象。
要想抽样后能够不失真地还原出原信号，则抽样频率必须大于两倍信号谱的最高频率（)2h s Ω>Ω，这就是奈奎斯特抽样定理。
2、抽样恢复如果满足奈奎斯特抽样定理，即信号谱的最高频率小于折叠频率，则抽样后不产生频谱混叠，有)2()(1)(sa a j X Tj X ΩΩΩΩ<=故将)(Ωj X a通过以下理想低通滤波器：就可得到原信号频谱)()()()(ΩΩΩΩj X j H j X j Y a a a =⋅=⎪⎪⎩⎪⎪⎨⎧≥<=)2(0)2()(ssTj H ΩΩΩΩΩ所以输出端即为原模拟信号)()(t x t y a a =【实验内容】用实验的方法验证抽样定理。
对下面连续信号：0()sin()()at a x t Ae t u t -=Ω进行采样，可得到采样序列)()sin()()(0n u nT Ae nT x n x nT a a Ωα-==其中A 为幅度因子，α为衰减因子，0Ω是模拟角频率，T 为采用间隔。
这些参数都要在实验过程中由键盘输入，产生不同的)(t x a 和)(n x a 。
产生采样信号序列)(n x a ，使A=144.128,πα250=，π250=0Ω ，图2—1给出了连续信号)(t x a 的幅频特性曲线，由此图可以确定对)(t x a 应采用的采样频率。
1、取采样频率KHz sf 1=，即T =1ms .观察所得采样()a x n 的幅频特性|)(|ωj e X 和图2-1中的|)(|jf X a ，在折叠频率附近有无明显差别。
应当注意，实验中所得频谱是用序列()a x n 的傅里叶变换公式求得的，所以在频率度量上存在关系：T Ωω=为数字频率，Ω为模拟频率。
2、改变采样频率，Hz sf 300= ,观察|)(|ωj e X 的变化，并作记录（打印曲线）；进一步降低采样频率，Hzs f 200= ，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的|)(|ωj e X 曲线。
3、变采样频率，KHz s f 2= ,观察|)(|ωl e X 的变化，并作记录（打印曲线）；观察过采样时的曲线，并记录（打印）这时的|)| j eX曲线。
(【思考题】采样频率不同时，相应理想采样序列的傅里叶变换频谱的数字频率刻度是否都相同。
它们所对应的模拟频率是否相同？为什么？实验三 离散时间系统【实验目的】1、熟悉时域离散系统的时域特性。
2、利用卷积方法观察分析系统的时域特性。
3、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。
【实验原理】离散信号和系统在时域均可用序列来表示。
序列图形给人以形象直观的印象，他可加深我们对信号和系统的时域特征的理解。
本实验将观察分析几种信号及系统的时域特性。
为了在数字计算机上观察分析各种序列的频域特性，通常对)(ωj e X 在]2,0[π上进行M 点采样来观察分析。
对长度为N 的有限长序列x(n),有∑-=-=1)()(N n nj j k kem x eX ωω （3-1）其中 1,1,02-==M k k Mk πω通常M 应该取得大一些，以便观察谱的细节变化。
取模|)(k j e X ω|可绘出幅频特性曲线。
一个时域离散线性非移变系统的输入/输出关系为y(n)=x(n)*h(n)=∑-∞-∞=m m n h m x )()( （3-2）这里，y(n)为系统的输出序列，x(n)为输入序列。
h(n)、x(n)可以无限长，也可以是有限长。
为了计算机绘图观察方便，主要讨论有限长情况。
如果h(n)和x(n)的长度分别为N 和M ，则y(n)的长度为L=N+M-1。
这样，（3-2）式所描述的卷积运算就是序列移位、相乘和累加的过程。
上述卷积运算也可以在频域实现)()()(ωωωj j j e H e X e Y ⋅= （3-3）式（3-3式右边的相乘是在各频域点上的频谱值相乘。
【实验内容】一、认真复习离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理。
二、编制实验用相关子程序。
①信号产生子程序，用于产生实验中要用到的下列信号序列： a.单位采样序列：)()(n n x a δ= b.矩形序列： 10),()(==N n R n x N b②系统单位脉冲响应序列产生子程序。
本实验要用到两种FIR 系统。
a. 10()()a h n R n =b. ()() 2.5(1) 2.5(2)(3)b h n n n n n δδδδ=+-+-+-③有限长序列线性卷积子程序，用于完成两个给定长度的序列的卷积。
可以直接调用MATLAB 语言中的卷积函数conv 。
conv 用于两个有限长序列的卷积，它假定两个序列都从n=0开始。
调用格式如下：y=(,)conv x h其中参数x 和y 是两个已赋值的行向量序列。
三、在完成编制上述子程序的基础上，完成时域离散信号、系统和系统响应分析。
1、观察信号)(n x a 和系统()b h n 的时域和频域特性；利用线性卷积求得信号()b x n 通过系统()b h n 的响应y(n),比较所求响应y(n)和()b h n 的时域及频域特性，注意它们之间有无差别，绘图说明，并用所学理论解释所得结果。
2、观察系统)(n h a 对信号)(n x b 的响应特性。
利用线性卷积求系统响应y(n),并判断y(n)图形及其非零值序列长度是否与理论结果一致。