用MATLAB实现常用地离散时间信号及其时域运算.docx

单 单 分 word 实验四 用 MATLAB实现常用的离散时间信号与其时域运算 —— 摘要： 在 MATLAB中，只能用向量来表示离散时间信号。与连续信号不同，离 散时间信号无法用符号运算来表示。 用适当的 MATLAB语句表示出信号后， 就可 以利用 MATLAB的绘图命令 stem来绘 出直观的信号波形图， stem 是专门用于绘制离散 时间信号的。 在 MATLAB 中离散序列的时域运算和变换不能用符号运算来实现， 而必须用向 量表示的方法，即在 MATLAB 中离散序列的相加、相乘需表示成两个向量的相加、相乘，因 而参加运算的两序列向量必须有一样的维数。 一、实验目的： 〔 1〕学习 MATLAB 语言与其常用指令； 〔 2〕学习和掌握用 MATLAB语言产生离散时间信号的编程方法； 〔 3〕通过编程绘制出离散时间信号的波形，加深理解信号的时域运算。 二、实验内容： 〔 1〕运用 MATLAB 的绘图指令绘制离散时间信号； 〔 2〕用 MATLAB语言实现离散时间信号的时域运算。 三、实验原理： 〔 1〕单位阶跃序列和单位样值序列。 离散时间信号只在某些离散的瞬时给出信号的值， 因此， 它是时间上不连续的序列。 位阶跃序列和单位样值序列在离散时间信号与系统的分析中是两个非常典型的序列， 别记为 u(n)和δ (n)。它们的定义分别如下： 1 n≥0 u(n)= 0 n0 1 n≥0 δ(n)= 0 n≠0 1 / 10 word 假如单位阶跃序列的起始点为 n0，单位样值序列出现在 n0 时刻，如此表达式分别为： 1 u(n-n0)= 0 n≥n0 1 n=n0 δ (n-n0)= nn0 0 n≠n0 应注意， 离散时间的单位阶跃序列与连续时间的单位阶跃信号的异同， 以与离散时间的 单位样值序列与连续时间的单位冲激信号的异同。 〔2〕离散时间信号的时域运算。 与连续时间系统的研究类似， 在离散系统分析中， 经常遇到离散时间信号的运算， 包括 两信号的相加、 相乘以与序列自身的移位、反褶、 尺度等等， 也需要了解在运算过程中序列 的表达式以与对应的波形的变化。 序列 x(n)的移位： x(n-n0) 序列 x(n)的反褶： x(-n) 序列 x(n)的尺度变换： x(an) 两序列 x1(n)与 x2(n)的相加减： x1(n) ±x2(n) 两序列与的相乘： x1(n) · x2(n) 〔3〕学习如何使用 MATLAB 语言产生离散时间信号并对离散时间信号进展时域运算。 四、实验任务： 〔1〕 编制用于产生如下信号的通用程序， 调试并运行这些通用的程序。 ①x(n)=A δ (n-n0) 要求对于任意给定的参数都能实现所要求的信号。 2 / 10 word 程序： function un(t1,t2,t0) t=t1:t2; n=length(t); tt=t1:t2; n1=length(tt); f=zeros(1,n); f(1,t0-t1+1)=3; stem(t,f),grid on title( μ￥ ??3?? ÷ D?o? ) axis([t1,t2 -0.2 4]) ②x(n)=A[u(n-m1)-u(n-m2)] 程序： function unn(t1,t2,A) t=t1:t2; n=length(t); 3 / 10 word f=[zeros(1,t1-1),A*ones(1,t2-t1),zeros(1,1)]; stem(t,f);grid on axis([t1 t2 -0.2 4]) 〔2〕离散序列波形。 ①f(k-2)u(k); u(k)可以直接调用 stepseq 函数，离散序列的平移编程如下： function[f,k]=lspy(ff,kk,k0) k=kk+k0; f=ff; 离散序列的相乘程序为：