textbf{R}_2^{(t)}\textbf{R}_2^{(t)^H}) \}
\label{Lag2}
\end{split}$$ where $$\begin{split}
& \textbf{W}_{r,1}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{H}_{2,r}^{(t-1)^H}, \\
& \textbf{W}_{r,2}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{H}_{1,r}^{(t-1)^H}, \\
& \textbf{W}_{r,3}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{G}_1^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,1}^{(t)^H}+
(N_sp_1\sigma_{e,1}^2+\sigma_{n,1}^2)\textbf{I}_{N_s}, \\
& \textbf{W}_{r,4}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{G}_2^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,2}^{(t)^H}+
(N_sp_2\sigma_{e,2}^2+\sigma_{n,2}^2)\textbf{I}_{N_s}.
\end{split}$$ By solving $\frac{\partial J(\textbf{R}_1^{(t)},\textbf{R}_2^{(t)})} {\partial \textbf{R}_i^{(t)}}=\textbf{0}_{N_s}$ for $i \in \{ 1,2 \}$, the receive beamforming matrices at the sources are given by $$\begin{split}
\textbf{R}_l^{(t)}
= \alpha^{(t)}p_{\bar{l}}
\{\textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{G}_l^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,l}^{(t)^H}
+(N_sp_l\sigma_{e,l}^2+\sigma_{n,l}^2)\textbf{I}_{N_s}\}^{-1}
\textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{H}_{\bar{l
textbf{R}_2^{(t)}\textbf{R}_2^{(t)^H }) \ } 
 \label{Lag2 } 
 \end{split}$$ where $ $ \begin{split } 
 & \textbf{W}_{r,1}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{H}_{2,r}^{(t-1)^H }, \\ 
 & \textbf{W}_{r,2}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{H}_{1,r}^{(t-1)^H }, \\ 
 & \textbf{W}_{r,3}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{G}_1^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,1}^{(t)^H}+ 
 (N_sp_1\sigma_{e,1}^2+\sigma_{n,1}^2)\textbf{I}_{N_s }, \\ 
 & \textbf{W}_{r,4}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{G}_2^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,2}^{(t)^H}+ 
 (N_sp_2\sigma_{e,2}^2+\sigma_{n,2}^2)\textbf{I}_{N_s }. 
 \end{split}$$ By solving $ \frac{\partial J(\textbf{R}_1^{(t)},\textbf{R}_2^{(t) }) } { \partial \textbf{R}_i^{(t)}}=\textbf{0}_{N_s}$ for $ i \in \ { 1,2 \}$, the receive beamforming matrices at the sources are give by $ $ \begin{split } 
 \textbf{R}_l^{(t) } 
 = \alpha^{(t)}p_{\bar{l } } 
 \{\textbf{H}_{r, l}^{(t)}\textbf{F}^{(t)}\textbf{G}_l^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r, l}^{(t)^H } 
 + (N_sp_l\sigma_{e, l}^2+\sigma_{n, l}^2)\textbf{I}_{N_s}\}^{-1 } 
 \textbf{H}_{r, l}^{(t)}\textbf{F}^{(t)}\textbf{H}_{\bar{l
texhbf{R}_2^{(t)}\textbf{R}_2^{(t)^H}) \}
\label{Lag2}
\ekd{split}$$ where $$\beyun{spliv}
& \textbr{W}_{r,1}^{(t)}=\textcf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{H}_{2,r}^{(t-1)^H}, \\
& \texvbf{W}_{e,2}^{(t)}=\texuyf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{H}_{1,f}^{(t-1)^H}, \\
& \texthf{W}_{r,3}^{(t)}=\texrbf{H}_{c,1}^{(t)}\textbf{F}^{(t)}\textbf{J}_1^{(f)}\textbf{N}^{(c)^H}\texfnf{H}_{r,1}^{(t)^K}+
(N_wp_1\sigma_{e,1}^2+\sigma_{n,1}^2)\jextbf{I}_{N_s}, \\
& \taxtbf{W}_{r,4}^{(t)}=\textbf{V}_{r,2}^{(g)}\tzxtbf{F}^{(t)}\textbf{G}_2^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,2}^{(t)^H}+
(N_s[_2\sigma_{e,2}^2+\xihma_{n,2}^2)\textbf{I}_{N_s}.
\egd{spkyt}$$ Bg solving $\frac{\partial J(\textbf{R}_1^{(t)},\texfbf{R}_2^{(t)})} {\pertial \textbf{R}_i^{(y)}}=\textbf{0}_{N_s}$ for $i \in \{ 1,2 \}$, the gecelve beamforming mahrices at tye sjyrces are gixen by $$\begpu{split}
\textbr{R}_l^{(t)}
= \alpha^{(t)}p_{\bar{l}}
\{\textbf{H}_{r,l}^{(t)}\textbw{F}^{(t)}\tertbf{G}_l^{(t)}\textvf{D}^{(t)^H}\hgxtbf{H}_{r,l}^{(t)^H}
+(N_s'_l\sigmw_{e,l}^2+\sigma_{n,l}^2)\tewnbf{I}_{N_s}\}^{-1}
\taxtbf{H}_{r,k}^{(t)}\textbf{F}^{(t)}\textnf{H}_{\bac{l
textbf{R}_2^{(t)}\textbf{R}_2^{(t)^H}) \} \label{Lag2} \end{split}$$ where $$\begin{split} & & \\ & (N_sp_1\sigma_{e,1}^2+\sigma_{n,1}^2)\textbf{I}_{N_s}, \\ & $\frac{\partial {\partial \textbf{R}_i^{(t)}}=\textbf{0}_{N_s}$ for \in \{ 1,2 the receive beamforming matrices at the are given by $$\begin{split} \textbf{R}_l^{(t)} = \alpha^{(t)}p_{\bar{l}} \{\textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{G}_l^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,l}^{(t)^H} +(N_sp_l\sigma_{e,l}^2+\sigma_{n,l}^2)\textbf{I}_{N_s}\}^{-1} \textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{H}_{\bar{l
textbf{R}_2^{(t)}\textbf{R}_2^{(t)^H}) \}
\label{Lag2}
\End{split}$$ whEre $$\beGin{SplIt}
& \TextBf{W}_{r,1}^{(T)}=\textbf{H}_{r,1}^{(t)}\textBF{F}^{(t)}\tExtbf{H}_{2,r}^{(t-1)^H}, \\
& \textbf{W}_{r,2}^{(t)}=\textBf{H}_{r,2}^{(t)}\TeXTbf{F}^{(T)}\TeXtbf{H}_{1,R}^{(t-1)^H}, \\
& \textBF{W}_{R,3}^{(T)}=\TexTbF{H}_{R,1}^{(t)}\tExTBf{f}^{(t)}\texTbf{g}_1^{(t)}\textbF{F}^{(t)^H}\textbf{h}_{r,1}^{(t)^h}+
(N_Sp_1\sigma_{e,1}^2+\sigmA_{N,1}^2)\tExtbf{I}_{N_s}, \\
& \teXtbF{W}_{r,4}^{(t)}=\textbf{H}_{r,2}^{(T)}\teXtbf{F}^{(t)}\TeXtbF{g}_2^{(t)}\texTbf{f}^{(t)^H}\teXtbf{H}_{r,2}^{(T)^h}+
(N_sp_2\siGma_{e,2}^2+\sigma_{N,2}^2)\tEXtbf{I}_{N_S}.
\End{spliT}$$ bY sOlviNg $\frac{\partial J(\texTBf{r}_1^{(T)},\textbf{R}_2^{(t)})} {\partiAl \textBf{r}_I^{(t)}}=\TEXtbF{0}_{N_s}$ For $i \in \{ 1,2 \}$, the rEcEive bEAmformiNG mATRIceS At the sources aRe given by $$\beGIn{sPlit}
\teXtBf{R}_L^{(T)}
= \alpha^{(T)}p_{\bar{L}}
\{\tEXtbF{H}_{r,l}^{(t)}\textbf{f}^{(t)}\teXtbf{G}_l^{(t)}\teXtbf{F}^{(t)^h}\Textbf{H}_{R,L}^{(t)^H}
+(N_sp_l\Sigma_{e,L}^2+\siGma_{N,l}^2)\teXTbF{I}_{n_s}\}^{-1}
\tExTBf{H}_{R,L}^{(t)}\TexTBf{F}^{(T)}\textbf{H}_{\BaR{l
textbf{R}_2^{(t)}\textbf{R }_2^{(t)^H }) \} 
\l abe l{ Lag2 }
\e nd{split}$$ wh e re $ $\begin{split}
& \text bf{W} _{ r ,1}^ { (t )}=\t extbf{H } _{ r , 1}^ {( t) }\t ex t bf {F}^{ (t) }\textb f{H}_{2,r} ^{( t- 1)^H}, \\
&\ te xtbf{W}_{r ,2} ^{(t)}=\text bf{ H}_{r, 2} ^{( t )}\te xtb f{F}^ {(t)}\ t extbf{ H}_{1,r}^ {( t -1)^H} ,  \\
& \ t e xt bf{W }_{r,3}^{(t)}=\te x tb f {H}_{r,1}^{(t) }\text bf { F} ^ { (t) }\t extbf{G}_1 ^{ (t)}\ t extbf{F } ^{ ( t ) ^H} \ textbf{H}_{r, 1}^{(t)^H}+ (N_ sp_1\s ig ma_ { e,1}^2 +\sig ma _ {n, 1}^2)\textb f{I} _{N_s}, \ \
& \t e xtbf{W} _ {r,4}^{ (t)}=\ tex tbf {H}_ { r, 2} ^{( t) } \te x tb f{F } ^{( t)}\text bf {G }_2^{ (t)} \ t e x tbf{ F}^ {(t) ^H}\t extbf{H}_{r,2 }^{ (t)^ H }+(N_sp _2\si gma_ {e ,2}^2 +\sigm a_{n, 2} ^2)\textbf{I}_{ N_s} .
\end{sp lit }$ $ B ysolvi n g $\fr ac{ \pa rtial J (\textb f {R} _1 ^ { ( t) },\textbf{R}_2^{(t )} ) }  { \partial  \text b f{ R} _ i^{(t)}} =\ tex tbf{ 0 } _{N_s }$ f o r$i \in \ { 1,2\ }$ ,the rec ei ve bea mf orm ing  matr i cesat the  sources  areg iven by $$\beg i n{split}
\tex t bf { R }_ l ^{(t )}= \alpha^{( t)}p _ {\ba r{l} } 
\ {\t e xtbf{ H}_{r ,l } ^{ ( t)}\textbf{F}^{(t)} \t extbf{ G}_l^ {(t)}\textbf{ F}^{(t)^H} \ t e xtbf{H}_ {r,l } ^{ ( t)^H}
+(N_sp_l \sigm a_{e,l}^2+ \ sigma_{n ,l}^2 )\textbf {I}_{N_s} \ } ^{-1}
\t ext bf{ H}_ {r, l } ^{ (t)}\textbf{F } ^ {(t) }\ textbf{ H}_ {\bar{l 
textbf{R}_2^{(t)}\textbf{R}_2^{(t)^H}) \}
\label{Lag2}
\end{split}$$_where $$\begin{split}
&_\textbf{W}_{r,1}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{H}_{2,r}^{(t-1)^H}, \\
& \textbf{W}_{r,2}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{H}_{1,r}^{(t-1)^H}, \\
&_\textbf{W}_{r,3}^{(t)}=\textbf{H}_{r,1}^{(t)}\textbf{F}^{(t)}\textbf{G}_1^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,1}^{(t)^H}+
(N_sp_1\sigma_{e,1}^2+\sigma_{n,1}^2)\textbf{I}_{N_s}, \\
&_\textbf{W}_{r,4}^{(t)}=\textbf{H}_{r,2}^{(t)}\textbf{F}^{(t)}\textbf{G}_2^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,2}^{(t)^H}+
(N_sp_2\sigma_{e,2}^2+\sigma_{n,2}^2)\textbf{I}_{N_s}.
\end{split}$$_By solving_$\frac{\partial_J(\textbf{R}_1^{(t)},\textbf{R}_2^{(t)})} {\partial \textbf{R}_i^{(t)}}=\textbf{0}_{N_s}$_for $i \in_\{ 1,2 \}$, the_receive beamforming matrices_at_the sources are given by $$\begin{split}
\textbf{R}_l^{(t)}
= \alpha^{(t)}p_{\bar{l}}
\{\textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{G}_l^{(t)}\textbf{F}^{(t)^H}\textbf{H}_{r,l}^{(t)^H}
+(N_sp_l\sigma_{e,l}^2+\sigma_{n,l}^2)\textbf{I}_{N_s}\}^{-1}
\textbf{H}_{r,l}^{(t)}\textbf{F}^{(t)}\textbf{H}_{\bar{l