lyapunov2.m

　システム
　　　
の安定解析問題に関する連立 LMI
　　　
を解く．
　なお，このシステムは不安定である．

実行結果 

>> lyapunov2

SeDuMi 1.3 by AdvOL, 2005-2008 and Jos F. Sturm, 1998-2003.
Alg = 2: xz-corrector, theta = 0.250, beta = 0.500
eqs m = 3, order n = 5, dim = 9, blocks = 3
nnz(A) = 9 + 0, nnz(ADA) = 9, nnz(L) = 6
 it :     b*y       gap    delta  rate   t/tP*  t/tD*   feas cg cg  prec
  0 :            2.89E+001 0.000
  1 :  0.00E+000 5.33E+000 0.000 0.1845 0.9000 0.9000  -0.90  1  1  1.4E+001
  2 :  0.00E+000 3.25E-001 0.000 0.0610 0.9900 0.9900  -0.55  1  1  4.1E+000
  3 :  0.00E+000 7.98E-004 0.000 0.0025 0.9990 0.9990  -0.99  1  1  5.5E+000
  4 :  0.00E+000 8.33E-005 0.000 0.1044 0.9450 0.9450  -0.91  1  1  4.7E+000
  5 :  0.00E+000 4.55E-006 0.030 0.0546 0.9900 0.9900  -0.90  2  2  4.4E+000
  6 :  0.00E+000 3.59E-007 0.000 0.0788 0.9900 0.9900  -0.93  3  3  3.5E+001
  7 :  0.00E+000 2.24E-008 0.313 0.0624 0.9900 0.9900  -1.16  4  4  5.3E+001
  8 :  0.00E+000 8.49E-009 0.000 0.3790 0.9000 0.9000  -1.22  5  5  3.4E+001
  9 :  0.00E+000 2.40E-009 0.000 0.2831 0.9000 0.9000  -0.71  5  5  2.2E+001
 10 :  0.00E+000 6.23E-010 0.000 0.2593 0.9000 0.9000  -0.90  5  5  1.9E+001

Dual infeasible, primal improving direction found.
iter seconds  |Ax|    [Ay]_+     |x|       |y|
 10      0.3  1.9e-010  3.9e-011  1.3e+000  5.4e-001

Detailed timing (sec)
   Pre          IPM          Post
1.099E-002    2.220E-001    1.099E-002    
Max-norms: ||b||=0, ||c|| = 1,
Cholesky |add|=0, |skip| = 1, ||L.L|| = 1.12307.
sol = 
    yalmiptime: 0.1520
    solvertime: 0.2450
          info: 'Infeasible problem (SeDuMi-1.3)'
       problem: 1
        dimacs: [1.8666e-010 0 0 0.5000 -0.5000 -0.5000]
前のページ (kawata_ebihara.zip) に戻る