5
$\begingroup$

Let $\boldsymbol{V}_{1},\dots,\boldsymbol{V}_{n}\in\mathbb{R}^{d\times m}$ be $n$ “tall” matrices (where $d\ge m$) with orthonormal columns.

And let $\boldsymbol{P}_{1},\dots,\boldsymbol{P}_{n}\in\mathbb{R}^{d\times d}$ be the orthogonal projection matrices defined as $\boldsymbol{P}_{i}=\boldsymbol{V}_{i}\boldsymbol{V}_{i}^{\top}$.

Finally, let $\boldsymbol{T}$ be an operator defined the concatenation of the projection matrices $\boldsymbol{T}=\boldsymbol{P}_{n}\cdots\boldsymbol{P}_{1}$.

Question: does there exists a sequence $\boldsymbol{V}_{1},\dots,\boldsymbol{V}_{n}$ such that the concatenation operator $\boldsymbol{T}$ is nilpotent with index $k\ge3$? That is,

$$\boldsymbol{T}^{k-1}=\left(\boldsymbol{P}_{n}\cdots\boldsymbol{P}_{1}\right)^{k-1}\neq\boldsymbol{0}_{d\times d},\,\,\,\,\,\,\,\boldsymbol{T}^{k}=\left(\boldsymbol{P}_{n}\cdots\boldsymbol{P}_{1}\right)^{k}=\boldsymbol{0}_{d\times d}$$


Special case (example): when $n=3$ and $d=2$, we can choose projections such that $\boldsymbol{T}$ is a nilpotent operator with an index of $k=2$.

Choose $\boldsymbol{v}_{1}=\left[0,1\right]^{\top}, \boldsymbol{v}_{2}=\frac{1}{\sqrt{2}}\left[1,1\right]^{\top},\boldsymbol{v}_{3}=\left[1,0\right]^{\top}$.

Then, the projection matrices are $$\boldsymbol{P}_{1}=\boldsymbol{v}_{1}\boldsymbol{v}_{1}^{\top}=\left[\begin{array}{cc} 0 & 0\\ 0 & 1 \end{array}\right],\,\,\,\boldsymbol{P}_{2}=\frac{1}{2}\left[\begin{array}{cc} 1 & 1\\ 1 & 1 \end{array}\right],\,\,\,\boldsymbol{P}_{3}=\left[\begin{array}{cc} 1 & 0\\ 0 & 0 \end{array}\right]$$ And the concatenation of these matrices is $$\boldsymbol{T}=\boldsymbol{P}_{3}\boldsymbol{P}_{2}\boldsymbol{P}_{1}=\frac{1}{2}\left[\begin{array}{cc} 0 & 1\\ 0 & 0 \end{array}\right]$$ which is a nilpotent matrix with an index of $k=2$ (i.e., $\boldsymbol{T}^{2}=\boldsymbol{0}_{d\times d}$).

$\endgroup$
4
  • 1
    $\begingroup$ Just to be sure, you don't assume that the columns of $V_i$ are orthonormal? Because if you just assume orthogonality and not normality $V_iV_i^T$ is not an orthogonal projection. $\endgroup$ Nov 8, 2021 at 12:54
  • $\begingroup$ Edited so the matrices are orthonormal and the projections are orthogonal. Thanks $\endgroup$
    – Itay
    Nov 8, 2021 at 12:57
  • 1
    $\begingroup$ The way you write things now, $P_3P_2P_1$ is not the matrix you write, since $P_3P_2P_1$ annihilates the first basis vector, and your matrix does not. $\endgroup$ Nov 8, 2021 at 12:59
  • 1
    $\begingroup$ @VladimirDotsenko OK I see it now. Thanks. Flipped the order accidently. Will fix. $\endgroup$
    – Itay
    Nov 8, 2021 at 13:03

2 Answers 2

5
$\begingroup$

I think your example is easily generalisable for any index. For example, let $$ Q_1=P_1\oplus(1), Q_2=P_2\oplus(1), Q_3=P_3\oplus(1)=(1)\oplus P_1, Q_4=(1)\oplus P_2, Q_5=(1)\oplus P_3. $$ Then $Q_5Q_4Q_3Q_2Q_1$ is nilpotent of index 3, and I think that a similar construction with $2n-1$ matrices of size $n$ gives nilpotence of index $n-1$.

$\endgroup$
2
  • $\begingroup$ Great! Many thanks. Just a comment: I think this construction requires $2^{n-1}+1$ matrices rather than $2n-1$. $\endgroup$
    – Itay
    Nov 8, 2021 at 14:26
  • $\begingroup$ @Itay the description in my answer does indeed suggest some sort of doubling, you are right. But what I was thinking about was product of the following $2n-1$ matrices: $n$ diagonal matrices with $n-1$ ones and one zero, and $n-1$ block-diagonal matrices with $n-2$ ones and one block equal to your $P_2$. I think that multiplying all of them in the most "natural" order gives an example. $\endgroup$ Nov 8, 2021 at 17:46
7
$\begingroup$

The paper Products of orthogonal projections shows in Theorem 1 that any singular matrix with norm less than $1$ is a product of projection matrices and gives an estimate on how many projection matrices you need. You can always scale a nilpotent matrix to have norm less than $1$ and so you can obtain a nilpotent matrix of any index that is a product of orthogonal projections. I'm not sure why the tall thing is needed here since you allow $d=m$ and you can just take $V_i=P_i$ for any $d\times d$ projection matrix $P_i$.

$\endgroup$
0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.