信号与系统仿真作业

1、ngdou-b 11112广 东 海 洋 大 学 学 生 实 验 报 告 书 （ 学 生 用 表 ）课 程 名 称课 程 号学 院 （ 系 )信 息 学 院专业班级学生姓名学号实验地点04002 实验日期实验一连时间信号的 matlab 表示和连续时间 lti 系统的时域分析一、实验目的1掌握matlab产生常用连续时间信号的编程方法，并熟悉常用连续时间信号的波形和特性；2运用 matlab符号求解连续系统的零输入响应和零状态响应；3运用 matlab数值求解连续系统的零状态响应; 4运用 matlab求解连续系统的冲激响应和阶跃响应；5运用 matlab卷积积分法求解系统的零状态响应。二、实

2、验原理1. 连续信号 matlab实现原理从严格意义上讲,ma tlab 数值计算的方法并不能处理连续时间信号.然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号,即当取样时间间隔足够小时，这些离散样值能够被matlab处理，并且能较好地近似表示连续信号. matlab提供了大量生成基本信号的函数.比如常用的指数信号、正余弦信号等都是matlab的内部函数。为了表示连续时间信号，需定义某一时间或自变量的范围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图. 三、实验内容1实例分析与验证根据以上典型信号的ma tlab 函数，分析与验证下列典型信号matlab程序，并实

3、现各信号波形图的显示，连续信号的图形显示使用连续二维图函数plot （). （1） 正弦信号：用matlab 命令产生正弦信号2sin(2/ 4) ，并会出时间0t 3的波形图。程序如下：k=2;w=2*pi ;phi=pi/4; t=0： 0.01：3；ft=k*sin （w*t+phi ） ；plot(t,ft ） ,grid on；axis（0,3,-2。2,2.2 ) title（正弦信号) （2） 抽样信号：用matlab中的 sinc(t)函数命令产生抽样信号sa(t),并会出时间为66t的波形图。程序如下：t=-6*pi ：pi/100:6*pi ；ft=sinc(t/pi);

4、plot(t，ft),grid on ；axis（-20,20，-0.5，1.2）title（抽样信号 ) (3) 矩形脉冲信号:用 ma tlab 中 y=rectpuls（t，width) 函数命令画出下列矩形脉冲信号的波形图2(01)( )0(0,1)tf ttt程序如下：t=-0.5：0.01:3；t0=0。5；width=1 ；ft=2*rectpuls （tt0，width) ；plot（t，ft）,grid on; axis（ -0.5，3,0.2，2.2 ）title(矩形脉冲信号) (4) 单位阶跃信号：用ma tlab 命令“ y=（ t=0)”绘出 1t 5 单位阶跃信号

5、 （t） 。程序如下：t=1：0.01:5; ft=（t=0）; plot（t，ft）,grid on ;网格axis( 1，5,0.5,1.5） ；title(单位阶跃信号) 也可通过定义阶跃函数来实现，function f=uct(t) ；f=（t=0); 保存为 uct.m 文件，上面实现阶跃信号的程序变为t=1:0。 01:5; ft=uct （t）; plot（t，ft),grid on; axis（ 1,5，0.5,1.5)；title（单位阶跃信号）(5) 连续时间系统零输入响应和零状态响应的符号求解试 用matlab命 令 求 解 微 分 方 程( )3( )2 ( )( )3

6、 ( )yty ty tx tx t， 当 输 入3( )( )tx tet，起始条件为(0 )1y、(0 )2y时系统的零输入响应、零状态响应及完全响应。源程序为：eq= d2y+3*dy+2 y=0; cond=y（0）=1，dy（0）=2 ；yzi=dsolve(eq，cond)；yzi=simplify （yzi）运行结果 : yzi = 3*exp( 2t）+4exp（t）eq1=d2y+3*dy+2*y=dx+3*x；eq2=x=exp(-3*t)*heaviside （t) ；cond=y(0.001)=0,dy （0。 001)=0；yzs=dsolve（eq1，eq2， co

7、nd） ；yzs=simplify （ yzs.y）运行结果：yzs = heaviside（ t) (exp（2*t） +exp(t） ）yt=simplify(yzi+yzs) 运行结果 : yt = -3*exp（2t)+4exp(t)exp(2t）*heaviside（t）+exp(t）*heaviside（t) subplot（311）ezplot(yzi ， 0，8)；grid on；xlabel(t）title(零输入响应 ) subplot(312）ezplot（yzs， 0， 8 ）;grid on xlabel（ t）title(零状态响应 ) subplot(313）ez

8、plot(yt ，0，8）;grid on xlabel(t）title（完全响应) (6） 连续时间系统零状态响应的数值求解试 用matlab数 值 求 解 微 分 方 程( )3( )2 ( )( )3 ( )yty ty tx tx t， 当 输 入3( )( )tx tet时系统的零状态响应。源程序为：ts=0；te=8;dt=0。01；sys=tf(1 ，3,1,3,2）; t=ts:dt:te; f=exp（-3*t) 。*uct(t ）; y=lsim （sys，f，t); plot（t，y）,grid on；axis(0 8 0。02 0。27）xlabel（t（s) ),yl

9、abel(y(t)title(零状态响应) (7) 连续时间系统冲激响应和阶跃响应的求解已知某 lti 系统的微分方程( )2 ( )32 ( )( )16 ( )y ty ty tftf t，试用 matlab 命令绘出04t范围内系统的冲激响应h(t)和阶跃响应s(t） 。源程序为 : t=0:0.001:4; sys=tf(1 ，16 ，1,2,32） ；h=impulse(sys,t）; s=step(sys，t); subplot(211） ；plot(t,h ） ，grid on xlabel（t（s）)，ylabel(h（t）title(冲激响应 ) subplot（212）;p

10、lot（t,s)，grid on xlabel(t（s） ） ,ylabel(s（t）) title（阶跃响应 ）(8） 利用卷积积分法求系统的零状态响应已知某 lti 系统的微分方程( )2 ( )32 ( )( ) 16 ( )y ty ty tftf t，其中，2( )tf te。试用 matlab 卷积积分方法绘出系统零状态响应y（t)的波形图。程序如下：dt=0.01;t1=0:dt:4 ；f1=exp（2t1)；t2=t1; sys=tf（ 1，16 ，1，2,32） ；f2=impulse(sys,t2) ；t,f=ctsconv （f1,f2 ，t1，t2，dt); 2。 实践

11、编程（1) 复指数信号：用matlab 命令画出 0t 3 复指数信号( 1.510)( )2jtf te的实部、虚部、模及相角随时间变化的曲线,并观察其时域特性. 程序： t=0：0。01：3；k=2;a=-1。5；b=10；ft=k*exp(a+i b）*t) ；subplot(2,2,1） ；plot(t,real(ft ）);title （ 实部） ；axis(0，3,-2,2 ）;grid on；subplot（2，2,2);plot(t,imag （ft） ）;title( 虚部 )； axis（ 0,3， 2，2)；grid on; subplot(2，2,3）;plot（t，

12、abs（ft） ）;title( 模 );axis（ 0，3， 0,2);grid on ；subplot(2，2,4） ；plot（t,angle(ft） );title（相角 );axis（0， 3， 4,4 ） ;grid on；（2） 用 matlab中 y=square（t，duty) 产生频率为10hz、占空比为30%的周期方波信号。程序：t=0： 0.001:3; y=square(2pi10*t ，30); plot(t， y） ；grid on; axis(0,0.3,1。2，1。2 ） ；title（周期方波信号) （3) 用上面定义的uct 阶跃函数实现幅度为1、宽度为1

13、的门函数( )gt。程序：t=-1:0.01：1; ft=uct （t+0.5)uct（ t-0。5)；plot（t,ft） ；grid on；axis( 1,1，0。2,1。2）title（门函数 ) （ 4) 已 知 系 统 的 微 分 方 程 和 激 励 信 号 为( )4( )4 ( )( )3 ( )y ty ty tftf t，( )( )tf tet,试用 matlab命令绘出系统零状态响应的时域仿真波形图。源程序 : ts=0;te=8;dt=0。 01; sys=tf(1，3,1,4，4 ）; t=ts:dt:te; f=exp(-t) 。*uct(t) ；y=lsim(sy

14、s,f ，t）; plot(t，y）,grid on ；xlabel(t(s）)，ylabel( y(t)title( 零状态响应 ）（5) 已知系统的微分方程为( )3( )2 ( )( )yty ty tf t，试用matlab命令求系统冲激响应和阶跃响应的数值解,并绘出冲激响应和阶跃响应的时域仿真波形图. 源程序：t=0：0。001:4; sys=tf(1,1,3 ，2); h=impulse（sys,t）; s=step（sys，t）; subplot（211）plot（t,h） ，grid on xlabel(t(s）),ylabel （ h(t ） ）title（冲激响应 ) su

15、bplot(212) plot（t，s)，grid on xlabel（ t(s）)，ylabel（ s（t)）title( 阶跃响应 ) 实验二 连续时间 lti 系统的频率特性及频域分析一、实验目的1运用 matlab分析连续系统的频率特性；2运用 matlab进行连续系统的频域分析. 二、实验原理1. 连续时间 lti 系统的频率特性一个连续时间lti 系统的数学模型通常用常系数线性微分方程来描述，即nn1010dddd( )( )ddddnnnnyyxxaaa y tbbb x ttttt对上式两边取傅里叶变换，并根据傅里叶变换的时域微分特性，得到系统的频率响应为1010()()( )

16、( )( )()()mmnnbjbjbyhxajajamatlab信号处理工具箱提供的freqs 函数可直接计算系统的频率响应的数值解，其语句格式为h=freqs（b,a,w）其中 ,b 和 a 分别表示h(w)的分子和分母多项式的系数向量；w 为系统频率响应的频率范围，其一般形式为w1:p: w2，w1 为频率起始值,w2 为频率终止值，p 为频率取样间隔。h 返回 w 所定义的频率点上系统响应频率响应的样值。注意，h 返回的样值可能为包含实部和虚部的复数。因此,如果想得到系统的幅频特性或相频特性，还需利用abs和 angle 函数来分别求得。2. 连续时间 lti 系统的频域分析连续 lt

17、i 系统的频域分析法，也称为傅里叶变换分析法。该方法是基于信号频谱分析的概念，讨论信号作用于线性系统时在频域中求解响应的方法.傅里叶分析法的关键是求系统的频率响应。 傅里叶分析法主要用来分析系统的频率响应特性，或分析输出信号的频谱，也可用来求解正弦信号作用下的稳态响应。对于周期信号激励而言，可首先将周期信号进行傅里叶级数展开，然后求系统在各傅里叶级数分解的频率分量作用下系统的稳态响应分量，再由系统的线性性质将这些稳态响应分量叠加 ,从而得出系统总的响应.该方法的理论基础是基于正弦信号作用下系统的正弦稳态响应。对于正弦激励信号0sin()at,当经过系统()h,其稳态响应为00000( )sin

18、()()() sin()zsytatha htangle h三、实验内容1。 实例分析与验证(1) 已知一个连续时间lti 系统的微分方程为)(7)( 13)(5)( 8)( 10)( txtxtytytyty求系统的频率响应，并用matlab绘出其幅频特性和相频特性图。解:对微分方程取傅里叶变换，得7)(13)(5)(8)(10)(23jwwxjwjwjwwy因此，频率响应为5)(8)(10)(7)(13)()()(23jwjwjwjwwxwywhmatlab源程序：w=-3*pi:0 。01:3pi; b=13,7 ； a=1， 10,8，5；h=freqs（b，a,w） ；subplot

19、（211）plot(w ，abs(h)）,grid on xlabel（ omega(rad/s)） ,ylabel(h(omega) ）title（h（w）的频率特性）subplot（212) plot（w，angle（h)）,grid on xlabel（ omega(rad/s）)， ylabel(phi(omega）| ）title（h（w）的相频特性) （2) 如图为 rc 低通滤波器电路,在输入端加入矩形脉冲u1(t).利用傅里叶分析法求输出端电压 u2(t)。u1(t)t11解： rc 低通滤波器的频率响应为ru2(t)( ) tt2)120ku (c10fujwwh)(其中15

20、rc激励信号的傅里叶变换为1( )(1)/jwuwejw因此，响应的傅里叶变换为2125(1)5(1)()()()(5)5jjeeuwh w uwjjjmatlab源程序：w=-6*pi:0.01 ：6*pi; b= 5;a=1,5; h1=freqs(b,a,w); plot(w,abs（h1）)，grid on xlabel（ omega（rad/s) ),ylabel（|h(omega）| )title(rc 低通滤波电路的幅频特性）u1t=sym（heaviside（t)-heaviside （t1） ；u1w=simplify(fourier（ u1t） ）; figure subplot(221)，ezplot(u1t, 0.2，2)， grid on title（矩形脉冲信号) xlabel（ time(sec)）,ylabel（u1(t）) subplot(222)，ezplot(abs(u1w), 6pi 6*pi ）,grid on title(矩形脉冲的频谱）xlabel(omega（rad/s）),ylabel（x （ omega) )u2w=sy