The set of the left essential polynomials could always be recovered with the help of the so-called conformation procedure of the right and left essential polynomials. Let us describe how for the given set of the right essential polynomials $R_1\left(z\right),$  $\dots ,R_{2(p+q)}\left(z\right)$  , $R_j\left(z\right)\in {\mathbb{C}}^{2q\times 1}\left[z\right]$  one can construct the conforming left essential polynomials $L_1\left(z\right),\dots ,L_{2(p+q)}\left(z\right)$  , $L_j\left(z\right)\in {\mathbb{C}}^{1\times 2p}\left[z\right]$  .

We introduce the matrix of the right essential polynomials $$\mathcal{ℛ}\left(z\right)=\left(R_1\left(z\right)\dots R_{2(p+q)}\left(z\right)\right)$$  and find the matrix polynomial ${\alpha }_{-}\left(z\right)$  from the next decomposition $$A\left(z\right)\mathcal{ℛ}\left(z\right)={\alpha }_{-}\left(z\right)d\left(z\right)-z^{n+1}{\beta }_{+}\left(z\right).$$  Here $d\left(z\right)=diag\left[z^{{\mu }_1},\dots ,z^{{\mu }_{2(p+q)}}\right],$  ${\beta }_{+}\left(z\right)$  is the matrix polynomial in $z$  , ${\alpha }_{-}\left(z\right)$  is the matrix polynomial in $z^{-1}$  . Both of them have sizes $2p\times 2(p+q)$  .

Let $U_{-}\left(z\right)$  be the matrix polynomial in $z^{-1}$  such that:

$$U_{-}\left(z\right)=\left(\begin{array}{c}{\mathcal{ℛ}}_{-}\left(z\right)\\ {\alpha }_{-}\left(z\right)\end{array}\right),$$  with ${\mathcal{ℛ}}_{-}\left(z\right)=z^{-m-1}\mathcal{ℛ}\left(z\right)d^{-1}\left(z\right).$  The matrix polynomial $U_{-}\left(z\right)$  is shown in [8] to be unimodular, i.e. its determinant is equal to a constant.

We pick the $2(p+q)\times 2p$  block $\mathcal{ℒ}\left(z\right)$  out $U_{-}^{-1}\left(z\right)$  :

$$U_{-}^{-1}\left(z\right)=\left(\begin{array}{cc}*& \mathcal{ℒ}\left(z\right)\end{array}\right).$$  The matrix polynomial $$\mathcal{ℒ}\left(z\right)=\left(\begin{array}{c}{L}_1\left(z\right)\\ ...\\ L_{2(p+q)}\left(z\right)\end{array}\right)$$  turns out to be the matrix of the conforming left essential polynomials.

The case when $\alpha =0$  or $p\ge q$  may be considered in a similar manner with help of the left essential polynomials.

Now we may present the formula (5.13) from [8] for the generalized inverse of $T_A$  :

$$\begin{array}{c}{T}_A^{†}=\left(\begin{array}{ccc}{\mathcal{ℛ}}_0& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m& ...& {\mathcal{ℛ}}_0\end{array}\right)\Pi \left(\begin{array}{ccc}{\mathcal{ℒ}}_0& ...& {\mathcal{ℒ}}_{-n}\\ ...& ...& ...\\ 0& ...& {\mathcal{ℒ}}_0\end{array}\right).\end{array}$$

(2.3)

Here ${\mathcal{ℛ}}_j\in {\mathbb{C}}^{2q\times 2(p+q)},$  ${\mathcal{ℒ}}_j\in {\mathbb{C}}^{2(p+q)\times 2p}$  are the coefficients of the matrix polynomials $\mathcal{ℛ}\left(z\right)$  , $\mathcal{ℒ}\left(z\right)$  , respectively, and $R_j\left(z\right),L_j\left(z\right)$  are the conforming right and left essential polynomials of the sequence $A_{-m},\dots ,A_0,$  $\dots ,A_n$  . The generalized inversion for matrix $A$  is meant to be the matrix $A^{†}$  such that $A{A}^{†}A=A.$  The matrix $\Pi$  is constructed as follows. Let ${\lambda }_1,\dots ,{\lambda }_r$  be the distinct essential indices of the sequence $A_{-m},\dots ,A_0,\dots ,A_n$  and let ${\nu }_1,\dots ,{\nu }_r$  be their multiplicities ( ${\nu }_1+\dots +{\nu }_r=2(p+q)$  ).

Then $$\Pi =\left(\begin{array}{cccc}{\Pi }_0& {\Pi }_{-1}& ...& {\Pi }_{-n}\\ {\Pi }_1& {\Pi }_0& ...& {\Pi }_{-n+1}\\ ...& ...& ...& ...\\ {\Pi }_m& {\Pi }_{m-1}& ...& {\Pi }_{m-n}\end{array}\right).$$  Here ${\Pi }_k=0$  for $-n\le k\le m,k\ne -{\lambda }_1,\dots ,-{\lambda }_r$  , ${\Pi }_{-{\lambda }_j}=\parallel {\varepsilon }_i^j{\delta }_{ik}{\parallel }_{i,k=1}^{2(p+q)}$  , $${\varepsilon }_i^j=\{\begin{array}{cc}1,& i={\nu }_1+\cdots +{\nu }_{j-1}+1,\dots ,{\nu }_1+\cdots +{\nu }_j,\\ 0,& otherwise.\end{array}$$  The following partition of the right essential polynomials $R_j\left(z\right)$  will be useful for the generalized inversion of the $T+H$  matrix:

$$R_j\left(z\right)=\left(\begin{array}{c}{R}_j^1\left(z\right)\\ \\ R_j^2\left(z\right)\end{array}\right).$$  Here $R_j^{1,2}\in {\mathbb{C}}^{q\times 1}\left[z\right].$  In similar way we partition the left essential polynomials:

$$L_j\left(z\right)=\left(\begin{array}{cc}{L}_j^1\left(z\right)& L_j^2\left(z\right)\end{array}\right),$$  with $L_j^{1,2}\in {\mathbb{C}}^{1\times p}\left[z\right].$  Then the matrix of these essential polynomials may be represented as:

$$\begin{array}{c}\mathcal{ℛ}\left(z\right)=\left(\begin{array}{c}{\mathcal{ℛ}}^1\left(z\right)\\ {\mathcal{ℛ}}^2\left(z\right)\end{array}\right),\mathcal{ℒ}\left(z\right)=\left(\begin{array}{cc}{\mathcal{ℒ}}^1\left(z\right)& {\mathcal{ℒ}}^2\left(z\right)\end{array}\right),\end{array}$$

(2.4)

with ${\mathcal{ℛ}}^{1,2}\left(z\right)\in {\mathbb{C}}^{q\times 2(p+q)},$  ${\mathcal{ℒ}}^{1,2}\left(z\right)\in {\mathbb{C}}^{2(p+q)\times p}.$ 

3. Generalized inversion

In the section we will present our main result.

Let us denote $$T_{{\mathcal{ℛ}}_j}=\left(\begin{array}{ccc}{\mathcal{ℛ}}_0^j& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m^j& ...& {\mathcal{ℛ}}_0^j\end{array}\right),T_{{\mathcal{ℒ}}_j}=\left(\begin{array}{ccc}{\mathcal{ℒ}}_0^j& ...& {\mathcal{ℒ}}_{-n}^j\\ ...& ...& ...\\ 0& ...& {\mathcal{ℒ}}_0^j\end{array}\right),j=1,2,$$  where ${\mathcal{ℛ}}_k^j\left({\mathcal{ℒ}}_k^j\right)$  are the coefficients of the polynomials ${\mathcal{ℛ}}^j\left({\mathcal{ℒ}}^j\right)$  . We also put $H_{{\mathcal{ℛ}}_2}=J{T}_{{\mathcal{ℛ}}_2},$  $H_{{\mathcal{ℒ}}_1}=T_{{\mathcal{ℒ}}_1}J$  .

Theorem 3.1. Generalized inverses of the $T+H$  and $T-H$  matrices are found by the formulas:

$$\begin{array}{c}{\left(T\pm H\right)}^{†}=\frac{1}{2}\left(T_{{\mathcal{ℛ}}_1}\pm H_{{\mathcal{ℛ}}_2}\right)\Pi \left(T_{{\mathcal{ℒ}}_2}\pm H_{{\mathcal{ℒ}}_1}\right).\end{array}$$

(3.1)

If $T\pm H$  is invertible (one-sided invertible), then ${\left(T\pm H\right)}^{†}$  is its inverse (one-sided inverse) matrix.

Proof. Let us construct a generalized inversion to $T_A\equiv T_0$  according to formula ( 2.3 ).

We are going to pass from block Toeplitz matrix $T_A$  to the mosaic matrix $$M_A=\left(\begin{array}{ccccccc}{b}_n& ...& b_{n+m}& & a_{n-m}& ...& a_n\\ b_{n-1}& ...& b_{n+m-1}& & a_{n-m-1}& ...& a_{n-1}\\ ...& ...& ...& & ...& ...& ...\\ b_0& ...& b_m& & a_{-m}& ...& a_0\\ a_0& ...& a_{-m}& & b_m& ...& b_0\\ a_1& ...& a_{-m+1}& & b_{m+1}& ...& b_1\\ ...& ...& ...& & ...& ...& ...\\ a_n& ...& a_{n-m}& & b_{n+m}& ...& b_n\end{array}\right).$$  At first, according to the block structure of $A_j$  ( 2.1 ), we partition each block column $X_j$  of the matrix $T_A$  into two block columns $X_j^1,X_j^2$  with sizes $2p(n+1)\times q$  :

$$X_j=\left(\begin{array}{cc}{X}_j^1& X_j^2\end{array}\right).$$  Then permute new block columns in $T_A$  and construct the matrix $$\left(\begin{array}{cccccc}{X}_1^1& ...& X_m^1& X_1^2& ...& X_m^2\end{array}\right)=$$  $$\left(\begin{array}{ccccccccc}{b}_n& b_{n+1}& ...& b_{n+m}& & a_{n-m}& a_{n-m+1}& ...& a_n\\ a_0& a_{-1}& ...& a_{-m}& & b_m& b_{m-1}& ...& b_0\\ b_{n-1}& b_n& ...& b_{n+m-1}& & a_{n-m-1}& a_{n-m}& ...& a_{n-1}\\ a_1& a_0& ...& a_{-m+1}& & b_{m+1}& b_m& ...& b_1\\ ...& ...& ...& ...& & ...& ...& ...& ...\\ b_0& b_1& ...& b_m& & a_{-m}& a_{-m+1}& ...& a_0\\ a_n& a_{n-1}& ...& a_{n-m}& & b_{n+m}& b_{n+m-1}& ...& b_n\end{array}\right).$$  This matrix is evident to be obtained by multiplying $T_A$  on a permutation matrix $P_2$  . Then we will do the analogous permutation with block rows in $T_A{P}_2$  .

As a result, we will get the matrix $P_1{T}_A{P}_2$  , where $P_1$  is a permutation matrix.

The matrix $P_1{T}_A{P}_2$  coincides with $M_A:$  $$M_A=P_1{T}_A{P}_2.$$  Thus we have passed from the block Toeplitz matrix $T_A$  to the mosaic matrix $M_A$  .

Taking into account that for a permutation matrix $P$  the equality $P^{-1}=P^t$  holds, we get the generalized inversion for $M_A:$  $$M_A^{†}=P_2^t{T}_A^{†}{P}_1^t.$$  Let us specify the structure of factors in $M_A^{†}=P_2^t{T}_A^{†}{P}_1^t$  . The operations which $P_2$  has done with the block columns of $T_A$  , the matrix $P_2^t$  now will carry out with the block rows of the matrix $$\left(\begin{array}{ccc}{\mathcal{ℛ}}_0& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m& ...& {\mathcal{ℛ}}_0\end{array}\right).$$  Thus $$P_2^t\left(\begin{array}{ccc}{\mathcal{ℛ}}_0& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m& ...& {\mathcal{ℛ}}_0\end{array}\right)=\left(\begin{array}{ccc}{\mathcal{ℛ}}_0^1& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m^1& ...& {\mathcal{ℛ}}_0^1\\ {\mathcal{ℛ}}_0^2& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m^2& ...& {\mathcal{ℛ}}_0^2\end{array}\right)\equiv \left(\begin{array}{c}{T}_{{\mathcal{ℛ}}_1}\\ T_{{\mathcal{ℛ}}_2}\end{array}\right),$$  where ${\mathcal{ℛ}}_j^{1,2}$  are the coefficients of the matrix polynomials ${\mathcal{ℛ}}^{1,2}\left(z\right)$  , presented in ( 2.4 ).

Similarly, we have $\left(\begin{array}{ccc}{\mathcal{ℒ}}_0& ...& {\mathcal{ℒ}}_{-n}\\ ...& ...& ...\\ 0& ...& {\mathcal{ℒ}}_0\end{array}\right)P_1^t=$  $$=\left(\begin{array}{cccccc}{\mathcal{ℒ}}_0^1& ...& {\mathcal{ℒ}}_{-n}^1& {\mathcal{ℒ}}_0^2& ...& {\mathcal{ℒ}}_{-n}^2\\ ...& ...& ...& ...& ...& ...\\ 0& ...& {\mathcal{ℒ}}_0^1& 0& ...& {\mathcal{ℒ}}_0^2\end{array}\right)\equiv \left(\begin{array}{cc}{T}_{{\mathcal{ℒ}}_1}& T_{{\mathcal{ℒ}}_2}\end{array}\right).$$  Then $$M_A^{†}=\left(\begin{array}{c}{T}_{{\mathcal{ℛ}}_1}\\ T_{{\mathcal{ℛ}}_2}\end{array}\right)\Pi \left(\begin{array}{cc}{T}_{{\mathcal{ℒ}}_1}& T_{{\mathcal{ℒ}}_2}\end{array}\right).$$  Let us apply now the well-known method [2] of reducing the mosaic matrix $M_A$  to the block-diagonal matrix formed from the Toeplitz-plus-Hankel and Toeplitz-minus-Hankel matrices:

$$M_A=\frac{1}{2}\left(\begin{array}{cc}J& J\\ I& -I\end{array}\right)\left(\begin{array}{cc}T+H& 0\\ 0& T-H\end{array}\right)\left(\begin{array}{cc}I& J\\ -I& J\end{array}\right).$$  We obtain that the matrix $G=\frac{1}{2}\left(\begin{array}{cc}I& J\\ -I& J\end{array}\right)M_A^{†}\left(\begin{array}{cc}J& J\\ I& -I\end{array}\right)=$  $$=\frac{1}{2}\left(\begin{array}{c}{T}_{{\mathcal{ℛ}}_1}+J{T}_{{\mathcal{ℛ}}_2}\\ -T_{{\mathcal{ℛ}}_1}+J{T}_{{\mathcal{ℛ}}_2}\end{array}\right)\Pi \left(\begin{array}{cc}{T}_{{\mathcal{ℒ}}_1}J+T_{{\mathcal{ℒ}}_2}& T_{{\mathcal{ℒ}}_1}J-T_{{\mathcal{ℒ}}_2}\end{array}\right)$$  is the generalized inversion for the matrix $\left(\begin{array}{cc}T+H& 0\\ 0& T-H\end{array}\right)$  .

Let us present $G$  in the block form $$G=\left(\begin{array}{cc}{G}_{11}& G_{12}\\ G_{21}& G_{22}\end{array}\right),$$  with $G_{ij}\in {\mathbb{C}}^{(m+1)q\times (n+1)p}.$  It is easy to get that $$G_{11}=\frac{1}{2}\left(T_{{\mathcal{ℛ}}_1}+H_{{\mathcal{ℛ}}_2}\right)\Pi \left(T_{{\mathcal{ℒ}}_2}+H_{{\mathcal{ℒ}}_1}\right),$$  $$G_{22}=\frac{1}{2}\left(T_{{\mathcal{ℛ}}_1}-H_{{\mathcal{ℛ}}_2}\right)\Pi \left(T_{{\mathcal{ℒ}}_2}-H_{{\mathcal{ℒ}}_1}\right)$$  are generalized inverses to $T+H$  $(T-H)$  .

The theorem statement concerning the invertibility (one-sided invertibility) is evident. The theorem has been proved.

Given $T\pm H$  matrices are block matrices with the sizes of their blocks $p\times q.$  The factors in the inverse formulas ( 3.1 ) have blocks with sizes $q\times 2(p+q)$  , $2(p+q)\times p$  . The compact form of the generalized inversion is in many respects because of such factors sizes. Sometimes it is convenient to have a formula for a generilized inversion where factors have blocks with sizes $q\times q,q\times p,p\times p$  .

In order to obtain it we partition $\mathcal{ℛ}\left(z\right)=\left(R_1\left(z\right)\dots R_{2(p+q)}\left(z\right)\right)$  into blocks:

$$\mathcal{ℛ}\left(z\right)=\left(\begin{array}{cccc}{\mathcal{ℛ}}_{11}& {\mathcal{ℛ}}_{12}& {\mathcal{ℛ}}_{13}& {\mathcal{ℛ}}_{14}\\ {\mathcal{ℛ}}_{21}& {\mathcal{ℛ}}_{22}& {\mathcal{ℛ}}_{23}& {\mathcal{ℛ}}_{24}\end{array}\right).$$  Here ${\mathcal{ℛ}}_{ij}$  have the sizes $q\times q$  for $i,j=1,2,$  and $q\times p$  for $i=1,2,$  $j=3,4.$  Then we will do the analogous partition with the matrix of the left essential polynomials:

$$\mathcal{ℒ}\left(z\right)=\left(\begin{array}{cc}{\mathcal{ℒ}}_{11}& {\mathcal{ℒ}}_{12}\\ {\mathcal{ℒ}}_{21}& {\mathcal{ℒ}}_{22}\\ {\mathcal{ℒ}}_{31}& {\mathcal{ℒ}}_{32}\\ {\mathcal{ℒ}}_{41}& {\mathcal{ℒ}}_{42}\end{array}\right),$$  here ${\mathcal{ℒ}}_{ij}$  have the sizes $q\times p$  for $i,j=1,2,$  and $p\times p$  for $i=3,4,j=1,2.$  Let us also partition $D=diag[z^{{\mu }_1}\dots z^{{\mu }_{2(p+q)}}]=\left(d_1{d}_2{d}_3{d}_4\right),$  where $d_{1,2}$  are diagonal matrices with the sizes $q\times q$  and $d_{3,4}$  are ones with the sizes $p\times p.$  We denote for $i,j=1,\dots ,4$  $$T_{{\mathcal{ℛ}}_{ij}}=\left(\begin{array}{ccc}{\mathcal{ℛ}}_0^{ij}& ...& 0\\ ...& ...& ...\\ {\mathcal{ℛ}}_m^{ij}& ...& {\mathcal{ℛ}}_0^j\end{array}\right),T_{{\mathcal{ℒ}}_{ij}}=\left(\begin{array}{ccc}{\mathcal{ℒ}}_0^{ij}& ...& {\mathcal{ℒ}}_{-n}^{ij}\\ ...& ...& ...\\ 0& ...& {\mathcal{ℒ}}_0^{ij}\end{array}\right).$$  Then it is easy to see that $${\left(T\pm H\right)}^{†}=\frac{1}{2}[{\sum }_{j=1}^4{T}_{{\mathcal{ℛ}}_{1j}}{\pi }_j{T}_{{\mathcal{ℒ}}_{j2}}+{\sum }_{j=1}^4{H}_{{\mathcal{ℛ}}_{2j}}{\pi }_j{H}_{{\mathcal{ℒ}}_{j1}}\pm$$  $$\pm \left({\sum }_{j=1}^4{T}_{{\mathcal{ℛ}}_{1j}}{\pi }_j{H}_{{\mathcal{ℒ}}_{j1}}+{\sum }_{j=1}^4{H}_{{\mathcal{ℛ}}_{2j}}{\pi }_j{T}_{{\mathcal{ℒ}}_{j2}}\right)],$$  where we denote $T_{\mathcal{ℒ},\mathcal{ℛ}}J=H_{\mathcal{ℒ},\mathcal{ℛ}}$  and ${\pi }_j$  are the matrices constructed by $d_j$  with the same manner as $\Pi$  by $d$  .

4. An example

Let us demonstrate the theorem with an example. We will find generalized inversions of the following $T+H$  and $T-H$  matrices:

$$T\pm H=\left(\begin{array}{cccc}1& 0& -1& 1\\ 1& 1& 0& -1\\ 1& 1& 1& 0\\ -1& 1& 1& 1\end{array}\right)\pm \left(\begin{array}{cccc}1& 0& -1& 1\\ 0& -1& 1& 0\\ -1& 1& 0& 0\\ 1& 0& 0& 1\end{array}\right).$$  Our calculations show that the indices of the sequence $A_{-3},\dots ,A_3$  are equal to ${\mu }_1=-1,{\mu }_2=0,{\mu }_3=0,{\mu }_4=1$  . The matrix $\mathcal{ℛ}\left(z\right)$  of the right essential polynomials and the matrix $\mathcal{ℒ}\left(z\right)$  of the conforming left essential polynomials are:

$$\mathcal{ℛ}\left(z\right)=\left(\begin{array}{c}{\mathcal{ℛ}}^1\left(z\right)\\ {\mathcal{ℛ}}^2\left(z\right)\end{array}\right),\mathcal{ℒ}\left(z\right)=\left(\begin{array}{cc}{\mathcal{ℒ}}^1\left(z\right)& {\mathcal{ℒ}}^2\left(z\right)\end{array}\right),$$  where $${\mathcal{ℛ}}^1\left(z\right)=-\left(\begin{array}{cccc}1& 2& 2& 0\end{array}\right)+\left(\begin{array}{cccc}-4& 5& 3& -1\end{array}\right)z+\left(\begin{array}{cccc}2& 1& 0& 0\end{array}\right)z^2-$$  $$-\left(\begin{array}{cccc}11& 4& 2& 0\end{array}\right)z^3+\left(\begin{array}{cccc}0& 4& 1& 0\end{array}\right)z^4+\left(\begin{array}{cccc}0& 0& 0& 1\end{array}\right)z^5,$$  $${\mathcal{ℛ}}^2\left(z\right)=\left(\begin{array}{cccc}11& -4& -2& 2\end{array}\right)-\left(\begin{array}{cccc}2& 11& 5& -2\end{array}\right)z+\left(\begin{array}{cccc}4& 0& 0& 0\end{array}\right)z^2+$$  $$+\left(\begin{array}{cccc}1& 0& 0& 0\end{array}\right)z^3+\left(\begin{array}{cccc}0& 0& 1& 0\end{array}\right)z^4,$$  $${\mathcal{ℒ}}^1\left(z\right)=\frac{1}{180}[\left(\begin{array}{c}-4\\ -12\\ 4\\ 0\end{array}\right)+\left(\begin{array}{c}4\\ 0\\ 12\\ 0\end{array}\right)z^{-1}+\left(\begin{array}{c}-1\\ 6\\ -23\\ 0\end{array}\right)z^{-2}-$$  $$-\left(\begin{array}{c}21\\ -6\\ 35\\ 180\end{array}\right)z^{-3}+\left(\begin{array}{c}-18\\ -84\\ 76\\ 0\end{array}\right)z^{-4}+\left(\begin{array}{c}16\\ 0\\ 0\\ 0\end{array}\right)z^{-5}],$$  $${\mathcal{ℒ}}^2\left(z\right)=\frac{1}{180}[\left(\begin{array}{c}25\\ 30\\ -25\\ 180\end{array}\right)+\left(\begin{array}{c}26\\ 108\\ -126\\ 0\end{array}\right)z^{-1}+\left(\begin{array}{c}-11\\ 30\\ -43\\ 0\end{array}\right)z^{-2}+$$  $$+\left(\begin{array}{c}-6\\ -24\\ 32\\ 0\end{array}\right)z^{-3}-\left(\begin{array}{c}10\\ 24\\ 26\\ 0\end{array}\right)z^{-4}-\left(\begin{array}{c}4\\ 0\\ 0\\ 0\end{array}\right)z^{-5}].$$  We construct the matrix $$\Pi =\left(\begin{array}{cccc}{\Pi }_0& {\Pi }_{-1}& {\Pi }_{-2}& {\Pi }_{-3}\\ {\Pi }_1& {\Pi }_0& {\Pi }_{-1}& {\Pi }_{-2}\\ {\Pi }_2& {\Pi }_1& {\Pi }_0& {\Pi }_{-1}\\ {\Pi }_3& {\Pi }_2& {\Pi }_1& {\Pi }_0\end{array}\right),$$  where $${\Pi }_0=\left(\begin{array}{cccc}0& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 0\end{array}\right),{\Pi }_{-1}=\left(\begin{array}{cccc}0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 1\end{array}\right),$$  $${\Pi }_1=\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right),$$  and the matrices ${\Pi }_{\pm 2},{\Pi }_{\pm 3}$  are zero matrices.

Then we construct the matrices $T_{{\mathcal{ℛ}}_1}$  , $T_{{\mathcal{ℛ}}_2}$  with help of the coefficients of ${\mathcal{ℛ}}^1,{\mathcal{ℛ}}^2$  and matrices $T_{{\mathcal{ℒ}}_1},$  $T_{{\mathcal{ℒ}}_2}$  with help of the coefficients of ${\mathcal{ℒ}}^1,{\mathcal{ℒ}}^2$  .

We use the formulas ( 3.1 ) for generalized inverses of the $T\pm H$  matrices and have $${\left(T+H\right)}^{†}=\frac{1}{20}\left(\begin{array}{cccc}5& 10& 0& 0\\ 2& -4& 12& -4\\ -4& 8& -4& 8\\ 1& -2& -4& 8\end{array}\right),$$  $${\left(T-H\right)}^{†}=\frac{1}{180}\left(\begin{array}{cccc}-113& 2& 50& -32\\ 88& 8& 20& 52\\ -134& -4& 80& 64\\ 17& -158& 10& 16\end{array}\right).$$  One may conclude that ${\left(T+H\right)}^{†}$  is the inverse matrix for the $T+H$  matrix and ${\left(T-H\right)}^{†}$  is the generalized inverse for $T-H.$  References

1. A. H. Sayed, H. Lev-Ari, and T. Kailath, Fast Triangular Factorization of the Sum of Quasi-Toeplitz and Quasi-Hankel Matrices, Linear Algebra Appl. 191 : 77 106, 1993.
2. G. A. Merchant and T. W. Parks, Efficient solution of a Toeplitz-plus-Hankel coefficient system of equations, IEEE Transactions on Acoustics, Speech and Signal Processing, 30(1): 40 44, February, 1982.
3. A. B. Nersesian, A. A. Papoyan, Construction of matrix inversed to the sum of Toeplitz and Hankel matrices (Russian),Izv. AN Arm. SSR. Matematika, 8, N2 (1983), 441-463.
4. G. Heinig, K. Rost, Algebraic Methods for Toeplitz-like Matrices and Operators, Academie-Verlag, Berlin, 1984.
5. G. Heinig, K. Rost, On the Inverses of Toeplitz-plus-Hankel matrices, Linear Algebra Appl. 106 : 39 52, 1988.
6. I. Gohberg, T. Shalom, On inversion of square matrices partitioned into non-square blocks, Integral Equations and Operator Theory, 12 (1989), 539–566.
7. V. M. Adukov, O. L. Ibryaeva, On the kernel structure of the Toeplitz-plus-Hankel matrices (Russian), Vestnik Jujno-Uralskogo gosudarstvennogo universiteta, 7, 3-12, 2001.
8. V. M. Adukov, Generalized inversion of block Toeplitz matrices, Linear Algebra Appl. 274 : 85 124, 1998.