线性系统理论大作业.docx
- 文档编号:25671712
- 上传时间:2023-06-11
- 格式:DOCX
- 页数:18
- 大小:57.48KB
线性系统理论大作业.docx
《线性系统理论大作业.docx》由会员分享,可在线阅读,更多相关《线性系统理论大作业.docx(18页珍藏版)》请在冰豆网上搜索。
线性系统理论大作业
题目一2
(1)状态反馈加积分器校正的输出反馈系统设计2
(1)建立被控对象的状态空间模型,并判断系统性质2
(2)状态反馈增益矩阵和积分增益常数的设计4
(3)全维观测器设计6
(4)如何在闭环调速系统中增加限流环节7
(2)二次型最优全状态反馈控制和按负载扰动前馈补偿的复合控制系统设计8
(1)线性二次型最优全状态反馈设计8
(2)降维观测器设计15
题目二17
(1)判断系统是否存在最优控制律17
(2)非零给定点的最优控制设计和仿真分析17
(3)权矩阵的各权值对动态性能影响分析19
题目一
(一)状态反馈加积分器校正的输出反馈系统设计
(1)建立被控对象的状态空间模型,并判断系统性质
1)画出与题目对应的模拟结构图,如图1所示:
图1原始系统结构图
取状态变量为x1=n,x2=ld,x3=ud,控制输入u=uc
将已知参数代人并设输出y=n=Xt,得被控对象的状态空间表达式为
0
0
375
GD2
-30.488
B
0=
=0,E
0:
=0,C100
Ks
23529.41
0
0
Ts
2)检查被控系统的结构性质
判断系统能控性、能观性、稳定性
程序如下:
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];C=[100];
Qc=ctrb(A,B);
Qo=obsv(A,C);
L=length(A);
ifrank(Qc)==L
disp('系统是状态完全能控');
else
disp('系统是状态不完全能控');
end
ifrank(Qo)==L
disp('系统是状态完全能观');
else
disp('系统是状态不完全能观');
end
disp(eig(A))%利用A的特征值判断系统稳定性
运行结果:
系统是状态完全能控
系统是状态完全能观
1.0e+02*
-0.0893+0.0820i
-0.0893-0.0820i
-5.8823+0.0000i
由于矩阵A全部特征值均具有负实部,因此系统渐近稳定。
原系统
设负载转矩为0,输入为阶跃信号,系统simulink仿真如下:
图2原始开环系统结构框图
图3原始开环系统仿真
分析:
由系统仿真图可以看出,调节时间大于0.5s,不满足性能指标。
(2)状态反馈增益矩阵和积分增益常数的设计
由于原系统能控,可以使用状态反馈。
为满足设计指标,采用状态反馈加积分器校正的输出反馈系统。
因增广系统能控,故可采用线性状态反馈控制律uK1xK2w。
将闭环系统极点配置到复平面左半开平面的任意期望位置且可消除阶跃扰动及阶跃参
考输入作用下的稳态误差
wvy,v为系统参考输入。
式中:
Ki=K11K12K13,x=XiX2X3
由经典控制理论,闭环极点为1,2njn的欠阻尼二阶线性定常系统的超
调量及调节时间为%=e1$100%,ts仝5。
n
系统需满足%<8%,tsW0.5s,计算可得>0.6266,n>11.1714,取=0.7,n=12,设计指标的期望闭环主导极点对为;28.4j8.57。
选择2个期望的闭环非主导极点离虚轴为主导极点5倍以上,取3,450,据期望闭环极点,采用MATLA极点配置函数可求出增广系统状态反馈增益阵,程序如下:
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];C=[100];
Az=[A[0;0;0];-C0];Bz=[B;0];Cz=[C0];
P=[-8.4+j*8.57;-8.4-j*8.57;-50;-50];
Km=acker(Az,Bz,P);
K=[Km(1,1),Km(1,2),Km(1,3),-Km(1,4)]
运行程序可得:
系统simulink仿真如下:
图4状态反馈加积分器校正系统结构框图
图5状态反馈加积分器校正系统仿真
由图可知,超调量%(1.0420-1)/14.2%8%,调节时间ts0.5s,满足要求
图6加负载扰动后系统状态反馈加积分器校正系统仿真
0时刻扰动,最终系统稳定在1,因此系统稳态误差为0。
(3)全维观测器设计
由于系统能观,可以使用状态观测器
00
B0=0,C100,
Ks23529.41
TT
新系统的特征根为:
-61,27.8072j12.1936,基于通常选择观测器的响应速度比
所考虑的状态反馈闭环系统快2-5倍这一经验规则,取观测器期望极点为:
-150,-60,-
70。
应用MATLA极点配置函数求解新系统全维观测器,程序如下:
A=[039.7680;-3.696-17.85727.056;-38.8235-88.9412-98.8233];
B=[0;0;23529.41];C=[100];
P=[-150;-60;-70];
Gt=acker(A',C',P);%求对偶系统的状态反馈增益阵G
G=Gt';%求系统的观测器偏差反馈增益矩阵G
运行程序可得:
163.3197
G=8.1911。
129.8454
带观测器的状态反馈加积分调节系统仿真结构图如下:
图7带观测器的状态反馈加积分器校正系统结构框图
图8系统加全维观测器波形图
图9全维观测器波形图
由仿真图可知,系统的稳态误差为0,动态误差满足超调量%<8%,调节时间ts<0.5s的要求。
?
状态估计误差收敛速度与状态观测器极点的配置有关。
一般而言状态观测器极点在复平面的左半开平面距离虚轴距离越远,则估计误差收敛速度越快。
但是,观测器响应速度过快会产生大量噪声,影响系统的正常工作故不宜取值过大。
综合工程实际出发,一般取为比状态反馈闭环系统快2—5倍。
(4)如何在闭环调速系统中增加限流环节
从加快启动电动机的角度来看,闭环调速系统应允许有较大的启动电流,而造成堵转的故障消失后,系统电流应能自动恢复正常。
所以常规的熔断器或过流继电器在这里均不能作为限流保护措施。
因为它们是通过切断电路来保护设备的,虽然能起到保护作用,但
故障消失后,系统无法自动恢复正常。
为了充分利用设备的过流能力,又保证设备的安全运行,电流截止负反馈则可以限制电流的大小。
?
电流截止负反馈的作用是:
当电枢电流大于某一截止值时,电流负反馈起作用,限制电流不能过大。
当电枢电流小于截止值时,电流反馈被截止,对系统的稳态运行不产生影响。
电动机启动时,因为电流截止负反馈作用,从而限制启动电流。
正常工作时,电流截止负反馈作用很小。
电动机发生堵转时,由于电流截止负反馈的作用,使Ud大大下降,因而
使la不致过大。
允许的堵转电流一般为电动机额定电流的2~2.5倍。
系统工作在额定值
时,由于电流截止负反馈起作用,从而保证系统设备的安全。
?
电流截止负反馈如图所示:
图10电流截止负反馈结构图
(二)二次型最优全状态反馈控制和按负载扰动前馈补偿的复合控制系统设计
(1)线性二次型最优全状态反馈设计
1)判断系统是否存在最优控制律
要使系统阶跃响应具有良好的动、静态特性,可按非零给定点的最优控制律设计,即
u(t)KXt)Ug,由于输入维数和输出维数相等,所以u(t)KXt)W1(O)yw。
由于系统完全能控,因此,最优控制u(t)存在。
1
最优控制性能指标为:
J-txT(t)QXt)uT(t)RUt)dt,其中Q为状态加权系数
2to
矩阵,R为控制加权系数矩阵。
2)非零给定点的最优控制设计和仿真分析
q11
0
0
100
0
0
Q
0
q22
0
0
1
0,R=1
0
0
q33
0
0
1
程序如下:
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];
C=[100];
D=0;
R=1;
Q=[10000;010;001];
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv((-C/(A-B*K))*B);
bc=B*W;
cc=C;
dc=D;
step(ac,bc,cc,dc);
grid
运行结果如下:
K=9.86354.87020.9809,Wc110.0009。
图11非零给定点最优控制系统单位阶跃响应
3)权矩阵的各权值对动态性能影响分析
a)固定控制加权系数矩阵R=1,且另q22、q33都为1,%取不同值时,研究非零给定点的最
优控制仿真曲线。
程序如下:
a_color=['r','g','b','y','c','m','k'];
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];
C=[100];
D=0;
R=1;
symsQq11;
N=[1100,200,500,1000,10000];
symsiK;
fori=1:
6
q11=N(i);
Q=[q1100;010;001]
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv(-C/(A-B*K)*B);
bc=B*W;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3),sys(4),a_color(
4),sys(5),a_color(5),sys(6),a_color(6));
grid
结果曲线如下:
图12q11取不同值时二次型最优全状态反馈曲线
1000,
红色:
qn1,绿色:
也100,蓝色:
qn200,黄色:
500,蓝绿色:
紫红q1110000。
由图可知,随着qn的增大,调节时间减少;如果qn过大,超调量会增大。
b)固定控制加权系数矩阵R=1,qn、q?
2、q33分别取相同值时,研究非零给定点的最优控
制仿真曲线。
程序如下:
a_color=['r','g','b'];
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];
C=[100];
D=0;
R=1;
symsQq11q22q33;
fori=1:
3
q11=1;
q22=1;
q33=1;
Q=[q1100;0q220;00q33];
Q(i,i)=100;
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv(-C/(A-B*K)*B);
bc=B*W;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3));
grid
结果曲线如下:
图13q11、q22、q33分别取相同值时二次型最优全状态反馈曲线
红色:
q11100,绿色:
q22100,蓝色:
q33100。
由图可知,q11100时,调节时间最小,系统动态性能比另外两个好。
c)当状态加权系数矩阵Q不变,控制加权系数矩阵R取不同值时。
研究非零给定点的最优
控制仿真曲线。
程序如下:
a_color=['r','g','b','y','c'];
A=[039.7680;-3.696-17.85727.056;00-588.235];
B=[0;0;23529.41];
C=[100];
Q=[10000;010;001];
D=0;
symsR;
N=[1,100,1000,1500,10000];
fori=1:
5
R=N(i);
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv(-C/(A-B*K)*B);
bc=B*W;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3),sys(4),a_color(
4),sys(5),a_color⑸);
grid
结果曲线如下:
图14R取不同值时二次型最优全状态反馈曲线
由图可知,R越大时,调节时间也越大
(2)降维观测器设计
由于Tl平稳变化,n、Id、Ud均可准确测量,对负载转矩进行估计的降维观测器的设
计需要对原系统结构进行变化参考文献[3],系统增广矩阵的状态方程可写成:
010
0
n
y001
0
。
1
000
1
1d
Ud
参考文献[1],将直接可测的与不能直接可测的状态变量分开:
X
X2
A11A12X1
j“
A21A22X2
B1
u,
B2
X1
y0
1
X20
X2
-30.488
式中,
A11:
=0,^2=0
00,A21=
0
0
0
39.768
0
A22'
-3.696-17.857
27.056,
0
0
-588.235
0
n
B1=0,
目2:
=0
X1Tl,X2
1d0
23529.41
ud
需要设计一维观测器重构。
设降维观测器的反馈阵Ggig2©3。
则降维观测器特征多项式为:
f()det[I(AnGA21)]30.488©!
。
选择降维观测器期望极点为闭环极点的2-5倍,取=-150,贝冋得©1=4920。
引入降
维观测器状态方程:
例调节的前馈补偿:
„RK3KsK2K.10.9809404.8702401门…厂
Kp3s2sWc0.425。
pKsCt401.304410.0009
降维观测器仿真如下:
图15降维观测器仿真结构图
图16降维观测器无扰动时仿真波形
图17降维观测器1s扰动时仿真波形
题目二
(1)判断系统是否存在最优控制律
要使系统阶跃响应具有良好的动、静态特性,可按非零给定点的最优控制律设计,即
u(t)KXt)Ug,由于输入维数和输出维数相等,所以u(t)KXt)W1(0)yw。
由于系统为能控标准型,所以系统完全能控,因此,最优控制u(t)存在。
1
最优控制性能指标为:
J-txT(t)QXt)uT(t)RUt)dt,其中Q为状态加权系数2t0
矩阵,R为控制加权系数矩阵。
(2)非零给定点的最优控制设计和仿真分析
原系统动态仿真模型如下:
图18原系统结构框图
图19原系统单位阶跃响应
由图可知,系统输出响应发散,可引入最优控制u(t)Kx(t)Wc1(0)yw。
选取设
q110010000
Q0q220=010,R=1。
00q33001
程序如下:
A=[010;001;0-18-8];
B=[0;0;1];
C=[100];
R=1;
Q=[10000;010;001];
K=lqr(A,B,Q,R);
W=inv((-C/(A-B*K))*B);
运行结果如下:
K=104.28320.5771,Wc110。
非零给定点最有控制系统动态仿真模型如下:
图20非零给定点最有控制系统结构框图
图21非零给定点最优控制系统单位阶跃响应
(3)权矩阵的各权值对动态性能影响分析
1)固定控制加权系数矩阵R=1,且另q22、q33都为1,%取不同值时,研究非零给定点的最优控制仿真曲线。
程序如下:
a_color=['r','g','b','y','c','k'];
A=[010;001;0-18-8];
B=[0;0;1];
C=[100];
D=0;
R=1;
symsQq11;
N=[100,200,500,1000,10000];
symsiK;
fori=1:
5
q11=N(i);
Q=[q1100;010;001]
K=lqr(A,B,Q,R);
ac=A-B*K;
k1=inv(-C/(A-B*K)*B);
bc=B*k1;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3),sys(4),a_color(
4),sys(5),a_color(5));
grid
结果曲线如下:
图22q11取不同值时二次型最优全状态反馈曲线
红色:
q11100,绿色:
q11200,蓝色:
%500,黄色:
q111000,蓝绿色:
也10000。
由图可知,随着q11的增大,调节时间减少;如果q11过大,超调量会增大。
2)固定控制加权系数矩阵R=1,%、q22、433分别取相同值时,研究非零给定点的最优控
制仿真曲线。
程序如下:
a_color=['r','g','b'];
A=[010;001;0-18-8];
B=[0;0;1];
C=[100];
D=0;
R=1;
symsQq11q22q33;
fori=1:
3
q11=1;
q22=1;
q33=1;
Q=[q1100;0q220;00q33];
Q(i,i)=100;
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv(-C/(A-B*K)*B);
bc=B*W;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3));
grid
结果曲线如下:
图23q11、q22、q33分别取相同值时二次型最优全状态反馈曲线
红色:
q11100,绿色:
q22100,蓝色:
q33100。
由图可知,q11100时,系统动态性能比另外两个好。
3)当状态加权系数矩阵Q不变,控制加权系数矩阵R取不同值时。
研究非零给定点的最优
控制仿真曲线。
程序如下:
a_color=['r','g','b','y','c'];
A=[010;001;0-18-8];
B=[0;0;1];
C=[100];
Q=[10000;010;001];
D=0;
symsR;
N=[1,100,1000,1500,10000];
fori=1:
5
R=N(i);
K=lqr(A,B,Q,R);
ac=A-B*K;
W=inv(-C/(A-B*K)*B);
bc=B*W;
cc=C;
dc=D;
sys(i)=ss(ac,bc,cc,dc);
end
figure
(1)
step(sys
(1),a_color
(1),sys
(2),a_color
(2),sys(3),a_color(3),sys(4),a_color(
4),sys(5),a_color(5));
grid
结果曲线如下:
图24R取不同值时二次型最优全状态反馈曲线
由图可知,R越大时,调节时间也越大。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性 系统 理论 作业