lyapunov_eps.m

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

実行結果

>> lyapunov_eps

SeDuMi 1.3 by AdvOL, 2005-2008 and Jos F. Sturm, 1998-2003.
Alg = 2: xz-corrector, theta = 0.250, beta = 0.500
Put 1 free variables in a quadratic cone
eqs m = 3, order n = 7, dim = 11, blocks = 4
nnz(A) = 11 + 0, nnz(ADA) = 9, nnz(L) = 6
 it :     b*y       gap    delta  rate   t/tP*  t/tD*   feas cg cg  prec
  0 :            2.53E+001 0.000
  1 :  0.00E+000 7.07E+000 0.000 0.2793 0.9000 0.9000   1.44  1  1  4.2E+000
  2 :  0.00E+000 1.37E+000 0.000 0.1938 0.9000 0.9000   1.90  1  1  9.9E-001
  3 :  0.00E+000 4.93E-003 0.000 0.0036 0.9990 0.9990   1.10  1  1  4.0E-002
  4 :  0.00E+000 5.77E-010 0.000 0.0000 1.0000 1.0000   1.00  1  1  5.7E-009
  5 :  0.00E+000 2.61E-014 0.000 0.0000 1.0000 1.0000   1.00  1  4  2.0E-013

iter seconds digits       c*x               b*y
  5      0.2  14.7  1.5228194498e-015  0.0000000000e+000
|Ax-b| =  2.2e-015, [Ay-c]_+ =  1.5E-016, |x|= 7.1e-001, |y|= 8.0e-001

Detailed timing (sec)
   Pre          IPM          Post
9.998E-003    1.440E-001    2.997E-003    
Max-norms: ||b||=0, ||c|| = 1,
Cholesky |add|=0, |skip| = 1, ||L.L|| = 1.
sol = 
    yalmiptime: 0.1530
    solvertime: 0.1580
          info: 'No problems detected (SeDuMi-1.3)'
       problem: 0
        dimacs: [2.2275e-015 0 0 1.1102e-016 1.5228e-015 2.0216e-015]
P_feas =
    0.7403    0.1461
    0.1461    0.2597
ans =
     1     1
ans =
     1     1
    
前のページ (kawata_ebihara.zip) に戻る