Suppose we are given a list of $N$ positive definite quadratic forms $X^TQ_k X$ (where $k\in[1,N]$ and $Q_k\in\mathbb{R}^{p\times p}$ $\forall k$), and a positive vector $V$ of same length $N$ i.e. $V=(v_k>0)_{k\in[1,N]}$.
My question is the following: how do we find the vector $X^*\in\mathbb{R}^p$ minimising the squared distance between the array of quadratic forms and vector $V$, i.e. solving: $$X^*=\underset{X}{\mathrm{argmin}} \sum_{k=1}^N \left(X^TQ_k X - v_k\right)^2 \quad ?$$ The scalar case ($p$=1) is easy to solve, giving $X^*=\pm\sqrt{\frac{\sum q_kv_k}{\sum q_k^2}}$, but I'm not sure how to proceed in the general case $p$>1.