<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2014.517253</article-id><article-id pub-id-type="publisher-id">AM-50346</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  On Two Problems for Matrix Polytopes
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>erife</surname><given-names>Yılmaz</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Taner</surname><given-names>Büyükköroğlu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mathematics, Faculty of Science, Anadolu University, Eskisehir, Turkey</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>serifeyilmaz@anadolu.edu.tr(EY)</email>;<email>tbuyukkoroglu@anadolu.edu.tr(TB)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>09</day><month>10</month><year>2014</year></pub-date><volume>05</volume><issue>17</issue><fpage>2650</fpage><lpage>2656</lpage><history><date date-type="received"><day>21</day>	<month>July</month>	<year>2014</year></date><date date-type="rev-recd"><day>20</day>	<month>August</month>	<year>2014</year>	</date><date date-type="accepted"><day>8</day>	<month>September</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  We consider two problems from stability theory of matrix polytopes: the existence of common quadratic Lyapunov functions and the existence of a stable member. We show the applicability of the gradient algorithm and give a new sufficient condition for the second problem. A number of examples are considered.
 
</p></abstract><kwd-group><kwd>Stable Matrix</kwd><kwd> Matrix Family</kwd><kwd> Common Quadratic Lyapunov Functions</kwd><kwd> Switched System</kwd><kwd>  Gradient Method</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Consider the switched system</p><disp-formula id="scirp.50346-formula201"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x5.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x6.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x7.png" xlink:type="simple"/></inline-formula>. In Equation (1), the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x8.png" xlink:type="simple"/></inline-formula> switches among <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x9.png" xlink:type="simple"/></inline-formula> matrices<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x10.png" xlink:type="simple"/></inline-formula>.</p><p>Switching signal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x11.png" xlink:type="simple"/></inline-formula> is piecewise continuous from the right function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x12.png" xlink:type="simple"/></inline-formula> and the switching times are arbitrary. For the switched system (1) with initial condition <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x13.png" xlink:type="simple"/></inline-formula> and with switching signal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x14.png" xlink:type="simple"/></inline-formula> denotes the solution by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x15.png" xlink:type="simple"/></inline-formula>.</p><p>Definition 1. The origin is uniformly asymptotically stable (UAS) for the system (1) if for every <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x16.png" xlink:type="simple"/></inline-formula> there exists <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x17.png" xlink:type="simple"/></inline-formula> such that for every signal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x18.png" xlink:type="simple"/></inline-formula> and initial state <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x19.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x20.png" xlink:type="simple"/></inline-formula>, the inequality <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x21.png" xlink:type="simple"/></inline-formula> is satisfied for all <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x22.png" xlink:type="simple"/></inline-formula> and uniformly on <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x23.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.50346-formula202"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x24.png"  xlink:type="simple"/></disp-formula><p>If all systems in (1) share a common quadratic Lyapunov function (CQLF) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x25.png" xlink:type="simple"/></inline-formula>then the switched</p><p>system is UAS (T denotes the transpose).</p><p>In this case there exists a common <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x27.png" xlink:type="simple"/></inline-formula> such that</p><disp-formula id="scirp.50346-formula203"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x28.png"  xlink:type="simple"/></disp-formula><p>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x29.png" xlink:type="simple"/></inline-formula> is called a common solution to the set of Lyapunov matrix inequalities (2).</p><p>The problem of existence of common positive definite solution <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x30.png" xlink:type="simple"/></inline-formula> of (2) has been studied in a lot of works (see [<xref ref-type="bibr" rid="scirp.50346-ref1">1</xref>] -[<xref ref-type="bibr" rid="scirp.50346-ref9">9</xref>] and references therein). Numerical solution for common <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x31.png" xlink:type="simple"/></inline-formula> via nondifferentiable convex optimization has been discussed in [<xref ref-type="bibr" rid="scirp.50346-ref10">10</xref>] .</p><p>In the first part of the paper, the problem of existence of CQLF is investigated by Kelley’s method. This method is applied when CQLF problem is treated as a convex optimization problem.</p><p>Second part of the paper is devoted to the following question:</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x32.png" xlink:type="simple"/></inline-formula> be a compact, for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x33.png" xlink:type="simple"/></inline-formula> the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x34.png" xlink:type="simple"/></inline-formula> is a real <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x35.png" xlink:type="simple"/></inline-formula> matrix. Is there a Hurwitz stable member (all eigenvalues lie in the open left half plane) in the family</p><disp-formula id="scirp.50346-formula204"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x36.png"  xlink:type="simple"/></disp-formula><p>or equivalently is there <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x37.png" xlink:type="simple"/></inline-formula> such that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x38.png" xlink:type="simple"/></inline-formula> is stable? This problem is one of the hard and important problems in control theory (see [<xref ref-type="bibr" rid="scirp.50346-ref11">11</xref>] ). Numerical solution of this problem is considered in [<xref ref-type="bibr" rid="scirp.50346-ref12">12</xref>] . In this paper we reduce this problem to a non-convex optimization problem.</p></sec><sec id="s2"><title>2. Common Quadratic Lyapunov Function</title><p>For the switched system</p><disp-formula id="scirp.50346-formula205"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x39.png"  xlink:type="simple"/></disp-formula><p>consider the problem of determination of CQLF <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x40.png" xlink:type="simple"/></inline-formula> where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x41.png" xlink:type="simple"/></inline-formula>. We are going to investigate it by Kelley’s cutting-plane method. This method gives new sufficient condition (Theorem 2) and new algorithm (Algorithm 1) which is more effective in comparison with the algorithm from [<xref ref-type="bibr" rid="scirp.50346-ref10">10</xref>] .</p><p>Consider the problem of existence of a common <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x42.png" xlink:type="simple"/></inline-formula> such that</p><disp-formula id="scirp.50346-formula206"><label>. (3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x43.png"  xlink:type="simple"/></disp-formula><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x44.png" xlink:type="simple"/></inline-formula> be <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x45.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x46.png" xlink:type="simple"/></inline-formula> be an <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x47.png" xlink:type="simple"/></inline-formula> symmetric matrix defined as</p><disp-formula id="scirp.50346-formula207"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x48.png"  xlink:type="simple"/></disp-formula><p>Define</p><disp-formula id="scirp.50346-formula208"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x49.png"  xlink:type="simple"/></disp-formula><p>If there exists <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x50.png" xlink:type="simple"/></inline-formula> such that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x51.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x52.png" xlink:type="simple"/></inline-formula> then the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x53.png" xlink:type="simple"/></inline-formula> is required solution. This problem can be reduced to the minimization of a convex function under convex constraints.</p><p>Consider the following convex minimization problem</p><disp-formula id="scirp.50346-formula209"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x54.png"  xlink:type="simple"/></disp-formula><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula> be a convex set and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x56.png" xlink:type="simple"/></inline-formula> be convex function. The vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x57.png" xlink:type="simple"/></inline-formula> is said to be a subgradient of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x58.png" xlink:type="simple"/></inline-formula> at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x59.png" xlink:type="simple"/></inline-formula> if for all <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x60.png" xlink:type="simple"/></inline-formula></p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x61.png" xlink:type="simple"/></inline-formula>.</p><p>The set of all subgradients of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula> at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x63.png" xlink:type="simple"/></inline-formula> is denoted by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x64.png" xlink:type="simple"/></inline-formula>. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x65.png" xlink:type="simple"/></inline-formula> is an interior point of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x66.png" xlink:type="simple"/></inline-formula> then the set <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x67.png" xlink:type="simple"/></inline-formula> is nonempty and convex. The following proposition follows from nondifferentiable optimization theory.</p><p>Proposition 1. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x68.png" xlink:type="simple"/></inline-formula> be defined as</p><disp-formula id="scirp.50346-formula210"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x69.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x70.png" xlink:type="simple"/></inline-formula> is compact, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x71.png" xlink:type="simple"/></inline-formula>is continuous and differentiable in<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x72.png" xlink:type="simple"/></inline-formula>. Then</p><disp-formula id="scirp.50346-formula211"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x73.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x74.png" xlink:type="simple"/></inline-formula> is the set of all maximizing elements <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x75.png" xlink:type="simple"/></inline-formula> in (6), i.e.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x76.png" xlink:type="simple"/></inline-formula>.</p><p>If for a given <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x77.png" xlink:type="simple"/></inline-formula> the maximizing element is unique, i.e. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x78.png" xlink:type="simple"/></inline-formula>then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x79.png" xlink:type="simple"/></inline-formula> is differentiable at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x77.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x78.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x80.png" xlink:type="simple"/></inline-formula> and its gradient is</p><disp-formula id="scirp.50346-formula212"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x81.png"  xlink:type="simple"/></disp-formula><p>In the case of the Function (4)</p><disp-formula id="scirp.50346-formula213"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x82.png"  xlink:type="simple"/></disp-formula><p>If for the given <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x83.png" xlink:type="simple"/></inline-formula> the maximizing <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x84.png" xlink:type="simple"/></inline-formula> is unique and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x85.png" xlink:type="simple"/></inline-formula> is a simple eigenvalues, the differentiability of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x86.png" xlink:type="simple"/></inline-formula> at the point <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x87.png" xlink:type="simple"/></inline-formula> is guaranteed [<xref ref-type="bibr" rid="scirp.50346-ref13">13</xref>] .</p><p>We investigate problem (5) by Kelley’s cutting-plane method.</p><p>This method converts the problem (5) to the problem</p><disp-formula id="scirp.50346-formula214"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x88.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x89.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x90.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x91.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x91.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x92.png" xlink:type="simple"/></inline-formula>.</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x93.png" xlink:type="simple"/></inline-formula> be a starting point and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x94.png" xlink:type="simple"/></inline-formula> be <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x95.png" xlink:type="simple"/></inline-formula> distinct points.</p><p>At the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x96.png" xlink:type="simple"/></inline-formula>th iteration, the cutting-plane algorithm solves the following LP problem</p><disp-formula id="scirp.50346-formula215"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-7402419x97.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x98.png" xlink:type="simple"/></inline-formula> denotes a subgradient of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x99.png" xlink:type="simple"/></inline-formula> at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x100.png" xlink:type="simple"/></inline-formula>.</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x101.png" xlink:type="simple"/></inline-formula> be the minimizer of the problem (8).</p><p>If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x102.png" xlink:type="simple"/></inline-formula> satisfies the inequality<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x103.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x104.png" xlink:type="simple"/></inline-formula> is a tolerance then <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x102.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x105.png" xlink:type="simple"/></inline-formula> is an approx-</p><p>imate solution of the problem (7).</p><p>Otherwise define <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x106.png" xlink:type="simple"/></inline-formula> as the index for the most negative<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x106.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x107.png" xlink:type="simple"/></inline-formula>, update the constraints in (8) by including the linear constraint</p><disp-formula id="scirp.50346-formula216"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x108.png"  xlink:type="simple"/></disp-formula><p>and repeat the procedure.</p><p>Recall that our aim is to find <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x109.png" xlink:type="simple"/></inline-formula> such that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x110.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x111.png" xlink:type="simple"/></inline-formula>, but not the solution of the minimization problem (5), (7).</p><p>Theorem 2. If there exists <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x112.png" xlink:type="simple"/></inline-formula> such that</p><disp-formula id="scirp.50346-formula217"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x113.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x114.png" xlink:type="simple"/></inline-formula> is the minimizer of the problem (8), then the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x114.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x115.png" xlink:type="simple"/></inline-formula> is a common solution to (3).</p><p>Proof:</p><disp-formula id="scirp.50346-formula218"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x116.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.50346-formula219"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x117.png"  xlink:type="simple"/></disp-formula><p>and by (5), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x118.png" xlink:type="simple"/></inline-formula>is a common solution to (3).</p><p>For the problem (5), (7) Kelley’s method gives the following</p><p>Algorithm 1.</p><p>Step 1. Take an initial point<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x119.png" xlink:type="simple"/></inline-formula>. Compute <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x120.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x121.png" xlink:type="simple"/></inline-formula>. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x121.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x122.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x121.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x122.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x123.png" xlink:type="simple"/></inline-formula> stop; otherwise continue.</p><p>Step 2. Determine <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x124.png" xlink:type="simple"/></inline-formula> by solving LP problem in (8). If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x125.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x126.png" xlink:type="simple"/></inline-formula> then stop; otherwise continue. Set<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x124.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x126.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x127.png" xlink:type="simple"/></inline-formula>, update the constraints in (8) and repeat the procedure.</p><p>Example 1. Consider the switched system</p><disp-formula id="scirp.50346-formula220"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x128.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.50346-formula221"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x129.png"  xlink:type="simple"/></disp-formula><p>are Hurwitz stable matrices.</p><p>Choose the initial point<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x130.png" xlink:type="simple"/></inline-formula>, then</p><disp-formula id="scirp.50346-formula222"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x131.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x132.png" xlink:type="simple"/></inline-formula>, and</p><p>We obtain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x135.png" xlink:type="simple"/></inline-formula> by solving LP problem in (8). Calculations give the following <xref ref-type="table" rid="table1">Table 1</xref>, and</p><disp-formula id="scirp.50346-formula223"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x136.png"  xlink:type="simple"/></disp-formula><p>Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x137.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x138.png" xlink:type="simple"/></inline-formula>,</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Kelley’s algorithm for Example 1</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >k</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x139.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x140.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x141.png" xlink:type="simple"/></inline-formula></th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x142.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >‒209.7383</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x143.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x144.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >‒127.1153</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x145.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x146.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >‒106.2473</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x147.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x148.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >‒63.4433</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x149.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x150.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x151.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x152.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x153.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x154.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x155.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.2694</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x156.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x157.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.2075</td></tr></tbody></table></table-wrap><disp-formula id="scirp.50346-formula224"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x158.png"  xlink:type="simple"/></disp-formula><p>is a common positive definite solution for</p><disp-formula id="scirp.50346-formula225"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x159.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Stable Member in a Polytope</title><p>This part is devoted to the following question: Given a matrix family <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x160.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x160.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x161.png" xlink:type="simple"/></inline-formula> is a compact, is there a stable matrix in this family?</p><p>In [<xref ref-type="bibr" rid="scirp.50346-ref12">12</xref>] , a numerical algorithm has been proposed for a stable member in the affine matrix family<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x162.png" xlink:type="simple"/></inline-formula>. In this algorithm the uncertainty vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x163.png" xlink:type="simple"/></inline-formula> varies in the whole space<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x164.png" xlink:type="simple"/></inline-formula>. On the other hand we consider the case where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x164.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x165.png" xlink:type="simple"/></inline-formula> varies in a box <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x162.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x164.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x165.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x166.png" xlink:type="simple"/></inline-formula> and use the gradient algorithm for minimization of the nonconvex maximum eigenvalue function. By choosing appropriate step-size, we obtain the convergence.</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x167.png" xlink:type="simple"/></inline-formula> be a basis for the subspace of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x167.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x168.png" xlink:type="simple"/></inline-formula> symmetric matrices and</p><disp-formula id="scirp.50346-formula226"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x169.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.50346-formula227"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x170.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x171.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x171.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x172.png" xlink:type="simple"/></inline-formula>.</p><p>Consider the problem</p><disp-formula id="scirp.50346-formula228"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x173.png"  xlink:type="simple"/></disp-formula><p>Theorem 3. There is a stable matrix in the family <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x174.png" xlink:type="simple"/></inline-formula> if and only if <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x174.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x175.png" xlink:type="simple"/></inline-formula></p><p>Proof:</p><disp-formula id="scirp.50346-formula229"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x176.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.50346-formula230"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x177.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.50346-formula231"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x178.png"  xlink:type="simple"/></disp-formula><p>By Lyapunov theorem, the matrix <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x179.png" xlink:type="simple"/></inline-formula> is stable.</p><p>Example 2. Consider the family of matrices</p><disp-formula id="scirp.50346-formula232"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x180.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.50346-formula233"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x181.png"  xlink:type="simple"/></disp-formula><p>For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x182.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x182.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x183.png" xlink:type="simple"/></inline-formula>is unstable. We apply the gradient algorithm to find a stable member in the family.</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x184.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x184.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x185.png" xlink:type="simple"/></inline-formula>. So</p><disp-formula id="scirp.50346-formula234"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x186.png"  xlink:type="simple"/></disp-formula><p>Then</p><disp-formula id="scirp.50346-formula235"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x187.png"  xlink:type="simple"/></disp-formula><p>Maximum eigenvalue of this matrix and its corresponding unit eigenvector are</p><disp-formula id="scirp.50346-formula236"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x188.png"  xlink:type="simple"/></disp-formula><p>respectively. Gradient of the function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x189.png" xlink:type="simple"/></inline-formula> at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x189.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x190.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.50346-formula237"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x191.png"  xlink:type="simple"/></disp-formula><p>The first tencomponent of the vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x192.png" xlink:type="simple"/></inline-formula> should be on the ten dimensional unit sphere. Therefore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x192.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x193.png" xlink:type="simple"/></inline-formula> and</p><disp-formula id="scirp.50346-formula238"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x194.png"  xlink:type="simple"/></disp-formula><p>After 4 steps, we get</p><disp-formula id="scirp.50346-formula239"><graphic  xlink:href="http://html.scirp.org/file/2-7402419x195.png"  xlink:type="simple"/></disp-formula><p>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x196.png" xlink:type="simple"/></inline-formula>. Therefore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x196.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-7402419x197.png" xlink:type="simple"/></inline-formula> is stable.</p></sec><sec id="s4"><title>4. Conclusion</title><p>Two important problems from control theory are considered: the existence of common quadratic Lyapunov functions for switched linear systems and the existence of a stable member in a matrix polytope. We obtain new conditions which give new effective computational algorithms.</p></sec><sec id="s5"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.50346-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Boyd, S. and Yang, Q. (1989) Structured and Simultaneous Lyapunov Functions for System Stability Problems. International Journal of Control, 49, 2215-2240. http://dx.doi.org/10.1080/00207178908559769</mixed-citation></ref><ref id="scirp.50346-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">B&amp;#252;y&amp;#252;kk&amp;#246;roglu, T., Esen, &amp;#214;. and Dzhafarov, V. (2011) Common Lyapunov Functions for Some Special Classes of Stable Systems. IEEE Transactions on Automatic Control, 56, 1963-1967. http://dx.doi.org/10.1109/tac.2011.2137510</mixed-citation></ref><ref id="scirp.50346-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Cheng, D., Guo, L. and Huang, J. (2003) On Quadratic Lyapunov Functions. IEEE Transactions on Automatic Control, 48, 885-890. http://dx.doi.org/10.1109/tac.2003.811274</mixed-citation></ref><ref id="scirp.50346-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Dayawansa, W.P. and Martin, C.F. (1999) A Converse Lyapunov Theorem for a Class of Dynamical Systems Which Undergo Switching. IEEE Transactions on Automatic Control, 44, 751-760. http://dx.doi.org/10.1109/9.754812</mixed-citation></ref><ref id="scirp.50346-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">King, C. and Shorten, R. (2004) A Singularity Test for the Existence of Common Quadratic Lyapunov Functions for Pairs of Stable LTI Systems. Proceedings of the American Control Conference, Boston, 30 June-2 July 2004, 3881-3884.</mixed-citation></ref><ref id="scirp.50346-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Mason, O. and Shorten, R. (2006) On the Simultaneous Diagonal Stability of a Pair of Positive Linear Systems. Linear Algebra and Its Applications, 413, 13-23. http://dx.doi.org/10.1016/j.laa.2005.07.019</mixed-citation></ref><ref id="scirp.50346-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Narendra, K.S. and Balakrishnan, J. (1994) A Common Lyapunov Function for Stable LTI Systems with Commuting A-Matrices. IEEE Transactions on Automatic Control, 39, 2469-2471. http://dx.doi.org/10.1109/9.362846</mixed-citation></ref><ref id="scirp.50346-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Shorten, R.N. and Narendra, K.S. (2002) Necessary and Sufficient Conditions for the Existence of a Common Quadratic Lyapunov Function for a Finite Number of Stable Second Order Linear Time-Invariant Systems. International Journal of Adaptive Control and Signal Processing, 16, 709-728. http://dx.doi.org/10.1002/acs.719</mixed-citation></ref><ref id="scirp.50346-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Shorten, R.N., Mason, O., Cairbre, F.O. and Curran, P. (2004) A Unifying Framework for the SISO Circle Criterion and Other Quadratic Stability Criteria. International Journal of Control, 77, 1-8. 
http://dx.doi.org/10.1080/00207170310001633321</mixed-citation></ref><ref id="scirp.50346-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Liberzon, D. and Tempo, R. (2004) Common Lyapunov Functions and Gradient Algorithms. IEEE Transactions on Automatic Control, 49, 990-994. http://dx.doi.org/10.1109/tac.2004.829632</mixed-citation></ref><ref id="scirp.50346-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Polyak, B.T. and Shcherbakov, P.S. (2005) Hard Problems in Linear Control Theory: Possible Approaches to Solution. Automation and Remote Control, 66, 681-718. http://dx.doi.org/10.1007/s10513-005-0115-0</mixed-citation></ref><ref id="scirp.50346-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Polyak, B.T. and Shcherbakov, P.S. (1999) Numerical Search of Stable or Unstable Element in Matrix or Polynomial Families: A Unified Approach to Robustness Analysis and Stabilization. Robustness in Identification and Control Lecture Notes in Control and Information Sciences, 245, 344-358. http://dx.doi.org/10.1007/bfb0109879</mixed-citation></ref><ref id="scirp.50346-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Horn, R.A. and Johnson, C.R. (1985) Matrix Analysis. Cambridge University Press, Cambridge.</mixed-citation></ref></ref-list></back></article>