| 查看: 647 | 回复: 2 | |||
[求助]
谁帮调整下线性时不变系统的输出曲线，只是调整极点位置，救急啊。。。
|
|
要交个小作业，很简单的，可是我调不好，大家帮帮忙啊，周日晚上前尽量搞出来，周一交，说白了就是怎么调出合适的极点阵P，简单吧，不过找到估计也不容易。所以任务就是给出一个合适的极点阵P，这样我自己就能得到输出曲线了。 线性系统dx/dt=Ax+Bu, y=Cx, A,B,C阵如下所示， A=[0 0 1 0 0;0 -0.154 -0.0042 1.54 0;0 0.249 -1 -5.2 0;0.0386 -0.996 -0.0003 -0.117 0;0 0.5 0 0 -0.5]; B=[0 0;-0.744 -0.032;0.337 -1.12;0.02 0;0 0];C=[0 1 0 0 -1;0 0 1 0 0;0 0 0 1 0;1 0 0 0 0];其中A,B可控，A,C可观，现在让镇定系统，选择状态反馈u=-Kx，从而系统化为dx/dt=(A-BK)x,可控，所以可以配置极点。我选择的极点是P=[-3,-0.5,-3,-9,-7]，输出y1=x2-x5,y2=x3,y3=x4,y4=x1，仿真步长0.001s,仿真时间5s.可是得到的输出曲线有2个小问题，第一个是y1=x2-x5收敛太慢了（其中x1收敛挺快的，就是由于x5收敛慢导致y1收敛慢），第二个是y2=x3仿真一开始后的值跑的稍微有点远到-0.8左右了，希望跟其他状态一样，在[-0.5,0.5]范围内波动。我自己调半天极点P也调整不好，感觉第2个极点-0.5基本不能动了，第4个极点跟第5个极点单独变化都输出影响不大，基本就是调整第1个跟第3个，配合第4个跟第5个来得到比较好的收敛曲线：一方面y1,y2,y3,y4收敛的快点，2s左右收敛到0都能接受；二方面不要开始震的那么厉害。如果震荡实在解决不了，只要没我下面的图震的厉害也行，但是收敛速度得解决啊。我用的P阵调整的图除了y1,其他3个收敛时间还凑合。我简单附下程序，找人帮写的，大家凑合看吧。 K=place(A,B,P);h=0.001; t_end=5; t_inter=h; t=0:t_inter:t_end; t=t'; nps=length(t); x(1,  =[0.4328 0.0657 0.5233 0.0698 0.3489];for i=2:nps x(i, =x(i-1,+h*((A-B*K)*x(i-1,')';end figure(1); plot(t,x(:,2)-x(:,5),'y-',t,x(:,3),'g-',t,x(:,4),'r-',t,x(:,1),'k-'); 最后附上我仿的图，看图大家就知道我说的问题了，就是下图中的黄线收敛太慢，绿线开始想下窜的太狠。  m1.jpg “ ”是:加右)，居然搞成微笑了，伤不起啊。[ Last edited by everfx on 2013-11-9 at 16:47 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有34人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2013-11-09 18:04:35
3楼2013-11-09 19:35:49
50