线性系统理论基础试验三

信控学院上机实验实验报告课程 线性系统理论基础 实验日期2016年_5_月25日专业班级 姓名学号同组人实验名称 状态反馈极点配置方法的研究 评分批阅教师签字 一、实验目的.掌握状态反馈系统的极点配置；.研究不同配置对系统动态特性的影响。二、实验环境MATLAB6.5三、实验内容、源程序代码、实验数据及结果分析原系统如图3-2所示。图中，XJDX2是可以测量的状态变量。Xz~~I~~Xi图3-2系统结构图试设计状态反馈矩阵K=匡】及口],使系统加入状态反馈后其动态性能指标满足给定的要求状态反馈后的系统，如图3-3所示：第25页共37页信控学院上机实验图3-3状态反馈后系统结构图分析：K开环传递函数为：Wk(s)=一K一s(Ts+1)闭环传递函数为：Wb(s)=亚(k)= - = - 1+W(k)s(Ts+1)+-Ts2+s+K设特征方程为：f(s)=s2+2:ws+w2(1)已知：K=10,T=1秒，要求加入状态反馈后系统的动态性能指标为：o%W20%,tsW1秒。(2)已知：K=1,T=0.05秒，要求加入状态反馈后系统的动态性能指标为：o%W5%,tsW0.5秒。分别观测状态反馈前后两个系统的阶跃响应曲线，并检验系统的动态性能指标是否满足设计要求。10 1(1)Wb(s)=——-一 经计算可取自=▼二0.707,w=10s2+s+10 七：2 n第26页共37页信控学院上机实验则加入状态反馈后的特征方程为:f(s)=s+14.14s+10其特征根为：-7.07sqrt(-50)则状态反馈矩阵:Ky=[k1k2]=[10-1001-14.14]=[-90-13.14]程序如下：num=10;den=[1110];sys=tf(num,den);figure(1);step(sys);gridon;title('原系统的阶跃响应曲线')；denf=[114.14100];k1=den(:,3)-denf(:,3);k2=den(:,2)-denf(:,2);disp('系统的状态反馈增益矩阵K:')Ky=[k1k2][A,B,C,D]=tf2ss(num,den);disp('原系统的极点为:‘)；p=eig(A)'P=[-7.07-sqrt(-50);-7.07+sqrt(-50)];K=place(A,B,P)disp('配置后系统的极点为:‘)p=eig(A-B*K)'disp('配置后的闭环系统为:‘)sys=ss(A-B*K,B,C,D)figure(2);step(sys/dcgain(sys))gridon;title('加入反馈后系统的阶跃响应曲线')；[y,t]=step(sys);C=dcgain(sys);[Y,K]=max(y);Tp=t(K)percentover=(Y-C)/Ci=length(t);while(y(i)>0.98*C)&(y(i)<1.02*C)i=i-1;第27页共37页信控学院上机实验endTs=t(i)运行结果：系统的状态反馈增益矩阵Ky:Ky=-90.0000-13.1400原系统的极点为：p=-0.5000-3.1225i-0.5000+3.1225iK=13.1400 89.9849配置后系统的极点为：p=-7.0700-7.0711i-7.0700+7.0711i配置后的闭环系统为：a=x1x2x1-14.14-99.98x21 0b=u1x11x20c=x1x2y1010d=u1y10Continuous-timemodel.Tp=0.4452percentover=0.0432Ts=0.5936第28页共37页信控学院上机实验原系统的阶跃响应曲编./—-///加入反慑后系统的阶跃啕应墀*原系统的阶跃响应曲编./—-///加入反慑后系统的阶跃啕应墀*0.1 32 0.3 04 0.5 06ft.7D8。由1Time(sec)可见加入状态反馈后，系统超调量减小，且系统平滑过渡到稳态，从而系统性能变好。1 20(2)Wb(s)= = 0 0.05s2+s+1s2+20s+20第29页共37页信控学院上机实验经计算可取自1经计算可取自1===0.707v2w=10则加入状态反馈后的特征方程为：f(s)=s2+14.14S+10其特征根为：-7.07+sqrt(-50)则状态反馈矩阵为：K=[k1k2]=[20-10020-14.14]=[-805.86]程序如下：num=20;den=[12020];sys=tf(num,den);figure(1);step(sys);gridon;title('原系统的阶跃响应曲线')；denf=[114.14100];k1=den(:,3)-denf(:,3);k2=den(:,2)-denf(:,2);disp('系统的状态反馈增益矩阵Ky:')Ky=[k1k2][A,B,C,D]=tf2ss(num,den);disp('原系统的极点为:‘);p=eig(A)'P=[-7.07-sqrt(-50);-7.07+sqrt(-50)];K=place(A,B,P)disp('配置后系统的极点为:‘)p=eig(A-B*K)'disp('配置后的闭环系统为:‘)sys=ss(A-B*K,B,C,D)figure(2);step(sys/dcgain(sys))gridon;title('加入反馈后系统的阶跃响应曲线')；[y,t]=step(sys);C=dcgain(sys);[Y,K]=max(y);Tp=t(K)第30页共37页信控学院上机实验percentover=(Y-C)/Ci=length(t);while(y(i)>0.98*C)&(y(i)<1.02*C)i=i-1;endTs=t(i)运行结果：系统的状态反馈增益矩阵Ky:Ky=-80.0000 5.8600原系统的极点为：p=-18.9443 -1.0557K=-5.8600 79.9849配置后系统的极点为：p=-7.0700-7.0711i-7.0700+7.0711i配置后的闭环系统为：a=x1x2x1-14.14-99.98x21 0b=u1x11x20c=x1x2y1020d=u1y10Continuous-timemode