实验2 系统的时域分析

《实验2 系统的时域分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、重庆三峡学院《信号与系统分析》实验实验2系统的时域分析1实验目的1）学会利用MATLAB求解连续系统的零状态响应、冲激响应和阶跃响应；2）学会利用MATLAB求解离散系统的单位取样响应；3）学会利用MATLAB求解离散系统的卷积和。2实验原理及实例分析（实验原理见教材的第二章和第三章。）2.1连续系统零状态响应的数值求解例1已知某LTI系统的微分方程为，其中，。试用MATLAB命令绘出范围内系统零状态响应的波形图。解：程序如下：clc;%命令窗口清屏closeall;clearall;t=0:0.01:5;sys=tf([6],[156]);%用传输函数形式表示系统f=10*sin(2*p

2、i*t).*uCT(t);y=lsim(sys,f,t);%对输入信号模拟仿真plot(t,y,'linewidth',2);gridon;xlabel('t(sec)');title('y(t)');产生的图形如图1所示。图1例1程序产生的图形5重庆三峡学院《信号与系统分析》实验2.2连续系统的冲激响应和阶跃响应的数值求解例2已知某LTI系统的微分方程为，试用MATLAB命令绘出范围内系统的冲激响应和阶跃响应。解：MATLAB程序如下：clc;closeall;clearall;t=0:0.01:5;sys=tf([116],[1232]);h=impulse(sys,t);%计算系统的

3、冲激响应g=step(sys,t);%计算系统的阶跃响应subplot(2,1,1);plot(t,h,'Linewidth',2);gridon;xlabel('t(sec)');title('Impulseresponse--h(t)');subplot(2,1,2);plot(t,g,'Linewidth',2);gridon;xlabel('t(sec)');title('Stepresponse--g(t)');图2例2程序产生的图形2.3离散系统的响应例3已知系统的差分方程为，试用MATLAB命令绘出当激励信号为时，该系统的零状态响应。解：MATLAB程序为：clc;close

4、all;clearall;a=[3-42];b=[12];n=0:30;x=(1/2).^n;y=filter(b,a,x);5重庆三峡学院《信号与系统分析》实验stem(n,y,'fill');gridon;xlabel('n');title('y(n)');程序产生的图形如图3所示。图3例3程序产生的结果2.4离散系统的单位取样响应例4已知系统的差分方程为，利用MATLAB命令绘出该系统的单位取样响应。解：MATLAB程序如下：clc;closeall;clearall;a=[3-42];b=[12];impz(b,a,30);gridon;程序产生的图形如图4所示。图4例4程序产生的

5、图形5重庆三峡学院《信号与系统分析》实验2.5离散卷积和的计算例5已知系统的单位取样响应为，试用MATLAB计算当输入信号为时，系统的零状态响应。解：MATLAB程序如下：clc;closeall;clearall;nx=-1:5;nh=-2:10;x=uCT(nx)-uCT(nx-4);h=0.8.^nh.*(uCT(nh)-uCT(nh-8));y=conv(x,h);ny=(min(nx)+min(nh)):(max(nx)+max(nh));subplot(3,1,1);stem(nx,x,'fill');gridon;xlabel('n');title('x(n)');axis(

6、[-41603]);subplot(3,1,2);stem(nh,h,'fill');gridon;xlabel('n');title('h(n)');axis([-41603]);subplot(3,1,3);stem(ny,y,'fill');gridon;xlabel('n');title('y(n)=x(n)*h(n)');axis([-41603]);程序产生的图形如图5所示。图5例5程序产生的图形3实验报告与要求5重庆三峡学院《信号与系统分析》实验在实验报告中，请简要说明连续系统的零状态响应、冲激响应和阶跃响应等概念，说明离散系统的单位取样响应、卷积和的计算方法。思考并完成下列

7、各题：1）已知系统的微分方程为，如何利用MATLAB命令绘出该系统的冲激响应和阶跃响应的时域波形图。clccloseallclearallt=0:0.01:10;sys=tf([10],[122]);h=impulse(sys,t);subplot(2,1,1);plot(t,h,'linewidth',2);grid;axis([010-0.21.2]);xlabel('t(sec)');title('Impulser