We usually consider $\dot{x} = f(x)$, where $x$ is a vector.
Now, I want to consider $$\dot{X}=f(X,U),$$ where $X$ is a square matrix $\mathbb{R}^{n\times n}$ state, $U$ is a square matrix variable $\mathbb{R}^{n\times n}$, like a control input. If I want to consider its stability, I have to use Lyapunov stability theorem, $$V(X)>0,$$ for $X\neq 0$, And $$\dot{V}<0$$along $\dot{X}=f(X)$. One way to do this i let $\dot{X}=-\nabla V$, so $$\dot{V} = \langle\nabla V, \dot{X}\rangle=-\langle\nabla V, \nabla V\rangle = -\|\nabla V\|_F^2.$$ So I have to let $$f(X,U) = -\nabla V$$ to solve for $U$. However, it is not the case of vector, we cannot easily "make them equal", for example, they may not have the same rank.
------------------------------new idea-----------------------------
To relax such condition, I would consider the set $$\mathcal{C} = \{W : \langle W,-\nabla V \rangle \geq 0\},$$ which implies $W\in \mathcal{C}$ points to the descent direction if $V$ is a convex function in $Q$. So instead of letting $f(X,U) = -\nabla V$, let $f(X,U) = W$, $W\in \mathcal{C}$. So now I probably have more flexible choice.
Any reference or paper about the matrix version of dynamic system / stability analysis / Lyapunov theory?
Sincerely appreciate this help.