zhedongxiao 发表于 2011-5-22 13:46

求教 ：用LMI方法求倒立摆H无穷状态反馈控制器程序

本帖最后由 zhedongxiao 于 2011-5-22 13:48 编辑

求教 ：用LMI方法求倒立摆H无穷状态反馈控制器程序
求教 ：用LMI方法求倒立摆H无穷状态反馈控制器程序用LMI方法求倒立摆H无穷状态反馈控制器程序 但运行出现没有合适的gama值
显示“could not establish feasibility nor infeasibility”

A=;
B1=;
B2=;
C1=;
D=;
E=;
F=;
e=4;
setlmis([]);               %%初始化LMI系统

X=lmivar(1,);
W=lmivar(2,);
gama=lmivar(2,);   %%定义决策变量X,W,r

NO=newlmi;
lmiterm(,A,1);
lmiterm(,B2,1);
lmiterm(,1,A');
lmiterm(,1,B2');
lmiterm(,e*E*E');
lmiterm(,1,F');
lmiterm(,B1);
lmiterm(,1,C1');
lmiterm(,1,D');
lmiterm(,-e);
lmiterm(, -1);   
lmiterm(,-gama^2*eye(size(C1,1)));      %%对矩阵不等式进行LMI描述
    %27行可能出错处
N0=newlmi
lmiterm([-NO 1 1 X],1,1);                  %%对X>0进行描述
               
lmil=getlmis;                                                                         %%获取LMI系统描述

n=decnbr(lmil);                        
c=zeros(n,1);
for j=1:n
=defcx(lmil,j,gama);
c(j)=gamaj*gamaj;
end                                  %%生成最小化r所需的向量c,这里 x即为决策
   %34-39可能出错处                            %% X中r所对应的变量

=mincx(lmil,c);         %%由求解器mincx()求解最小化                     
X=dec2mat(lmil,xopt,X);
W=dec2mat(lmil,xopt,W);            %%从决策变量xopt中得到变量矩阵X,W的解
W*inv(X)                           %%求得K
求教 ：用LMI方法求倒立摆H无穷状态反馈控制器程序用LMI方法求倒立摆H无穷状态反馈控制器程序 但运行出现没有合适的gama值
显示“could not establish feasibility nor infeasibility”

A=;
B1=;
B2=;
C1=;
D=;
E=;
F=;
e=4;
setlmis([]);               %%初始化LMI系统

X=lmivar(1,);
W=lmivar(2,);
gama=lmivar(2,);   %%定义决策变量X,W,r

NO=newlmi;
lmiterm(,A,1);
lmiterm(,B2,1);
lmiterm(,1,A');
lmiterm(,1,B2');
lmiterm(,e*E*E');
lmiterm(,1,F');
lmiterm(,B1);
lmiterm(,1,C1');
lmiterm(,1,D');
lmiterm(,-e);
lmiterm(, -1);   
lmiterm(,-gama^2*eye(size(C1,1)));      %%对矩阵不等式进行LMI描述
    %27行可能出错处
N0=newlmi
lmiterm([-NO 1 1 X],1,1);                  %%对X>0进行描述
               
lmil=getlmis;                                                                         %%获取LMI系统描述

n=decnbr(lmil);                        
c=zeros(n,1);
for j=1:n
=defcx(lmil,j,gama);
c(j)=gamaj*gamaj;
end                                  %%生成最小化r所需的向量c,这里 x即为决策
   %34-39可能出错处                            %% X中r所对应的变量

=mincx(lmil,c);         %%由求解器mincx()求解最小化                     
X=dec2mat(lmil,xopt,X);
W=dec2mat(lmil,xopt,W);            %%从决策变量xopt中得到变量矩阵X,W的解
W*inv(X)                           %%求得K

yufeng 发表于 2011-5-23 18:10

lmiterm用法有问题

zhedongxiao 发表于 2011-5-25 12:29

回复 2 # yufeng 的帖子

什么问题 指教

xiaoga 发表于 2011-6-2 10:54

我也遇到同样的问题，mincx无解， The LMI constraints were found INFEASIBLE 。
??? Error using ==> dec2mat
DECVARS must be a vector of length 12

Error in ==> shiji at 66
Popt=dec2mat(lmisys,xopt,P);
我在网上也看到很多这个问题，但现在还没找到解，我也急求！

zhedongxiao 发表于 2011-6-8 15:25

回复 4 # xiaoga 的帖子

用 feasp 求次优解可以的

xiaoga 发表于 2011-6-17 08:55

回复 5 # zhedongxiao 的帖子

我用feasp求解还是显示 Result:could not establish feasibility nor infeasibility
          f-radius saturation:2.445% of R = 1.00e+009
Termination due to SLOW PROGRESS:
          t was decreased by less than 10.000% during
          the last 10 iterations.


These LMI constraints were found infeasible

??? Error using ==> dec2mat
DECVARS must be a vector of length 12

Error in ==> shiji at 66
Popt=dec2mat(lmisys,xopt,P);
是不是我的已知矩阵有问题，你qq多少可以加我吗？我的694227449
页: [1]
查看完整版本: 求教 ：用LMI方法求倒立摆H无穷状态反馈控制器程序