实验1-多种离散信号产生

电 子 科 技 大 学

实 验 报 告

学生姓名： 学 号： 指导教师：

一、实验室名称：数字信号处理实验室 二、实验项目名称：多种离散时间信号的产生 三、实验原理：

1、基本离散时间信号

利用MATLAB 强大的数值处理工具来实现信号的分析和处理，首先就是要学会应用MATLAB 函数来构成信号。常见的基本信号可以简要归纳如下： （1）．单位采样序列

在MATLAB 中可以利用zeros()函数实现。

如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： （2）．单位阶跃序列

在MATLAB 中可以利用ones()函数实现。 （3）．正弦序列

采用MATLAB 的实现方法，如： （4）．实指数序列

其中，A 、a 为实数。采用MATLAB 的实现方法，如： （5）．复指数序列

采用MATLAB 的实现方法，如：

为了画出复数信号x [n ]，必须要分别画出实部和虚部，或者幅值和相角。MATLAB 函数real 、

imag 、abs 和angle 可以逐次计算出一个复数向量的这些函数。 2、基本数字调制信号

（1）．二进制振幅键控（2ASK ）

最简单的数字调制技术是振幅键控（ASK ），即二进制信息信号直接调制模拟载波的振幅。二进制幅度键控信号的时域表达式：

其中，a n 为要调制的二进制信号，g (t)是单极性脉冲信号的时间波形，Ts 表示调制的信号间隔。典型波形如下：

图 1 – 1 二进制振幅键控信号时间波形

（2）．二进制频移键控（2FSK ）

在二进制数字调制中，若正弦载波的频率随二进制基带信号在f 1和f 2两个频率点间变化，则产

生二进制移频键控信号（2FSK 信号）。二进制频域键控已调信号的时域表达式为:

这里，22112,2f f πωπω==，的反码是n n a a 。典型波形如下：

（3） 1 和 0。

Y =

为载波频率(Hz)；Fd 为

METHOD 是可

以选择的调制方式：ask 、psk、qask、fsk、msk等。

3、双音多频DTMF信号

DTMF（Double Tone MulitiFrequency,双音多频）作为实现电话号码快速可靠传输的一种技术，它具有很强的抗干扰能力和较高的传输速度，因此，可广泛用于电话通信系统中。但绝大部分是用作电话的音频拨号。另外，它也可以在数据通信系统中广泛地用来实现各种数据流和语音等信息的远程传输。

DTMF是用两个特定的单音频组合信号来代表数字信号以实现其功能的一种编码技术。两个单音频的频率不同，代表的数字或实现的功能也不同。这种电话机中通常有16个按键，其中有10个数字键0～9和6个功能键*、#、A、B、C、D。由于按照组合原理，一般应有8种不同的单音频信号。因此可采用的频率也有8种，故称之为多频，又因它采用从8种频率中任意抽出2种进行组合来进行编码，所以又称之为“8中取2”的编码技术。

根据CCITT的建议，国际上采用的多种频率为697Hz、770Hz、852Hz、941Hz、1209Hz、1336Hz、1477Hz和1633Hz等8种。用这8种频率可形成16种不同的组合，从而代表16种不同的数字或功能键，具体组合见表1-1。

表1-1 双音多频的组合功能

因此，DTMF信号可以看作两个有限长度的正弦序列相加，正弦信号的频率由按键数字或字母符号对应的频率决定。如，数字“8”由行频852Hz和列频1336Hz决定。

四、实验目的：

1、 掌握几种基本的离散时间信号（包括单位采样序列，单位阶跃序列，单频正弦序列，单频复指数序列，实指数序列等）。

2、 能够熟练利用MATLAB 产生这些基本的离散时间信号。

3、 理解双音多频DTMF 信号、ASK 、FSK 、BPSK 等信号的产生原理。

4、 学习并运用MATLAB 产生各种通信中的调制信号及双音多频信号。

五、实验内容：

1、对几种基本离散时间信号（包括单位采样序列，单位阶跃序列，正弦序列，复指数序列，实指数序列等）在MATLAB 中编程产生。

2、（拓展要求）利用MATLAB 编程产生2ASK ，2FSK ，2PSK 等数字调制信号。

3、（拓展要求）利用MATLAB 编程产生理解双音多频DTFM 信号。

4、（拓展要求）利用MATLAB 编程产生高斯白噪声序列。

5、（拓展要求）利用MATLAB 中的谱分析函数对正弦信号的频谱进行分析。

6、通过硬件（DSP ）实验箱演示上述信号的时域（示波器）波形与频域波形（计算结果）。

六、实验器材（设备、元器件）：

安装MATLAB 软件的PC 机一台，DSP 实验演示系统一套。

七、实验步骤：

1、 在2020≤≤-n 内，画出单位下列信号：

(a ).单位采样序列][][1n n x δ=和单位阶跃序列][][2n u n x =的时域波形图。

(b ).]5[][11+=n x n y 、]8[][22-=n x n y 的波形。说明][1n x 与][1n y 、][2n x 与][2n y 之间的关系。

2、画出下列信号在1000≤≤n 内的波形。

观察][][][543n x n x n x 、、是否周期信号。如果是周期信号，信号的基波周期是什么？如果不是周期信号，说明原因。

3、在300≤≤n 内，画出下列信号：

对于复数序列，要求分别画出实部和虚部；幅值和相角。若把][6n x 中的底数0.8分别改为1.2、-0.8，讨论产生的时域波形有何变化。总结指数序列的底数对序列变化的影响。 4、（拓展要求）设计产生数字二进制序列：1 0 1 0 1 0 的2ASK 、2FSK 、2PSK 调制信号。

已知符号速率Fd ＝10Hz （即时间间隔Ts 为0.1），输出信号的采样频率为20Hz 。 (a ).2ASK 信号的载波频率Fc ＝5Hz ，

(b ).2FSK 信号载波1频率F1＝5Hz ，载波2频率F2＝1Hz 。 (c ).2PSK 载波频率Fc ＝1Hz 。

分别画出以上信号调制前后的时域波形图。

5、 （拓展要求）利用MATLAB 产生DTMF 双音多频信号。画出数字“0”的时域波形图。

6、 （拓展要求）MATLAB 函数randn(1,N)可以产生均值为0，方差为1的高斯随机序列，

也就是白噪声序列。试利用randn 函数产生均值为0.15，方差为0.1的高斯白噪声序列

][8n x ，要求序列时域范围为1000≤≤n 。画出时域波形图。同时将实验步骤2中产生的信

号][2n x 与][8n x 相加，将得到的波形与][2n x 的波形做比较。

7、 （拓展要求）利用MATLAB 中的谱分析函数画出][][][543n x n x n x 、、的频谱。与理论上根

据傅立叶变换的定义计算出的][][][543n x n x n x 、、的频谱进行比较。

8、通过硬件（DSP ）实验箱演示上述信号的时域（示波器）波形与频域波形（计算结果）。

八、实验数据及结果分析：

注：本次实验在寝室电脑上完成，所用MATLAB 版本为MATLAB R2010b 程序：

（1）产生][][][][][][][][][765432121n x n x n x n x n x n y n y n x n x 、、、、、、、、序列的程序 ? 产生][1n x ，][2n x ，][1n y ，][2n y 的程序：

产生][][][543n x n x n x 、、的程序：

（2）产生2ASK、2FSK、2PSK调制信号的程序（拓展要求）

（3）产生DTMF信号的程序（拓展要求）

（4）高斯白噪声序列的产生程序（扩展要求）

根据随机变量函数的分布性质生成均值为0.15，方差为0.1的高斯白噪声序列：

如果正态随机变量),(~2σμN X ，则其线性函数)0(≠+=a b a X Y 也服从正态分布且

))(,(~2σμa b a N Y +.

注：为了便于观察，这里用信号x3[n]而没用x2[n]。

（5）正弦信号频谱分析的程序（扩展要求）

结果：

（1）][][][][][][][][][765432121n x n x n x n x n x n y n y n x n x 、、、、、、、、的时域波形

][1n x ： ][2n x ： ][1n y ： ][2n y ： ][3n x ：

][4n x ： ][5n x ： ][6n x ： ][7n x ：

实部： 虚部： 幅值： 相角：

（2）信号的时移：][1n x 与][1n y 、][2n x 与][2n y 之间的关系。

][1n x 与][1n y ：

][1n y 由][1n x 向左平移5个单位得到。 ][2n x 与][2n y ：

][2n y 由][2n x 向右平移8个单位得到。

（2）正弦序列][][][543n x n x n x 、、周期的判断

][3n x 是周期信号，信号的基波周期是32。

][4n x 是周期信号，信号的基波周期是π4。 ][5n x 是周期信号，信号的基波周期是48。

（3）指数序列底数与序列变化的关系总结

若把][6n x 中的底数0.8改为1.2，时域波形就呈现指数增长；若改为-0.8，就呈现振荡衰减，最后趋于0.

总结：对形如n A n x α=][的序列，

若底数α的绝对值小于1，

? 若底数为正，时域波形呈指数衰减，幅值为正； ? 若底数为负，时域波形呈指数衰减，但是振荡的。 若底数α的绝对值大于1，

? 若底数为正，时域波形呈指数增长，幅值为正； ? 若底数为负，时域波形呈指数增长，但是振荡的。 （4）2ASK 、2FSK 、2PSK 调制信号时域波形（拓展要求） 注：学生觉得要求中载波频率太小，就把载波频率都扩大了10倍。 2ASK ： 2FSK ： 2PSK ：

（5）数字“0” DTMF 时域波形（拓展要求） 取样频率为100Hz

（6）高斯白噪声序列的时域波形；正弦序列加上高斯白噪声后的时域波形（扩展要求） 注：为了观察方便，这里用信号x3[n]而没用x2[n]

（6）正弦序列的频谱图（扩展要求）

九、实验结论：

1. 利用MATLAB的zeros()和ones(),等函数可以方便地产生各种基本的离散时间信号。

2. 从MATLAB模拟调制的波形图看，2ASK调制中，码元1的已调波形即为载波幅度，而码元0无幅度，2FSK调制中码元1和0的已调波形频率不同，2PSK调制中码元1和0的已调波形初始相位不同，相差 。

3. 双音多频信号可以看作是两个有限长度的正弦信号的叠加。如数字“0”由行频941Hz和列频1336Hz决定。

十、总结及心得体会：

通过实验，练习了MATLAB的使用，练习了一些基本离散时间信号的生成，实现了3种模拟调制，结合实际运用了通信原理的一些知识，并了解双音多频DTMF信号的产生原理，既巩固了已学到的知识，又有不少新的收获，拓展了动手能力。

十一、对本实验过程及方法、手段的改进建议：

1.建议将模拟调制实验中的载波频率设置得大一些；

2.建议将实验内容4的原始信号设为正弦信号x3[n]，这样更便于观察。

3.建议适当调整实验时间段，不要和考试时间段重合，让同学们有机会认真研究实验内容，并

有时间和老师交流。

报告评分：

指导教师签字：

实验1离散时间信号的产生与运算

数字信号处理 实验报告 班级： 学号： 姓名：word文档可自由复制编辑

实验1离散时间信号的产生与运算 一、实验目的 (1)了解离散时间信号的特点。 (2)掌握在计算机中生成及绘制各种常用离散时间信号序列的方法。 (3)掌握序列的加、减、乘、除和平移、反转、尺度变换等基本运算及计算机的 实现方法。 二、实验原理 信号是随时间变化的物理量，而计算机只能处理离散信号。离散信号是在某些不连续的时间上有信号值，而在其它时间点上没有定义的一类信号。离散信号一般可以由连续信号通过模数转换得到。 常用的离散信号有单位脉冲序列、单位阶跃序列、复指数序列、正弦信号序列、随机序列等。 离散信号的基本运算包括信号的加、减、乘、除。离散信号的时域变换包括信号的平移、反转、尺度变换等。 三、实验内容与方法 1、编写程序，生成如下数字信号：sqrt(2*k)u(k错误！未找到引用源。3)， δ(k+5)。 (1) f(k)=sqrt(2*k)u(k错误！未找到引用源。3) 代码： k=(1:10); n=3; u=[(k-n)>=0]; a=sqrt(2*k); stem(k,a.*u); title('sqrt(2*k)u(k 3)的图像'); xlabel('时间(k)');ylabel('幅值f(k)'); 运行图： word文档可自由复制编辑

(2) f(k)= δ(k+5) 代码： k1=-10;k2=0;k=k1:k2; n=-5; %单位脉冲出现的位置 f=[(k-n)==0]; stem(k,f,'filled');title('δ(k+5)序列的图像') xlabel('时间(k)');ylabel('幅值f(k)'); 运行图： word文档可自由复制编辑

信号与系统实验

《信号与系统及MATLAB实现》实验指导书

前言 长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。

信号与测试技术实验一

实验一基本信号分析实验报告 一实验目的 1掌握基本信号的时域和频域分析方法； 2掌握信号的自相关和互相关分析，了解其应用。 二实验内容与图像结果分析 （1）产生不同的周期信号，包括正弦信号、方波信号、锯齿波，在时域分析这些波形特征（幅值、频率（周期））。 （2）在Matlab中产生不同的非周期信号，包括随机噪声、阶跃信号、矩形脉冲。（3）对产生的信号进行Fourier变换，从频率域分析信号的特征，并说明方波信号和锯齿波信号的信号带宽； 从图中可以看到，正弦信号基频为10rad/s，因此其Fourier变换在w=10处出现了峰值，而方波信号依据佛利叶级数展开可知是由一系列不同频率的正弦波构成，基频是w=10，基频的幅值最大，同时其他频率为基频的整数倍（不含20,40…）,且幅值依次减少。

锯齿波信号的基频为w=10，因此傅里叶级数展开同样在10出出现了峰值，而其他出现的依次是基频的整数倍，且幅值依次减少。由于随机噪声信号是随机信号，不具有规律性，因此在傅里叶变换后我们可以看到它含有各个频率的谐波。 阶跃信号的傅里叶变换为冲击函数。矩形信号为非周期信号，因此它的傅里叶变换为连续函数，频率在各处均有分布。 （4）产生复合信号：由3个不同频率、幅值的正弦信号叠加的信号，从图形上判断信号的特征； 产生由正弦信号和随机信号叠加的混合信号，从图形上判断信号的特征；产生由正弦信号和方波叠加的信号，从图形上判断信号的特征。 （5）对（4）中的3种复合信号进行FFT计算，从图上判断信号的特征。

三种不同幅值、频率的正弦信号叠加后，在时域图上我们看不出很有规律性的东西，然而进行傅里叶变换后，放到频域图之后，我们可以很清楚的看到叠加信号的组成规律，在三个频率出现了峰值。正弦信号叠加随机噪声，我们在时域图上也看不到很明显的规律特征，进行傅里叶变换后，我们看到时域图上在一处出现了峰值，则这个频率处实际就是正弦信号的频率。正弦信号叠加方波信号在时域图中同样规律不明显，在进行傅里叶变换后，在频域图上我们看到有两处峰值，这两个频率实际就是正弦波的频率和方波的基频信号，其余较小的为方波的谐波信号。 由此可以看出，通过傅里叶变换，将时域波形变换到频域波形，更加有助于我们分析信号的本质特征，也有利于从噪声信号出提取有用的信号。

matlab实现：常见的离散时间信号

1． 单位抽样序列，或称为离散时间冲激，单位冲激： ? ??=01)(n δ 00≠=n n 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-01)(k n δ 0≠=n k n 2．单位阶跃序列 ? ??01)(n u 00<≥n n 在MATLAB 中可以利用ones( )函数实现。 );,1(N ones x = 3．正弦序列 )(cos )(0φω+=n A n x 这里， ,,0ωA 和φ都是实数，它们分别称为本正弦信号)(n x 的振幅，角频率和初始相位。 πω200=f 为频率。 4．复正弦序列 n j e n x ω=)( 5．实指数序列 n A n x α=)( 6. 随机序列 长度为N 的随机序列 基本数学函数参考教材P69页以及随后的使用说明。 注意使用行向量，特别是冒号运算符。 举例，长度为N 的实指数序列在MATLAB 中实现： n a x N n .^1 :0=-= 1. 单位采样 长度为N 的单位采样序列u(n)可以通过下面的MATLAB 命令获得：

u=[1 )1,1(-N zeros ]； 延迟M 个采样点的长度为N 的单位采样序列ud(n)(M

信号与系统实验报告一

1. 实验原理 2. 设描述连续时间系统的微分方程为： ) ()()()()()()()(01) 1(1) (01)1(1)(t f b t f b t f b t f b t y a t y a t y a t y a m m m m n n n n +'+++=+'+++---- 则可用向量a 和b 表示该系统，即 ],,,,[011a a a a a n n -= ],,,,[011b b b b b m m -= 注意，向量a 和b 的元素一定要以微分方程时间求导的降幂次序排列，且缺项要用0补齐。 如微分方程 )()()(2)(3)(t f t f t y t y t y +''=+'+'' 表示该系统的向量为 ]2 3 1[=a ]1 0 1[=b （1）求解冲激响应：impulse()函数 impulse()函数有以下四种调用格式： ① impulse(b,a) 该调用格式以默认方式绘制由向量a 和b 定义的连续时间系统的冲激响应的时域波形。 ② impulse(b,a,t) 该调用格式绘制由向量a 和b 定义的连续时间系统在t ~0时间范围内的冲激响应的时域波形。 ③ impulse(b,a, t1:p:t2) 该调用格式绘制由向量a 和b 定义的连续时间系统在21~t t 时间范围内，且以时间间隔 p 均匀抽样的冲激响应的时域波形。 ④ y=impulse(b,a,t1:p:t2) 该调用格式并不绘制系统冲激响应的波形，而是求出由向量a 和b 定义的连续时间系统在21~t t 时间范围内以时间间隔p 均匀抽样的系统冲激响应的数值解。 （2）求解阶跃响应：step()函数 step()函数也有四种调用格式： ① step(b,a) ② step(b,a,t) ③ step(b,a, t1:p:t2) ④ y=step(b,a,t1:p:t2) 上述调用格式的功能与impulse()函数完全相同。 （3）求解零状态响应：lsim()函数 lsim()函数有以下二种调用格式：

实验1-多种离散时间信号产生(答案)

实 验 报 告 学生姓名： 学 号： 指导教师： 一、实验室名称：数字信号处理实验室 二、实验项目名称：多种离散时间信号的产生 三、实验原理： 1、基本离散时间信号 利用MATLAB 强大的数值处理工具来实现信号的分析和处理，首先就是要学会应用MATLAB 函数来构成信号。常见的基本信号可以简要归纳如下： （1）．单位采样序列 ???=0 1 )(n δ ≠=n n 在MATLAB 中可以利用zeros()函数实现。 ; 1)1(); ,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-0 1 )(k n δ ≠=n k n （2）．单位阶跃序列 ???=0 1 )(n u 00<≥n n 在MATLAB 中可以利用ones()函数实现。 );,1(N ones x = （3）．正弦序列

)2sin()(?π+=fn A n x 采用MATLAB 的实现方法，如： ) ***2sin(*1 :0?+=-=n f pi A x N n （4）．实指数序列 n a A n x ?=)( 其中，A 、a 为实数。采用MATLAB 的实现方法，如： n a x N n .^1 :0=-= （5）．复指数序列 n j e A n x )(0)(ωσ+?= 采用MATLAB 的实现方法，如： ) *)*exp((*1 :00n j A x N n ωσ+=-= 为了画出复数信号x [n ]，必须要分别画出实部和虚部，或者幅值和相角。 MATLAB 函数real 、imag 、abs 和angle 可以逐次计算出一个复数向量的这些函数。 2、基本数字调制信号 （1）．二进制振幅键控（2ASK ） 最简单的数字调制技术是振幅键控（ASK ），即二进制信息信号直接调制模拟载波的振幅。二进制幅度键控信号的时域表达式： ∑-=n c s n ASK t nT t g a t S ωcos )]([)( 其中，a n 为要调制的二进制信号，g (t)是单极性脉冲信号的时间波形，Ts 表示调制的信号间隔。典型波形如下：

信号与系统实验报告_1(常用信号的分类与观察)

实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：

图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)

离散时间信号的表示及运算

第2章 离散时间信号的表示及运算 2.1 实验目的 ● 学会运用MATLAB 表示的常用离散时间信号； ● 学会运用MATLAB 实现离散时间信号的基本运算。 2.2 实验原理及实例分析 2.2.1 离散时间信号在MATLAB 中的表示 离散时间信号是指在离散时刻才有定义的信号，简称离散信号，或者序列。离散序列通常用)(n x 来表示，自变量必须是整数。 离散时间信号的波形绘制在MATLAB 中一般用stem 函数。stem 函数的基本用法和plot 函数一样，它绘制的波形图的每个样本点上有一个小圆圈，默认是空心的。如果要实心，需使用参数“fill”、“filled”，或者参数“.”。由于MA TLAB 中矩阵元素的个数有限，所以MA TLAB 只能表示一定时间范围内有限长度的序列；而对于无限序列，也只能在一定时间范围内表示出来。类似于连续时间信号，离散时间信号也有一些典型的离散时间信号。 1. 单位取样序列 单位取样序列)(n δ，也称为单位冲激序列，定义为 ) 0()0(0 1)(≠=?? ?=n n n δ （12-1） 要注意，单位冲激序列不是单位冲激函数的简单离散抽样，它在n =0处是取确定的值1。在MATLAB 中，冲激序列可以通过编写以下的impDT .m 文件来实现，即 function y=impDT(n) y=(n==0); %当参数为0时冲激为1,否则为0 调用该函数时n 必须为整数或整数向量。 【实例2-1】 利用MATLAB 的impDT 函数绘出单位冲激序列的波形图。 解：MATLAB 源程序为 >>n=-3:3; >>x=impDT(n); >>stem(n,x,'fill'),xlabel('n'),grid on >>title('单位冲激序列') >>axis([-3 3 -0.1 1.1])

信号与系统实验三

信号与系统实验实验三：信号的卷积 小组成员： 黄涛13084220 胡焰焰13084219 洪燕东13084217

一、实验目的 1. 理解卷积的物理意义； 2. 掌握运用计算机进行卷积运算的原理和方法； 3. 熟悉卷积运算函数conv的应用； 二、预习内容 1. 卷积的定义及物理意义； 2. 卷积计算的图解法； 3. 卷积的应用 三、实验原理说明 1.卷积的定义 连续时间和离散时间卷积的定义分别如下所示： 2.卷积的计算 由于计算机技术的发展，通过编程的方法来计算卷积积分和卷积和已经不再是冗繁的工作，并可以获得足够的精度，因此信号的时域卷积分析法在系统分析中得到了广泛的应用。 卷积积分的数值运算可以应用信号的分段求和来实现，即： 数值运算只求当时的信号值，则由上式可以得到： 上式中实际上就是连续信号等间隔均匀抽样的离散序列的卷积和，当足够小的时候就是信号卷积积分的数值近似。因此，在利用计算机计算两信号卷积积分时，实质上是先将其转化为离散序列，再利用离散卷积和计算原理来计算。 3.卷积的应用 3.1 求解系统响应 卷积是信号与系统时域分析的基本手段，主要应用于求解系统响应，已知一LTI系统的单位冲激响应和系统激励信号则系统响应为激励与单位冲激响应的卷积。 需要注意的是利用卷积分析方法求得的系统响应为零状态响应。 3.2 相关性分析 相关函数是描述两个信号相似程度的量。两信号之间的相关函数一般称之为互相关函数或者互关函数，定义如下： 若是同一信号，此时相关函数称为自相关函数或者自关函数： 对于相关函数与卷积运算有着密切的联系，由卷积公式与相关函数比较得： 可见，由第二个信号反转再与第一个信号卷积即得到两信号的相关函数。 4.涉及的Matlab函数 4.1 conv函数 格式w = conv(u,v)，可以实现两个有限长输入序列u，v的卷积运算，得到有限冲激响应系统的输出序列。输出序列长度为两个输入序列长度和减一。 四、实验内容 给定如下因果线性时不变系统： y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3] (1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本; clc; N = 0:19;

信号与系统实验二

实验二 常用信号分类与观察 一、实验目的 1、观察常用信号的波形特点及产生方法。 2、学会使用示波器对常用波形参数的测量。 二、实验内容 1、信号的种类相当的多，这里列出了几种典型的信号，便于观察。 2、这些信号可以应用到后面的“基本运算单元”和“无失真传输系统分析”中。 三、实验仪器 1、信号与系统实验箱一台（主板）。 2、20MHz 双踪示波器一台。 四、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。 1、正弦信号：其表达式为)sin()(θω+=t K t f ，其信号的参数：振幅K 、角频率ω、与初始相位θ。其波形如下图所示： 图 1-5-1 正弦信号 2、指数信号：指数信号可表示为at Ke t f =)(。对于不同的a 取值，其波形表现为不同的形式，如下图所示：

图 1-5-2 指数信号 3、指数衰减正弦信号：其表达式为 ?? ???><=-)0()sin()0(0)(t t Ke t t f at ω 其波形如下图： 图 1-5-3 指数衰减正弦信号 4、抽样信号：其表达式为： sin ()t Sa t t = 。)(t Sa 是一个偶函数，t = ±π,±2π,…,±n π时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

信号与系统实验DOC

信号与系统实验讲义 雷明东编 重庆文理学院 电子电气学院 2014年10月

实验注意事项 1、不准迟到早退，开始做实验前需要签字； 2、在离开实验室前，要整理好实验设备、桌椅、收拾好垃圾后，待老师检查完毕，方可离开实验室； 3、做实验期间不准大声喧哗，如有问题需举手示意； 4、不准在无老师授权的情况下随意拆卸实验设备； 5、在每次做新实验前，需交前个实验的实验报告。

实验一 常用信号的分类和观察 一 实验目的： 1、观察和了解常见信号的波形和特点。 2、理解相关信号参数的作用和意义。 3、掌握信号的FFT 变换。 3、熟练掌握示波器的使用。 二 实验原理： 描述信号的基本方法是写出它的数学表达式，此表达式是时间的函数，绘出函数的图像称为信号的波形。 对于各种信号，可以从不同的角度分类。如分成确定性信号与随机信号；周期信号与非周期信号；连续时间信号与离散时间信号等。 常见信号除了包括正弦波)sin()(0φω+=t A t x 、指数函数信号t Ke t x α=)(、抽样函数信号t t A t x /)(sin )(=、高斯函数信号τ/)(t Ke t x -=、方波、三角波、锯齿波，还包括一些直流信号。 三 预习练习： 1、预习有关信号的分类和描述。 2、理解信号的函数表达式和相关参数的意义。 四 实验内容及步骤： 1、 根据实验箱上函数信号发生器模块的提示选择相应的信号波形代码。 01：正弦波 02：方波 03：锯齿波 04：三角波

05：阶梯波 06：衰减指数信号 07：高斯函数信号 08：抽样函数信号 09：抽样脉冲 10：调幅信号 11：扫频信号 2、用示波器测量信号，读取信号的幅度和频率，并用坐标纸记录信号波形； 在信号与系统实验箱上的电源模块用电压表（或万用表）与示波器来观 测电源信号的特点，并测量电源的幅度。 3、在示波器上观测扫频信号的波形特征，大致画出扫频信号的波形。 4、利用示波器中的FFT函数，来观看信号的FFT变换形式。 5、用频谱分析仪观测各个信号的频谱（选做）。 五实验仪器： 1、信号系统实验箱（函数信号发生器模块） 2、双踪示波器 六实验报告内容： 1、根据实验测量所得数据，绘制各个信号的波形图。 2、绘制各个波形的FFT变换波形。 3、写出相应的函数表达式与频域变换表达式。 4、用示波器直流档观测函数信号的波形特点，并说明原因（提示：本函数发生器所产生的信号均由单片机AT89C51产生）。

信号与测试实验1时率与频率

基本信号分析 一、实验目的 1.掌握基本信号的时域和频域分析方法 2.掌握信号的自相关和互相关分析，了解其应用 二、数据处理与分析 （1）幅值为1，频率为100Hz的正弦信号，上图为时域图，下图为利用快速傅里叶变换获得的频谱图。从频谱图上看出，f=100Hz时频域的幅值最大。 （2）频域为100Hz，幅值为1的方波信号，上图为时域图，下图为借助快速傅立叶变换获得的频域图。从频谱图上看出，f=100Hz时频域的幅值最大，随着频域增大，频域的幅值逐渐衰减。

（3）频率为100Hz，幅值为1的锯齿波信号图，上图为时域图，下图为借助傅立叶变换而获得的频域图。从频域图看出，在100Hz的整数倍频率上，频域幅值都出现了峰值，随着频率的增大，峰值逐渐收敛至0. （4）平均振幅为1的噪声信号，上图为时域图，下图为通过快速傅立叶变

换得出的频谱图，从频谱图可以看出，白噪声信号的频谱杂乱无章，无明显规律。 （5）由频率为50Hz、100Hz、150Hz的正弦信号组成的复合信号，上图为时域图，下图为频域图，从图中可以看出，频谱图在50、100、150Hz处出现了峰值。 （6）频率为100Hz 的正弦信号叠加噪声信号：上图为时域信号图，下图为

通过快速傅立叶变换获得的频谱图。与没有叠加噪声信号的正弦波相比，时域波形出现了毛刺，而频谱图中除了在100Hz处有峰值外，在其他频率点处也出现了一些较低的峰值。 （7）频率为100Hz的正弦信号和频率为100Hz的方波信号进行叠加，上图为时域信号，下图为频谱图。从时域图上可以看出，正弦波形叠加方波后有了明显的畸变。从频谱图上可以看出，除了100Hz处出现峰值以外，在其他频率点也出现了一些峰值。

离散时间信号表与运算

离散时间信号表与运算

————————————————————————————————作者：————————————————————————————————日期：

实验一 离散时间信号的表示与运算 一 实验目的 1、熟悉MATLAB 的绘图函数； 2、掌握单位取样序列、单位阶跃序列、矩形序列和正余弦序列的产生方法； 3、掌握离散时间信号基本运算的MATLAB 实现； 4、掌握离散时间信号线性卷积和运算的MATLAB 实现。 二 实验设备 1、计算机 2、MA TLAB R2007a 仿真软件 三 实验原理 1）序列相加和相乘 设有序列)(1n x 和)(2n x ，它们相加和相乘如下： ) ()()()()()(2121n x n x n x n x n x n x ?=+= 注意，序列相加（相乘）是对应序列值之间的相加（相乘），因此参加运算的两个序列必须具有相同的长度，并且保证位置相对应。如果不相同，在运算前应采用zeros 函数将序列左右补零使其长度相等并且位置相对应。在MATLAB 中，设序列用x1和x2表示，序列相加的语句为：x=x1+x2；然而要注意，序列相乘不能直接用x=x1*x2,该式表示两个矩阵的相乘，而不是对应项的相乘。对应项之间相乘的实现形式是点乘“.*”，实现语句为：x=x1.*x2。 2）序列翻转 设有序列：)()(n x n y -=，在翻转运算中，序列的每个值以n=0为中心进行翻转，需要注意的是翻转过程中序列的样值向量翻转的同时，位置向量翻转并取反。MATLAB 中，翻转运算用fliplr 函数实现。设序列)(n x 用样值向量x 和位置向量nx 表述，翻转后的序列 )(n y 用样值向量y 和位置向量ny 描述。 3）序列的移位 移位序列)(n x 的移位序列可表示为：)()(0n n x n y -=，其中，00>n 时代表序列右移 0n 个单位；00

测试信号分析与处理作业实验一二

王锋 实验一：利用FFT 作快速相关估计 一、实验目的 a.掌握信号处理的一般方法，了解相关估计在信号分析与处理中的作用。 b.熟悉FFT算法程序；熟练掌握用FFT作快速相关估计的算法。 c.了解快速相关估计的谱分布的情况。 二、实验内容 a.读入实验数据[1]。 b.编写一利用FFT作相关估计的程序[2]。 c.将计算结果表示成图形的形式，给出相关谱的分布情况图。 注[1]：实验数据文件名为“Qjt.dat”。 实验数据来源：三峡前期工程 “覃家沱大桥” 实测桥梁振动数据。 实验数据采样频率：50Hz。 可从数据文件中任意截取几段数据进行分析，数据长度N 自定。 注[2]：采用Matlab 编程。 三、算法讨论及分析 算法为有偏估计，利用FFT计算相关函数 Step 1: 对原序列补N个零，得新序列x2N(n) Step2: 作FFT[x2N(n)]得到X2N(k) Step 3: 取X2N(k)的共轭，得 Step 4: 作 Step 5: 调整与的错位。 四、实验结果分析 1. 该信号可以近似为平稳信号么？ 可以近似为平稳信号，随机过程的统计特性不随样本的采样时刻而发生变化。取N=8192，分别取间隔m=500,m=700,m=1000,所得到的均值均为0.5366，方差为47369，与时间无关。

图1-1 自相关函数图 (上图表示的R0，下图为调整后的R0) 2. 该信号是否具有周期性，信噪比如何？ >> load Qjt.dat; %加载数据 N=32768; %数据长度 i=1:1:N; %提取数据 plot(i,Qjt(i)); 抛去几个极值点，从图1-2可以看出，数据具有一定的周期性，杂音比较少，说明信噪比较高。 图1-2 数据图

离散时间信号分析

离散时间信号分析 实验目的：利用MA TLAB进行离散时间序列的基本运算，掌握基本的MA TLAB函数的编写和调试方法。 实验内容： （1）信号相加 x(n)=x1(n)+x2(n) 当两个相加的序列长度不同时或位置不对应时，首先必须调整二者的位置对齐，然后通过zeros函数左右补零使其长度相等后再相加。下面的参考代码利用函数sigadd说明了这些运算，其验证将在后续实验中进行。 MATLAB参考代码 function[y,n]=sigadd(x1,n1,x2,n2) %implements y(n)=x1(n)+x2(n) %--------------------------------------------- %[y,n]=sigadd(x1,n1,x2,n2) %y=sum sequence over n,which includes n1 and n2 %x1=first sequence over n1 %x2=second sequence over n2(n2 can be different from n1) % n=min(min(n1),min(n2)):max(max(n1),max(n2));%duration of y(n) y1=zeros(1,length(n)); y2=y1; y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;%x1 with duration of y y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;%x2 with duration of y y=y1+y2;%sequence addition （2）信号相乘 信号相乘，即两个序列的乘积（或称“点乘”），表达式为： x(n)=x1(n)?x2(n) 在MA TLAB中，用运算符“.*”实现。

北京理工大学信号与系统实验实验报告

实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号

t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。

电子科大实验1-多种离散时间信号产生答案

实 验 报 告 学生姓名： 学 号： 指导教师： 一、实验室名称：数字信号处理实验室 二、实验项目名称：多种离散时间信号的产生 三、实验原理： 1、基本离散时间信号 利用MATLAB 强大的数值处理工具来实现信号的分析和处理，首先就是要学会应用MATLAB 函数来构成信号。常见的基本信号可以简要归纳如下： （1）．单位采样序列 ? ??=01)(n δ 00 ≠=n n 在MATLAB 中可以利用zeros()函数实现。 ; 1)1();,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即：

???=-0 1)(k n δ ≠=n k n （2）．单位阶跃序列 ? ??=01)(n u 00 <≥n n 在MATLAB 中可以利用ones()函数实现。 );,1(N ones x = （3）．正弦序列 )2sin()(?π+=fn A n x 采用MATLAB 的实现方法，如： ) ***2sin(*1 :0?+=-=n f pi A x N n （4）．实指数序列 n a A n x ?=)( 其中，A 、a 为实数。采用MATLAB 的实现方法，如： n a x N n .^1:0=-= （5）．复指数序列 n j e A n x )(0)(ωσ+?= 采用MATLAB 的实现方法，如： ) *)*exp((*1 :00n j A x N n ωσ+=-= 为了画出复数信号x [n ]，必须要分别画出实部和虚部，或者幅值和相角。MATLAB 函数real 、imag 、abs 和angle 可以逐次计算出一个复数向量的这些函数。

信号与系统实验总结及心得体会

信号与系统实验总结及心得体会 2011211204 刘梦颉2011210960 信号与系统是电子信息类专业的一门重要的专业核心基础课程，该课程核心的基本概念、基本理论和分析方法都非常重要，而且系统性、理论性很强，是将学生从电路分析领域引入信号处理与传输领域的关键性课程，为此开设必要的实验对我们加强理解深入掌握基本理论和分析方法，以及对抽象的概念具体化有极大的好处，而且为后续专业课程的学习提供了理论和大量实验知识储备，对以后的学术科研和创新工作都是十分重要的。下面我将从实验总结、心得体会、意见与建议等三方面作以总结。 一．实验总结 本学期我们一共做了四次实验，分别为：信号的分类与观察、非正弦周期信号的频谱分析、信号的抽样与恢复（PAM）和模拟滤波器实验。 1.信号的分类与观察 主要目的是：观察常用信号的波形特点以及产生方法，学会用示波器对常用波形参数进行测量。主要内容是：利用实验箱中的S8模块分别产生正弦信号、指数信号和指数衰减正弦信号，并用示波器观察输出信号的波形，测量信号的各项参数，根据测量值计算信号的表达式，并且与理论值进行比较。 2.非正弦信号的频谱分析 主要目的是：掌握频谱仪的基本工作原理和正确使用方法，掌握非正弦周期信好的测试方法，理解非正弦周期信号频谱的离散性、谐波性欲收敛性。主要内

容是：通过频谱仪观察占空比为50%的方波脉冲的频谱，和占空比为20%的矩形波的频谱，并用坐标纸画图。 3.信号的抽样与恢复 主要目的是：验证抽样定理，观察了解PAM信号的形成过程。主要内容是：通过矩形脉冲对正弦信号进行抽样，再把它恢复还原过来，最后用还原后的图形与原图形进行对比，分析实验并总结。 4.模拟滤波器实验 主要目的是：了解RC无源和有源滤波器的种类、基本结构及其特性，比较无源和有源滤波器的滤波特性，比较不同阶数的滤波器的滤波效果。主要内容：利用点频法通过测试无源低通、高通、带通和有源带阻，以及有源带通滤波器的幅频特性，通过描点画图形象地把它们的特点表现出来。 通过对信号与实验课程的学习，我掌握了一些基本仪器的使用方法，DDS 信号源、实验箱、示波器、频谱仪等四种实验仪器。初步了解了对信号的测试与分析方法对以前在书本上看到的常见信号有了更加具体的认识，使得书本上的知识不再那么抽象。 DDS信号源，也就是函数发生器，可以产生固定波形，如正弦波、方波或三角波，频率和幅度可以调节。实验箱是很多个信号实验装置的集合，可谓集多种功能于一身，其中包括函数发生器、模拟滤波器、函数信号的产生与测量、信号的抽样与恢复等模块。示波器能把抽象的电信号转换成具体的图像，便于人们研究各种电现象的变化过程。利用示波器能观察各种不同的信号幅度随时间变化的波形曲线，还可以用它测试各种不同的电量，如电压、电流、频率、相位差、

信号与系统实验报告

实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10()0 t u t t >?=? 0); %定义函数体，即函数所执行指令

离散时间信号的产生及信号的卷积和运算实验报告2

离散时间信号的产生及信号的卷积和运算 实验报告 班级：___________ 姓名：__________ 学号:____________ 一、实验目的和原理 实验原理： （一）DTFT 和DFT 的定义及其相互关系： 序列x[n] 的DTFT 定义：∑=∞ -∞ =-n jn ωj ω x[n]e )X(e 它是关于自变量ω的复函数，且是以π2为周期的连续函数。)X(e j ω 可以表示为： )(e jX )(e X )X(e j ωim j ωre j ω+= 其中，)(e X j ω re 和)(e X j ωim 分别是)X(e j ω的实部和虚部；还可以表示为： )(ωj j ωj ωe )X(e )X(e θ= 其中，)X(e j ω 和}arg{)()X(e j ω=ωθ分别是)X(e j ω的幅度函数和相位函数； 它们都是ω的实函数，也是以π2为周期的周期函数。 序列x[n]的N 点DFT 定义： ∑∑-=-=-===10 1 22][][)(][N n kn N N n kn N j k N j W n x e n x e X k X ππ ][k X 是周期为N 的序列。 )X(e j ω与][k X 的关系：][k X 是对)X(e j ω在一个周期中的谱的等间隔N 点采样，即： k N j ω)X(e k X πω2| ][= =， 而)X(e j ω 可以通过对][k X 内插获得，即：

]2/)1)][(/2([1 ) 22sin() 22sin( ] [1----=?--= ∑N N k j N k j ω e N k N k N k X N )X(e πωπωπω (二) 线性时不变离散时间系统的变换域表示： LTI 离散时间系统的时域差分方程为： ∑∑==-=-M k k N k k k n x p k n y d )()( (1) 传递函数： 对上面的差分方程两边求z 变换，得： ∑∑∑∑=-=-=-=-=? =N k k k M k k k M k k k N k k k z d z p z X z Y z p z X z d z Y 0 00 ) () ()()( 我们定义LTI 离散时间系统的输出的Z 变换Y(z)与输入的Z 变换X(z)的比值为该系统的传递函数，即) () ()(z X z Y z H = 为系统的传递函数。 N N M M z d z d d z p z p p z D z p z H ----++++++= =......)()()(110110 分解因式 ∏-∏-=∑∑= =-=-=-=-N i i M i i N i i k M i i k z z K z d z p z H 11 11 0)1()1()(λξ ，其中i ξ和i λ称为零、极点。 利用系统的传递函数)(z H ，我们可以分析系统的零极点，稳定性及实现结构等特点。 (2) 频率响应： 因为大多数离散时间信号都可以分解为n j e ω的线性组合，所以研究输入n j e ω-的响应具有极大的意义，即当输入为n j e n x ω=][时，输出为： )()()(][) (ωωωωωj n j m m j n j m n j m e H e e m h e e m h n y === ∑∑∞ -∞ =--∞ -∞ = 这里，∑∞-∞ =-= n n j j e n h e H ωω )()(是h(n)的DTFT ，称为LTI 离散时间系统的频率