Skip to main content
added 20 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17

System of matrix equations with trace

Problem definition: Let $x_i \in \mathbb{R}^d$ and $a_i \in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Attempted solution: Note that, since $\text{trace}(Ayx^T) = \text{trace}(x^TAy) = x^TAy$, the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l, $$

andwith $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$$X: \text{col}(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$$A = \text{diag}(a)$.

Hence we have

   $$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n$n$-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

System of matrix equations with trace

Problem definition: Let $x_i \in \mathbb{R}^d$ and $a_i \in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l, $$

and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

 $$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

System of matrix equations

Problem definition: Let $x_i \in \mathbb{R}^d$ and $a_i \in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Attempted solution: Note that, since $\text{trace}(Ayx^T) = \text{trace}(x^TAy) = x^TAy$, the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l, $$

with $M(a) = X^TAX$, where $X: \text{col}(X) = \{x_i\}_{i\in[k]}$, and $A = \text{diag}(a)$.

Hence we have  $$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its $n$-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

added 20 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17

Problem definition: Let $x_i\in\mathbb{R}^d$$x_i \in \mathbb{R}^d$ and $a_i\in [0,1]$$a_i \in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$$$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Solution for the simple case $d=1$ (by Carlo Beenakker): $a_i=\frac{x_i^2}{\sum_{j = 1}^kx_j^2}, \forall i \in[k].$

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l , $$ and$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l, $$

and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.Any suggestion on the correctness of this solution would be appreciated.

Problem: Let $x_i\in\mathbb{R}^d$ and $a_i\in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Solution for the simple case $d=1$ (by Carlo Beenakker): $a_i=\frac{x_i^2}{\sum_{j = 1}^kx_j^2}, \forall i \in[k].$

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l , $$ and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

Problem definition: Let $x_i \in \mathbb{R}^d$ and $a_i \in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l, $$

and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

deleted 3 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17

Problem: Let $x_i\in\mathbb{R}^d$ and $a_i\in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Solution for the simple case $d=1$ (by Carlo Beenakker): $a_i=\frac{x_i^2}{\sum_{j = 1}^kx_j^2}, \forall i \in[k].$

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l , $$ and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any hintsuggestion on the correctness of this proceduresolution would be greatly appreciated.

Problem: Let $x_i\in\mathbb{R}^d$ and $a_i\in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Solution for the simple case $d=1$ (by Carlo Beenakker): $a_i=\frac{x_i^2}{\sum_{j = 1}^kx_j^2}, \forall i \in[k].$

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l , $$ and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any hint on the correctness of this procedure would be greatly appreciated.

Problem: Let $x_i\in\mathbb{R}^d$ and $a_i\in [0,1]$, for all $i = 1,\dots, k$ (with $k\geq d$). Define $$M(a) = \sum_{i = 1}^k a_i x_ix_i^T,$$ and assume $M(a) \succ 0.$

Question: Is there any closed-form solution (for $a$) to this set of equations? If not, there is at least an efficient way to solve it?

$$\begin{cases}\text{trace}\left(M(a)^{-1}x_jx_j^T\right) = \text{trace}\left(M(a)^{-1}x_lx_l^T\right), \forall j \neq l,\\ \sum_{i =1}^k a_i= 1.\end{cases}$$

Related question: here.

Solution for the simple case $d=1$ (by Carlo Beenakker): $a_i=\frac{x_i^2}{\sum_{j = 1}^kx_j^2}, \forall i \in[k].$

Attempted solution: Note that the first equation can be rewritten as

$$ x_j^T M(a)^{-1}x_j =x_l^T M(a)^{-1}x_l , $$ and $M(a) = X^TAX$, where $X: col(X) = \{x_i\}_{i\in[k]}$, and $A = diag(a)$.

Hence we have

$$ (x_j^TX^{-T})A^{-1}(X^{-1}x_j) = (x_l^TX^{-T})A^{-1}(X^{-1}x_l) $$

By denoting $\tilde{x}_i = x_i^TX^{-T}$, for all $i\in[k]$, and for a vector $x$, its n-th component $x(n)$, we can write

$$ \sum_{i = 1}^k \frac{\tilde{x_j}(i)^2}{a_i} = \sum_{i = 1}^k \frac{\tilde{x_l}(i)^2}{a_i}. $$

Any suggestion on the correctness of this solution would be appreciated.

Added attempted solution
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
deleted 22 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 4 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 8 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 50 characters in body; edited title
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
deleted 15 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 93 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 93 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 3 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 14 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
added 1 character in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
deleted 54 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
deleted 10 characters in body
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading
Source Link
Apprentice
  • 245
  • 3
  • 17
Loading