电子科技大学信号处理实验1_信号的基本表示及时域分析

电 子 科 技 大 学

实 验 报 告

一、实验室名称：数字信号处理实验室

二、实验项目名称：信号的基本表示及时域分析 三、实验原理：

1、MATLAB基础

MATLAB（Matrix+Laboratory）是MathWorks公司1982年推出的数值计算和可视化数学软件。其基本元素是矩阵，基本方法是程序实现。使用MATLAB可以方便地解决复杂的数值计算问题，并且具有丰富的工具箱（ToolBox)专用于辅助各类与数学相关问题的分析。在高校中，MATLAB已成为数学，信息，控制等诸多学科有关课程的有效教学工具。 1.1 MATLAB功能特点

MATLAB命令可控、可编程，有上百个预先定义好的命令和函数，这些函数能通过用户自定义函数进一步扩展。MATLAB本身带有许多强有力的命令，例如，能够用一个单一的命令求解线性系统，能完成大量的高级矩阵处理，能绘制色彩绚丽的二维、三维图形，能与其他程序一起使用。

MATLAB语言的特点包括：

? 具有编程效率高，允许用数学的语言来编写程序

? 用于使用方便，把程序的编辑、编译、连接和执行融为一体 ? 高效方便的矩阵和数组运算 ? 语句简单，内涵丰富 ? 扩充功能强，交互性，开放性 ? 方便的绘图功能

? 该软件由c语言编写，移植性好

MATLAB环境如图1.1所示，其中包括：

图1.1 MATLAB的环境窗口

? 菜单栏 ? 工具栏

? [Command Window]命令窗口 ? [Launch Pad]分类帮助窗口

? [Command History]指令历史记录窗口 ? [Current Directory]当前目录选择窗口 1.2 MATLAB的数组与矩阵

MATLAB中最常见的操作就是对数组与矩阵的操作。MATLAB中的变量一般都默认为向量或矩阵，对于变量的定义不需要任何类型及维度的声明，变量名的第一个字符必须是字母。

常用的常数包括：pi（π），i或j（虚数单位），Inf（无穷），NaN（空值），e（以10为底的幂次，例如1.602e12）等。

常用运算符包括：+ 加，- 减，* 乘，/ 除 , \\ 左除，^ 幂， ’复共轭转置，. 点运算。

常用函数包括sqrt（）开平方；abs（）绝对值；exp（）以e为底的指数；log（x）以e为底的对数；round（）取整；syms定义符号变量等。需要了解更多函数的使用方法请利用帮助函数help。

构造数组的方法有：直接构造（用空格或逗号间隔数组元素），用增量法构

造（数组起值：数组终值或数组起值：步长：数组终值），用linspace函数构造数组（x=linspace(起值，终值，个数）。

构造矩阵的方法有：简单创建方法（直接进行赋值型定义）；构造特殊矩阵（ones创建一个所有元素都为1 的元素；zeros创建一个所有元素都为0的元素；eye创建对角元素为1，其他元素为0的元素；rand创建一个矩阵或数组，其中的元素服从均匀分布；randn创建一个矩阵或数组，其中的元素服从正态分布；diag创建对角矩阵）；聚合矩阵。

获取矩阵元素方法有：取单个元素（矩阵变量名（行，列））；取多个元素（矩阵变量名（：，列），矩阵变量名（行，：））；获取所有元素（矩阵变量名（：））。获取与矩阵有关信息的函数为length（返回最长维长度），ndims（返回维数），numel（返回元素个数），size（返回每一维的长度）。

矩阵的基本运算有*（矩阵相乘），rank（求矩阵的秩），det（矩阵的行列式），poly（求特征多项式），sym（将矩阵转换为符号矩阵）。 1.3 MATLAB绘图

常见绘图函数有plot（x,y）（绘制二维连续图形，其中x为横坐标，y为纵坐标），三维图形绘制plot3（x,y,z）(绘制三维连续图形，其中x为横坐标，y为纵坐标，z为与x，y所在平面垂直的坐标），stem（n，y）（绘制二维离散图形，其中n为离散的横坐标，y为纵坐标），subplot（n，m，i）（在一幅大图中绘制n行m列共计nm幅小图，当前的图是第i幅图）。 1.4 MATLAB中的程序流控 （1）顺序结构

按指令顺序一条条执行下去。 （2）循环结构

for -end ：特点是循环判断条件通常是对循环次数的判断，即循环次数是预先设定的。调用格式为： for 变量=表达式 执行语句 …… 执行语句 end

表达式是一个向量，可表示为m：s：n

while-end：判断控制一般是逻辑判断语句，使用范围更大。调用格式为： while 表达式 执行语句 …… 执行语句 end

表达式一般是关系运算式或逻辑运算式组成的逻辑判断语句，以确定循环是否继续。通常表达式的值非零，即为逻辑真，程序继续循环，否则停止循环。 三. 选择结构

if-else-end 又称为条件语句，根据表达式的情况判断是否满足条件来确定程序下一步的运行。大致可分为三步进行： 首先计算紧跟if后面表达式。

接着判断表达式计算结果，若结果为0，判断值为假；若结果为1，判断值为真。

然后若判断值为真，则执行其后的执行语句组；否则跳过，执行下一个条件表达式或者结束该选择语句。 调用格式为： 1. if 表达式 执行语句组 end 2. if 表达式

执行语句组1 else

执行语句组2 end 3. if 表达式1 执行语句组1 elseif 表达式2 执行语句组2 elseif 表达式3

执行语句组3 …… else

执行语句组n end

关键字if 或elseif后面的条件表达式为条件，通常是由关系运算或与逻辑运算式组成的逻辑判断语句，如果if 或elseif后面的表达式的值为真时，执行紧跟其后的语句内容，否则跳过去，并根据选择语句的表达形式执行后面的elseif表达式语句、跟在else后的执行语句或end语句。 四. 分支语句

switch-case-end 又称为开关语句，使程序在不同的情况下进行相应的操作。调用格式为： switch 表达式

case 常量表达式1 执行语句组1 case 常量表达式2 执行语句组2 ……

case 常量表达式n 执行语句组n otherwise 执行语句组n+1 end

在switch后面的表达式为开关条件，它可以是数字或字符串。当表达式的值与某个case后面的常量表达式的值相等时，就执行相应的语句组，若没有值与所有常量表达式的值相等，执行othwewise后面的语句组。与C语言的switch不同的是，在MATLAB中，当程序执行完某个case语句组后便会立即直接跳出switch语句，执行后续的语句。 五. 其他控制语句

1. 中断语句break：通常用于循环控制中，如for、while等循环，通过if语句判断调用条件，程序在满足条件下调用break语句后，在循环未自然终止之前跳出

当前循环体。在多层循环嵌套中，break只是终止包含break指令的最内层的循环体。

2. 继续语句continue：通常用于循环控制中，可采用if语句判断调用条件，调用continue语句后，程序不再执行循环体内剩余部分的语句而是直接转到循环的终点，继续下一次的循环运行。

3. 返回语句return：使当前正在运行的M文件函数正常结束并返回调用它的函数或程序继续运行，或返回到调用它的环境，如命令窗口。return 语句通常用在函数M文件里面，对某些输入参数或执行结果进行判断，如果出现问题，便调用return语句终止当前程序的运行并返回。 2、信号处理Toolbox介绍

在MATLAB中，针对不同的应用设计了相应的Toolbox，为该领域中的常见应用提供了方便快捷的函数，扩充了MATLAB的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。信号处理Toolbox是其中专用于信号处理学科的领域型工具箱，可方便地应用于信号处理相关问题的仿真与设计。

信号处理Toolbox的内容包括数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT、DCT等常用变换，参数化模型等。具体函数包括如下： 2.1 波形产生

sin，cos，exp 产生正余弦及复指数函数（掌握） chirp 产生扫描频率余弦

diric 产生Dirichlet函数或周期Sinc函数 gauspuls 产生高斯调制正弦脉冲 pulstran 产生脉冲串 rectpuls 产生非周期矩形信号

sawtooth 产生锯齿波或三角波（掌握） sinc 产生sinc函数（掌握） square 产生方波（掌握） strips 产生条图

tripuls 产生非周期三角波 2.2 变换

czt Chirp z-变换 dct 离散余弦变换

dftmtx 离散傅立叶变换矩阵 fft 一维FFT（掌握） fft2 二维FFT

fftshift 函数fft和fft2输出的重新排列 hilbert 希尔伯特（Hilbert）变换 idct 离散余弦逆变换 ifft 一维逆FFT（掌握） ifft2 二维逆FFT 2.3线性系统变换 convmtx 卷积矩阵

latc2tf 格型滤波器转换为传递函数形式 poly2rc 多项式系数转换为反射系数 rc2poly 反射系数转换为多项式系数

residuez z-传递函数的部分分式展开（掌握） sos2ss 二阶级联转换为状态空间

sos2tf 二阶级联转换为传递函数（了解） sos2zp 二阶级联转换为零极点增益形式（了解） ss2sos 状态空间转换为二阶级联形式 ss2tf 状态空间转换为传递函数 ss2zp 状态空间转换为零极点增益 tf2latc 传递函数转换为格型滤波器 tf2ss 传递函数转换为状态空间

tf2zp 传递函数转换为零极点增益（掌握）

zp2sos 零极点增益形式转换为二阶级联形式（掌握） zp2ss 零极点增益形式转换为状态空间 zp2tf 零极点增益转换为传递函数（掌握） 2.4 统计信号处理

cohere 两个信号相干函数估计 corrcoef 相关系数矩阵 cov 协方差矩阵

csd 互功率谱密度估计（CSD） pmem 最大熵功率谱估计

pmtm 多窗口功率谱估计（MTM） pmusic 特征值向量功率谱估计（MUSIC） psd 自功率谱密度估计 tfe 传递函数估计 xcorr 互相关函数估计 xcorr2 二维互相关函数估计 xcov 互协方差函数估计 2.5 滤波器分析与实现 abs 绝对值（幅值）（掌握） angle 相位角（掌握）

conv 卷积和多项式乘法（掌握） conv2 二维卷积

fftfilt 基于FFT重叠加法的数据滤波（掌握） filter IIR或FIR滤波器的数据滤波（掌握） firter2 二维数字滤波

filtfilt 零相位数字滤波（了解） filtic 函数filter初始条件确定 freqs 模拟滤波器频率响应 freqspace 频率响应的频率空间设置 freqz 数字滤波器频率响应（掌握） grpdelay 群延迟（掌握）

impz 数字滤波器的脉冲响应（掌握） latcfilt 格型梯形滤波器实现 unwrap 相位角展开（掌握）

zplane 零极点图（掌握） 2.6 IIR滤波器设计

buttord Butterworth型滤波器阶数的选择（掌握） cheb1ord Chebyshev1型滤波器阶数的选择（掌握） cheb2ord Chebyshev2型滤波器阶数的选择（掌握） ellipord 椭圆滤波器阶次选择（掌握） besself Bessel模拟滤波器设计 butter Butterworth滤波器设计（掌握）

cheby1 Chebyshev I型滤波器设计 （ 通带波纹）（掌握） cheby2 chebyshevII型滤波器设计（阻带波纹）（掌握） ellip 椭圆滤波器设计（掌握）

maxflat 通用数字Butterworth滤波器设计 yulewalk 递归数字滤波器设计 2.7 FIR滤波器设计

cremez 复响应和非线性相位等波纹FIR滤波器设计

fir1 基于窗函数的有限冲激响应滤波器设计----标准响应（掌握） fir2 基于窗函数的有限冲激响应滤波器设计----任意响应（掌握） fircls 多频带滤波的最小方差FIR滤波器设计

fircls1 低通和高通线性相位FIR滤波器的最小方差设计 firs 最小线性相位滤波器设计 firrcos 升余弦FIR滤波器设计 intfilt 插值FIR滤波器设计

kaiserord 用Kaiser窗估计函数fir1参数（掌握） remez Parks-McClellan优化滤波器设计 remezord Parks-McCllan优化滤波器阶估计 bartlett Bartlett窗（掌握） blackman Blackman窗（掌握） boxcar 矩形窗 chebwin Chebyshev窗

hamming Hamming窗（掌握） hanning Hanning窗（掌握） kaiser Kaiser窗（掌握） triang 三角窗（掌握） 2.8 其他 参数建模

invfreqs 由频率响应辨识连续时间（模拟）滤波器 invfreqz 由频率响应辨识离散时间滤波器 levinson Levinson-Durbin递归算法 lpc 线性预测系数

prony Prong法的时域IIR滤波器设计

stmcb 利用Steiglitz-McBride迭代法求线性模型 模拟原型设计

besselap Bessel模拟低通滤波器原型设计 buttap Butterworth模拟低通滤波器原型设计 cheblap Chevbyshev1型模拟低通滤波器原型设计 cheb2ap Chevbyshev2型模拟低通滤波器原型设计 ellipap 椭圆低通滤波器原型设计 频率变换

lp2bp 低通至带通模拟滤波器变换 lp2bs 低通至带阻模拟滤波器变换 lp2hp 低通至高通模拟滤波器变换 lp2lp 低通至低通模拟滤波器变换 滤波器离散变换

bilinear 双线性变换（掌握）

impinvar 冲激不变法的模拟至数字滤波器变换（掌握） 交互式工具

sptool 交互式信号、滤波器和频谱分析工具 特殊运算

cceps 复时谱分析

cplxpair 重新排列组合复数 decimate 降低序列的采样频率 deconv 解卷积和多项式除法 demod 通信仿真中的解调制 detrend 去除线性趋势 dpss Slepain序列

dpssclear 去除数据库Slepain序列 dpssdir 从数据库目录消去Slepain序列 dpssload 从数据库调入Slepain序列 dpsssave Slepain序列存入数据库 icceps 倒复时谱

interp 整数倍提高采样频率 medfilt1 一维中值滤波 modulate 通讯仿真调制 polystap 稳定多项式 rceps 实时谱和最小相位重构 resample 任意倍数改变采样速率 specgram 频谱分析

upfirdn 利用fir滤波器转换采样频率 vco 电压控制振荡器 3、时域信号发生 3.1基本离散时间信号

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

（1）单位采样序列

?1n?0 （1.1） ?[n]??

0n?0?

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

x?zeros(1,N);

x[0]?1;如果?[n]在时间轴上延迟了k个单位，得到?[n?k]即：

?1n?k （1.2） ?[n?k]??

0n?0?

（2）单位阶跃序列

n?0?1 （1.3） u[n]??

n?0?0

在MATLAB中可以利用函数ones实现。

x?ones(1,N);

（3）正弦序列

x[n]?Asin(2?fn??) （1.4）

采用MATLAB的实现方法，如：

n?0:N?1 x?A*sin(2*pi*f*n??)（4）实指数序列

nx[n]?A?a （1.5）

其中，A、a为实数。采用MATLAB的实现方法，如：

n?0:N?1 x?a.^n（5）复指数序列

x[n]?A?e(??j?0)n （1.6）

采用MATLAB的实现方法，如：

n?0:N?1x?A*exp((??j*?0)*n)

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

数real、imag、abs和angle可以逐次计算出一个复数向量的这些函数。

3.2 基本数字调制信号

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

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

SASK(t)?[?ag(t?nT)]cos?tnscn（1.7）

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

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

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

在二进制数字调制中，若正弦载波的频率随二进制基带信号在f1和f2两个频率点间变

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

????S2FSK(t)???ang(t?nTS)?cos?1t???ang(t?nTS)?cos?2t（1.8）

?n??n?

这里，?1?2?f1,?2?2?f2，an是an的反码。典型波形如下：

an

载波信号1

1 0 1 1 0 0 1 t

载波信号2

t

图 1.3 二进制频移键控信号时间波形

（3）二进制相移键控（2PSK或BPSK）

在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二

进制移相键控(2PSK)信号。通常用已调信号载波的0°和 180°分别表示二进制数字基带信号的 1 和 0。二进制移相键控信号的时域表达式为：

??S2PSK(t)???ang(t?nTS)?cos(?ct??i),?i?0或? （1.9）

?n?

典型波形如下：

图 1. 4 二进制相移键控信号时间波形

理解了这些调制信号的产生原理，借助于MATLAB中的循环操作，正弦函数及绘图命令，即可编程产生并显示这类信号。

2PSK信号 an 1 0 0 1 T s t 3.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 双音多频的组合功能

高频群Hz 功能 低频群Hz 1209 1 4 7 * 1336 2 5 8 0 1477 3 6 9 # 1633 A B C D 697 770 852 941

因此，DTMF信号可以看作两个有限长度的正弦序列相加，正弦信号的频率

由按键数字或字母符号对应的频率决定。如，数字“8”由行频852Hz和列频1336Hz决定。理解了DTMF信号的原理，利用MATLAB中的正弦函数，即可产生该信号。

四、实验目的：

1、 掌握MATLAB的基本编程方法，包括MATLAB中输入输出，赋值，简单向量或矩阵运算，绘图等基本功能。

2、 掌握信号处理Toolbox中基本函数使用方法（参见前面信号处理Toolbox函数列表中标注了“掌握”字样的函数），为后续实验打好基础。 3、 掌握几种基本的离散时间信号（包括单位采样序列，单位阶跃序列，单频正弦序列，单频复指数序列，实指数序列等），并能够熟练利用MATLAB产生这些信号。

4、 理解ASK、FSK、BPSK信号，及DTMF信号的产生原理，并能够利用MATLAB产生它们。

五、实验内容：

1、设计并编写MATLAB程序，要求包含输入输出操作，赋值，简单程序流控（循环，逻辑判断，选择结构，分支语句等），绘图操作。

2、查阅MATLAB中help，学习信号处理Toolbox中要求掌握及了解的函数使用方法。

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

4、（拓展要求）利用MATLAB编程产生并绘图显示2ASK，2FSK，2PSK

等数字调制信号。

5、（拓展要求）利用MATLAB编程产生并绘图显示DTFM信号。 6、（拓展要求）利用MATLAB编程产生并绘图显示高斯白噪声序列。 7、（拓展要求）利用MATLAB中的谱分析函数对正弦信号的频谱进行分析

并绘图显示。

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

安装了MATLAB软件的PC机一台。

七、实验步骤：

1、设计一段MATLAB程序，要求包括从键盘输入，赋值，基本矩阵和向量计算，程序流控（循环，逻辑判断，选择结构，分支语句等）和绘图操作。

2、查阅MATLAB中的help，将信号处理Toolbox中要求掌握及了解的函数做成列表，分别列出它们的输入输出。（参照表1-2例子） 表1-2 信号处理Toolbox中要求掌握及了解的函数的列表

函数名 sin,cos,exp 输入 2pi*f0*n/fT（数字角频率*n） 输出 每个n对应的函数值 功能备注 计算正余弦及复指数函数的值

3、在?20?n?20内，画出单位下列信号：

(a)单位采样序列x1[n]??[n]和单位阶跃序列x2[n]?u[n]的时域波形图。 (b)y1[n]?x1[n?5]、y2[n]?x2[n?8]的波形。说明x1[n]与y1[n]、x2[n]与

y2[n]之间的关系。

4、画出下列信号在0?n?100内的波形。

??n?x3[n]?sin???16??n?x4[n]?sin??

2????n??3?n?x5[n]?cos???cos??128????

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

5、在0?n?30内，画出下列信号：

x6[n]?0.2(0.8)n

x7[n]?e??1/12?j?/6?n

对于复数序列，要求分别画出实部和虚部；幅值和相角。若把x6[n]中的底数0.8分别改为1.2、-0.8，讨论产生的时域波形有何变化。总结指数序列的底数对序列变化的影响。

6、（拓展要求）设计产生数字二进制序列：1 0 1 0 1 0 的2ASK、2FSK、2PSK调制信号。已知符号速率Fd＝10Hz（即时间间隔Ts为0.1），输出信号的采样频率为200Hz。

(a)2ASK信号的载波频率Fc＝40Hz，

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

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

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

8、（拓展要求）MATLAB函数randn(1,N)可以产生均值为0，方差为1的高斯随机序列，也就是白噪声序列。试利用randn函数产生均值为0.15，方差为0.1的高斯白噪声序列x8[n]，要求序列时域范围为0?n?100。画出时域波形图。同时将实验步骤2中产生的信号x2[n]与x8[n]相加，将得到的波形与x2[n]的波形做比较。

9、（拓展要求）利用MATLAB中的谱分析函数画出x3[n]、x4[n]、x5[n]的频谱。与理论上根据傅立叶变换的定义计算出的x3[n]、x4[n]、x5[n]的频谱进行比较。

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

需要的程序代码：

（1）一段MATLAB包括从键盘输入，赋值，基本矩阵和向量计算，程序流控（循环，逻辑判断，选择结构，分支语句等）和绘图操作的MATLAB简单程序。

实验1代码如下：

disp('请选择本函数所要执行的功能');

disp('k=1：绘出零极图；k=2：绘出相位延迟；k=3：绘出群延迟'); k=input('k='); num1=[-0.2 0.8 0.4 1]; num2=[2.4 -0.9 1.2 2]; num=num1+num2; den1=[1 0.4 0.18 -0.2]; den2=[ 2 3.5 -2.1 1.05]; den=den1.*den2; figure(1); switch k

case 1;

zplane(num,den); case 2

phasedelay(num,den); case 3

grpdelay(num,den); otherwise

disp('错误！'); end disp(num); disp(den);

（2）产生x1[n]、x2[n]、y1[n]、y2[n]、x3[n]、x4[n]、x5[n]、x6[n]、x7[n]序列的程序

实验2代码如下： N=41; n=-20:20; figure(1); %x1[n] x1=zeros(1,N); x1(21)=1; subplot(2,2,1); stem(n,x1,'fill');

xlabel('n'),ylabel('x1[n]'),title('x1[n]的时域波形图'); %x2[n]

x2=[zeros(1,20) ones(1,21)]; subplot(2,2,2);

stem(n,x2,'fill')

xlabel('n'),ylabel('x2[n]'),title('x2[n]的时域波形图'); %y1[n]

y1=[x1(1+5:length(x1)),zeros(1,5)]; subplot(2,2,3); stem(n,y1,'fill');

xlabel('n'),ylabel('y1[n]'),title('y1[n]的时域波形图'); %y2[n]

y2=[zeros(1,8),x2(1:length(x2)-8)]; subplot(2,2,4); stem(n,y2,'fill')

xlabel('n'),ylabel('y2[n]'),title('y2[n]的时域波形图'); n=0:100; figure(1); %x3[n] x3=sin(pi*n/16); subplot(3,1,1) stem(n,x3,'fill')

xlabel('n'),ylabel('x3'),title('x3[n]的时域波形图'); %x4[n] x4=sin(n/2); subplot(3,1,2) stem(n,x4,'fill')

xlabel('n'),ylabel('x4'),title('x4[n]的时域波形图'); %x5[n]

（3）信号的时移：x1[n]与y1[n]、x2[n]与y2[n]之间的关系。

通过比较x1[n]与y1[n]、x2[n]与y2[n]的波形可知，y1[n]相当于将x1[n]向左搬移5个单位，y2[n]相当于x2[n]向右搬移8个单位。

（4）正弦序列x3[n]、x4[n]、x5[n]周期的判断

通过观察x3[n]、x4[n]、x5[n]的波形可知，x3[n]是周期信号，x4[n]不是周期信号，x5[n]是周期信号。

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

通过比较他们的图形，当指数序列的底数大于0小于1时，该序列能在正数范围内单调递减；当指数序列的底数大于1时，该序列在正数范围内能单调递增；当指数序列的底数小于0且绝对值小于1时，该序列在正数范围内出现正值与负值的交替且绝对值单调递减。

（6）2ASK、2FSK、2PSK调制信号时域波形（拓展要求） （7）数字“0” DTMF时域波形（拓展要求）

（8）高斯白噪声序列的时域波形；正弦序列加上高斯白噪声后的时域波形

（扩展要求）

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

九、实验结论：

通过本次实验，我掌握了MATLAB的基本编程方法，包括MATLAB中输入输出，赋值，简单向量或矩阵运算，绘图等基本功能。掌握了信号处理Toolbox中基本函数使用方法，为后续实验打好基础。掌握了几种基本的离散时间信号（包括单位采样序列，单位阶跃序列，单频正弦序列，单频复指数序列，实指数序列等），并能够熟练利用MATLAB产生这些信号。理解了ASK、FSK、BPSK信号，及DTMF信号的产生原理，并能够利用MATLAB产生它们。

十、总结及心得体会：

我了解了多种离散时间信号的产生方法，加深了matlab对该课程的的重要作用。

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

本次实验可以使用更加多种多样的信号处理函数。

报告评分：

指导教师签字：

百度搜索“yundocx”或“云文档网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，云文档网，提供经典综合文库电子科技大学信号处理实验1_信号的基本表示及时域分析在线全文阅读。