Digital Beamforming - Channel Matrix Singular Value Decomposition

Consider a point-to-point MIMO beamforming system, assuming that the transmitter has N t N_tNtAntennas, N r N_rNrThe root receiving antenna, the transmitting end sends the modulated complex number information vector as xxx , after the transmit beamforming matrixFFAfter the function of F , it is loaded on each array element of the transmitting antenna, and passes through the MIMO channelHHThe received signal vector after H
is: y = HF x + ny=HFx+ny=HFx+nwhere
, nnn is complex Gaussian white noise, the mean isθ \thetaθ , the variance isσ 2 \sigma^2p2 .
At the receiver using the receive beamforming matrixWWW performs a weighted combination operation on the received signals:
WH y = ​​WHHF x + WH n W^Hy=W^HHFx+W^HnWHy=WHHFx+WH nFor
the convenience of analysis, the beamforming matrixWWW andFFF normalization, ieWHW = FHF = IW^HW=F^HF=IWHW=FHF=I , using the criterion of maximizing the received signal-to-noise ratio to obtain the beamforming matrix. In order to maximize the receiving SNR, the transmitting beamforming matrix should be obtained according to the optimization criterion of the following formula:
F t = argmax ⁡ F [ ( HF ) HHF ] \mathbf{F}_{t}=\underset{\boldsymbol{ F}}{\operatorname{argmax}}\left[(\mathbf{H} \mathbf{F})^{H} \mathbf{H} \mathbf{F}\right]Ft=Fargmax[(HF)H HF]
According to the KKT condition in convex optimization theory,F t F_tFtThe columns are represented by HHHH^HHHH HN s N_sNsThe eigenvectors corresponding to the largest eigenvalues ​​are formed, where N s N_sNsis the number of data streams. Similarly, at the receiver, the receive beamforming matrix W r W_rWrThe columns are represented by HHHH^HHHH HN s N_sNsThe eigenvectors corresponding to the largest eigenvalues ​​are formed. According to matrix theory, F t F_tFtGive W r W_rWrThey are the largest N s N_s in the right singular matrix and the left singular matrix obtained after Singular Value Decomposition (SVD) of the channel matrixNsSingular vectors corresponding to singular values.
For channel matrix HHH is decomposed by SVD to get:
H = U Λ VHH=U\Lambda V^HH=UΛVH
使用,U = [ u 1 , u 2 , . . . , u N r ] U=[u_1,u_2,...,u_{N_r}]U=[u1,u2,...,uNr] isN r N_rNrOrder unitary matrix, V = [ v 1 , v 2 , . . . , v N t ] V=[v_1,v_2,...,v_{N_t}]V=[v1,v2,...,vNt] isN t N_tNtOrder unitary matrix, Λ \LambdaΛ N r × N t N_r\times N_t Nr×NtOrder diagonal matrix, the elements on the diagonal are H's p = min ( N t , N r ) p=min(N_t,N_r)p=m i n ( Nt,Nr) singular values​​σ 1 , σ 2 , . . . , σ p \sigma_1,\sigma_2,...,\sigma_pp1,p2,...,pp, arrange them in descending order, that is, σ 1 > σ 2 > . . . > σ p \sigma_1>\sigma_2>...>\sigma_pp1>p2>...>pp. Then, the known channel matrix HHH , the optimal beamforming matrixW r W_rWrand F t F_tFtare the right singular vector and left singular vector corresponding to the largest singular value of H respectively: W r = U [ 1 : N s ] W_r=U[1:N_s]Wr=U[1:Ns] andF t = V [ 1 : N s ] F_t=V[1:N_s]Ft=V[1:Ns] , if the transmitted signal power is normalized, the maximum SNR isσ i 2 / σ 2 , i = 1 , . . . , N s \sigma_i^2/\sigma^2,i=1,. ..,N_spi2/ p2,i=1,...,Ns


References:
[1] Fu Zigang. Research on Beamforming in Millimeter Wave MIMO System[D]. 2016.

Guess you like

Origin blog.csdn.net/han_xj/article/details/111410019
Recommended