实验5_连续时间系统的复频域分析报告.docx

一，实验目的 针对拉普拉斯变换及其反变换，了解定义、并掌握matlab实现 方法；掌握连续时间系统函数的定义和复频域分析方法 ；利用 MATLAB加深掌握系统零极点和系统分布。 二，实验原理 1?拉普拉斯变换 调用laplace和ilaplace函数表示拉氏变换和拉氏反变换 ： L=laplace(F)符号表达式F的拉氏变换，F中时间变量为t，返回变量为s的结果表 达式。 L=laplace(F,t)用t替换结果中的变量 s。 F=ilaplace(L)以s为变量的符号表达式 L的拉氏反变换，返回时间变量为t的结果 表达式。 F=ilaplace(L,x)用x替换结果中的变量 t。 2?连续时间系统的系统函数 3?连续时间系统的零极点分析 求多项式的根可以通过 roots来实现： r=roots(c) c为多项式的系数向量，返回值r为多项式的根向量。 绘制系统函数的零极点分布图 ，可调用pzmap函数： Pzmap(sys)绘出由系统模型sys描述的系统的零极点分布图。 [p，z]=pzmap(sys)返回极点和零点，不绘出分布图。 三，实验内容 ，试采用复频域的方法求(1 )已知系统的冲激响应 h(t)=u(t)-u(t-2),输入信号x(t)=u(t) ，试采用复频域的方法求 解系统的响应，编写MATLAB程序实现。 MATLAB程序如下： syms t h x y H X h = heaviside(t) - heaviside(t - 2) x = heaviside(t) H = laplace(h) X = laplace(x) Y = X*H y = ilaplace(Y) disp(y) ezplot(y,[-5,4]) title(h(t)) 程序执行结果如下： 冲) 所以解得 (2 )已知因果连续时间系统的系统函数分别如下 ： ① 一 ② 试采用matlab画出其零极点分布图，求解系统的冲激响应 h(t)和频率响应H(w)，并判断 系统是否稳定。 ① —— MATLAB程序如下： syms H s a = [1,2,2,1] H = tf(b,a) pzmap(H) axis([-2,2,-2,2]) figure impulse(H) 程序执行结果如下： Pole-Zero Map 5 0.5 0 -0 5 -1 *1.5 -2 -2 -15 / Q5 0 Q.5 1 1.5 Real Axis (seconds1) Impule Response- dJDn 三 cul12 14 dJDn 三 cul 12 14 2 4 6 ft 10 Time (swonds) 0,45 0 4 0.i5 0 3 0.25 0.2 0.15 0.1 0.05 D -005 该因果系统所有极点位于 s面左半平面，所以是稳定系统 ② MATLAB程序如下： b = [1,0,1] a=[1,2,-3,3,3,2] H = tf(b,a) figure pzmap(H) axis([-3.5,3.5,-3.5,3.5]) figure impulse(H) 程序执行结果如下： (cfipuauas) 5一 XV AJecu 一 *2-1012 Real Axis [seconds1) ajpn三d￡ ■2.50 Impute Response 1 0 5 0 -0 5 -1 -1.5 -2 60 70 20 30 40 50 Tims (seconds) 该因果系统的极点不全位于 S平面的左半平面，所以系统是不稳定系统 (3)已知连续时间系统函数的极点位置分别如下所示 ： 试用MATLAB绘制下述6种不同情况下，系统函数的零极点分布图，并绘制响应冲激 响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响 。 p=0 z =[] p = [0] k = [1] [b,a] = zp2tf(z,p,k) sys = tf(b,a) pzmap(sys) impulse(sys) impulse Response- cupn 七QLU 2 ta 16 14 12 1 os 0 02 0 04 O.OS 0 00 0.1 DJ2 O ld 0.1B 0 tfl CL2 Time (swonds) p=-2 z =[] P = [-2] k = [1] [b,a] = zp2tf( z, p,k) sys = tf(b,a) pzmap(sys) impulse(sys) cupn七-duuimpulse- Response0 9 0 8 0 7 0 6 0.5 0 4 0 3 cupn七-duu impulse- Response 0 9 0 8 0 7 0 6 0.5 0 4 0 3 0.2 0.1