实验1 信号的ATLAB 表示及系统的时域分析.doc

广东技术师范学院实验报告 学院： 自动化学院 专业： 班级： 成绩： 姓名： 学号： 计算机编号： 实验地点： 实验日期： 指导教师签名： 预习情况 良好 操作情况 良好 考勤情况 全勤 数据处理情况 良好 实验 （一） 项目名称： 信号S=symadd(f1,f2) Ezplot(s) 以上是用Matlab的符号运算命令来表示两连续信号的相加，然后用ezplot命令绘制出其结果波形图。 相乘 w=symmul(f1,f2) ezplot(w) 以上是用Matlab的符号运算命令来表示两连续信号的相乘，然后用ezplot命令绘制出其结果波形图。 时移 y=subs(f,t,t-t0) ezplot(y) 以上是用Matlab的符号运算命令来实现连续信号的平移及其结果的可视化，subs命令将连续时间信号中的时间变量t用t-t0替换。 反褶 y=subs(f,t,-t) ezplot(y) 以上是用Matlab的符号运算命令来实现连续信号的反褶及其结果的可视化。 尺度变换 y=subs(f,t,a*t) ezplot(y) 以上是用Matlab的符号运算命令来实现连续信号的尺度变换及其结果的可视化。 连续时间信号的卷积 Matlab的库函数conv()只能计算离散序列f1与f2的卷积和，构造函数sconv()实现两连续函数的卷积积分。 conv()的调用形式：conv(f1,f2)。 f1与f2的卷积过程如下： 将f1(t)与f2(t)以时间间隔λ进行抽样，得到离散序列f1(kλ)和f2(kλ)； 构造与f1(kλ) 和f2(kλ)相对应的时间向量k1和k2； 调用conv()函数计算卷积积分f(t)的近似向量f(kλ)； 构造f(kλ)对应的时间向量K。 连续时间系统的响应 连续时间系统零状态响应的MATLAB实现 Matlab库函数中的lsim()能对微分方程描述的LIT连续时间系统响应进行仿真。该函数能绘制连续时间系统在指定的任意时间范围内系统响应的时域波形图，还能求出连续时间系统在指定的任意时间范围内系统响应的数值解。 函数lsim()的调用格式： lsim(b,a,x,t) 该调用格式中，a和b是由描述系统的微分方程左边和右边系数构成的两个行向量；x和t是表示输入信号的行向量，其中t为输入信号时间范围的向量，x为输入信号在向量t定义的时间点上的抽样值。 图描述线性移不变离散时间系统的数学模型是常系统差分方程，它与系统的结构流图之间可以互相推导。求解差分方程最常用的方法是迭代解法，也是实现数字滤波器的一种基本方法。 用与分别表示系统的激励和响应，差分方程通式为： ( – 1 ) 离散系统的响应可分为零输入响应分量和零状态响应分量，零输入分量仅由系统的初始状态引起，外激励；而零状态分量仅由外激励引起，初始状态均为零。 对于线性移不变离散时间系统，如果已经知道了系统的单位冲激响应 ，就可以利用离散卷积求出系统在任何激励作用下的零状态响应 。在本实验中通过输入序列、，得到的直观图形表示。 离散卷积的计算公式如下： ( – 2 ) 在离散卷积中，多讨论有限长序列。若和长度分别为 M 和 N，则卷积结果即响应序列也是有限长序列，长度为 L=M+N-1。上式形象描述了离散卷积中两个有限长序列反褶、移位、相乘、累加的过程。 本实验差分方程解法中只限于激励是单位阶跃信号，即的情况，通过给定系统阶数 N 和系数向量和以及初始状态的值可以求出系统在单位阶跃信号激励下的响应，包括单位冲激响应以及激励下的全响应和零输入响应、零状态响应。至于其它激励下的零状态响应，可以用它的单位冲激响应与输入信号的离散卷积求出。 x=conv(x1,x2) 上式中，输入参量x1、 x2分别为包含序列x1(n)和x2(n)的所有非零样值点的行向量；输出参量x为返回序列x(n)=x1(n) ＊ x2(n) 的所有非零样值点的行向量。 三．实验内容 1.设f(t)=(t+2)(u(t+2)-u(t))+2*(u(t)-u(t-2))。请用Matlab绘出下列信号的时域波形。 （1）f(t) (2) f(3t)（3）f(t-3)u(t-3) （4）f(3-t) 头文件：function f=heaviside(t) f=(t0) syms t; f=sym((t+2)*(heaviside(t+2)+2*(heaviside(t)-heaviside(t-2)))); subplot(2,2,1),ezplot(f,[-3,3]); title(f(t)); grid; y1=subs(f,t,3^