  + \sum_{k=1}^{K_P} \bar{\lambda}_k^\star\left( \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \otimes\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\right)\mathbf{w}_f^\star \right.   \nonumber\\
  &\left. + (\sigma_I^2+\sigma_n^2) \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T\otimes\mathbf{I}_{N_r}\right)\mathbf{w}_f^\star\right) = \bar{\mathbf{0}}_{N_tN_r}.\end{aligned}$$

By applying the properties of Kronecker product [@meyer2000matrix] and after some manipulations, we can obtain the optimal $\mathbf{w}_f^\star$ as $$\begin{aligned}
 \label{E:optimalwf2}
  & \mathbf{w}_f^\star  \!=\! \!-\!e^{j\phi} \big( \big(\sum_{k=1}^{K_P} \bar{\lambda}_k^\star(\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \!+\! \mathbf{I}_{N_t}\big)^{-\!1} \big(\sum_{k=1}^{K_P} \bar{\lambda}_k^\star\bar{\mathbf{h}}_{Mk} \mathbf{h}_{Pk}^H\big)^T  \nonumber\\
  & \otimes \!\big(\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\!+\! (\sigma_I^2 \!+\! \sigma_n^2)\mathbf{I}_{N_r}\big)^{-\!1}\bar{\mathbf{H}}_{MP}^H\big)  \mathtt{vec}(\mathbf{I}_{N_tN_r}).\end{aligned}$$

Then we can recover the optimal $\mathbf{W}_f^\star$ from $\mathbf{w}_f^\star$ as $$\begin{aligned}
 \label{E:Wf}
  &{\mathbf{W}_f^\star = } -e^{j\phi}\Big(\sum_{k=1}^{
+ \sum_{k=1}^{K_P } \bar{\lambda}_k^\star\left (\left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \otimes\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\right)\mathbf{w}_f^\star \right.    \nonumber\\ 
   & \left. + (\sigma_I^2+\sigma_n^2) \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T\otimes\mathbf{I}_{N_r}\right)\mathbf{w}_f^\star\right) = \bar{\mathbf{0}}_{N_tN_r}.\end{aligned}$$ 

 By applying the properties of Kronecker product   [ @meyer2000matrix ] and after some manipulations, we can receive the optimum $ \mathbf{w}_f^\star$ as $ $ \begin{aligned } 
  \label{E: optimalwf2 } 
   & \mathbf{w}_f^\star   \!=\! \!-\!e^{j\phi } \big (\big(\sum_{k=1}^{K_P } \bar{\lambda}_k^\star(\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \!+\! \mathbf{I}_{N_t}\big)^{-\!1 } \big(\sum_{k=1}^{K_P } \bar{\lambda}_k^\star\bar{\mathbf{h}}_{Mk } \mathbf{h}_{Pk}^H\big)^T   \nonumber\\ 
   & \otimes \!\big(\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\!+\! (\sigma_I^2 \!+\! \sigma_n^2)\mathbf{I}_{N_r}\big)^{-\!1}\bar{\mathbf{H}}_{MP}^H\big)   \mathtt{vec}(\mathbf{I}_{N_tN_r}).\end{aligned}$$ 

 Then we can recover the optimal $ \mathbf{W}_f^\star$ from $ \mathbf{w}_f^\star$ as $ $ \begin{aligned } 
  \label{E: Wf } 
   & { \mathbf{W}_f^\star = } -e^{j\phi}\Big(\sum_{k=1}^ {
  + \dum_{k=1}^{K_P} \bar{\lambda}_k^\star\lefu( \left((\mathbf{h}_{Pk}\majhvf{h}_{Pk}^H)^V \otimea\bar{\mathcf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\right)\mathbh{w}_f^\srar \rught.   \nonumber\\
  &\left. + (\sigma_I^2+\sihma_n^2) \lefr((\matibf{h}_{Pk}\mathbf{h}_{Pk}^H)^V\ktimes\mathbf{I}_{H_v}\righc)\methbf{w}_f^\star\righj) = \bar{\mathbf{0}}_{T_tN_r}.\end{aligned}$$

Ty a'plying the properties of Kronecker [roduct [@kejer2000matrix] and wfteg fome manipulations, we can obtain the oltimal $\kathbf{w}_f^\star$ ss $$\begin{aligned}
 \label{E:optilalwv2}
  & \mathbf{w}_f^\star  \!=\! \!-\!e^{j\phi} \big( \vig(\stn_{k=1}^{K_P} \bar{\lambaa}_k^\star(\mathbf{h}_{Pk}\mathbf{g}_{Pk}^H)^T \!+\! \mathbf{I}_{N_t}\big)^{-\!1} \big(\sum_{k=1}^{K_P} \bxr{\lamyda}_k^\star\bar{\narhbv{v}}_{Mk} \mathbf{h}_{'k}^H\big)^N  \nonumber\\
  & \otimes \!\tig(\bar{\msthbf{H}}_{MP}^H\bar{\mabhbf{H}}_{KP}\!+\! (\wigma_I^2 \!+\! \sigma_n^2)\mathbf{I}_{I_r}\big)^{-\!1}\bar{\mathbf{H}}_{MP}^H\bid)  \mathtt{eee}(\mathbf{I}_{N_tN_r}).\end{alignee}$$

Tyen wg can recucer tht o'tijal $\mahhbh{W}_f^\star$ froj $\mathbf{w}_f^\srar$ as $$\begin{aligned}
 \lwvel{E:Wf}
  &{\mathbr{W}_f^\staw = } -e^{j\phi}\Big(\sum_{k=1}^{
+ \sum_{k=1}^{K_P} \bar{\lambda}_k^\star\left( \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \otimes\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\right)\mathbf{w}_f^\star \right. \nonumber\\ (\sigma_I^2+\sigma_n^2) = \bar{\mathbf{0}}_{N_tN_r}.\end{aligned}$$ applying the properties after manipulations, we can the optimal $\mathbf{w}_f^\star$ $$\begin{aligned} \label{E:optimalwf2} & \mathbf{w}_f^\star \!=\! \!-\!e^{j\phi} \big(\sum_{k=1}^{K_P} \bar{\lambda}_k^\star(\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \!+\! \mathbf{I}_{N_t}\big)^{-\!1} \big(\sum_{k=1}^{K_P} \bar{\lambda}_k^\star\bar{\mathbf{h}}_{Mk} \mathbf{h}_{Pk}^H\big)^T \nonumber\\ & \otimes \!\big(\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\!+\! (\sigma_I^2 \!+\! \mathtt{vec}(\mathbf{I}_{N_tN_r}).\end{aligned}$$ Then we can recover the optimal $\mathbf{W}_f^\star$ from $\mathbf{w}_f^\star$ as $$\begin{aligned} \label{E:Wf} = -e^{j\phi}\Big(\sum_{k=1}^{
  + \sum_{k=1}^{K_P} \bar{\lambda}_k^\star\left( \lEft((\mathbf{h}_{pk}\matHbf{H}_{Pk}^h)^T \OtimEs\baR{\mathbf{H}}_{MP}^H\bar{\MAthbF{H}}_{MP}\right)\mathbf{w}_f^\star \rIght.   \nOnUMber\\
  &\LEfT. + (\sigmA_I^2+\sigma_N^2) \LeFT((\MatHbF{h}_{pk}\mAtHBf{H}_{Pk}^H)^T\OtiMes\mathBf{I}_{N_r}\right)\MatHbF{w}_f^\star\right) = \BAr{\Mathbf{0}}_{N_tN_r}.\End{Aligned}$$

By appLyiNg the pRoPerTIes of kroNeckeR produCT [@meyer2000Matrix] and AfTEr some MAnipulaTIOnS, we cAn obtain the optimaL $\MaTHbf{w}_f^\star$ as $$\begIn{aligNeD}
 \LaBEL{E:oPtiMalwf2}
  & \mathbF{w}_F^\star  \!=\! \!-\!E^{J\phi} \big( \BIg(\SUM_{K=1}^{K_P} \BAr{\lambda}_k^\star(\Mathbf{h}_{Pk}\maTHbf{H}_{Pk}^H)^T \!+\! \mAtHbf{i}_{n_t}\big)^{-\!1} \bIg(\sum_{K=1}^{K_p} \Bar{\Lambda}_k^\star\Bar{\mAthbf{h}}_{Mk} \mAthbf{h}_{pK}^H\big)^T  \nONumber\\
  & \oTimes \!\bIg(\bAr{\mAthbF{h}}_{Mp}^H\Bar{\MaTHbf{h}}_{mP}\!+\! (\SigMA_I^2 \!+\! \sIgma_n^2)\matHbF{I}_{n_r}\big)^{-\!1}\Bar{\mATHBF{H}}_{MP}^h\biG)  \matHtt{veC}(\mathbf{I}_{N_tN_r}).\eNd{aLignED}$$

ThEn we cAn recOver ThE optiMal $\matHbf{W}_f^\StAr$ from $\mathbf{w}_f^\sTar$ aS $$\begin{aliGneD}
 \lAbeL{E:wf}
  &{\matHBf{W}_f^\stAr = } -e^{J\phI}\Big(\sum_{K=1}^{
  + \sum_{k=1}^{K_P} \bar{ \lambda}_k ^\sta r\l eft (\lef t((\ mathbf{h}_{Pk} \ math bf{h}_{Pk}^H)^T \otime s\bar {\ m athb f {H }}_{M P}^H\ba r {\ m a thb f{ H} }_{ MP } \r ight) \ma thbf{w} _f^\star \ rig ht .   \nonumbe r \\ 
  &\left.  +(\sigma_I^2+ \si gma_n^ 2)  \l e ft((\ mat hbf{h }_{Pk} \ mathbf {h}_{Pk}^ H) ^ T\otim e s\mathb f { I} _{N_ r}\right)\mathbf{ w }_ f ^\star\right)= \bar {\ m at h b f{0 }}_ {N_tN_r}.\ en d{ali g ned}$$By a p ply i ng the proper ties of Kro n eck er pro du ct[ @meyer 2000m at r ix]  and aftersome  manipula tions, we cano btain t he opt ima l $ \mat h bf {w }_f ^\ s tar $  a s $ $ \be gin{alig ne d} 
 \la bel{ E : o p tima lwf 2}
 & \m athbf{w}_f^\s tar   \! = \!\!-\! e^{j\ phi}  \ big(\big(\ sum_{ k= 1}^{K_P} \bar{\ lamb da}_k^\st ar( \m ath bf {h}_{ P k}\mat hbf {h} _{Pk}^H )^T \!+ \ ! \ ma t h b f{ I}_{N_t}\big)^{-\! 1} \ bi g(\sum_{ k=1}^{ K _P }\ bar{\lam bd a}_ k^\s t a r\bar {\ma t hb f{h}}_{M k} \ma t hb f{ h}_{Pk} ^H \big)^ T \n onu mber\ \ 
  &  \otim es \!\bi g(\ba r {\mathbf{H}}_{ M P}^H\bar{\mat h bf { H }} _ {MP} \!+ \! (\sigma_ I^2\ !+\!  \si g ma _n^ 2 )\mat hbf{I }_ { N_ r }\big)^{-\!1}\bar{\ ma thbf{H }}_{M P}^H\big)  \m athtt{vec} ( \ m athbf{I} _{N_ t N_ r }).\end{aligne d}$$
Then we c a n recove r the  optimal  $\mathbf { W }_f^\sta r$fro m $ \ma t h bf {w}_f^\star$a s  $$\ be gin{ali gne d}
 \la bel {E: Wf} 
 &{ \mathbf{W }_f^\sta r=}-e ^{j \phi} \ Big(\sum _{ k=1 }^ {
 _+ \sum_{k=1}^{K_P}_\bar{\lambda}_k^\star\left( \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T \otimes\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\right)\mathbf{w}_f^\star \right._ _\nonumber\\
_ &\left._+_(\sigma_I^2+\sigma_n^2) \left((\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T\otimes\mathbf{I}_{N_r}\right)\mathbf{w}_f^\star\right) =_\bar{\mathbf{0}}_{N_tN_r}.\end{aligned}$$

By applying the_properties of Kronecker product [@meyer2000matrix]_and after some_manipulations,_we can obtain the optimal $\mathbf{w}_f^\star$ as $$\begin{aligned}
 \label{E:optimalwf2}
  & \mathbf{w}_f^\star  \!=\!_\!-\!e^{j\phi}_\big( \big(\sum_{k=1}^{K_P}_\bar{\lambda}_k^\star(\mathbf{h}_{Pk}\mathbf{h}_{Pk}^H)^T_\!+\!_\mathbf{I}_{N_t}\big)^{-\!1} \big(\sum_{k=1}^{K_P} \bar{\lambda}_k^\star\bar{\mathbf{h}}_{Mk} \mathbf{h}_{Pk}^H\big)^T _\nonumber\\
  & \otimes \!\big(\bar{\mathbf{H}}_{MP}^H\bar{\mathbf{H}}_{MP}\!+\!_(\sigma_I^2 \!+\!_\sigma_n^2)\mathbf{I}_{N_r}\big)^{-\!1}\bar{\mathbf{H}}_{MP}^H\big)  \mathtt{vec}(\mathbf{I}_{N_tN_r}).\end{aligned}$$

Then we can recover the optimal_$\mathbf{W}_f^\star$_from $\mathbf{w}_f^\star$ as_$$\begin{aligned}
 \label{E:Wf}
  &{\mathbf{W}_f^\star = } -e^{j\phi}\Big(\sum_{k=1}^{