《信号与系统》实验3(连续系统的s域分析).docVIP

实验三 连续系统的s域分析 学号： 姓名： 成绩： 一、实验目的 （1）熟悉拉氏变换。 （2）掌握系统响应s域求法。 （3）熟悉系统的频率响应。 二、实验原理 连续LTI系统，在s域可以用系统函数H(s)描述，其实质是系统冲激响应h(t)的拉氏变换。 拉氏逆变换 若H(s)的极点分别为p1,…,pn，则H(s)可表示为：部分分式+多项式 由此可以方便的求出其拉氏逆变换（即对应的时间域信号）。 （2）s域求响应 变换到s域，系统响应等于激励信号与系统函数相乘 （3）系统的频率响应 如果系统函数H(s)的收敛域包含虚轴，则令s=jω，得到系统的频率响应H(jω)。 三、验证性实验 已知系统，其系统函数为。 求零、极点。 程序： 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); 问题：该系统的零点能够抵销什么形式的激励信号？ 求冲激响应h(t)（系统函数的逆变换） 程序： clear; b=[3,9]; %分子多项式系数 a=[1,6,8]; %分母多项式系数 [r,p,k]=residue(b,a) 运行结果： r = 1.5000 1.5000 p = -4 -2 k = [] 则 问题：该系统是不是稳定系统？ e(t)=u(t)时，求零状态响应 程序： clear; b=[3,9]; %分子多项式系数 a=[1,6,8,0]; %分母多项式系数 [r,p,k]=residue(b,a); %求留数、极点 t=0:0.1:10; f=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t); plot(t,f); 问题：响应的极点有哪些，与激励相同的极点是哪一个，对应着响应的什么分量？ 求频率响应H(jω)。 程序： b=[3 9]; %分子多项式的系数向量 a=[1 6 8]; %分母多项式的系数向量 w=0:0.01:100; %生成角频率w的矢量 H =freqs(b,a,w); subplot(2,1,1),plot(w,abs(H)); %画幅频特性 title(abs(H(jw)));grid on subplot(2,1,2),plot(w,angle(H)); %画相频特性 title(angle(H(jw)));grid on 问题：考察一个系统的频响特性，对其系统函数有什么要求？ 根据频响判断，这是一个什么类型的滤波器？ 四、设计性实验 已知系统，当e(t)=cos(t)u(t)时，写出其系统函数，利用拉氏变换求系统的冲激响应和零状态响应。零状态响应是否有界？该系统是否稳定？ 答：零状态有界，振荡变化。 由题可知：系统函数H（s）=1/(s*s+s); 源函数： 1，冲激响应 clear; b=[1]; a=[1,1,0]; [r,p,k]=residue(b,a) 运行结果： r = -1 1 p = -1 0 k = [] H(s)=-1/(s+1)+1/s; H(t)=-exp(-t)+u(t) 由：e(t)=cos(t)u(t); E(s)=L[cos(t)u(t)]=s/(s*s+1); R(s)=E(s)*H(s)=1/(s^3+s*s+s+1); 2，零状态响应 源函数： clear; b=[1]; %分子多项式系数 a=[1,1,1,1]; %分母多项式系数 [r,p,k]=residue(b,a); %求留数、极点 t=0:0.1:100; f=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t); plot(t,f); 运行结果： 频率响应： 3，,源函数： b=[1]; %分子多项式的系数向量 a=[1,1,1,1]; %分母多项式的系数向量 w=0:0.01:100; %生成角频率w的矢量 H =freqs(b,a,w); subplot(2,1,1),plot(w,abs(H)); %画幅频特性 title(abs(H(jw)));grid on subplot(2,1,2),plot(w,angle(H)); %画相频特性 title(angle(H(jw)));grid on 运行结果：