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

实验四 连续系统分析 姓名： 闫寒旭 学号： 1004020227 班级： 电信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.01660301743799 2.23604806301553 1.61922679061303 1.13434261071363 0.755564109062196 0.461924575103335 0.236422836217017 0.0652887336109231 -0.0626167616706915 单位阶跃响应： 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.348954618334942 0.610072098109208 0.801613935716061 0.938308763696642 1.03201393503785 1.09225516658937 1.12666653757861 1.14134907937037 1.14116287044624 （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.348954618334942 0.610072098109208 0.801613935716061 0.938308763696642 1.03201393503785 1.09225516658937 1.12666653757861 1.14134907937037 当输入为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.0182361347241379 0.0665742212393625 0.136706266755931 0.221747200266776 0.315985111668556 0.414679525767692 0.513898703167481 0.610388586551587 0.701467343894638 当输入为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.331571684607128 0.548093157685628 0.677874493159808 0.743172734015535 0.761396952468469 0.746082689758419 0.707678921600608 0.654182729227792 0.591650349325202 （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 r