实验一典型离散信号的生成及运算.docVIP

实验一 典型离散信号的生成及运算 一、实验目的： 1）熟悉MATLAB环境。 2）理解离散时间信号。 3）掌握几个典型离散信号类型及其运算。 4）学会利用MATLAB表示和实现典型离散信号类型及其运算。 二、实验原理： 广义讲，信号可以分为模拟和离散信号【1】： 一个模拟信号：，t可以代表任何物理量，现假定t代表时间，连续的。 一个离散信号：，n是整数值，现假定n代表时间上离散的时刻。 因此称为离散时间信号，是一个数值的序列。 在MATLAB中，只需要两一个行向量来表示这个有限长序列： 例如： 图1 命令窗口(Command Window) 图2 工作空间（WorkSpace） 注意，要准确的表示一个离散信号，还需要另外一个向量表示n（样本位置信息），但是一般来说，当序列从位置n = 1 开始时（MATLAB中数组下标从1开始而不是从0开始），只用x向量表示该信号。 下面说明几种典型离散信号类型和运算，及其MATLAB的表示： 1、单位抽样序列 又称为单位脉冲序列/单位样本序列，其特点是在n=0时取值为1，其它取值为0。 MATLAB中，使用函数zeros（1，N）可以产生一个N个零的行向量。然后对其中 n = 0 这个位置进行单独赋值为1即可。如： 图3 没有明确样本的位置信息及其stem图 图4 明确样本的位置信息及其stem图 这里还有另外一种便利的表示方法——利用逻辑关系式表示方法： 图5 使用逻辑关系式表示及其stem图 2、单位阶跃序列 MATLAB中，使用函数ones（1，N）可以产生一个N个1的行向量。下面通过N个0序列合并上N+1个1序列产生单位阶跃序列。 注意1和0的样本位置。 图6 单位阶跃序列表示及其stem图 请同学们思考，如果利用逻辑关系式，程序又该怎么写？ 3、实值指数序列 MATLAB中，使用运算符“.^”实现实指数序列。 图7 单位阶跃序列表示及其stem图 请同学们思考： 1）当a = 1 时是什么序列？ 2）图7横坐标为0时，纵坐标是多少？ 3）试画出a = -0.5 、a = 3 时的指数序列图，判断收敛和发散。 4、正弦序列 注意：是初始相位，正弦函数与余弦函数仅初相位不同。 1）模拟正弦信号 t是时间上的自变量，连续变化 从单位圆来看，f=1时，旋转一圈表示信号变化一个周期，那么f可以用来描述变化快慢即模拟频率，表示每秒信号旋转的圈数（即每秒可以变化f次/（圈））。单位Hz，1/s。 2）从单位圆的角度变化来看，旋转一圈对应的角度变化范围是，旋转f圈则角度变化范围是，所以用角频率描述变化快慢。单位rad/s，每秒经过多少弧度。 所以，模拟正弦信号 t还是时间上的自变量，连续变化 3）数字正弦信号是从模拟正弦信号采样而来。所以讲数字频率就离不开采样周期Ts（采样频率Fs = 1/Ts）。 在单位圆上以N点等间隔采样(点与点之间经过周期Ts这么长时间），那么相邻两个采样点之间的弧度值，表示数字角频率也称数字频率，它是一个相对频率，仅一个称呼，又称归一化数字角频率。此时自变量t=n*Ts=n/Fs了，离散值。 所以，数字正弦信号 N是数字周期，经过需多少个点，即一个周期内有N个点，单位1，。 一般地，表示正弦信号， 需要确定幅值A，初始相位，频率f或， 确定采样点数N及采样频率Fs，则样本位置序列t=n*Ts=[0:N-1]*Ts, 最后利用公式得到正弦序列。 如产生一个频率为50Hz，幅值为2，初始相位为的正弦信号： 图8 正弦序列表示及其stem图 注意: 1）采样频率Fs取值一定要满足奈奎斯特定理。 2）采样频率和信号频率的比值：Fs/f = T/Ts，也就是一个信号周期内有多少个采样周期，即数字周期，也就确定了一个信号周期的采样点数。（上例是20）。 3）总采样点数N和数字周期的区别。（上例是100） 4）绘图时，注意样本位置序列的选取。 5、复正弦序列 设，那么复指数函数定义为： 复指数序列可以表示任意序列 当时，复正弦序列： 如产生一个频率为50Hz的复正弦信号，绘出它的实部，虚部，模，相角： clear; clc; Fs = 1000; N = 100; f = 50; Ts=1/Fs; t=[0:N-1]*Ts; x = cos(2*pi*f*t)+j*sin(2*pi*f*t); x_real = real(x); x_imag = imag(x); x_abs = abs(x); x_angle = angle(x); subplot(2,2,1),stem(t,x_real,filled),title(实部),xlabel(t); subplot(2,2,3),