《计算机控制实验报告.docx》由会员分享，可在线阅读，更多相关《计算机控制实验报告.docx（28页珍藏版）》请在金锄头文库上搜索。
1、实验一：MATLAB控制系统设计工具箱的使用一、实验目的：了解和掌握MATLAB的工具箱中用于离散控制系统的各个函数的用法。二、实验设备及软件计算机 MATLAB软件三、实验内容1、当T等于1秒时用matlab验证Gs=的广义Z传递函数（P37 例2.14） 2、求下列函数的Z变换 （P37 第1题）（1） （2） k=0 （3） （4） 单位阶跃信号的Z变换 （5）单位速度信号的Z变换 （6）单位加速度信号的Z变换(7) 广义Z变换 延迟0.25 T 的速度信号的Z变换 3、求下列函数的初值和终值 （P37 第2题）（1）：（2）： 4、求下列各函数的Z反变换。（P38 第3题）（1）：（2
2、）：5、 试求如图所示的采样控制系统在单位阶跃信号作用下的输出响应y（t）。设G(s)=，采样周期T=0.1s。 （P56 第1题）6求如图所示的采样控制系统在单位速度作用下的稳态误差。设G(s)=，采样周期T=0.1s。 （P56 第2题）7、分析系统稳定性 （P56 第5、6题）8、设某系统的Z传递函数为，求状态空间表达式。（P150 第2题）9. 求离散化状态空间方程 (P151 第3题)10. 求传递函数和特征值 (P151 第4题)11、.设离散系统的系数矩阵为A= ，试根据系统稳定的充要条件确定该系统的稳定性。 (P151 第6题)12、设离散系统的系数矩阵为 A= 试用Liapu
3、nov法确定该系统的稳定性。(P151 第7题)13.试确定下列离散系统的可控性 (P151 第7题)（1）A= ，B=14.试确定下列离散系统状态的可测性。(P151 第7题)（1）A= ，C= 15.设被控对象的状态空间方程为 （P208 第8题） X(k+1)= x(k)+ u(k) y(k)= 1x(k) 试用极点配置法确定状态反馈矩阵K，使状态反馈闭环系统的特征值为0.4和0.7，并画出状态反馈系统方块图四、 实验结果及分析1、当T等于1秒时用matlab验证Gs=的广义Z传递函数（P37 例2.14） num=1; den=1 1; Gs=tf(num,den,inputdelay
4、,0.75); Gz=c2d(Gs,1) Transfer function:0.2212 z + 0.4109- z2 - 0.3679 z Sampling time: 12、求下列函数的Z变换 （P37 第1题）（1） （2） k=0 （3） （4） 单位阶跃信号的Z变换 （5）单位速度信号的Z变换 （6）单位加速度信号的Z变换(7) 广义Z变换 延迟0.25 T 的速度信号的Z变换（1） syms a n T;FZ=(ztrans(1-exp(-a*n*T) FZ = z/(z-1)-z/exp(-a*T)/(z/exp(-a*T)-1)（2） syms k;FZ=ztrans(1/4
5、)k) FZ =4*z/(4*z-1) FZ=ztrans(sys,z) FZ = z/exp(-log(4)*t)/(z/exp(-log(4)*t)-1)（3） syms sft=ilaplace(6/(s*(s+2) ft = 6*exp(-t)*sinh(t) FZ=ztrans(ft) FZ = 6*(-z+z/exp(-1)*exp(1)/(2*z2/exp(-1)2-2*z-2*z/exp(-1)*exp(1)+2*exp(1)*exp(-1)（4） syms n fz=ztrans(n/n) fz = z/(z-1)（5） syms n syms T fz=ztrans(n*T
6、) fz = T*z/(z-1)2（6） syms n T fz=ztrans(1/2)*(n*T)2) fz = 1/2*T2*z*(z+1)/(z-1)3（7） syms n T fz=ztrans(n*T+0.75*T) fz = T*z/(z-1)2+3/4*T*z/(z-1) 3、求下列函数的初值和终值 （P37 第2题）（1）：syms z FZ=(10*z(-1)/(1-z(-1)2) FZ = 10/z/(1-1/z)2 limit(FZ,z,inf) ans = 0 limit(FZ*(z-1),z,1) ans = NaN（2）： syms z FZ=z2/(z-0.8)/
7、(z-0.1) FZ = z2/(z-4/5)/(z-1/10) limit(FZ,z,inf) ans = 1 limit(FZ*(z-1),z,1) ans = 04、 求下列各函数的Z反变换。（P38 第3题）（1）： syms z fz=z/(z-0.5) fz = z/(z-1/2) iztrans(fz) ans = (1/2)n（2）： syms z fz=z2/(z-0.8)*(z-0.1) fz = z2/(z-4/5)/(z-1/10) iztrans(fz) ans = 8/7*(4/5)n-1/7*(1/10)n 5、 试求如图所示的采样控制系统在单位阶跃信号作用下的输
8、出响应y（t）。设G(s)=，采样周期T=0.1s。 （P56 第1题） syms z fz=z2/(z-0.8)*(z-0.1) fz = z2/(z-4/5)/(z-1/10) iztrans(fz) ans = 8/7*(4/5)n-1/7*(1/10)n Gs=tf(20,1 10 0); Gz=c2d(Gs,0.1) Transfer function: 0.07358 z + 0.05285-z2 - 1.368 z + 0.3679 Sampling time: 0.1 ft=feedback(GZ,1) Transfer function: 0.07358 z + 0.0528
9、5-z2 - 1.294 z + 0.4207 Sampling time: 0.1 step(ft)6求如图所示的采样控制系统在单位速度作用下的稳态误差。设G(s)=，采样周期T=0.1s。 （P56 第2题）gs=tf(1,0.1 1 0);T=0.1;gz=c2d(gs,T,imp);gzb=feedback(gz,1); rz = tf(0.1 0,1 -2 1,T); rz1 = zpk(0,1 1,T,T); yz=rz*gzb;impulse(yz);t=0:0.1:10; ramp=t;lsim(gzb,ramp,t)y,t1 = lsim(gzb,ramp,t);ER = r
10、amp - yplot(ER,t),grid gs=tf(1,0.1 1 0); gsb=feedback(gs,1);rs = tf(1,1 0 0); ys=rs*gsb;t1=0:0.01:10;impulse(ys,t1);t=0:0.01:10; ramp=t;lsim(gsb,ramp,t)7、分析系统稳定性 （P56 第5、6题）gs=tf(1,0.1 1 0);T=0.1;gz=c2d(gs,T,imp);gzb=feedback(gz,1); rz = tf(0.1 0,1 -2 1,T); rz1 = zpk(0,1 1,T,T); yz=rz*gzb;impulse(yz
11、);t=0:0.1:10; ramp=t;lsim(gzb,ramp,t)y,t1 = lsim(gzb,ramp,t);ER = ramp - yplot(ER,t),grid gs=tf(1,0.1 1 0); gsb=feedback(gs,1);rs = tf(1,1 0 0); ys=rs*gsb;t1=0:0.01:10;impulse(ys,t1);t=0:0.01:10; ramp=t;lsim(gsb,ramp,t)gs=tf(1,1 1 0);T=1;gz=c2d(gs,T,imp);gzb=feedback(gz,1);pzmap(gzb)gz1=tf(1,45 -117 -119 -39,1);pzmap(gz1)8、设某系统的Z传递函数为，求状态空间表达式。（P150 第2题） syms z num=1,0.4; dun=1,0.7,0.06; FZ=tf(num,dun,1); sys=ss(FZ) a = x1 x2 x1 -0.7 -0.24 x2 0.25 0 b = u1 x1 2 x2 0 c = x1 x2 y1 0.5 0.8 d = u1 y1 0 Sampling time: 1Discrete-time model.9，10. 求离散化状态空间方程及传递函数和特征值 (P151 第3题) gz=ss(0 1;0 -2,0;1,1
