实验报告基于Matlab的控制系统仿真1

1、实验一 基于Matlab的控制系统模型姓名 学号 班级机械一、实验目的1) 熟悉Matlab的使用环境，学习Matlab软件的使用方法和简单编程方法。2) 学习使用Matlab软件进行拉氏变换和拉式反变换的方法。3) 学习使用Matlab软件建立、转换连续系统数学模型的方法。4) 学习使用Matlab软件分析控制系统稳定性的方法。二、实验原理1. 拉氏变换和反拉氏变换 (1) 拉氏变换(2) 拉氏反变换2. 控制系统模型的建立和转化传递函数模型：零极点增益模型：(1) 建立系统传递函数模型(2) 建立系统的零极点模型(3) 传递函数模型转化为零极点模型(4) 零极点模型转化为传递函数模型3.

2、用Matlab进行传递函数部分分式展开4. 连续系统稳定性分析已知传递函数，试求该系统的闭环极点并判断系统的稳定性。(1)计算结果：-1.2888 + 0.4477i-1.2888 - 0.4477i-0.7244 + 1.1370i-0.7244 - 1.1370i0.1364 + 1.3050i0.1364 - 1.3050i0.8767 + 0.8814i0.8767 - 0.8814i由计算结果可知，该系统有4个极点具有正实部，因此系统不稳定。(2) 零极点分布图：三、实验内容1. 用Matlab求下两式的拉氏变换。syms tg1=1-0.5*t2*exp(-t)-t*exp(-t)

3、-exp(-t)laplace(g1)g1 =1-1/2*t2*exp(-t)-t*exp(-t)-exp(-t) ans =1/s-1/(1+s)3-1/(1+s)2-1/(1+s)syms tg2=0.6-0.3*exp(-t)*sin(-2*t)-0.6*exp(-t)*cos(2*t)laplace(g2)g2 =3/5+3/10*exp(-t)*sin(2*t)-3/5*exp(-t)*cos(2*t)ans =3/5/s+3/20/(1/4*(1+s)2+1)-3/20*(1+s)/(1/4*(1+s)2+1)2. 用Matlab求解拉式反变换。syms s g1=(s+2)/(s

4、*(s+1)2*(s+3)ilaplace(g1)g1 =(s + 2)/(s*(s + 1)2*(s + 3)ans =exp(-3*t)/12 - (3*exp(-t)/4 - (t*exp(-t)/2 + 2/3syms sg2=(s2+3*s+5)/(s2+s+2)ilaplace(g2)g2 =(s2 + 3*s + 5)/(s2 + s + 2)ans =dirac(t) + 2*exp(-t/2)*(cos(7(1/2)*t)/2) + (2*7(1/2)*sin(7(1/2)*t)/2)/7)3. 用Matlab将以下零极点模型转换为一般多项式传递函数模型。z=-1p=0,0,

5、-2,-5k=5Gs1=zpk(z,p,k)num,den=zp2tf(z',p',k)Gs2=tf(num,den)z = -1p = 0 0 -2 -5k = 5Gs1 = 5 (s+1) - s2 (s+2) (s+5)Continuous-time zero/pole/gain modelnum = 0 0 0 5 5den = 1 7 10 0 0Gs2 = 5 s + 5 - s4 + 7 s3 + 10 s24. 用Matlab将以下传递函数转换为零极点模型。 num=1,7,24,24den=1,2,3,4,5,6,7,8,9Gs1=tf(num,den)z,p

6、,k=tf2zp(num,den)Gs2=zpk(z,p,k)num = 1 7 24 24den = 1 2 3 4 5 6 7 8 9Gs1 = s3 + 7 s2 + 24 s + 24 - s8 + 2 s7 + 3 s6 + 4 s5 + 5 s4 + 6 s3 + 7 s2 + 8 s + 9Continuous-time transfer function.z = -2.7306 + 2.8531i -2.7306 - 2.8531i -1.5388 + 0.0000ip = -1.2888 + 0.4477i -1.2888 - 0.4477i -0.7244 + 1.1370

7、i -0.7244 - 1.1370i 0.1364 + 1.3050i 0.1364 - 1.3050i 0.8767 + 0.8814i 0.8767 - 0.8814ik =1Gs2 = (s+1.539) (s2 + 5.461s + 15.6) - (s2 + 2.578s + 1.861) (s2 - 1.753s + 1.545) (s2 + 1.449s + 1.817) (s2 - 0.2728s + 1.722) 5. 用Matlab对以下传递函数按部分分式展开。num=1,5,9,7den=conv(1, 1,1 ,2)Gs1=tf(num,den)r,p,k=resid

8、ue(num,den)num = 1 5 9 7den = 1 3 2Gs1 = s3 + 5 s2 + 9 s + 7 - s2 + 3 s + 2Continuous-time transfer function.r = -1 2p = -2 -1k = 1 2Gs=2s+1-1s+2+s+26. 已知传递函数，试用两种方法求该系统的闭环极点并判断系统的稳定性。方法一：根据稳定的充分必要条件判别线性系统的稳定性， 最简单的方法是求出系统所有极点，并观察是否含有实部大于0的极点，如果有，系统则不稳定。（1）den=1 1 2 2 3 5roots(den)den = 1 1 2 2 3 5ans = 0.7207 + 1.1656i 0.7207 - 1.1656i -0.6018 + 1.3375i -0.6018 - 1.3375i -1.2378 + 0.0000i由计算结果可知，该系统有2个极点具有正实部，因此系统不稳定。方法2：（2）num=1 -15 126den=1 1 2 2 3 5pzmap(num,den)title('Pole-Zero Map')零极点分布图：从计算机结果以及零极点图可以看出,该系统的极点并