24小时热门版块排行榜

>论坛更新日志 (2079)
>虫友互识 (121)
>导师招生 (88)
>招聘信息布告栏 (31)
>硕博家园 (22)
>考博 (21)
>公派出国 (21)
>博后之家 (19)
>考研 (14)
>休闲灌水 (13)
>文献求助 (12)
>论文道贺祈福 (11)
>找工作 (11)
>药学 (11)
>教师之家 (10)
>论文投稿 (10)

返回列表

zjc1987

金虫 (著名写手)

应助: 5 (幼儿园)
金币: 2059.2
散金: 1871
红花: 12
帖子: 1661
在线: 649.5小时
虫号: 2374941
注册: 2013-03-22
性别: GG
专业: 常微分方程与动力系统

[求助] 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该如何写？谢谢！

回复此楼