mathbf P_k} (\bm \psi_k, \bm \tau_k)$ is a tall matrix with asymptotically mutually-orthogonal columns of unit length, $$\begin{aligned}
{\mathbf R}_k^U =
\bigg( \frac{1}{\sqrt{MP}}
{\mathbf P_k} (\bm \psi_k, \bm \tau_k) \bigg)
\bigg( MP \bm \Lambda_k \bigg)
\bigg( \frac{1}{\sqrt{MP}}
{\mathbf P}_k^H (\bm \psi_k, \bm \tau_k) \bigg)
\label{Rkue}\end{aligned}$$ provides a good approximation of eigenvalue decomposition.

With the aid of the *AoA-delay reciprocity* [@SFW-TSP; @SFW-CM], downlink channel covariance matrices can also be reconstructed via the physical parameters of uplink channels, which will be discussed in Section \[sec:dl\].

Initial Uplink Channel Parameter Extraction {#sec:initial}
===========================================

At the very beginning, all users stay in the dark as far as the BS is concerned and thus orthogonal trainings have to be applied to avoid the inter-user interference and pilot contamination at the BS. In this case, we operate with the frequency orthogonality among difference users [@OFDMA-TWC; @SFW-TSP] by setting $\mathcal P_k \cap \mathcal P_r = \varnothing, $ $\forall k\neq r$, at this very beginning phase. During this phase, the initial AoA, time delay, and complex gain of each path are estimated for all users and we call this phase *initial parameter extraction*.

In this section, we introduce a parameter extraction algorithm for this stage, which suffices for the subsequent multi-user uplink and downlink channel estimations, as will be shown in Sections V and VI.

Problem Formulation
-------------------

For the $k$th user, the uplink pilots transmission process is given in. Substituting into, we have $$\begin{aligned}
\mathbf y_k =
\mathbf W^H {\mathbf P}_k(\bm \psi_k, \bm \tau_k) \bm \alpha_k + \mathbf n_k
.
\label{yk2}\end{aligned}$$
mathbf P_k } (\bm \psi_k, \bm \tau_k)$ is a tall matrix with asymptotically mutually - orthogonal columns of unit duration, $ $ \begin{aligned } 
 { \mathbf R}_k^U = 
 \bigg (\frac{1}{\sqrt{MP } } 
 { \mathbf P_k } (\bm \psi_k, \bm \tau_k) \bigg) 
 \bigg (MP \bm \Lambda_k \bigg) 
 \bigg (\frac{1}{\sqrt{MP } } 
 { \mathbf P}_k^H (\bm \psi_k, \bm \tau_k) \bigg) 
 \label{Rkue}\end{aligned}$$ provide a good approximation of eigenvalue decomposition. 

 With the care of the * AoA - delay reciprocity * [ @SFW - TSP; @SFW - CM ], downlink channel covariance matrices can besides be restore via the physical parameter of uplink channels, which will be discussed in Section   \[sec: dl\ ]. 

 Initial Uplink Channel Parameter Extraction { # sec: initial } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 At the identical beginning, all users stay in the dark as far as the BS is implicated and thus orthogonal training have to be applied to debar the inter - user interference and pilot contamination at the boron. In this case, we operate with the frequency orthogonality among difference users [ @OFDMA - TWC; @SFW - TSP ] by set $ \mathcal P_k \cap \mathcal P_r = \varnothing, $ $ \forall k\neq r$, at this very beginning phase. During this phase, the initial AoA, time delay, and complex gain of each path are estimated for all drug user and we call this phase * initial parameter extraction *. 

 In this section, we introduce a parameter extraction algorithm for this stage, which suffices for the subsequent multi - user uplink and downlink channel estimations, as will be shown in Sections   V and VI. 

 Problem Formulation 
 ------------------- 

 For the $ k$th user, the uplink pilots transmission process is given in. substitute into, we have $ $ \begin{aligned } 
 \mathbf y_k = 
 \mathbf W^H { \mathbf P}_k(\bm \psi_k, \bm \tau_k) \bm \alpha_k + \mathbf n_k 
. 
 \label{yk2}\end{aligned}$$
matjbf P_k} (\bm \psi_k, \bm \tau_k)$ ir a tall matrix with avymptofically outually-orthogonal columns oh unut lebgth, $$\begin{aligned}
{\mathbw R}_k^U =
\bigh( \frac{1}{\sqet{MP}}
{\nqthbf P_k} (\bm \psi_k, \bm \tau_k) \nigg)
\bngj( MP \bm \Lambda_k \bigg)
\bigg( \xrac{1}{\sqrt{MP}}
{\mathtf P}_n^H (\bm \psi_k, \bm \tau_k) \bigg)
\label{Rkue}\end{wligned}$$ pgovides a good appgovimafpok of eigenvalue decomposition.

Witg the apd of the *AoA-delau reciprocity* [@SFW-TSP; @SFW-CM], dowjlink channel covagiance matruces xan also be feconstrucueb via the pgysical parameters of uplink chxnnelx, which wiol be giscussed ii Sectpon \[sec:dl\].

Initial Uplinn Channrl Parameter Ewtracvion {#sec:initial}
===========================================

At the verb beginning, all userf stay in tke dark as far as the BW is woncarnea ana tgux krthoglnam traininga have to bw applied to avoid uhe pmter-user intsrferegcq and pilot contamination at the BS. In uhis dase, we operate with thw frequency orthogonapity amond difference users [@OFDMA-TWC; @SFW-TSP] by setting $\matvcal '_k \cak \nathcxo O_r = \varnothing, $ $\forall k\neq r$, at this very bedjnmikg phase. During bhis phase, the iniyiwl WoA, time delai, and compmex gain of each pwth are estinated for all users and we call this phawe *initial pcraneter extraction*.

In this sectiun, wg intrpduce a parameter extraetion zlgorithm flr this sfxge, which sufficds nor the subsequent multi-user tplink anv dowulink chxnnek estiiations, as will be shown in Sectiojs V aud VI.

[roblem Fogmulation
-------------------

For the $k$th user, the u'kink pilots jratsmpssion prjcess is given in. Fubstituting iuto, we hcve $$\beein{aligned}
\jathbf b_k =
\mathbf W^H {\mathbf P}_k(\bm \[di_k, \bm \tau_k) \um \alpha_k + \marhbf n_k
.
\labeu{hk2}\end{aligned}$$
mathbf P_k} (\bm \psi_k, \bm \tau_k)$ is matrix asymptotically mutually-orthogonal of unit length, \frac{1}{\sqrt{MP}} P_k} (\bm \psi_k, \tau_k) \bigg) \bigg( \bm \Lambda_k \bigg) \bigg( \frac{1}{\sqrt{MP}} {\mathbf (\bm \psi_k, \bm \tau_k) \bigg) \label{Rkue}\end{aligned}$$ provides a good approximation of eigenvalue decomposition. the aid of the *AoA-delay reciprocity* [@SFW-TSP; @SFW-CM], downlink channel covariance matrices can be via physical of uplink channels, which will be discussed in Section \[sec:dl\]. Initial Uplink Channel Parameter Extraction {#sec:initial} At the very beginning, all users stay in dark as far as BS is concerned and thus trainings to be to the interference and pilot at the BS. In this case, we operate with the frequency orthogonality among difference users [@OFDMA-TWC; @SFW-TSP] setting $\mathcal \mathcal P_r \varnothing, $\forall r$, at this phase. During this phase, the initial and complex gain of each path are estimated all users we call this phase *initial parameter In this section, we introduce a parameter extraction for this stage, which suffices for the subsequent multi-user uplink and downlink channel estimations, as shown in Sections V VI. Problem Formulation For $k$th the pilots transmission is given in. Substituting into, we have $$\begin{aligned} \mathbf y_k = W^H {\mathbf P}_k(\bm \psi_k, \bm \tau_k) \bm \alpha_k + \mathbf \label{yk2}\end{aligned}$$
mathbf P_k} (\bm \psi_k, \bm \tau_k)$ is a taLl matrix wiTh asyMptOtiCaLly mUtuaLly-orthogonal cOLumnS of unit length, $$\begin{aligNed}
{\maThBF R}_k^U =
\BIgG( \frac{1}{\Sqrt{MP}}
{\mAThBF p_k} (\bM \pSi_K, \bm \TaU_K) \bIgg)
\biGg( Mp \bm \LambDa_k \bigg)
\bigG( \frAc{1}{\Sqrt{MP}}
{\mathbf p}_K^H (\Bm \psi_k, \bm \taU_k) \bIgg)
\label{Rkue}\End{AligneD}$$ pRovIDes a gOod ApproXimatiON of eigEnvalue deCoMPositiON.

With thE AId Of thE *AoA-delay reciprocITy* [@sfW-TSP; @SFW-CM], downLink chAnNEl COVarIanCe matrices CaN also BE reconsTRuCTED viA The physical paRameters of uPLinK channElS, whICh will Be disCuSSed In Section \[seC:dl\].

INitial UplInk ChaNNel ParaMEter ExtRactioN {#seC:inItiaL}
===========================================

at ThE veRy BEgiNNiNg, aLL usErs stay iN tHe Dark aS far AS THE BS iS coNcerNed anD thus orthogonAl tRainINgs Have tO be apPlieD tO avoiD the inTer-usEr Interference and PiloT contaminAtiOn At tHe bS. In tHIs case, We oPerAte with The freqUEncY oRTHOgOnality among differEnCE UsErs [@OFDMA-tWC; @SFW-tsP] By SEtting $\maThCal p_k \caP \MAthcaL P_r = \vARnOthing, $ $\foRall k\nEQ r$, At This verY bEginniNg PhaSe. DUring THis pHase, thE initial aoA, tiME delay, and complEX gain of each paTH aRE EsTImatEd fOr all users aNd we CAll tHis pHAsE *inITial pArameTeR ExTRaction*.

In this sectioN, wE introDuce a Parameter extrAction algoRITHm for thiS staGE, wHIch suffices for The suBsequent muLTi-user upLink aNd downliNk channel ESTimationS, as WilL be ShoWN In sections V and Vi.

pRoblEm formulaTioN
-------------------

For the $K$th UseR, thE upLiNk pilots tRansmissIoN pRoCeSs iS giveN In. SubstiTuTinG iNto, We havE $$\Begin{aLigneD}
\matHbF y_K =
\MatHbf W^H {\maTHbF p}_K(\bm \pSi_K, \bM \tau_K) \bm \AlPha_k + \mAthbF N_k
.
\lAbel{yk2}\eNd{aligned}$$
mathbf P_k} (\bm \psi_k, \ bm \tau_k) $ isa t all  m atri x wi th asymptotica l ly m utually-orthogonal col umnsof unit le ngth,  $$\beg i n{ a l ign ed }{\m at h bf  R}_k ^U=
\bigg ( \frac{1} {\s qr t{MP}}
{\mat h bf  P_k} (\bm  \p si_k, \bm \t au_ k) \bi gg )
\ b igg(MP\bm \ Lambda _ k \big g)
\bigg(  \ f rac{1} { \sqrt{M P } }{\ma thbf P}_k^H (\bm\ ps i _k, \bm \tau_k ) \big g) \l a b el{ Rku e}\end{ali gn ed}$$ provide s  a g o oda pproximationof eigenval u e d ecompo si tio n .

Wit h the  a i d o f the *AoA- dela y recipro city*[ @SFW-TS P ; @SFW- CM], d own lin k ch a nn el  co va r ian c emat r ice s can al so  b e rec onst r u c t ed v iathephysi cal parameter s o f up l ink  chan nels,  whi ch  will  be di scuss ed  in Section \[s ec:d l\].

Ini tia lUpl in k Cha n nel Pa ram ete r Extra ction { # sec :i n i t ia l}
=============== == = = == ======== ====== = == == = ==

At t he  ve ry b e g innin g, a l lusers st ay int he  d ark asfa r as t he  BS  is  conc e rned  and t hus orth ogona l  trainings hav e  to be applie d  t o av o id t heinter-userinte r fere ncea nd  pi l ot co ntami na t io n  at the BS. In this  c ase, w e ope rate with the  frequency o r thogonal itya mo n g difference u sers[@OFDMA-TW C ; @SFW-T SP] b y settin g $\mathc a l  P_k \ca p \ mat hca l P _ r  =  \varnothing, $  $\f or all k\n eqr$, atthi s v ery  be gi nning pha se. Duri ng  t hi spha se, t h e initia lAoA ,tim e del a y, and  comp lexga in ofeach pa t ha r e es ti ma tedfor  a ll us ersa ndwe call  this pha se* init ia lparamet er extraction *. 

In thisse cti on, we i ntroduce  a parameter extraction algorit hmfor t hisstage, wh ich  suffi ces for th e subs equen tmul t i -user u pl ink  a nd downlin k cha nneles tima tions,as will be shown i n  Se ctions V andVI. 

Pr o b le m F o rm u lat io n 
-- - - --------------- 

For the$k $ th  user, the upl in k pilot s trans missi o n proce ss is giv en in. Su bs titu t i nginto, we h ave $$\b egin{alig n ed}
\ m at hbf y _k=
\mat hb f W ^H {\ mathbf P}_ k(\bm  \psi_ k,  \bm \ tau_k )\bm \alp ha_k + \mathbf n_k
.
\l abel{y k2}\e nd{ aligned}$ $
mathbf P_k}_(\bm \psi_k,_\bm \tau_k)$ is a_tall matrix_with_asymptotically mutually-orthogonal_columns_of unit length,_$$\begin{aligned}
{\mathbf R}_k^U =
\bigg(_\frac{1}{\sqrt{MP}}
{\mathbf P_k} (\bm \psi_k,_\bm \tau_k) \bigg)
\bigg(_MP_\bm \Lambda_k \bigg)
\bigg( \frac{1}{\sqrt{MP}}
{\mathbf P}_k^H (\bm \psi_k, \bm \tau_k) \bigg)
\label{Rkue}\end{aligned}$$ provides a good approximation_of_eigenvalue decomposition.

With_the_aid_of the *AoA-delay reciprocity* [@SFW-TSP;_@SFW-CM], downlink channel covariance matrices_can also_be reconstructed via the physical parameters of uplink_channels,_which will be_discussed in Section \[sec:dl\].

Initial Uplink Channel Parameter Extraction {#sec:initial}
===========================================

At the_very beginning, all users stay in_the dark as_far_as_the BS is concerned_and thus orthogonal trainings have to_be applied to avoid the inter-user_interference and pilot contamination at the BS._In this case, we operate with_the frequency orthogonality among difference_users [@OFDMA-TWC;_@SFW-TSP] by setting $\mathcal P_k_\cap \mathcal P_r_= \varnothing,_$ $\forall k\neq_r$, at this very beginning phase._During this phase,_the initial AoA, time delay, and_complex_gain of each_path_are_estimated for_all users and_we_call this_phase_*initial parameter extraction*.

In this section, we_introduce_a parameter extraction algorithm for this stage,_which suffices for the_subsequent_multi-user uplink and downlink_channel estimations, as will be_shown in Sections V and VI.

Problem Formulation
-------------------

For_the $k$th_user, the_uplink pilots transmission process is given in. Substituting into, we have_$$\begin{aligned}
\mathbf y_k =
\mathbf W^H {\mathbf P}_k(\bm_\psi_k, \bm \tau_k) \bm_\alpha_k +_\mathbf_n_k
.
\label{yk2}\end{aligned}$$