matlab gevp,LMI mincx 和gevp 问题 我得到的矩阵的数量级为啥那么大?

在使用MATLAB的GEVP求解问题时,作者遇到了矩阵P、F数值过大的问题,与他人结果相差甚远。代码中定义了矩阵并设置了LMI,但经过迭代优化后,得到的矩阵数量级达到10^7。作者疑惑是否程序存在错误,并询问如何使用mincx来解决此问题。
摘要由CSDN通过智能技术生成

C_a=[1 0 0 0 0; 0 1 0 0 0;0 0 1 0 0; 0 0 0 1 0];

D=[0 0; 0 0; 0 -1; 1 0; 0 0];

A=[0 1 0 0 0 ; 0 0 1 0 0; 0 0 0 1 0; 0 0 0 0 1; -1 -5 -10 -10 -5];

setlmis([])

P=lmivar(1,[5 1]);%1代表方块阵

%[P,n,sP]=lmivar(1,[4 1]);

F=lmivar(2,[2 4]);%2代表长方形矩阵

K=lmivar(2,[5 4]);%长方形矩阵

%lmiterm([1 1 1 0],1) ;        % P > I :: I

lmiterm([-1 1 1 P],1,1);         % P > I :: P

lmiterm([2 1 1 P],1,A,'s');%A_T*P+P_T*A

lmiterm([2 1 1 K],1,C_a,'s');%

lmiterm([3 1 2 F],1,C_a);

lmiterm([3 1 2 P],D',-1);%D_T*P-F*C_a

lmiterm([-3 1 1 0],1);

lmiterm([-3 2 2 0],1);

lmisys=getlmis;

[alpha,popt]=gevp(lmisys,1);

%[deta,xopt]=(lmisys,1);[alpha,popt]=feasp(lmisys);

P=dec2mat(lmisys,popt,P)

K=dec2mat(lmisys,popt,K)

F=dec2ma

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值