1、实验二:1、函数 initial(A,B,C,D,x0)可以得到系统输出对初始状态 x0 的时间响应;2、函数 step(A,B,C,D)给出了系统的单位阶跃响应曲线;3、函数 impulse(A,B,C,D) 给出了系统的单位脉冲响应曲线;4、函数 y,T,x=lsim(sys,u,t,x0) 给出了一个状态空间模型对任意输入u 的响应,其中的 sys 表示贮存在计算机内的状态空间模型,它可以由函数 sys=ss(A,B,C,D)得到,x0 是初始状态实验三:MATLAB 提供了计算离散化状态空间模型中状态矩阵和输入矩阵的函数:G,H=c2d(A,B,T) 其中的 T 是离散化模型的采样周期
2、。实验四:MATLAB软件提供了两个函数acker和place来确定极点配置状态反馈控制器的增益矩阵K。函数acker是基于求解极点配置问题的爱克曼公式,它只能应用到单输入系统,要配置的闭环极点中可以包括多重极点。函数acker和place的一般形式是:K=acker(A,B,J)K=place(A,B,J)其中的J是一个向量,是n个期望的闭环极点。得到了所要求的反馈增益矩阵后,可以用命令eig(A-B*K)来检验闭环极点。实验五:应用MATLAB中极点配置的函数来确定所需要的观测器增益矩阵。例如,对于单输入单输出系统,观测器的增益矩阵可以由函数L=(acker(A,C,V)得到。其中的V是由
3、期望的观测器极点所构成的向量。类似的,也可以用L=(place(A,C,V)来确定一般系统的观测器矩阵,但这里要求V不包含相同的极点。对于降阶观测器的设计,使用MATLAB软件中的函数L=(acker(Abb,Aab,V)或L=(place(Abb,Aab,V) 可以得到观测器的增益矩阵L。其中的V是由降阶观测器的期望极点所组成的向量。三、仪器设备PC计算机1台(要求P4-1.8G以上),MATLAB6.X或MATLAB7.X软件1套。四、线路示图五、内容步骤1、根据所给系统的传递函数或(A、B、C、D),依据系统的传递函数阵和状态空间模型之间的关系,采用MATLAB的相关函数编写m-文件。2
4、、求取系统的状态和输出响应;3、在MATLA界面下调试程序,并检查是否运行正确。六、数据处理实验一程序代码:num=0 0 1 2;0 1 5 3den=1 2 3 4A,B,C,D=tf2ss(num,den)结果:A = -2 -3 -4 1 0 0 0 1 0B = 1 0C = 0 1 2 1 5 3D =实验二A=0 1;-1 0;B=0;0;D=B;C=1 0;0 1;x0=0;1;y,x,t=initial(A,B,C,D,x0);plot(t,x(:,1),t,x(:,2)gridtitle(Response to Initial Condition)xlabel(Time(s
5、ec)ylabel(x1,x2text(0.55,1.15,x1text(0.4,-2.9,x2syms t;F=expm(A*t)F = cos(t), sin(t) -sin(t), cos(t)实验三A=-1 0;0 -2;B=1;G,H=c2d(A,B,1)G = 0.3679 0 0 0.1353H = 0.6321G,H=c2d(A,B,0.5) 0.6065 0 0 0.3679 0.3935实验四4.1A=0 1 0;0 0 1;-1 -5 -6;0;J=-2+j*4 -2-j*4 -10;sys=ss(A-B*K,0;0,eye(3),0);t=0:0.01:4;x=init
6、ial(sys,1;0,t);x1=1 0 0*x;x2=0 1 0*xx3=0 0 1*xsubplot(3,1,1);plot(t,x1),gridsubplot(3,1,2);plot(t,x2),gridsubplot(3,1,3);plot(t,x3),gridt (sec)x34.5-3 -4;C=3 2;D=0;J=-4 -5;x0=1;0,eye(2),0);x1=1 0*xx2=0 1*x实验五5.1a=0 1 0;1.244 0.3965 -3.145;b=0;1.244;c=1 0 0;v=-5+j*5*sqrt(3) -5-j*5*sqrt(3) -10;l=(acke
7、r(a,c,v)Warning: Pole locations are more than 10% in error.l = 16.8550 147.3875 544.39325.2Aaa=0;Aab=1 0;Aba=0;Abb=0 1;0.3965 -3.145;Ba=0;Bb=0;v=-5+j*5*sqrt(3) -5-j*5*sqrt(3);l=(acker(Abb,AabAhat=Abb-l*AabBhat=Ahat*l+Aba-l*AaaFhat=Bb-l*Ba 6.8550 78.8375Ahat = -6.8550 1.0000 -78.4410 -3.1450Bhat = 31
8、.8465 -784.4132Fhat = 1.24405.5num=4;den=1 2 0;A,B,C,D=tf2ss(num,den);J=-2+j*2*sqrt(3) -2-j*2*sqrt(3);K=acker(A,B,J);V=-8 -8;L=(acker(A,C,V)A0=A-B*K B*K;zeros(2,2) A-L*C;sys=ss(A0,eye(4),eye(4),eye(4);5;1;x1=1 0 0 0*xx2=0 1 0 0*xe1=0 0 1 0*xe2=0 0 0 1*xfigure(1)subplot(2,1,1);Response to initial con
9、ditionsubplot(2,1,2);figure(2)plot(t,e1),gride1plot(t,e2),gride2L = 9.00003.5000七、分析讨论 通过本次试验,我学习了系统状态空间模型的建立方法、了解了状态空间模型与传递函数相互转换的方法,通过编程、上机调试,掌握系统状态空间模型与传递函数相互转换的方法。同时我学会了根据状态空间模型分析系统由初始状态和外部激励所引起的响应,通过编程、上机调试,掌握系统运动的分析方法。掌握了建立基于对象的一个连续时间状态空间模型,并导出其相应的离散化状态空间模型,并通过编程、上机调试,掌握离散系统运动分析方法。配置极点后,系统的稳态性能变好了,但是稳态误差变大了。为了消除几点配置对稳态性能的负面影响,对于一个n维系统,必须制定n个实极点或共轭复极点,在确定极点位置的时候,要充分考虑他们对于系统性能的主导影响及其余系统零点分布状况关系,还要兼顾系统抗干扰能力和对参数飘逸低敏感性的要求。通过实验五我学习掌握了观测器设计计算方法,通过编程上机调试,掌握了全阶观测器的输出反馈控制系统设计方法。
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2