| 查看: 1864 | 回复: 3 | ||
zjc1987金虫 (著名写手)
|
[求助]
LMI mincx 和gevp 问题 我得到的矩阵的数量级为啥那么大？
|
|
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=dec2mat(lmisys,popt,F) P_inv=eye(5)/P; eig(P)%P为正定对称阵 L=-P_inv*K; 以上是我的代码，下面是结果 Solver for generalized eigenvalue minimization Iterations : Best objective value so far * switching to QR 1 309.375000 2 146.808105 3 101.297593 4 3.147201 5 2.171569 6 0.725681 7 0.725681 8 0.500720 9 0.112361 10 0.080799 11 0.080799 12 0.055751 13 0.017428 14 0.012025 15 8.297436e-003 16 7.467693e-003 17 6.720923e-003 18 6.048831e-003 19 5.443948e-003 20 4.899553e-003 21 4.409598e-003 22 3.968638e-003 23 3.571774e-003 24 3.214597e-003 25 2.893137e-003 26 2.603823e-003 *** new lower bound: -0.006076 27 2.441084e-005 *** new lower bound: -0.005879 28 2.196976e-005 29 1.977278e-005 *** new lower bound: -0.001453 30 1.779551e-005 31 1.601596e-005 *** new lower bound: -7.414335e-005 32 5.198251e-007 33 5.198251e-007 34 5.198251e-007 *** new lower bound: -8.813072e-006 35 8.234556e-008 *** new lower bound: -4.010205e-006 36 1.839946e-008 37 1.839946e-008 *** new lower bound: -1.995903e-006 Result: feasible solution best value of t: 1.839946e-008 f-radius saturation: 91.573% of R = 1.00e+008 Termination due to SLOW PROGRESS: the gen. eigenvalue t decreased by less than 1.000% during the last 5 iterations. P = 1.0e+007 * 2.7203 -0.6492 0.0771 0.1632 0.0438 -0.6492 3.0753 -0.5240 0.1544 0.3525 0.0771 -0.5240 1.9095 -0.1469 -0.0000 0.1632 0.1544 -0.1469 2.4675 -0.0000 0.0438 0.3525 -0.0000 -0.0000 0.5666 K = 1.0e+007 * -2.9074 -0.4745 0.1665 0.1722 -0.4745 -2.3822 -0.1845 0.5828 0.1665 -0.1845 -3.3104 -0.8516 0.1722 0.5828 -0.8516 -3.1560 0.2846 1.2298 2.2264 1.2542 F = 1.0e+007 * 0.1632 0.1544 -0.1469 2.4675 -0.0771 0.5240 -1.9095 0.1469 ans = 1.0e+007 * 0.5048 1.6938 2.1521 2.6687 3.7197 >> 我的问题是得到的矩阵P F为什么那么大的数量级？别人做出来的结果很小。是不是我的程序有问题？如果用mincx该如何写？谢谢！ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有4人回复
-
让我中一个面上吧！
已经有17人回复
-
基金不中，共勉
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LMI问题
已经有0人回复
- LMI问题求助
已经有1人回复
- 关于matlab的LMI求助
已经有8人回复
- LMI TOOLBOX中用defcx遇到参数太多的问题
已经有1人回复
- 【讨论】LMI到底好用不好用？这是什么问题，不可思议！
已经有0人回复
- 【求助】一个lmi问题。。求高手解答。。
已经有3人回复
2楼2014-07-26 16:19:06
3楼2015-09-11 10:28:33
4楼2016-12-23 15:28:26
投票:
