信号与系统信号分析与处理实验报告2012.docVIP

实验四 连续系统分析 姓名： 闫寒旭 学号： 班级： 电信102 一 实验目的 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。 二 实验内容 1、已知连续系统的系统函数为 (1)计算并画出该系统的单位冲激响应h(t),单位阶跃响应g(t) (2)求输入分别为 时，计算并画出系统的零状态输出响应y(t)， (3)计算该系统的零极点，并作出零极点分布图； (4)画出并分析系统的幅频特性。 %代码（1）单位冲击响应 a=[1,3,2]; b=[4,1]; sys=tf(b,a); t=0:0.1:10; h=impulse(sys, t); plot(h); xlabel(t); title(h(t)) h(t)的值： 4 3.799 2.553 1.303 1.363 0.2196 0.3335 0.7017 0.09231 -0.06915 单位阶跃响应： a=[1,3,2]; b=[4,1]; sys=tf(b,a); t=0:0.1:10; g=step(sys, t); plot(g); xlabel(t); title(g(t)) g(t)的值 0 0.4942 0.9208 0.6061 0.6642 1.785 1.937 1.861 1.037 1.624 （2）当输入为u（t）时 a=[1,3,2]; b=[4,1]; sys=tf(b,a); t=0:0.1:10; x=(t0); y=lsim(sys, x, t) ; plot(t,y); xlabel(t); title(y(t)) 输出值： 0 0 0.4942 0.9208 0.6061 0.6642 1.785 1.937 1.861 1.037 当输入为sin(t)u(t)时： a=[1,3,2]; b=[4,1]; sys=tf(b,a); t=0:0.1:10; x=sin(t); y=lsim(sys, x, t) ; plot(t,y); xlabel(t); title(y(t)) 输出值： 0 0.41379 0.93625 0.5931 0.6776 0.8556 0.7692 0.7481 0.1587 0.4638 当输入为exp（-t）u（t）时 a=[1,3,2]; b=[4,1]; sys=tf(b,a); t=0:0.1:10; x=exp(-t); y=lsim(sys, x, t) ; plot(t,y); xlabel(t); title(y(t)) 输出值： 0 0.7128 0.5628 0.9808 0.5535 0.8469 0.8419 0.0608 0.7792 0.5202 （3） b=[4,1];a=[1,3,2]; z=roots(b); sys=tf(b,a); pzmap(sys) 零极点有：-2，-1，-0.25 （4）。b=[4,1]; a=conv([1,1],[1,2]); [H,w]=freqs(b,a) ; plot(w,abs(H)); xlabel(Frequency(rad/s)); ylabel(Amplitude); title(Magnitude response);