数字信号处理实验指导手册【模板】
数字语音信号处理实验指导书(学生用)
数字语音信号处理实验指导书(学生用）前言语音信号处理是研究用数字信号处理技术和语音学知识对语音信号进行处理的新兴的学科，是目前发展最为迅速的信息科学研究领域的核心技术之一。
通过语音传递信息是人类最重要、最有效、最常用和最方便的交换信息形式。
同时，语言也是人与机器之间进行通信的重要工具，它是一种理想的人机通信方式，因而可为信息处理系统建立良好的人机交互环境，进一步推动计算机和其他智能机器的应用，提高社会的信息化程度。
语音信号处理是一门新兴的学科，同时又是综合性的多学科领域和涉及面很广的交叉学科。
虽然从事这一领域研究的人员主要来自信号与信息处理及计算机应用等学科，但是它与语音学、语言学、声学、认知科学、生理学、心理学等许多学科也有非常密切的联系。
20世纪60年代中期形成的一系列数字信号处理的理论和算法，如数字滤波器、快速傅立叶变换（FFT）等是语音信号数字处理的理论和技术基础。
随着信息科学技术的飞速发展，语音信号处理取得了重大的进展：进入70年代之后，提出了用于语音信号的信息压缩和特征提取的线性预测技术（LPC），并已成为语音信号处理最强有力的工具，广泛应用于语音信号的分析、合成及各个应用领域，以及用于输入语音与参考样本之间时间匹配的动态规划方法；80年代初一种新的基于聚类分析的高效数据压缩技术—矢量量化（VQ）应用于语音信号处理中；而用隐马尔可夫模型（HMM）描述语音信号过程的产生是80年代语音信号处理技术的重大发展，目前HMM已构成了现代语音识别研究的重要基石。
近年来人工神经网络(ANN)的研究取得了迅速发展，语音信号处理的各项课题是促进其发展的重要动力之一，同时，它的许多成果也体现在有关语音信号处理的各项技术之中。
要求：1、用自己的手机录音“我是贵州大学科学院XX级XX专业的学生，名叫XXX，男（女），毕业于XX省XX市（县）XX中学。
2、在实验报告中要注明手机品牌，录音的采样频率：单（双）声道3、每个人实验报告只能处理自己的录音作为输入，实验报告采用电子文档提交，不接收纸质报告，附录音记录——切记！！！！！。
数字信号处理实验指导书2016-通信
数字信号处理实验指导书2016-通信数字信号处理实验徐俊2016年8⽉⽬录实验⼀离散时间系统及系统响应 (2)实验⼆离散傅⽴叶变换及其特性验证 (8)实验三时域采样与频域采样 (17)实验四冲激响应不变法IIR数字滤波器设计 (24)实验⼀离散时间系统及系统响应⼀、实验⽬的1、掌握求解离散时间系统脉冲响应和阶跃响应的⽅法；2、掌握⽤线性卷积求解离散时间系统响应的基本⽅法。
⼆、实验原理与设计⽅法1、⽤impz和dstep函数求解离散系统的单位脉冲响应和阶跃响应【例1-1】已知某因果系统的差分⽅程为y n+0.5y n?1=x n+2x(n?2)系统为零状态，求系统的脉冲响应和阶跃响应。
解：该系统是⼀个2阶系统，列出b m和a k系数为a0=1，a1=0.5，a2=0，b0=1，b1=0，b2=2MALAB程序如下（取16点作图）：a=[1,0.5,0];b=[1,0,2];n=16;hn=impz(b,a,n); %脉冲响应gn=dstep(b,a,n); %阶跃响应subplot(1,2,1),stem(hn,'k');title('系统的单位脉冲响应');ylabel('h(n)');xlabel('n');axis([0,n,1.1*min(hn),1.1*max(hn)]);subplot(1,2,2),stem(gn,'k');title('系统的单位阶跃响应');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);结果如下图所⽰：2、⽤conv函数进⾏卷积计算求系统响应【例1-2】某离散时间系统的脉冲响应为h b(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)激励信号为x t=Ae?αnT sinΩ0nT 0≤n<50设A=444.128，α=502π，Ω0=502π。
数字信号处理实验指导书
注意此书用的时候N要先付值数字信号处理实验指导书目录前言 (1)第一章MATLAB基础知识 (1)第二章MATLAB基本数值运算 (4)第三章MATLAB的图形处理功能 (8)第四章MATLAB的程序设计 (11)第五章常用数字信号处理函数 (16)第六章MATLAB在数字信号处理中的应用 (23)实验一常见离散信号的MATLAB产生和图形显示 (33)实验二离散系统的频率响应分析和零、极点分布 (37)实验三序列线性卷积、圆周卷积的计算及其关系的研究 (39)实验四利用DFT分析信号的频谱 (41)实验五信号时间尺度变换的研究 (43)实验六快速傅里叶变换及其应用 (47)实验七IIR滤波器的实现与应用 (56)实验八FIR滤波器的实现与应用 (61)第一章MATLAB基础知识§1-1 MA TLAB软件简介MATLAB，Matrix Laboratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。
它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个界面友好的用户环境，在这个环境中，问题与求解都能方便地以数学的语言（主要是矩阵形式）或图形方式表示出来。
与大家常用的Fortran 和C等高级语言相比，MA TLAB的语法规则更简单，更贴近人的思维方式，被称为“草稿纸式的语言”。
§1-2 MA TLAB应用入门1．MATLAB的安装与卸载MATLAB软件在用户接口时具有较强的亲和力，其安装过程比较典型，直接运行光盘中的安装向导支撑程序SETUP.exe，按其提示一步步选择即可。
MATLAB自身带有卸载程序，在其安装目录下有uninstall子目录，运行该目录下uninstall.exe的即可；也可以通过Windows系统的安装卸载程序进行卸载。
2．MATLAB的启动与退出MATLAB安装完成后，会自动在Windows桌面上生成一个MA TLAB图标，它是指向安装目录下\bin\win32\matlab.exe的链接，双击这个图标即可来到MATLAB集成环境的基本窗口；也可以在开始菜单的程序选项中选择MATLAB 快捷方式；还可以在MA TLAB的安装路径的bin子目录中双击可执行文件matlab.exe。
数字信号处理实验指导书(定稿新)
数字信号处理实验指导书前言数字信号处理是是电子信息工程专业的一门专业基础课。
本课程主要研究如何对信号进行分析、变换、综合、估计与识别等加工处理的基本理论和方法，数字信号处理实验是验证、巩固和补充课堂讲授的理论知识的必要环节。
通过实验，使学生巩固所学基本理论，掌握最基本的数字信号处理的理论和方法，提高综合运用所学知识，提高计算机编程的能力。
进一步加强学生独立分析问题、解决问题的能力、综合设计及创新能力的培养，同时注意培养学生实事求是、严肃认真的科学作风和良好的实验习惯，为今后的工作打下良好的基础。
数字信号处理实验指导书针对每个实验介绍了MATLAB语言数字信号处理工具箱中的相应函数，举例并附有相应的程序。
为配合课堂理论教学，实验内容安排仍从认识性和验证性入手，逐步增加设计性和工程应用性内容，达到训练实验技能和积累工程实际应用经验之目的。
数字信号处理实验成绩按百分制核定。
预习占20%，实验过程占40%，实验报告占40%。
电气电子信息工程系电工电子基础教研室2007.03实验要求在实验过程中，要求学生做到：（1）预习实验指导书有关部分，认真做好实验内容的准备工作，就实验可能出现的情况提前作出思考和分析，需要计算的参数提前完成计算工作，并认真写出预习报告。
（2）仔细观察实验过程中图形随参数的变化，记录图形变化的主要情况，作出必要说明和分析。
（3）认真书写实验报告并在规定的时间内把实验报告交给辅导教师。
实验报告包括实验目的和要求，实验情况及其分析。
对需要编程的实验，写出程序设计说明，给出源程序框图和清单。
（4）遵守机房纪律，服从辅导教师指挥，爱护实验设备。
（5）实验课程不迟到。
如有事不能出席，所缺实验一般不补。
实验验收分为两个部分。
第一部分是上机操作，包括检查程序运行和即时提问。
第二部分是提交书面的实验报告。
每个实验都应当在规定的时间内完成并检查通过，过期视为未完成该实验，扣该实验操作成绩。
为避免期末集中检查方式产生的诸多不良问题，希望同学们抓紧时间，合理安排，认真完成。
《数字信号处理》实验指导书学生版
2015－2016学年第1学期学院物理与电子信息学院教研室信息与通信工程教研室课程名称数字信号处理授课班级13电信本、13电信卓越主讲教师黄隆胜职称副教授2014年03月目录前言 (1)实验一熟悉MATLAB环境 (4)实验二用MATLAB进行离散系统的Z域分析 (6)实验三傅立叶变换 (8)实验四IIR及FIR滤波器的MATLAB实现 (11)前言MATLAB是由美国Math Works公司推出的软件产品。
MATLAB是“Matrix Laboratory”的缩写，意及“矩阵实验室”。
MATLAB是一完整的并可扩展的计算机环境，是一种进行科学和工程计算的交互式程序语言。
它的基本数据单元是不需要指定维数的矩阵，它可直接用于表达数学的算式和技术概念，而普通的高级语言只能对一个个具体的数据单元进行操作。
因此，解决同样的数值计算问题，使用MATLAB要比使用Basic、Fortran和C语言等提高效率许多倍。
许多人赞誉它为万能的数学“演算纸”。
MATLAB采用开放式的环境，你可以读到它的算法，并能改变当前的函数或增添你自己编写的函数。
在欧美的大学和研究机构中，MATLAB是一种非常流行的计算机语言，许多重要的学术刊物上发表的论文均是用MATLAB来分析计算以及绘制出各种图形。
它还是一种有利的教学工具，它在大学的线性代数课程以及其它领域的高一级课程的教学中，已成为标准的教学工具。
最初的MATLAB是用FORTRAN编写的，在DOS环境下运行。
新版的MATLAB 是C语言编写的高度集成系统。
它在几乎所有流行的计算机机种，诸如PC、MACINTOSH、SUN、VAX上都有相应的MATLAB版本。
新版的MATLAB增强了图形处理功能，并在WINDOWS环境下运行。
现今，MATLAB的发展已大大超出了“矩阵实验室”的范围，在许多国际一流专家学者的支持下，Maths Works公司还为MATLAB 配备了涉及到自动控制、信息处理、计算机仿真等种类繁多的工具箱（Tool Box），这些工具箱有数理统计、信号处理、系统辨识、最优化、稳健等等。
数字信号处理实验指导书(2010)v1.001
数字信号处理 实验指导书（V1.001版）成都理工大学信息工程学院2010年9月目录第一部分实验基础 (3)第一章序列的产生 (3)1.1 单位采样序列 (3)1.2 单位阶跃序列 (4)1.3 指数序列 (5)1.4 实正弦序列 (6)1.5 斜坡序列 (7)1.6 复正弦序列 (8)1.7 随机序列 (9)第二章序列的基本运算 (10)2.1 信号的延迟 (10)2.2 信号相加 (11)2.3 信号相乘 (12)2.4 信号乘以标量值 (14)2.5 信号翻转 (14)2.6 信号和 (15)2.7 信号积 (16)2.8 信号能量 (16)第三章序列的变换 (18)3.1 Z变换 (18)3.2 Chirp Z变换 (19)3.3 DFT (21)3.4 FFT (22)3.5 DCT (23)3.6 Hilbert变换 (25)第四章 IIR数字滤波器设计 (28)4.1 Butterworth低通滤波器 (29)4.2 Chebyshev I型低通滤波器 (31)4.3 Chebyshev I型高通滤波器 (32)4.4 Chebyshev I型带通滤波器 (33)4.5 Chebyshev I型带阻滤波器 (34)4.6 基于冲激响应不变法的IIR数字滤波器 (35)第五章 FIR数字滤波器设计 (38)5.1 基于频率抽样法的FIR数字滤波器设计 (38)5.2 基于Chebyshev逼近法的FIR数字滤波器设计 (40)第二部分实验内容 (42)实验一信号、系统及系统响应 (42)实验二使用FFT作谱分析 (42)实验三使用双线性Z变换法设计IIR滤波器 (42)实验四使用窗函数法设计FIR滤波器 (43)第一部分实验基础第一章序列的产生1.1 单位采样序列（1）例程% 程序1.1-1% 一个单位采样序列的产生clf; % 清屏n = 0:20; % 产生从0到20的一个向量u = [1, zeros(1,20)]; % 产生单位采样序列stem(n,u); % 绘制单位采样序列xlabel(‘时间序号(n)’);ylabel(‘振幅’);title(‘单位采样序列’);axis([0 20 0 1.2]);% 程序1.1-2% 一个移位单位采样序列的产生clf; % 清屏n = -10:20; % 产生从-10到20的一个向量u = [zeros(1,10), 1, zeros(1,20)]; % 产生单位采样序列stem(n,u); % 绘制单位采样序列xlabel(‘时间序号(n)’);ylabel(‘振幅’);title(‘单位采样序列’);axis([-10 20 0 1.2]);（2）程序运行结果时间序号(n)振幅单位采样序列时间序号(n)振幅单位采样序列图1.1-1 例程1.1-1运行结果 图1.1-2 例程1.1-2运行结果1.2 单位阶跃序列（1）例程% 程序1.2-1% 一个单位阶跃序列的产生clf; % 清屏n = 0:20; % 产生从0到20的一个向量 s = [ones(1,20)]; % 产生单位阶跃序列stem(n,u);% 绘制单位采样序列 xlabel(‘时间序号(n)’); ylabel(‘振幅’);title(‘单位采样序列’); axis([0 20 0 1.2]);（2）程序运行结果时间序号(n)振幅单位阶跃序列图1.2-1 例程1.2-1的运行结果1.3 指数序列（1）例程% 程序1.3-1% 一个复指数序列的产生clf; % 清屏c = -(1/12) + (pi/6)*I;K = 2;n = 0:40;x = K*exp(c*n);subplot(2,1,1);stem(n, real(x)); % 绘制复指数序列实部xlabel(‘时间序号(n)’);ylabel(‘振幅’);title(‘实部’);subplot(2,1,2);stem(n, imag(x)); % 绘制复指数序列虚部xlabel(‘时间序号(n)’);ylabel(‘振幅’);title(‘虚部’);% 程序1.3-2% 一个实指数序列的产生clf; % 清屏n = 0:35;a = 1.2;K = 0.2;x = K*a.^n;stem(n,x); % 绘制实指数序列xlabel(‘时间序号(n)’);ylabel(‘振幅’);（2）程序运行结果时间序号(n)振幅时间序号(n)振幅虚部时间序号(n)振幅图1.3-1 例程1.3-1的运行结果 图1.3-2 例程1.3-2的运行结果1.4 实正弦序列（1）例程% 程序1.4-1% 一个实正弦序列的产生clf; % 清屏 n = 0:40; f = 0.1; phase = 0; A = 1.5;arg = 2*pi*f*n – phase; x = A * cos(arg);stem(n, x); % 绘制正弦序列 xlabel(‘时间序号(n)’); ylabel(‘振幅’); title(‘正弦序列’); grid;axis([0 40 -2 2]);（2）程序运行结果时间序号(n)振幅正弦序列图1.4-1 例程1.4-1的运行结果1.5 斜坡序列（1）例程% 程序1.5-1% 一个斜坡序列的产生clear all ; N = 32; k = 4; B = 3;x = [zeros(1,k),ones(1,N-k)]; for i=1:N x(i) = B*x(i)*(i-k); endxn = 0:N-1; stem(xn,x);axis([-1 32 0 90]);（2）程序运行结果图1.5-1 例程1.5-1的运行结果1.6 复正弦序列（1）例程% 程序1.6-1% 一个复正弦序列的产生clear all;N = 32;A = 3;w = 314;xn = 0:N-1;x = A*exp(j*w*xn);stem(xn,x);axis([-1 32 -3.2 3.2]);（2）程序运行结果图1.6-1 例程1.6-1的运行结果1.7 随机序列（1）例程% 程序1.7-1% 一个随机序列的产生lear all;N = 8;x = zeros(1,N);x(1) = 8.0;x(2) = 3.4;x(3) = 1.8;x(4) = 5.6;x(5) = 2.9;x(6) = 0.7;xn = 0:N-1;stem(xn,x);axis([-1 8 0 8.2]);（2）程序运行结果图1.7-1 例程1.7-1的运行结果第二章 序列的基本运算在数字信号处理中，信号的基本运算包括移位、相加、相乘和翻转等。
《数字信号处理》实验指导书v0.3
i
实验一 时域采样与频域采样
1. 实验目的
(1) 进一步加深 DFT 算法原理和基本性质的理解。 (2) 熟悉 FFT 算法原理和 FFT 子程序的应用。 (3) 学习用 FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误 差及其原因。
2.实验原理
(1) 复习 DFT 的定义、性质和用 DFT 作谱分析的有关内容。 (2) 复习 FFT 算法原理与编程思想。 (3) 查询 Matlab 中 fft 函数的用法，学会通过改变 fft 的参数计算不同 DFT。 (4) 编制信号产生子程序，产生以下典型信号供谱分析用：
画图时，用 20 lg | H (k ) | 打印幅度特性。第 k 点对应的频率
ωk =
2π k N 。为使曲线包络
jω DFT 变换区间要选大些。 例如窗口长度 N = 33 时， 线更接近 H (e ) 的幅度特性曲线，
可通过在 h(n) 末尾补零的方法，使长度变为 64，再进行 64 点 DFT，则可得到更精 确的幅度衰减特性曲线。
《数字信号处理》实验指导
浙江工业大学信息工程学院通信工程系 2010.4
目录
目录 ........................................................................................................................ i 实验一 时域采样与频域采样 ............................................................................. 1 实验二 用双线性变换法设计IIR数字滤波器 ................................................... 4 实验三 用窗函数法设计FIR数字滤波器........................................................... 6 附录：部分Matlab函数........................................................................................ 9
《数字信号处理》实验指导书
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
《数字信号处理》实验指导书-167
《数字信号处理》实验指导书实验1 数字滤波器的设计实验序号：1 实验名称：数字滤波器的设计 适用专业：通信工程、电子信息工程 学 时 数：4学时一、实验目的1．掌握双线性变换法设计IIR 数字滤波器的具体设计方法及其原理和窗函数设计FIR 滤波器的设计原理与基本方法。
2．观察双线性变换的频域特性。
熟悉Butterworth 滤波器的频率特性。
3．了解各种不同窗函数对滤波器性能的影响。
4．熟悉Matlab 计算机编程。
二、实验原理1．用双线性变换法设计IIR 数字滤波器方法（1）设计思想：将模拟滤波器转换成数字滤波器的实质是，用一种从s 平面到z 平面的映射函数将Ha(s)转换成H(z)。
对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。
(2)数字滤波器的频率响应模仿模拟滤波器的频响，s 平面的虚轴映射z 平面的单位圆，相应的频率之间成线性关系。
脉冲响应不变法和双线性变换法都满足如上要求。
s 平面与z 平面之间满足以下映射关系：1111--+-=zz ss 平面的虚轴单值地映射于z 平面的单位圆上，s 平面的左半平面完全映射到z 平面的单位圆内。
双线性变换不存在混叠问题。
双线性变换时一种非线性变换)/ω(tg 2=Ω，这种非线性引起的幅频特性畸变可通过预畸而得到校正。
（2）以低通数字滤波器为例，将设计步骤归纳如下：·确定数字滤波器的性能指标：通带临界频率f p 、阻带临界频率f s ；通带内的最大衰减A p ；阻带内的最小衰减A s ；· 确定相应的数字角频率，ωp=2πfp ；ωs=2πfs ；·计算经过预畸的相应模拟低通原型的频率，)2/(ωtg =Ω；·根据Ωp 和Ωs 计算模拟低通原型滤波器的阶数N ，并求得低通原型的传递函数Ha(s)； ·用上面的双线性变换公式代入Ha(s)，求出所设计的传递函数H(z)； ·分析滤波器特性，检查其指标是否满足要求。
数字信号处理实验指导书
% Program P1_2
% Generation of a sinusoidal sequence
n = 0:40;
பைடு நூலகம்f = 0.1;
phase = 0;
A = 1.5;
arg = 2*pi*f*n - phase;
x = A*cos(arg);
clf;
% Clear old graph
附录A MATLAB系统的常用概念 .........................................................................28
附录B
信号处理工具箱函数...........................................................................33
分析，从而进一步研究它们的性质。 2．熟悉离散时间序列的 3 种表示方法：离散时间傅立叶变换（DTFT）,离
散傅立叶变换（DFT）和 Z 变换。
二．实验相关知识准备 1．用到的 MATLAB 命令 运算符和特殊字符： < > .* ^ .^ 语言构造与调试： error function pause 基本函数： angle conj rem 数据分析和傅立叶变换函数： fft ifft max min 工具箱： freqz impz residuez zplane
数字信号处理 实验指导书
长沙理工大学电气与信息工程学院
目
录
实验一：离散时间信号的时域分析........................................................................3 实验二：离散时间系统的时域分析........................................................................6 实验三：离散时间信号的频域分析........................................................................9 实验四：线性时不变离散时间系统的频域分析..................................................13 实验五： IIR数字滤波器的设计...........................................................................17 实验六： FIR数字滤波器的设计..........................................................................24
数字信号处理实验指导书
《数字信号处理》实验指导书安阳工学院电子信息与电气工程学院目录实验一离散时间信号的表示及运算 (1)实验二离散时间LTI系统的时域分析 (10)实验三 z变换及离散时间LTI系统的z域分析 (17)实验四离散傅立叶变换及其快速算法 .................... ....... .. (29)实验五 IIR数字滤波器的MATLAB实现 (46)实验六 FIR数字滤波器的MATLAB实现.............................. .. (54)附录 MATLAB主要命令函数表....................... . (61)实验一 离散时间信号的表示及运算一、实验目的1、学会运用MATLAB 表示的常用离散时间信号；2、学会运用MATLAB 实现离散时间信号的基本运算。
二、实验条件装有MATLAB 软件的PC 机 三、 实验原理1、离散时间信号在MATLAB 中的表示离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。
离散序列通常用)(n x 来表示，自变量必须是整数。
离散时间信号的波形绘制在MATLAB 中一般用stem 函数。
stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。
如果要实心，需使用参数“fill”、“filled”，或者参数“.”。
由于MATLAB 中矩阵元素的个数有限，所以MATLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。
类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。
①单位取样序列单位取样序列)(n δ，也称为单位冲激序列，定义为)0()0(01)(≠=⎩⎨⎧=n n n δ （1-1） 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。
在MATLAB 中，冲激序列可以通过编写以下的impDT.m 文件来实现，即function y=impDT(n)y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。
《数字信号处理》实验指导书-新
实验一：离散时间序列卷积和MATLAB实现实验学时：2实验类型：验证实验要求：必修（一）实验目的：学会用MATLAB对信号与系统分析的方法，理解离散序列卷积和的计算对进行离散信号与系统分析的重要性。
（二）实验原理：1、离散时间序列f1(k)和f2(k)的卷积和定义：f(k)=f1(k)*f2(k)=∑∞-∞=-∙iikfif)(2)(12、在离散信号与系统分析中有两个与卷积和相关的重要结论：a、f(k)= ∑∞-∞=-∙iikif)()(δ=f(k)* δ(k)即离散序列可分解为一系列幅度由f(k)决定的单位序列δ(k)及其平移序列之积。
b、对线性时不变系统，设其输入序列为f(k)，单位响应为h(k)，其零状态响应为y(k)，则有：y(k)= ∑∞-∞=-∙iikhif)()(（三）实验内容：conv.m用来实现两个离散序列的线性卷积。
其调用格式是：y=conv(x,h)若x的长度为N，h的长度为M，则y的长度L=N+M-1。
题一：令x(n)= {}5,4,3,2,1，h(n)＝{}246326，，，，，，y(n)=x(n)*h(n)，求y(n)。
要求用subplot和stem画出x(n),h(n),y(n)与n的离散序列图形。
题二：已知序列f1(k)=⎩⎨⎧≤≤其它0201k f2(k)=⎪⎪⎩⎪⎪⎨⎧===其它332211k k k调用conv()函数求上述两序列的卷积和题三：编写计算两离散序列卷积和f(k)=f1(k)*f2(k)的实用函数dconv().要求该程序在计算出卷积和f(k)的同时，还绘出序列f1(k),f2(k)和f(k)的时域波形图，并返回f(k)的非零样值点的对应向量。
function[f,k]=dconv(f1,f2,k1,k2)%f1(k),f2(k)及f(k)的对应序号向量分别为k1，k2和k 。
题四：试用MATLAB 计算如下所示序列f1(k)与f2(k)的卷积和f(k)，绘出它们的时域波形，并说明序列f1(k)与f2(k)的时域宽度与序列f(k)的时域宽度的关系。
数字信号处理指导书
实验一 信号、系统及系统响应一、 实验目的：1、 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定的理解。
2、 熟悉时域离散系统的时域特性。
3、 利用卷积方法观察分析系统的时域特性。
4、 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。
二、实验内容及步骤：1、 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有内容，阅读本实验原理与方法。
2、 了解要使用到的MA TLAB 命令：1）基于DTFT 离散时间信号分析函数：freqz ，real ，imag ，abs ，angle ，unwrap 。
函数freqz 可以用来计算一个以ωj e的有理分式形式给出的序列的DTFT 值。
freqz 的形式多样，常见的有H=freqz(num,den,w),其中num 表示序列有理分式DTFT 的分子多项式系数，den 表示分母多项式系数（均按z 的降幂排列），矢量w 表示在0～π2中给定的一系列频率点集合。
freqz 函数的其他形式参见帮助。
在求出DTFT 值后，可以使用函数real, imag, abs 和angle 分别求出并绘出其实部，虚部，幅度和相位谱。
如果需要，还可以用unwrap 函数消除相位中的跳变。
2）函数fft(x)可以计算R 点序列的R 点DFT 值；而fft(x,N)则计算R 点序列的N 点DFT ，若R>N ，则直接截取R 点DFT 的前N 点，若R<N ，则x 先进行补零扩展为N 点序列再求N 点DFT 。
函数ifft(X)可以计算R 点的谱序列的R 点IDFT 值；而ifft(X,N)同fft(x,N)的情况。
卷积函数2）conv 用求于两个有限长度序列的卷积。
3、 实验准备工作：1）产生实验中要用到的下列信号序列： a ．采样信号序列：对下面连续信号： )()s i n ()(0t u t Aet x ata Ω=-进行采样，可得到采样序列)()sin()()(0n u nT AenT x n x anta a Ω==-，500<≤n其中A 为幅度因子，a 为衰减因子，0Ω是模拟角频幸，T 为采样间隔。
数字信号处理实验指导书_学生
The experiment instruction Digital Signal ProcessingSchool of AutomationContents1.1 Discrete-time Signal Analysis1.2 Discrete Time Domain System Analysis1.3 Frequency domain system analysis2.1 FFT2.2 Digital filter design2.3 TEXT BOOK3.1 Signal analysis and processing (Gas-solid flow signal)3.2 Signal analysis and processing (Magnetic resonance signal)4 Additional experimentsExperiment 1.1 Discrete-time Signal Analysis1. PurposeBe familiar with the characteristics of discrete-time sinusoidal signal.2. Contents(1) Consider x M (n)=sin(2πMn/N), assume N=12, when M=4,5,7 and 10, draw x M (n)respectively with stem in the range 0≤n ≤2N-1. What are the minimum periods of above signals? How to determine the minimum period of signal with arbitrary positive integers M and N? The case of M>N must be considered.(2) Consider x k (n)=sin(ωk n), where ωk =2πk/5. Let k=1,2,4 and 6,draw the curves ofx k (n) with stem in the range 0≤n ≤9. And all the curves should be put in one graph with subplot . How many different signals in the graph? Please give the reason why the signals are same with different ωk ?3. Experiment report requirements(1) Give the M files and the figures.(2) Answer the questions.Experiment 1.2 Discrete Time Domain System Analysis1. Purpose(1) Calculation of discrete-time convolution(2) Determine the linearity and shift-invariance of the systems.2. Content2.1 Calculation of discrete-time convolution(1) Calculate the convolutions of the following signals with conv. Draw the curves of x(n), h(n) and y(n) ( x(n)* h(n)) respectively.140)5/4()(,01401)()1≤≤=⎩⎨⎧≤≤=n n h others n n x n 190)5.0sin(5.0)(,0901)()2≤≤=⎩⎨⎧≤≤=n n n h others n n x(2) Consider x(n) = 0.5n u(n)h(n)= u(n)If calculate y(n)=x(n)*h(n) with conv , the problem of infinite lengths of x(n) and h(n) should be solved. Save the x(n) in the range of 0≤n≤24 in vector x, the h(n) in the range of 0≤n≤14 in vector h, the result conv(h,x) in vector y. Draw the y(n) with stem, comparing with 0.5n u(n)* u(n), please point out which values are true and which values are false.2.2 Determine the linearity and shift-invariance of the systemsConsider the following three systemsSystem1: w(n)=x(n)-x(n-1)-x(n-2)System2: y(n)=cos(x(n))System3: z(n)=nx(n)where x(n) is the input, w(n) ,y(n) and z(n) are the outputs.(1)Consider the three input signalsx1(n)=δ(n)x2(n)=δ(n-1)x3(n)=δ(n)+2δ(n-1)the output of system1 are stored respectively in vectors w1, w2 and w3(0≤n≤5). Draw the curves of w1, w2, and w1+2w2 in one graph.(2)Repeat (1) with system2.(3)Repeat (1) with system3.(4)Which of the systems is the linear and shift invariant system?3. Experiment report requirements(1) Give the M files and the figures.(2) Answer the question.Experiment 1.3 Frequency domain system analysis1. Purposes(1)Be familiar with the frequency characteristics of discrete-time systems.(2)Learn to analyze systems in analog domain and digital domain by a computer.2. Contents[1] Consider the system H(Z)=Z/（Z-A ） , analyzing the amplitude spectrum and the phase spectrum of )(ωj e H in the range 0≤ω≤4π when(1) A= 0.5; (2) A= -0.5; (3) A= 0[2] Analyzing the amplitude spectrum of H （j Ω）in the range 0≤Ω≤2π.2432)(SS S H ++= [3] Design the DF and analyze the amplitude spectrum in the range 0≤ω≤2π.(1) The prototype is as above, assume the sampling interval Ts is 0.5s, determine the digitalfilter H 1(z) with impulse invariance method.(2) Assume the sampling interval Ts is 0.1s. Repeat (1) and obtain the digital filter H 2(z).(3) Plot H 1(z) and H 2(z) in the same graph with subplot .3. Experiment report requirements[1] Calculate H 1(z) and H 2(z);[2] Give the M files and the figures.Experiment 2.1 FFT1. Purposes(1) Learn to analyze the signal in frequency domain with FFT.(2) Learn to calculate the convolution with FFT.2. Contents2.1 Assume T=0.0125s ，N=16，calculate the amplitude spectrum and the phase spectrum of signal x(t)=e -t (t ≥0). And draw the graphs with stem .2.2 Calculate the convolution with FFT and IFFT.90................1.0)(90)........1.01()(≤≤=≤≤-=n n n h n n n x 2.3 Give the amplitude spectrums of the following squences. (1) 1()sin(),074x n n n π=≤≤(2) 2()sin(),0154x n n n π=≤≤ (3) 3()sin(),084x n n n π=≤≤ (4) 42()sin(),045x n n n π=≤≤ 3. Experiment report requirements（1）Give the M files and the figures.（2）Discuss the exprimental results of the third experiment.Experiment 2.2 Digital filter design1. PurposeDigital filter design and digital signal processing.2. ContentsSave the data in vector y 0. (数据文件，后缀名为DAT)[1] High pass filter design(1) Exact correspondence at 0.006π rad in the behavior of the digital filter with respect to 10 rad/s of the reference prototype analog filter. The prototype is10)(1+=s s s H Determine the digital filter H 1(z) using the bilinear transformation.(2) Filter the signal y0 by the above high pass filter and save the filtered signal in vector y1.[2] Low pass filter design (1)The prototype is11)(2+=s s H Assume the sampling interval is 0.001s, determine the digital filter H 2(z) using the impulse invariance transformation.(2) Filter the signal y0 by the above low pass filter and save the filtered signal in vector y2. 3 Experiment report requirements[1] Calculate H 1(z) and H 2(z);[2] Give the M files;[3] Give the figures of y 0, y 1,y 2, and the amplitude spectrums of )(1ωj e H and )(2ωj e H .Experiment 2.3TEXT BOOK: 2.41，2.43，4.25，5.24Experiment 3.1 Signal analysis and processing1. PurposeAnalyze the characteristics of the gas-solid signal and calculate the mass of solid passing through the pipe(管道输运粉体质量的计算).2. Contents[1] Read the books about measuring the gas-solid mass flow;[2] Analyze the gas-solid mass flow signal;[3] Calculate the solid mass.3 Experiment report requirements[1] Describe the method of measuring the gas-solid mass flow;[2] Give the signal processing steps;[3] Give the M files;Experiment 3.2 Signal analysis and processing磁共振测量数据的分析和处理1. 实验目的从实验装置的输出测量信号中恢复出原始的给定输入信号。
数字信号处理实验指导书(更新)
.《数字信号处理》实验指导书吴启宏编华北电力大学2007 年 3 月前言1．实验总体目标通过实验，学生应能深入了解“数字信号处理”这门课程的知识点，加深对基本概念的理解，同时基本掌握MATLAB语言，为今后进一步学习和研究打下一个良好的基础。
⒉适用专业全校工科各专业⒊先修课程信号与系统⒋实验课时分配⒌实验环境普通配置计算机，安装MATLAB6.5或以上版本。
⒍实验总体要求自行编写MATLAB程序，完成上机实验内容，并提交实验报告。
⒎本实验的重点、难点及教学方法建议本实验的重点为“数字信号处理”的基本理论及应用；难点在MATLAB语言的熟练使用。
教学方法建议：学生在实验前先自学MATLAB语言，指导教师将实验中用到的主要的MATLAB函数列出，提前给学生，供参考。
目录实验一、熟悉MATLAB环境 (3)实验二、快速傅里叶变换（FFT）及其应用 (5)实验三、IIR数字滤波器的设计 (10)实验四、FIR数字滤波器的设计 (13)实验一 熟悉MATLAB 环境一、实验目的1. 熟悉MATLAB 环境的主要操作命令。
2. 学会简单的矩阵输入和数据读写。
3. 掌握简单的绘图命令4. 用MA TLAB 编程并学会创建函数。
5. 观察离散系统的频率响应。
二、实验类型验证型三、实验仪器装有MATLAB 语言的计算机四、实验原理MATLAB 的基本命令和函数五、实验内容和要求1. 数组的加、减、乘、除和乘方运算。
输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B ，D=A-B ，E=A.*B ，F=A./B ，G=A.^B ，并用stem 语句画出A 、B 、C 、D 、E 、F 、G 。
2. 用MATLAB 实现下列序列： (1) 1508.0)(≤≤=n n x n(2) 150)()32.0(≤≤=+n en x nj(3) 150)1.025.0sin(2)2.0125.0cos(3)(≤≤+++=n n n n x ππππ(4) 将（3）中的)(n x 扩展为以16为周期的函数)16()(16+=n x n x ，绘出四个周期。
《数字信号处理》实验指导书
的相角, Ai 就是极点 pi 到单位圆上的点 e jω 的矢量长度(距离),而θ i 就是该矢量 的相角,因此有:
M
∏ B e j(ψ1 +ψ 2 +⋅⋅⋅⋅+ψ M ) j
H (e jω ) =
j =1 N
= H (e jω ) e jϕ (ω )
∏ A e j(θ1+θ2 +⋅⋅⋅⋅+θ N ) i
(1) 设有直流信号 g(t)=1,现对它进行均匀取样，形成序列 g(n)=1。试讨 论若对该序列分别作加窗、补零，信号频谱结构有何变化。 四、实验过程及结果（含程序）
12
13
14
15
16
实验三 IIR 数字滤波器的设计
一、实验目的 (1)掌握双线性变换法及脉冲相应不变法设计 IIR 数字滤波器的具体设计 方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和 带通 IIR 数字滤波器的计算机编程。 (2)观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双 线性变换法及脉冲响应不变法的特点。 (3)熟悉 Butterworth 滤波器、Chebyshev 滤波器和椭圆滤波器的频率特 性
《数字信号处理》
实验指导书
班级： 学号： 姓名： 苏州科技学院 电子教研室
实验一 信号、系统及系统响应
一、实验目的
(1) 熟悉 MATLAB 平台的使用，掌握离散信号、离散系统的 MATLAB 实现。 (2)掌握根据系统函数绘制系统零极点分布图的基本原理和方法。 (3)理解离散系统频率特性分析的基本原理，掌握根据系统函数零极点分布来分 析离散系统频率响应的几何矢量法。
17
变换类型 低通
Байду номын сангаас
数字信号处理实验指导书
1.5
1
0.5
0
-0.5
-1
-1.5
-2
0
2
4
6
8
10
12
14
16
18
20
（3）用 impz 函数 a1=[1,0.75,0.125]; b1=[1,-1]; impz(b1,a1,21);
Impulse Response 1.5
1
0.5
Amplitude
0
-0.5
-1
-1.5
-2
0
2
4
6
8
10 12 n (samples)
5
……
程序计算结果： I.
y[n] + 0.75 y[n − 1] + 0.125 y[n − 2] = x[n] − x[n − 1]
a. 单位冲激响应： (1) 用 filter 函数 a1=[1,0.75,0.125]; b1=[1,-1]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title('y1filter'); xlabel('x'); ylabel('y');
3. 编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分
别用 filter、conv、impz 三种函数完成。
y[n] + 0.75 y[n − 1] + 0.125 y[n − 2] = x[n] − x[n − 1] y[n] = 0.25{x[n − 1] + x[n − 2] + x[n − 3] + x[n − 4]}
数字信号处理实验指导书
实验一DES 综合外设实验1.1实验目的和要求DES320E 提供了键盘，液晶，数码管，直流电机，步进电机，交通灯等外设。
本实验学习这些外设的控制原理。
本实验为大型综合性实验，要求学生掌握DSP编程的基本方法。
通过实验，学生能编写外设控制程序。
例如，使用交通灯和定时器实现十字路口红绿灯的控制，直流电机的调速控制，使用液晶数码管显示和键盘实现计算器等。
1.2实验原理1) C54XX 的I O 空间读写C54XX 提供64K 字的I O 空间访问能力。
在汇编指令中分别提供了读和写命令：portr和p ortw。
你也可以在C中实现该I O 操作，方法如下：首先定义I O 空间变量，如：ioport unsigned portXXXX；/* 其中，XXXX 代表具体I O 口地址*/然后，就可以象访问普通变量一样访问I O 口。
如portXXXX=0x55；/* 将0x55 写到X XXX 指定的I O 口*/2) 交通灯的控制DES320E 提供了 12 个 LED，其控制地址为 IO 空间的 0x0c000h。
该地址的 D0-11比特位分别对应这12 个L ED。
将1写入可以点亮L ED，0 则关闭。
3) 直流电机控制DES320E 实验系统配有一个小型直流电机，可以 DSP 编程完成直流电机的调速控制。
其控制方法为：当向0x0e000h（…VC5402 的I O 空间）的D0 比特位写入1时，电机正向转动；当写入 0 时，电机反向转动。
用户可以通过 D0 位为 1 或 0 的持续时间（即D0 输出方波的占空比）控制电机的转速。
注意，使用直流电机时，应该先接通电机的电源，方法如下：向I O 空间的0x8000 地址的D0 比特位写入1。
若要关闭电源，请写入0。
当写入1或0时，你可以听到继电器动作的声音。
4) 步进电机的控制DES320E 实验系统还配有一个步进电机。
IO 空间的0x0f000h 的D0，D1，D2，D3 四个比特位分别对应步进电机的四相驱动端。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字信号处理实验指导手册西安文理学院机械电子工程系目录实验一离散时间信号 (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实现信号的相加、移位、反转运算。
2、用MATLAB计算信号的能量。
实验二时域采样定理【实验目的】1、掌握奈奎斯特抽样定理的含义，以及在实际应用中需要注意的问题。
2、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。
【实验原理】1、理想抽样信号连续信号经过时间间隔为T的均匀采样后，得到的理想抽样信号可表示为其傅里叶变换为可见，一个连续时间信号经过理想抽样后，其频谱将以抽样频率为间隔而重复，这就是频谱产生周期延拓。
如果信号的最高频谱超过，则各周期延拓分量产生频谱的交叠，称为混叠现象。
要想抽样后能够不失真地还原出原信号，则抽样频率必须大于两倍信号谱的最高频率（，这就是奈奎斯特抽样定理。
2、抽样恢复如果满足奈奎斯特抽样定理，即信号谱的最高频率小于折叠频率，则抽样后不产生频谱混叠，有故将通过以下理想低通滤波器：就可得到原信号频谱所以输出端即为原模拟信号【实验内容】用实验的方法验证抽样定理。
对下面连续信号：进行采样，可得到采样序列其中A为幅度因子，为衰减因子，是模拟角频率，T为采用间隔。
这些参数都要在实验过程中由键盘输入，产生不同的和。
产生采样信号序列，使A=144.128,，，图2—1给出了连续信号的幅频特性曲线，由此图可以确定对应采用的采样频率。
1、取采样频率，即T=1ms.观察所得采样的幅频特性和图2-1中的，在折叠频率附近有无明显差别。
应当注意，实验中所得频谱是用序列的傅里叶变换公式求得的，所以在频率度量上存在关系：为数字频率，为模拟频率。
2、改变采样频率，,观察的变化，并作记录（打印曲线）；进一步降低采样频率，，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的曲线。
3、变采样频率，,观察的变化，并作记录（打印曲线）；观察过采样时的曲线，并记录（打印）这时的曲线。
【思考题】采样频率不同时，相应理想采样序列的傅里叶变换频谱的数字频率刻度是否都相同。
它们所对应的模拟频率是否相同？为什么？实验三离散时间系统【实验目的】1、熟悉时域离散系统的时域特性。
2、利用卷积方法观察分析系统的时域特性。
3、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。
【实验原理】离散信号和系统在时域均可用序列来表示。
序列图形给人以形象直观的印象，他可加深我们对信号和系统的时域特征的理解。
本实验将观察分析几种信号及系统的时域特性。
为了在数字计算机上观察分析各种序列的频域特性，通常对在上进行M点采样来观察分析。
对长度为N的有限长序列x(n),有（3-1）其中通常M应该取得大一些，以便观察谱的细节变化。
取模||可绘出幅频特性曲线。
一个时域离散线性非移变系统的输入/输出关系为y(n)=x(n)*h(n)= （3-2）这里，y(n)为系统的输出序列，x(n)为输入序列。
h(n)、x(n)可以无限长，也可以是有限长。
为了计算机绘图观察方便，主要讨论有限长情况。
如果h(n)和x(n)的长度分别为N和M，则y(n)的长度为L=N+M-1。
这样，（3-2）式所描述的卷积运算就是序列移位、相乘和累加的过程。
上述卷积运算也可以在频域实现（3-3）式（3-3式右边的相乘是在各频域点上的频谱值相乘。
【实验内容】一、认真复习离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理。
二、编制实验用相关子程序。
①信号产生子程序，用于产生实验中要用到的下列信号序列：a.单位采样序列：b.矩形序列：②系统单位脉冲响应序列产生子程序。
本实验要用到两种FIR系统。
a.b.③有限长序列线性卷积子程序，用于完成两个给定长度的序列的卷积。
可以直接调用MATLAB语言中的卷积函数。
用于两个有限长序列的卷积，它假定两个序列都从n=0开始。
调用格式如下：y=其中参数x和y是两个已赋值的行向量序列。
三、在完成编制上述子程序的基础上，完成时域离散信号、系统和系统响应分析。
1、观察信号和系统的时域和频域特性；利用线性卷积求得信号通过系统的响应y(n),比较所求响应y(n)和的时域及频域特性，注意它们之间有无差别，绘图说明，并用所学理论解释所得结果。
2、观察系统对信号的响应特性。
利用线性卷积求系统响应y(n),并判断y(n)图形及其非零值序列长度是否与理论结果一致。
调用序列傅里叶变换数值计算子程序，求得，观察特性曲线，定性判断结果的正确性。
改变的长度，取N=5，重复该实验。
注意参数变化的影响，说明变化前后的差异，并解释所得结果。
3、卷积定理的验证。
对信号和分别求出其频谱，按（3-3）式计算：并绘出曲线与直接对y(n)进行傅里叶变换所得幅频特性进行比较，验证时域卷积定理。
【思考题】1、试用MATLAB编写卷积函数y(n)=conv(x,n)。
2、在卷积定理验证的实验中，如果选用不同的频率采样点数M值，例如，选M=10和M=20，分别做序列的傅里叶变换，求得0，1，…，M-1所得结果之间有无差异？为什么？实验四线性卷积与圆周卷积【实验目的】理解离散序列的线性卷积与圆周卷积的原理，比较其相同点和不同点，掌握线性卷积与圆周卷积的计算步骤和计算方法，能熟练使用MATLAB的相关命令。
【实验原理】1．线性卷积设两序列为和，则和的线性卷积定义为卷积和的运算在图形表示上可分为4步（图4-1）：（1）翻褶：先在哑变量坐标m上做出x(m)和h(m),将h(m)以m=0的垂直轴为对称轴翻褶成h(-m)。
（2）移位：将h(-m)移位，即得h(n-m)。
当n为正整数时，右移n位。
当n为负整数时，左移n位。
（3）相乘：再将h(n-m)和x(m)的相同m值的对应点值相乘。
（4）相加：把以上所有对应点的乘积叠加起来，即得y(n)值。
注意：对于得到的结果仍然是一个序列，若x(n)的长度是N，h(n)的长度是M，则y(n)的长度是N+M-1。
2．圆周卷积设两序列为x(n)和h(n)，则x(n)和h(n)的圆周卷积和定义为=圆周卷积过程（图4-2）：（1）补零：若x(n)的长度是N，h(n)的长度是M，取H>=max(N,M)，对序列补零至点H。
（2）周期延拓：先在哑变量坐标m上做出x(m)和h(m)，将h(m)周期延拓。
（3）翻褶，取主值序列：对h(m)以m=0的垂直轴为对称轴翻褶成h(-m)，然后取主值序列。
（4）圆周移位：对得到的序列进行圆周移位。
（5）相乘相加：与x(m)对应项相乘，并累加，得到y(n)。
3．线性卷积与圆周卷积的关系时域的圆周卷积在频域上相当于两序列的DFT相乘，因而可以采用DFT的快速算法。
它与线性卷积相比，计算速度可以大大加快。
但是，一般实际问题都是线性卷积运算。
4.结论若两序列为x(n)和h(n),长度分别为N、M；则其线性卷积的长度为N+M-1，而圆周卷积的长度为K=max(N,M)。
只有L>=N1+N2-1时，则L点圆周卷积和线性卷积相等。
【实验内容】（1）计算下面信号的线性卷积、圆周卷积，比较什么情况下二者相等。