Skip to main content
added 699 characters in body
Source Link
loup blanc
  • 3.7k
  • 17
  • 32

OvrOver the real numbers:

Ok, assume that $\alpha<L$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

EDIT: Now we are over $\mathbb{C}$. Remark 1. The case 1 is impossible. Proof: $A^HA$ is hermit $\geq 0$. Then we may assume $A^HA=diag(t_1,\cdots,t_L,0_{K-L})$ with $t_i>0$. Let $P$ s.t. $A^HAP=0$. Then the $L$ first rows of $P$ are $0$. Then $||I-P||^2\geq L>\alpha$.

Remark 2. For the case 2., we have a good candidate $P=(I+\dfrac{1}{\lambda} A^HA)^{-1}$. It can happen that there is a better complex solution, but it is unlikely (except if you are unlucky!). Here $P=diag(\dfrac{\lambda}{\lambda+t_1},\cdots,\dfrac{\lambda}{\lambda+t_L},I_{K-L})$ and $\lambda$ is a solution of $L+\sum_{i=1}^L \dfrac{\lambda^2}{(\lambda+t_i)^2}-2\sum_{i=1}^L \dfrac{\lambda}{\lambda+t_i}=\alpha$.

Ovr the real numbers:

Ok, assume that $\alpha<L$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

Over the real numbers:

Ok, assume that $\alpha<L$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

EDIT: Now we are over $\mathbb{C}$. Remark 1. The case 1 is impossible. Proof: $A^HA$ is hermit $\geq 0$. Then we may assume $A^HA=diag(t_1,\cdots,t_L,0_{K-L})$ with $t_i>0$. Let $P$ s.t. $A^HAP=0$. Then the $L$ first rows of $P$ are $0$. Then $||I-P||^2\geq L>\alpha$.

Remark 2. For the case 2., we have a good candidate $P=(I+\dfrac{1}{\lambda} A^HA)^{-1}$. It can happen that there is a better complex solution, but it is unlikely (except if you are unlucky!). Here $P=diag(\dfrac{\lambda}{\lambda+t_1},\cdots,\dfrac{\lambda}{\lambda+t_L},I_{K-L})$ and $\lambda$ is a solution of $L+\sum_{i=1}^L \dfrac{\lambda^2}{(\lambda+t_i)^2}-2\sum_{i=1}^L \dfrac{\lambda}{\lambda+t_i}=\alpha$.

deleted 6 characters in body
Source Link
loup blanc
  • 3.7k
  • 17
  • 32

if $AP=0$, then $\alpha=K$.Ovr the real numbers:

Ok, assume that $\alpha<K$$\alpha<L$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

if $AP=0$, then $\alpha=K$.

Ok, assume that $\alpha<K$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

Ovr the real numbers:

Ok, assume that $\alpha<L$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.

Source Link
loup blanc
  • 3.7k
  • 17
  • 32

if $AP=0$, then $\alpha=K$.

Ok, assume that $\alpha<K$ and the matrices are real. Let $f(P)=tr(P^TA^TAP), g(P)=K-2tr(P)+tr(PP^T)-\alpha$. By the Lagrange's method, we seek $\lambda\in\mathbb{R}$ s.t. $Df_P+\lambda Dg_P=0$, that is, for any matrix $H$,

$tr(H^TA^TAP+P^TA^TAH)+\lambda(-2tr(H)+tr(PH^T+HP^T))=0$, that is

$tr(P^TA^TAH)=\lambda tr(H-P^TH)$. This is equivalent to $P^TA^TA=\lambda(I-P^T).$

Case 1. $\lambda=0$. Then $P^TA^TA=0$ and $f(P)=0$. It remains to find $P$ s.t.one has also $g(P)=\alpha$.

Case 2. $\lambda\not=0$. Then $P^T=(I+\dfrac{1}{\lambda} A^TA)^{-1}$ and $P$ is symmetric. It remains to find $\lambda\not= 0$ s.t. $(I+\dfrac{1}{\lambda} A^TA)$ is invertible and $K-2tr(P)+tr(P^2)=\alpha$.