Line 17: Line 17:
 
'''Definition: Feasible Direction'''  
 
'''Definition: Feasible Direction'''  
  
<math>\text{A vector } d\in\Re^{n}, d\neq0, \text{ is a feasible direction at } x\in\Omega</math><br>  
+
&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<math>\text{A vector } d\in\Re^{n}, d\neq0, \text{ is a feasible direction at } x\in\Omega</math><br>  
  
<math>\text{if there exists } \alpha_{0}>0 \text{ such that } x+\alpha d\in\Omega \text{ for all } \alpha\in\left[ 0,\alpha_{0}\right]</math>  
+
&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<math>\text{if there exists } \alpha_{0}>0 \text{ such that } x+\alpha d\in\Omega \text{ for all } \alpha\in\left[ 0,\alpha_{0}\right]</math>  
  
 
'''FONC:'''  
 
'''FONC:'''  
  
If x<span style="font-size: 11px;">*</span>&nbsp;is a local minimizer of f over <span class="texhtml">Ω</span>, then for any feasible direction d at x*, we have&nbsp;<sup></sup><sup></sup>  
+
&nbsp; &nbsp; &nbsp; &nbsp; If x<span style="font-size: 11px;">*</span>&nbsp;is a local minimizer of f over <span class="texhtml">Ω</span>, then for any feasible direction d at x*, we have&nbsp;<sup></sup><sup></sup>  
  
<math>d^{T} \nabla f\left ( x^{*} \right )\geq0</math>  
+
&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<math>d^{T} \nabla f\left ( x^{*} \right )\geq0</math>
 +
 
 +
'''FONC Interior Case:'''
 +
 
 +
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<math>\nabla f\left ( x^{*} \right )=0</math>
  
 
'''SONC:&nbsp;'''  
 
'''SONC:&nbsp;'''  
  
Let x* a local minimizer of f and d a feasible direction at x*, then  
+
&nbsp; &nbsp; &nbsp; &nbsp; Let x* a local minimizer of f and d a feasible direction at x*,
 +
 
 +
&nbsp; &nbsp; &nbsp; &nbsp; If <math>d^{T} \nabla f\left ( x^{*} \right )=0</math>&nbsp;, then &nbsp;<math>d^{T} F\left ( x^{*} \right )d\geq 0</math>
 +
 
 +
'''SONC Interior Case:&nbsp;'''
  
<math>d^{T} F\left ( x^{*} \right )d\geq 0</math>  
+
&nbsp; &nbsp; &nbsp; &nbsp; If&nbsp;<math>\nabla f\left ( x^{*} \right )=0</math>&nbsp; , then&nbsp;<math>d^{T} F\left ( x^{*} \right )d\geq 0</math>
  
 
----
 
----
Line 39: Line 47:
 
===== <math>\color{blue}\text{Solution 1:}</math>  =====
 
===== <math>\color{blue}\text{Solution 1:}</math>  =====
  
<math>\text{We need to find a direction }d\text{, such that } \exists\alpha_{0}>0,</math>  
+
<math>\text{We need to find a direction }d\text{, such that } \exists\alpha_{0}>0,</math>&nbsp;
  
<math>\left( \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right) + \alpha d \text{ for all } \alpha\in\Omega \left[0,\alpha_{0}\right]</math><br>  
+
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<math>\left( \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right) + \alpha d \text{ for all } \alpha\in \left[0,\alpha_{0}\right]</math><br>  
  
 
<math>\text{As } x_{1}\geq0, x_{2}\geq0, d= \left( \begin{array}{c} x \\ y \end{array} \right)\text{where } x\in\Re, \text{ and } y\geq0.</math>  
 
<math>\text{As } x_{1}\geq0, x_{2}\geq0, d= \left( \begin{array}{c} x \\ y \end{array} \right)\text{where } x\in\Re, \text{ and } y\geq0.</math>  
Line 49: Line 57:
 
<math>\color{blue}\text{Solution 2:}</math>  
 
<math>\color{blue}\text{Solution 2:}</math>  
  
<math>d\in\Re_{2}, d\neq0 \text{ is a feasible direction at } x^{*} \text{, if } \exists \alpha_{0} \text{ that } \left[ \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right] + \alpha\left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right] \in\Omega \text{ for all } 0\leq\alpha\leq\alpha_{0}</math>&nbsp;<br>  
+
<math>d\in\Re^{2}, d\neq0 \text{ is a feasible direction at } x^{*}</math>&nbsp;<br>
 +
 
 +
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<math>\text{ if } \exists \alpha_{0} \text{ that } \left[ \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right] + \alpha\left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right] \in\Omega \text{ for all } 0\leq\alpha\leq\alpha_{0}</math>
  
 
'''<math>\because \begin{Bmatrix}x\in\Omega: x_{1}\geq0, x_{2}\geq0\end{Bmatrix}</math>'''  
 
'''<math>\because \begin{Bmatrix}x\in\Omega: x_{1}\geq0, x_{2}\geq0\end{Bmatrix}</math>'''  
  
 
<br> <math>\therefore d=
 
<br> <math>\therefore d=
\left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right], d_{1}\in\Re^{2}, d_{2}\geq0</math><br>  
+
\left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right], d_{1}\in\Re, d_{2}\geq0</math>
  
 
----
 
----
Line 65: Line 75:
 
</span></font>  
 
</span></font>  
  
'''<font face="serif"><math>\text{It is equivalent to minimize } f\left(x\right) \text{,  }</math>&nbsp;&nbsp;</font><math>\text{  subject to } g_{1}(x)\leq0, g_{2}(x)\leq0</math>'''  
+
'''<font face="serif"><math>\text{It is equivalent to minimize } f\left(x\right) \text{,  }</math>&nbsp;&nbsp;</font>'''
 +
 
 +
'''<font face="serif"></font>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<math>\text{  subject to } g_{1}(x)\leq0, g_{2}(x)\leq0</math>'''  
  
 
<font color="#ff0000" style="font-size: 17px;">'''<math>\left\{\begin{matrix}
 
<font color="#ff0000" style="font-size: 17px;">'''<math>\left\{\begin{matrix}
l\left(x,\mu \right) = \nabla f(x)+\mu_{1}\nabla g_{1}(x)+ \mu_{2}\nabla g_{2}(x)=\left( \begin{array}{c} 2x_{1}-1+x_{2} \\ 1+x_{1} \end{array} \right) + \left( \begin{array}{c} -\mu_{1} \\ 0 \end{array} \right) +\left( \begin{array}{c} 0 \\ -\mu_{2} \end{array} \right) =0\\  
+
l\left(x,\mu \right) = \nabla f(x)+\mu_{1}\nabla g_{1}(x)+ \mu_{2}\nabla g_{2}(x) \\ =\left( \begin{array}{c} 2x_{1}-1+x_{2} \\ 1+x_{1} \end{array} \right) + \left( \begin{array}{c} -\mu_{1} \\ 0 \end{array} \right) +\left( \begin{array}{c} 0 \\ -\mu_{2} \end{array} \right) =0\\  
 
-\mu_{1}x_{1}-\mu_{2}x_{2} = 0 \\  
 
-\mu_{1}x_{1}-\mu_{2}x_{2} = 0 \\  
 
x_{1} = \frac{1}{2},x_{2} = 0
 
x_{1} = \frac{1}{2},x_{2} = 0
\end{matrix}\right.</math>'''</font><br> <math>\Rightarrow \mu_{1}=0 , \mu_{2}=3/2</math>&nbsp; &nbsp;  
+
\end{matrix}\right.</math>'''</font><br><math>\Rightarrow \mu_{1}=0 , \mu_{2}=3/2</math>&nbsp; &nbsp;  
 +
 
 +
<math>\therefore x^{*} \text{ satisfies FONC}</math>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;
 +
 
 +
<math>\color{green} \text{There exist } \mu \text{ which make point } x^{*} \text{ satisfies FONC.}</math>
  
<math>\therefore x^{*} \text{ satisfies FONC}</math>  
+
<math>\text{SONC: } L(x^{*},\mu^{*}) = \nabla l(x^{*},\mu^{*})=\left( \begin{array}{cc} 2 & 1 \\ 1 & 0 \end{array} \right)</math>  
  
<math>\text{SONC: } L(x^{*},\mu^{*}) = \nabla l(x,\mu)=\left( \begin{array}{cc} 2 & 1 \\ 1 & 0 \end{array} \right)</math>  
+
<font color="#ff0000" size="5">'''<math>T(x^{*},\mu^{*}): \begin{cases} y^{T}\nabla g_{1}(x)=0 \\ y^{T}\nabla g_{2}(x)=0 \end{cases} : \begin{cases} y^{T}\left( \begin{array}{c} -1 \\ 0 \end{array} \right)=0 \\ y^{T}\left( \begin{array}{c} 0 \\-1 \end{array} \right)=0 \end{cases} \Rightarrow y=\left( \begin{array}{c} 0 \\0 \end{array} \right)</math><br>'''</font>  
  
<font color="#ff0000" size="5">'''<math>T(x^{*},\mu^{*}): \begin{cases} y^{T}\nabla g_{1}(x)=0 \\ y^{T}\nabla g_{2}(x)=0 \end{cases} : \begin{cases} y^{T}\left( \begin{array}{c} -1 \\ 0 \end{array} \right)=0 \\ y^{T}\left( \begin{array}{c} 0 \\-1 \end{array} \right)=0 \end{cases} \Rightarrow y=\left( \begin{array}{c} 0 \\0 \end{array} \right)</math><br>'''</font>
+
<font color="#ff0000" size="5"</font><math>\color{green} \text {Here did not use formal set expression.} T\left( x^{* },\mu^{* } \right) \text{ here should be } T\left( x^{* } \right)</math>
  
 
<math>\text{The SONC condition is for all } y\in T \left(x^{*},\mu^{*} \right) , y^{T}L\left(x^{*},\mu^{*} \right)y \geq 0</math>  
 
<math>\text{The SONC condition is for all } y\in T \left(x^{*},\mu^{*} \right) , y^{T}L\left(x^{*},\mu^{*} \right)y \geq 0</math>  
Line 100: Line 116:
 
<math>\text{SONC for local minimizer } x^{*}=\begin{bmatrix} \frac{1}{2}\\0 \end{bmatrix}</math>  
 
<math>\text{SONC for local minimizer } x^{*}=\begin{bmatrix} \frac{1}{2}\\0 \end{bmatrix}</math>  
  
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<math>d^{T} \nabla f\left ( x^{*} \right )=0  \cdots \left ( 1 \right )</math>  
+
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<math>d^{T} \nabla f\left ( x^{*} \right )=0  \cdots \left ( 1 \right )</math>&nbsp; &nbsp; &nbsp;
  
 
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>d^{T} F\left ( x^{*} \right )d\geq 0  \cdots \left ( 2\right )</math><br>  
 
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>d^{T} F\left ( x^{*} \right )d\geq 0  \cdots \left ( 2\right )</math><br>  

Revision as of 12:11, 29 June 2012

ECE Ph.D. Qualifying Exam in "Automatic Control" (AC)

Question 3, Part 1, August 2011

Part 1,2,3,4,5

 $ \color{blue}\text{1. } \left( \text{20 pts} \right) \text{ Consider the optimization problem, } $

               $ \text{maximize} -x_{1}^{2}+x_{1}-x_{2}-x_{1}x_{2} $

               $ \text{subject to } x_{1}\geq0, x_{2}\geq0 $


Definition: Feasible Direction

        $ \text{A vector } d\in\Re^{n}, d\neq0, \text{ is a feasible direction at } x\in\Omega $

        $ \text{if there exists } \alpha_{0}>0 \text{ such that } x+\alpha d\in\Omega \text{ for all } \alpha\in\left[ 0,\alpha_{0}\right] $

FONC:

        If x* is a local minimizer of f over Ω, then for any feasible direction d at x*, we have 

        $ d^{T} \nabla f\left ( x^{*} \right )\geq0 $

FONC Interior Case:

         $ \nabla f\left ( x^{*} \right )=0 $

SONC: 

        Let x* a local minimizer of f and d a feasible direction at x*,

        If $ d^{T} \nabla f\left ( x^{*} \right )=0 $ , then  $ d^{T} F\left ( x^{*} \right )d\geq 0 $

SONC Interior Case: 

        If $ \nabla f\left ( x^{*} \right )=0 $  , then $ d^{T} F\left ( x^{*} \right )d\geq 0 $


$ \color{blue}\left( \text{i} \right) \text{ Characterize feasible directions at the point } x^{*}=\left[ \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right] $

$ \color{blue}\text{Solution 1:} $

$ \text{We need to find a direction }d\text{, such that } \exists\alpha_{0}>0, $ 

         $ \left( \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right) + \alpha d \text{ for all } \alpha\in \left[0,\alpha_{0}\right] $

$ \text{As } x_{1}\geq0, x_{2}\geq0, d= \left( \begin{array}{c} x \\ y \end{array} \right)\text{where } x\in\Re, \text{ and } y\geq0. $


$ \color{blue}\text{Solution 2:} $

$ d\in\Re^{2}, d\neq0 \text{ is a feasible direction at } x^{*} $ 

         $ \text{ if } \exists \alpha_{0} \text{ that } \left[ \begin{array}{c} \frac{1}{2} \\ 0 \end{array} \right] + \alpha\left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right] \in\Omega \text{ for all } 0\leq\alpha\leq\alpha_{0} $

$ \because \begin{Bmatrix}x\in\Omega: x_{1}\geq0, x_{2}\geq0\end{Bmatrix} $


$ \therefore d= \left[ \begin{array}{c} d_{1} \\ d_{2} \end{array} \right], d_{1}\in\Re, d_{2}\geq0 $


$ \color{blue}\left( \text{ii} \right) \text{Write down the second-order necessary condition for } x^{*} \text{. Does the point } x^{*} \text{ satisfy this condition?} $

$ \color{blue}\text{Solution 1:} $

$ \text{Let } f\left(x\right)=x_{1}^{2}-x_{1}+x_{2}+x_{1}x_{2} \text{ , } g_{1}\left(x\right)=-x_{1} \text{ , } g_{2}\left(x\right)=-x_{2} $

$ \text{It is equivalent to minimize } f\left(x\right) \text{, } $  

                                           $ \text{ subject to } g_{1}(x)\leq0, g_{2}(x)\leq0 $

$ \left\{\begin{matrix} l\left(x,\mu \right) = \nabla f(x)+\mu_{1}\nabla g_{1}(x)+ \mu_{2}\nabla g_{2}(x) \\ =\left( \begin{array}{c} 2x_{1}-1+x_{2} \\ 1+x_{1} \end{array} \right) + \left( \begin{array}{c} -\mu_{1} \\ 0 \end{array} \right) +\left( \begin{array}{c} 0 \\ -\mu_{2} \end{array} \right) =0\\ -\mu_{1}x_{1}-\mu_{2}x_{2} = 0 \\ x_{1} = \frac{1}{2},x_{2} = 0 \end{matrix}\right. $
$ \Rightarrow \mu_{1}=0 , \mu_{2}=3/2 $   

$ \therefore x^{*} \text{ satisfies FONC} $              

$ \color{green} \text{There exist } \mu \text{ which make point } x^{*} \text{ satisfies FONC.} $

$ \text{SONC: } L(x^{*},\mu^{*}) = \nabla l(x^{*},\mu^{*})=\left( \begin{array}{cc} 2 & 1 \\ 1 & 0 \end{array} \right) $

$ T(x^{*},\mu^{*}): \begin{cases} y^{T}\nabla g_{1}(x)=0 \\ y^{T}\nabla g_{2}(x)=0 \end{cases} : \begin{cases} y^{T}\left( \begin{array}{c} -1 \\ 0 \end{array} \right)=0 \\ y^{T}\left( \begin{array}{c} 0 \\-1 \end{array} \right)=0 \end{cases} \Rightarrow y=\left( \begin{array}{c} 0 \\0 \end{array} \right) $

<font color="#ff0000" size="5"</font>$ \color{green} \text {Here did not use formal set expression.} T\left( x^{* },\mu^{* } \right) \text{ here should be } T\left( x^{* } \right) $

$ \text{The SONC condition is for all } y\in T \left(x^{*},\mu^{*} \right) , y^{T}L\left(x^{*},\mu^{*} \right)y \geq 0 $

$ y^{T}L\left(x^{*},\mu^{*} \right)y =0 \geq 0 \text{. So } x^{*} \text{satisfies SONC.} $


$ \color{blue}\text{Solution 2:} $

$ \text{The problem is equivalent to min} f\left(x_{1},x_{2}\right) = x_{1}^{2}-x_{1}+x_{2}+x_{1}x_{2} $  

                                                                      $ \text{subject to } x_{1}\leq0, x_{2}\leq0 $

$ Df\left ( x \right )=\left ( \nabla f\left ( x \right ) \right )^{T} = \left [ \frac{\partial f}{\partial x_{1}}\left ( x \right ),\frac{\partial f}{\partial x_{2}}\left ( x \right ) \right ]=\left [ 2x_{1}-1+x_{2},1+x_{1} \right ] $

$ F\left ( x \right ) =D^{2}f\left ( x \right )=\begin{bmatrix} \frac{\partial^{2} f}{\partial x_{1}^{2}}\left ( x \right ) & \frac{\partial^{2} f}{\partial x_{2}\partial x_{1}}\left ( x \right )\\ \frac{\partial^{2} f}{\partial x_{1}\partial x_{2}}\left ( x \right ) & \frac{\partial^{2} f}{\partial x_{2}^{2}}\left ( x \right ) \end{bmatrix}=\left [ \begin{array}{cc} 2 & 1 \\ 1 & 0 \end{array} \right ] $

$ \text{SONC for local minimizer } x^{*}=\begin{bmatrix} \frac{1}{2}\\0 \end{bmatrix} $

                  $ d^{T} \nabla f\left ( x^{*} \right )=0 \cdots \left ( 1 \right ) $     

                  $ d^{T} F\left ( x^{*} \right )d\geq 0 \cdots \left ( 2\right ) $

$ \text{For (1), } \begin{bmatrix} d_{1} & d_{2} \end{bmatrix}\begin{bmatrix} 0\\ \frac{3}{2}\end{bmatrix} =0 \Rightarrow d_{2}=0, d_{1}\in\Re $

$ \text{For (2), } F\left ( x \right ) = \begin{bmatrix} 2 &1 \\ 1 &0\end{bmatrix}>0 $       $ \color{green} A=\begin{bmatrix} a &b \\ c &d\end{bmatrix} \text{ is positive definite when } a>0 \text{ and } ac-b^{2}>0 $

$ \therefore \text{ for all } d\in\Re^{n}, d^{T}F\left ( x^{*} \right )d\geq 0 $

$ \text{The point } x^{*}=\begin{bmatrix} \frac{1}{2}\\0 \end{bmatrix} \text{ satisfies SONC for local minimizer.} $


$ \color{blue}\text{Related Problem:} $


Automatic Control (AC)- Question 3, August 2011

Go to


Back to ECE Qualifying Exams (QE) page

Alumni Liaison

Abstract algebra continues the conceptual developments of linear algebra, on an even grander scale.

Dr. Paul Garrett