信号系统实验3.docVIP

实验三 连续LTI系统的复频域分析 集成111 王茂林 1108040102 一. 实验目的 掌握基于MATLAB的拉普拉斯变换和反变换分析应用。进一步了解用MATLAB计算复杂系统的方法。　　 掌握用MATLAB分析并绘制连续系统零极点图以判断因果系统稳定的方法。 掌握用MATLAB实现连续系统的频率特性及其幅度特性、相位特性。 二. 实验原理 1.拉普拉斯变换和反变换的符号运算 在MATLAB符号运算工具箱中，提供了拉普拉斯正变换和反变换的函数。 正变换的调用格式为 F = laplace(f) 其中f为时间函数的符号表达式，F为拉普拉斯变换式，也是符号表达式。 反变换的调用格式为 f = ilaplace(F) 其中F为拉普拉斯变换式的符号表达式，f为时间函数，是符号形式。 2.连续系统的零极点分析 系统函数H(s)通常是一个有理分式，其分子和分母均为多项式。计算H(s)的零极点可以应用MATLAB的roots函数，求出分子和分母多项式的根即可。绘制系统的零极点分布图可以根据已求出的零极点，利用plot语句画图，还可以由H(s)直接应用pzmap函数画图。 pzmap函数的调用形式 pzmap(sys) 表示出sys所描述系统的零极点图。LTI系统模型要借助tf函数获得，其调用方式为 sys=tf(b,a) 式中b和a分别表示系统函数H(s)分子多项式和分母多项式的系数向量。 MATLAB提供了专门对连续系统频率响应H(jw)进行分析的函数freqs()。该函数可以求出系统频率响应的数值解，并可绘制出系统的幅频及相频响应曲线。 freqs()函数调用格式: freqs(b,a) 该格式并不返回系统频率响应的样值，而是以对数坐标的方式绘出系统的幅频响应和相频响应曲线。= abs(h)：求模数，求得系统的幅度频率响应，返回值存于hm之中。 real(h)：求h的实部； imag(h)：求h的虚部； phase = atan(-imag(h)./(real(h)+eps))：求相位频率响应特性，atan()用来计算反正切值；或 phase = angle(h)：求相位频率响应特性； 三. 实验内容 1.求函数的拉普拉斯变换。 解：实现1所求运算的m文件如下， syms t ; e = sym(Heaviside(t)); F=laplace (exp(-2*t)*cos(10*pi*t)*e) 运行结果： F = 1/100/pi^2*(s+2)/(1/100*(s+2)^2/pi^2+1) 2.求的原函数f(t)。 实现2所求运算的m文件如下， syms s f=ilaplace((s+３)/((s^2))) 运行结果： f = exp(-3*t)*cos(t) 3.已知系统函数为，画出该系统的零极点分布图，并判断系统是否稳定。若稳定，再求解系统的幅频、相频特性以及幅频的实部和虚部。 解：实现3所求运算的m文件如下， b=[1]; a=[1 2 2 1]; figure(1) pzmap(b,a); p=roots(a); pxm=max(real(p)); if pxm=0 系统不稳定 else figure(2) freqs(b,a) end 运行结果： 零极点分布图 幅频、相频特性 四. 实验报告 1. 报告组成部分：实验目的、实验原理、实验内容、实验结果及分析。 2. 在实验结果中列出全部m文件。 3. 简述心得体会。 信号与系统实验讲义-青岛科技大学信息学院信息工程教研室编 1