《实验 信号的matlab 表示及系统的时域分析》由会员分享，可在线阅读，更多相关《实验 信号的matlab 表示及系统的时域分析（21页珍藏版）》请在金锄头文库上搜索。
1、广东技术师范学院实验报告学院：自动化学院专业：班级：成绩：姓名：学号：计算机编号：实验地点：实验日期：指导教师签名：预习情况良好操作情况良好考勤情况全勤数据处理情况良好实验 （一） 项目名称： 信号的MATLAB 表示及系统的时域分析 一实验目的和要求1掌握函数的向量表示法2掌握函数的符号运算表示法，实现连续信号的时域运算和变换3.掌握连续时间信号卷积的原理及其matlab实现4.掌握连续时间系统响应的原理及其matlab实现5.掌握离散系统单位序列响应； 6.掌握离散序列卷积和的计算方法 。二实验原理1向量表示法Matlab的信号处理工具箱有大量的函数可用于产生信号，函数中大部分要求用矢量来
2、表示时间。如：他0：0.005：1；表明该矢量抽样频率为200Hz。Matlab可以精确表示离散时间信号，对连续时间信号只能近似表示。因此应选择足够小的抽样间隔Ts，以保证反映信号的全部细节。2符号运算表示法用符号运算可以较快地实现连续信号地时域运算和变换。例：（1） 相加S=symadd(f1,f2)Ezplot(s)以上是用Matlab的符号运算命令来表示两连续信号的相加，然后用ezplot命令绘制出其结果波形图。（2） 相乘w=symmul(f1,f2)ezplot(w)以上是用Matlab的符号运算命令来表示两连续信号的相乘，然后用ezplot命令绘制出其结果波形图。（3） 时移y=s
3、ubs(f,t,t-t0)ezplot(y)以上是用Matlab的符号运算命令来实现连续信号的平移及其结果的可视化，subs命令将连续时间信号中的时间变量t用t-t0替换。（4） 反褶y=subs(f,t,-t)ezplot(y) 以上是用Matlab的符号运算命令来实现连续信号的反褶及其结果的可视化。（5） 尺度变换y=subs(f,t,a*t)ezplot(y)以上是用Matlab的符号运算命令来实现连续信号的尺度变换及其结果的可视化。3. 连续时间信号的卷积Matlab的库函数conv()只能计算离散序列f1与f2的卷积和，构造函数sconv()实现两连续函数的卷积积分。conv()的调
4、用形式：conv(f1,f2)。f1与f2的卷积过程如下：（1） 将f1(t)与f2(t)以时间间隔进行抽样，得到离散序列f1(k)和f2(k)；（2） 构造与f1(k) 和f2(k)相对应的时间向量k1和k2；（3） 调用conv()函数计算卷积积分f(t)的近似向量f(k)；（4） 构造f(k)对应的时间向量K。4. 连续时间系统的响应A 连续时间系统零状态响应的MATLAB实现Matlab库函数中的lsim()能对微分方程描述的LIT连续时间系统响应进行仿真。该函数能绘制连续时间系统在指定的任意时间范围内系统响应的时域波形图，还能求出连续时间系统在指定的任意时间范围内系统响应的数值解。函
5、数lsim()的调用格式：lsim(b,a,x,t)该调用格式中，a和b是由描述系统的微分方程左边和右边系数构成的两个行向量；x和t是表示输入信号的行向量，其中t为输入信号时间范围的向量，x为输入信号在向量t定义的时间点上的抽样值。 图2.3B连续时间系统冲激响应的MATLAB实现MATLAB为用户提供了专门用于连续系统冲激响应并绘制其时域波形的impulse函数。函数impulse()的调用格式：impulse(b,a)C连续时间系统冲激响应的MATLAB实现MATLAB为用户提供了专门用于连续系统冲激响应并绘制其时域波形的step函数。函数step()的调用格式：step(b,a)5. 描
6、述线性移不变离散时间系统的数学模型是常系统差分方程，它与系统的结构流图之间可以互相推导。求解差分方程最常用的方法是迭代解法，也是实现数字滤波器的一种基本方法。用与分别表示系统的激励和响应，差分方程通式为： ( 3 1 )离散系统的响应可分为零输入响应分量和零状态响应分量，零输入分量仅由系统的初始状态引起，外激励；而零状态分量仅由外激励引起，初始状态均为零。对于线性移不变离散时间系统，如果已经知道了系统的单位冲激响应 ，就可以利用离散卷积求出系统在任何激励作用下的零状态响应 。在本实验中通过输入序列、，得到的直观图形表示。离散卷积的计算公式如下： ( 3 2 )在离散卷积中，多讨论有限长序列。若
7、和长度分别为 M 和 N，则卷积结果即响应序列也是有限长序列，长度为 L=M+N-1。上式形象描述了离散卷积中两个有限长序列反褶、移位、相乘、累加的过程。6. 本实验差分方程解法中只限于激励是单位阶跃信号，即的情况，通过给定系统阶数 N 和系数向量和以及初始状态的值可以求出系统在单位阶跃信号激励下的响应，包括单位冲激响应以及激励下的全响应和零输入响应、零状态响应。至于其它激励下的零状态响应，可以用它的单位冲激响应与输入信号的离散卷积求出。Matlab为用户提供了专门用于求离散系统单位序列响应数值解并绘制其时域波形的impz函数。其调用形式之一为：Impz（b,a,N）其中，a和b分别为描述离散
8、系统的差分方程系数所决定的行向量，N为正整数，用于指定绘制系统单位序列响应的时间范围，即该调用格式将绘出由向量a和b定义的离散系统在0N个时间样点区间单位序列响应的时域波形。7. Matlab提供了用于求两个有限时间区间非零的离散时间序列卷积和的专用函数conv。其调用格式为：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绘出下列信号的时域
9、波形。（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,3t);subplot(2,2,2),ezplot(y1,-5,1);title(f(3t);grid;b=subs(f,t,t-3);y2=b* heaviside(t-3);subplot(2,
10、2,3),ezplot(y2,1,7);title(f(t-3)u(t-3);grid;a=subs(f,t,-t);y3=subs(a,t,t+3);subplot(2,2,4),ezplot(y3,-5,1);title(f(3-t);grid;2. 已知，试用Matlab计算，并绘出的函数波形。syms t;f=sym( (heaviside(t)-heaviside(t-2);subplot(2,2,1),ezplot(f,-3,3);title(f1(t);y1=(0.5*t)*subs(f,t,t+2);subplot(2,2,2),ezplot(y1,-3,3);title(f2
11、(t);grid;3. 已知描述连续时间系统的微分方程和激励信号分别为试用Matlab求该系统的零状态响应，并绘制其时域仿真波形。4.已知描述连续时间系统的微分方程为试求系统的冲激响应和阶跃响应，并绘出时域仿真波形。a=1 6 8;b=1 0subplot(2,2,1)y=impulse(b,a)plot(t,y)subplot(2,2,2)y=step(b,a)plot(t,y)%program3ts=0;te=5;dt=0.01sys=tf(1,0,1,6,8)T=ts:dt:tesubplot(2,2,1)y=impulse(sys,t)plot(t,y)subplot(2,2,2)y=
12、step(sys,t)plot(t,y)xlabel(Time(sec)ylabel(h(t)xlabel(Time(sec)ylabel(h(t)5. 已知x1k=1,1,0,1,0,1，x2k=1,1,0,1,0,1,1,0,1,1用Matlab计算离散卷积和y(k)=x1(k)*x2(k)，并绘出x(k)和y(k)的时域波形。6. 已知描述离散系统的差分方程为y(n)-0.6y(n-1)+0.8y(n-2)=x(n)-3x(n-1)试利用Matlab绘出该系统单位序列响应的时域波形4 实验总结functionf,k=sconv(f1,f2,k1,k2,p) 计算连续信号卷积积分f(t)=
13、f1(t)*f2(t) f: 卷积积分f(t)对应的非零样值向量 k: f(t)的对应时间向量 f1: f1(t)非零样值向量 f2: f2(t)非零样值向量 k1: f1(t)的对应时间向量 k2: f2(t)的对应时间向量 p: 取样时间间隔f=conv(f1,f2); %计算序列f1和f2的卷积和f f=f*p; k0=k1(1)+k2(1); %计算序列f非零样值的起点位置 k3=length(f1)+length(f2)-2 %计算卷积和f非零样值的宽度 k=k0:p:k3*p %确定卷积和f非零样值的时间量 subplot(2,2,1) plot(k1,f1) %在子图1绘出f1(t)的时域波形图 title(f1(t) xlabel(t)ylabel(f1(t)subplot(2,2,2) plot(k2,f2) %在子图2绘出f2(t)的时域波形图 title(f2(t) xlabel(t) ylabel(f2(t) subplot(2,2,3) plot(k,f) %画卷积f(t)的时域波形图 h=get(gca, position) h(3)=2.5*h(3) set(gca, position，h) %将第三个子图的横坐标范围扩为原来的2.5倍 title(f(t)=f1(t)*f2(t) xlabel(t)ylabel(f(t)
