| 查看: 1851 | 回复: 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该如何写？谢谢！ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LMI问题
已经有0人回复
- LMI问题求助
已经有1人回复
- 关于matlab的LMI求助
已经有8人回复
- LMI TOOLBOX中用defcx遇到参数太多的问题
已经有1人回复
- 【讨论】LMI到底好用不好用？这是什么问题，不可思议！
已经有0人回复
- 【求助】一个lmi问题。。求高手解答。。
已经有3人回复
3楼2015-09-11 10:28:33
2楼2014-07-26 16:19:06
4楼2016-12-23 15:28:26
25