声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1905|回复: 6

[控制理论] 求助,关于Hinf控制器

[复制链接]
发表于 2009-4-15 09:52 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
num=400;
den=[1 2 400]
g=tf(num,den);
nw1=100*conv([0.005,1],[0.005,1]);
dw1=conv([0.2,1],[0.2,1]);
nw3=[1,0,0];
dw3=[40000];

[a,b,c,d]=tf2ss(num,den);
s=mksys(a,b,c,d);
w1=[nw1;dw1];
w3=[nw3;0,0,dw3];
tss=augtf(s,w1,[],w3);
[cf1,ccl]=hinf(tss);
[af,bf,cf,df]=branch(cf1);
gc=zpk(ss(af,bf,cf,df))

刚涉入鲁棒控制,有基本教材上都有这个例子,可是却实现不了,显示[cf1,ccl]=hinf(tss);这里有错误,请高手指点,谢谢。
回复
分享到:

使用道具 举报

发表于 2009-4-15 15:39 | 显示全部楼层
可以运行,运行结果如下:
<< H-inf Optimal Control Synthesis >>

            Computing the 4-block H-inf optimal controller
          using the S-L-C loop-shifting/descriptor formulae
  
     Solving for the H-inf controller F(s) using U(s) = 0 (default)
     Solving Riccati equations and performing H-infinity
     existence tests:
        1.  Is D11 small enough?                       OK
        2.  Solving state-feedback (P) Riccati ...
             a.  No Hamiltonian jw-axis roots?         OK
             b.  A-B2*F stable (P >= 0)?               OK
        3.  Solving output-injection (S) Riccati ...
             a.  No Hamiltonian jw-axis roots?         OK
             b.  A-G*C2 stable (S >= 0)?               OK
        4.  max eig(P*S) < 1 ?                         OK
    -------------------------------------------------------
       all tests passed -- computing H-inf controller ...
                         DONE!!!
    -------------------------------------------------------

Zero/pole/gain:
8597.8554 (s+42.84) (s^2  + 2s + 400)
-------------------------------------
(s+5)^2 (s^2  + 308.5s + 4.388e004)
 楼主| 发表于 2009-4-15 16:22 | 显示全部楼层
哦,谢谢楼上,你倒是提醒了我,也许是我的matlab没安装完全,我只装了一个盘,怪不得help  hinf 都找不到。下载个再安装试试
 楼主| 发表于 2009-4-16 08:45 | 显示全部楼层
请问2楼,为什么我完全安装了matlab后还是不能运行,错误如下:
??? Attempt to execute SCRIPT hinf as a function.

Error in ==> hinf at 18
[cf1,ccl]=hinf(tss);
 楼主| 发表于 2009-4-16 08:46 | 显示全部楼层
我用的是7.1版本,按理来说这个不是问题
 楼主| 发表于 2009-4-16 09:27 | 显示全部楼层
不好意思,原来是文件名不对,改个文件名就好了。
发表于 2009-4-16 09:52 | 显示全部楼层
有时候一个不注意的小问题就能让人憋半天的。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-17 15:55 , Processed in 0.058792 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表