连续系统的s域分析实验报告.docxVIP

连续系统的s域分析实验报告 　　实验四连续系统的S域分析实验　　一实验目的　　学习MATLAB语言的编程方法及熟悉MATLAB指令掌握符号运算函数，虚数单位变量以及结构的概念。用软件函数laplace和ilaplace拉氏变换和逆变换。二实验内容　　信号与系统MATLAB综合实验教材第五章练习题三。三实验编程　　symss;　　h=ilaplace(+1/(s+2)-4/(s+3));　　b=[,,2];　　a=[1,5,6];　　figure;　　subplot(2,1,1);　　pzmap(b,a);　　subplot(2,1,2);　　impulse(b,a);　　figure(2);　　freqs(b,a);　　h(t)=exp(-2*t)-4*exp(-3*t)+(3*dirac(t))/2;　　四实验结果　　五实验总结　　通过本实验学会了软件学会了用MATLAB届微分方程，计算出自由项以及各个分式的系数和极点，另外能绘制零、极点分布图和频率响应曲线。在今后的实验中，要多了练习使用pzmap和freqs。　　实验三连续系统的s域分析　　学号：姓名：成绩：　　一、实验目的　　熟悉拉氏变换。　　掌握系统响应s域求法。熟悉系统的频率响应。　　二、实验原理　　连续LTI系统，在s域可以用系统函数H(s)描述，其实质是系统冲激响应h(t)的拉氏变换。　　H(s)?　　B(s)　　A(s)　　拉氏逆变换　　若H(s)的极点分别为p1,…,pn，则H(s)可表示为：部分分式+多项式　　M　　r1r2　　H(s)?????????cmsm　　s?p1s?p2s?pnm?0　　由此可以方便的求出其拉氏逆变换。　　s域求响应　　变换到s域，系统响应等于激励信号与系统函数相乘　　R(s)?E(s)H(s)　　系统的频率响应　　如果系统函数H(s)的收敛域包含虚轴，则令s=jω，得到系统的频率响应H(jω)。　　三、验证性实验已知系统r　　(2)　　(t)?6r(1)(t)?8r(t)?3e(1)(t)?9e(t)，其系统函数为H(s)?　　3s?9　　。2　　s?6s?8　　求零、极点。程序：clear;　　b=[3,9];%分子多项式系数a=[1,6,8];%分母多项式系数zs=roots(b);ps=roots(a);　　plot(real(zs),imag(zs),go,real(ps),imag(ps),rx);grid;　　legend(zero,pole);　　-4　　-问题：该系统的零点能够抵销什么形式的激励信号？　　求冲激响应h(t)程序：clear;　　b=[3,9];%分子多项式系数a=[1,6,8];%分母多项式系数[r,p,k]=residue(b,a)运行结果：r=　　p=-4-2k=[]则　　H(s)?　　s?4?　　s?2　　h(t)??4t??2t问题：该系统是不是稳定系统？　　e(t)=u(t)时，求零状态响应　　E(s)?L[u(t)]?　　1s　　R(s)?H(s)E(s)?3s?9　　s3　　?6s2?8s　　程序：clear;　　b=[3,9];%分子多项式系数a=[1,6,8,0];%分母多项式系数　　-2　　[r,p,k]=residue(b,a);%求留数、极点t=0::10;　　f=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t);plot(t,f);　　1　　问题：响应的极点有哪些，与激励相同的极点是哪一个，对应着响应的什么分量？　　求频率响应H(jω)。程序：　　b=[39];%分子多项式的系数向量a=[168];%分母多项式的系数向量w=0::100;%生成角频率w的矢量H=freqs(b,a,w);　　subplot(2,1,1),plot(w,abs(H));%画幅频特性title(abs(H(jw)));gridon　　subplot(2,1,2),plot(w,angle(H));%画相频特性title(angle(H(jw)));gridon　　abs(H(jw))　　　　20　　40　　60　　80　　100　　　　angle(H(jw))　　-1　　-2　　0　　问题：考察一个系统的频响特性，对其系统函数有什么要求？根据频响判断，这是一个什么类型的滤波器？　　四、设计性实验　　已知系统r(2)(t)?r(1)(t)?e(t)，当e(t)=cos(t)u(t)时，写出其系统函数，利用拉氏变换求系统的冲激响应和零状态响应。零状态响应是否有界