KR101599189B1 - Channel Estimation Method for Large Scale MIMO Systems - Google Patents

Channel Estimation Method for Large Scale MIMO Systems Download PDF

Info

Publication number
KR101599189B1
KR101599189B1 KR1020150037990A KR20150037990A KR101599189B1 KR 101599189 B1 KR101599189 B1 KR 101599189B1 KR 1020150037990 A KR1020150037990 A KR 1020150037990A KR 20150037990 A KR20150037990 A KR 20150037990A KR 101599189 B1 KR101599189 B1 KR 101599189B1
Authority
KR
South Korea
Prior art keywords
channel
matrix
toeplitz
receiver
present
Prior art date
Application number
KR1020150037990A
Other languages
Korean (ko)
Inventor
이문호
모하마드하니프
Original Assignee
전북대학교산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 전북대학교산학협력단 filed Critical 전북대학교산학협력단
Priority to KR1020150037990A priority Critical patent/KR101599189B1/en
Application granted granted Critical
Publication of KR101599189B1 publication Critical patent/KR101599189B1/en

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Radio Transmission System (AREA)

Abstract

A problem of a memory space for big data is on the rise because of too many users and a limited memory space. In a big MIMO system, Toeplitz channel can help the performance be improved, as well as higher electricity efficiency. The present invention relates to a method for channel estimation, and more specifically, to a method for channel estimation in a big MIMO system. According to an embodiment of the present invention, provided are a Toeplitz channel division based on matrix vectorization while using Toeplitz matrix in a channel for a big MIMO system, and a Toeplitz Jacket matrix divided into Cooley-Tukey sparse matrix like a fast Fourier transform (FFT). According to another embodiment of the present invention, a multiple-input multiple-output (MIMO) communication system comprises a transmitter which transmits a signal by using multiple antennas; and a receiver which receives a signal by using multiple antennas and estimates a channel by using the received signal.

Description

[0001] The present invention relates to a channel estimation method for large scale MIMO systems,

The present invention relates to a channel estimation method, and more particularly, to a channel estimation method in a large scale MIMO system.

Multiple-input multiple-output (MIMO) antenna technology is emerging as an effective technology for increasing the capacity and robustness of wireless communication systems. Due to the very high spectral efficiency and reliability improvement in free space and the potential of actually achieving the theoretical predicted merits of MIMO in terms of power efficiency, large scale MIMO systems with tens to hundreds of antennas are of great interest have.

A large-scale MIMO system can be thought of as a system with several thousand terminals simultaneously using an antenna array of several hundred antennas on the same time-frequency resource. The basic premise of large-scale MIMO since then not only accommodates all the advantages of existing MIMO, but also has a much larger advantage. Overall, large-scale MIMO is an energy-efficient, secure, robust, and optimal means of enabling future fixed and mobile broadband networks to use spectrum efficiently.

In wireless MIMO communication, the channel between the transmitter and the receiver can be designed according to Kronecker-model in some cases. This model is possible when there is a certain relationship between different antennas. When this model is designed, the radio channel covariance matrix can be modeled as a two-matrix Kronecker product of a smaller dimension. When these conditions are met, it is important to have an algorithm that fits this particular problem to increase accuracy and reduce computational costs.

Korean Patent No. 10-1573001 (Nov. 21, 2015) Korean Patent No. 10-1263257 (2013.05.06)

SUMMARY OF THE INVENTION It is an object of the present invention to provide a channel estimation method in a MIMO communication system using vectorization based on Toeplitz channel matrix decomposition.

According to an aspect of the present invention, there is provided a channel estimation method including: receiving a signal; And estimating the channel using the received signal, wherein the estimating step estimates the channel using the following equation.

Figure 112015026784704-pat00001

y (m) is the received signal vector of the m-th channel

x (m) is the transmission signal vector of the m-th channel

z (m) is the noise

Figure 112015026784704-pat00002
≪ RTI ID = 0.0 > m <

Figure 112015026784704-pat00003

I N is a unit matrix

The Kronecker product is the Kronecker product.

Then, the channel matrix can be decomposed as shown in the following equation.

Figure 112015026784704-pat00004

Meanwhile, a multiple-input multiple-output (MIMO) communication system according to another embodiment of the present invention includes: a transmitter for transmitting a signal using a plurality of antennas; And a receiver that receives a signal using a plurality of antennas and estimates a channel using the received signal, and the receiver estimates the channel using the above equation.

As described above, according to the embodiments of the present invention, channel estimation in a MIMO communication system becomes possible by using a vectorization technique based on Toeplitz channel matrix decomposition. This can improve memory efficiency, data processing speed, performance improvement, and power efficiency in the channel estimation process through channel matrix decomposition.

1 illustrates a forward fast algorithm,
FIG. 2 illustrates a reverse fast algorithm,
3 is a graph showing a large-scale MIMO capacity vs. CDF,
4 is a graph showing a large-scale MIMO capacity for the SNR difference,
5 is a diagram illustrating a MIMO communication system to which the present invention is applicable.

Hereinafter, the present invention will be described in detail with reference to the drawings.

In the embodiment of the present invention, a vectorization based on Toeplitz channel matrix decomposition is proposed using a Kronecker product. Also, we propose the possibility of switching between antennas to maximize system throughput. Also, we propose an optimal switching between SB (statistical beamforming) with high correlation and SM (spatial multiplexing) with low correlation. This switching method improves performance and requires minimum feedback information. This is because the switching method depends only on the two channel statistics, the average SNR and the spatial correlation.

1. Toeplitz matrix structure

The Toeplitz matrix has a constant element along the diagonal and is a structure that transitions from the upper left to the lower right

Figure 112015026784704-pat00005
procession
Figure 112015026784704-pat00006
to be.

Figure 112015026784704-pat00007

The Toeplitz matrix mentioned here is classified into two types as follows. The first classification is formed by a circulation matrix, where each row vector is rotated from left to right in relation to the previous column vector. In particular, the rotation as in equation (1)

Figure 112015026784704-pat00008
About
Figure 112015026784704-pat00009
to be. The second Toeplitz matrix is a negacycle matrix. Here,
Figure 112015026784704-pat00010
About
Figure 112015026784704-pat00011
to be.

This matrix shown in equation (1) is applied to various fields. For example, in the following equation

Figure 112015026784704-pat00012
Represents the input to the column vector,
Figure 112015026784704-pat00013
when
Figure 112015026784704-pat00014
Is zero.

Figure 112015026784704-pat00015

Therefore, the vector can be expressed as follows.

Figure 112015026784704-pat00016

Figure 112015026784704-pat00017

At this time, the matrix elements are as follows.

Figure 112015026784704-pat00018

The vector

Figure 112015026784704-pat00019
Lt; RTI ID = 0.0 >
Figure 112015026784704-pat00020
Time-invariant filter of the discrete-time causal.

Figure 112015026784704-pat00021
Toeplitz Jacket (TJ) matrix
Figure 112015026784704-pat00022
Can be decomposed into sparse matrices so that a fast algorithm can be implemented as follows. Figures 1 and 2 show the forward fast algorithm and the reverse fast algorithm.

Figure 112015026784704-pat00023

Figure 112015026784704-pat00024

therefore,

Figure 112015026784704-pat00025
Wow
Figure 112015026784704-pat00026
, We can see that it becomes a unit matrix as follows.

Figure 112015026784704-pat00027

Figure 112015026784704-pat00028
The fast algorithm of Fig.
Figure 112015026784704-pat00029
The high-speed algorithm of Fig.

A special case of the Toeplitz matrix is that each row of the matrix moves to the right one cycle from the top row

Figure 112015026784704-pat00030
About
Figure 112015026784704-pat00031
(5) and becomes a circulant matrix.

Figure 112015026784704-pat00032

This matrix is used for cyclic coding for error correction including DFT (Discrete Fourier Transform).

In the Toeplitz procession,

Figure 112015026784704-pat00033
The Toeplitz matrix, which deals with the role of the eigenvalue as it infinitely increases
Figure 112015026784704-pat00034
For the sequence of
Figure 112015026784704-pat00035
. A nonzero vector such as:
Figure 112015026784704-pat00036
If there is a complex scalar
Figure 112015026784704-pat00037
The matrix
Figure 112015026784704-pat00038
.

Figure 112015026784704-pat00039

in this case

Figure 112015026784704-pat00040
Is an eigenvector of.

Eigenvalue

Figure 112015026784704-pat00041
end
Figure 112015026784704-pat00042
, It can be approximated using integration, and the Hermitian Toeplitz matrix
Figure 112015026784704-pat00043
Eigenvalues for sequences
Figure 112015026784704-pat00044
Dealing with the asymptotic role of
Figure 112015026784704-pat00045
Which means that generality is not lost. This theorem has some technical requirements to be satisfied. For example, the coefficients associated with each other by
Figure 112015026784704-pat00046
Fourier series is a condition that must exist.

Figure 112015026784704-pat00047

So,

Figure 112015026784704-pat00048
Function
Figure 112015026784704-pat00049
, And vice versa. So the sequence of the matrix is sometimes
Figure 112015026784704-pat00050
. if,
Figure 112015026784704-pat00051
If it is Hermitian,
Figure 112015026784704-pat00052
If so,
Figure 112015026784704-pat00053
ego
Figure 112015026784704-pat00054
Is a real number value.

If you make a reasonable assumption

Figure 112015026784704-pat00055
Functions that are continuous in the range of
Figure 112015026784704-pat00056
The Fourier series and the Toeplitz sequence can be represented by the same face.

Figure 112015026784704-pat00057

2. System and Channel Model

The transmitting antenna

Figure 112015026784704-pat00058
And the receive antenna is
Figure 112015026784704-pat00059
And these antennas assume a large, spatially multiplexed point-to-multipoint MIMO system with thousands of antennas.
Figure 112015026784704-pat00060
Represents a channel gain matrix using an m-th channel. At this time, it is assumed that these elements have an average of 0, a variance of 1 and follow the iid Gaussian distribution. The reception vector y (m) in the case of using the m-th channel is expressed by Equation (13).

Figure 112015026784704-pat00061

here,

Figure 112015026784704-pat00062

Figure 112015026784704-pat00063

ego,

Figure 112015026784704-pat00064
The elements i
Figure 112015026784704-pat00065
. ≪ / RTI > In equation (14)
Figure 112015026784704-pat00066
Has the Toeplitz structure as shown in Eq. (15).

Figure 112015026784704-pat00067

Therefore, equation (13) can be rewritten as in equation (16).

Figure 112015026784704-pat00068

Matrix vectorization and Kronecker product identity

Figure 112015026784704-pat00069
(16) can be rewritten as Equation (17).

Figure 112015026784704-pat00070

Several channel vector decomposition examples are presented here. For example, the following cases are assumed.

Figure 112015026784704-pat00071

here,

Figure 112015026784704-pat00072
Assuming that

Figure 112015026784704-pat00073

, Where < RTI ID = 0.0 >

Figure 112015026784704-pat00074
Lt; / RTI &
Figure 112015026784704-pat00075
And receiving antenna
Figure 112015026784704-pat00076
Lt; RTI ID = 0.0 > impulse < / RTI > Equation (18) can now be rewritten as Equation (19).

Figure 112015026784704-pat00077

From Eq. (16) it can be written as

Figure 112015026784704-pat00078

Therefore, the Kronecker product used in Eq. (17) can be expressed as Eq. (21).

Figure 112015026784704-pat00079

Figure 112015026784704-pat00080
sign
Figure 112015026784704-pat00081
and
Figure 112015026784704-pat00082
Assuming

Figure 112015026784704-pat00083

And the received signal can be obtained by the matched filter as follows.

Figure 112015026784704-pat00084

Toeplitz

Figure 112015026784704-pat00085
An example of a channel is shown below.

Figure 112015026784704-pat00086

As a result, the following equation (24) can be obtained.

Figure 112015026784704-pat00087

In the receiver, the reception capacity can be expressed by Equation (25).

Figure 112015026784704-pat00088

In a similar way, Eq. (23) can be written as Eq. (26) and its capacity can be obtained as Eq. (27).

Figure 112015026784704-pat00089

3. Simulation

To analyze the meaning of the Toeplitz channel decomposition, a monte-carlo simulation was performed. 3 and 4 show results of a more realistic scenario, where the channel coefficients

Figure 112015026784704-pat00090
Is a real number with a Gaussian distribution with an average of 0 and a variance of 0.5 and a complex number of imaginary parts. channel
Figure 112015026784704-pat00091
Is random, so the capacity is also a random variable with a special distribution. The CCDF (Complimentary Cumulative Distribution Function) is used as an important unit for measuring the capacity of such a channel. This curve basically suggests the probability that the MIMO capacity will be above a certain threshold.

4. MIMO communication system

5 is a diagram illustrating a MIMO communication system to which the present invention is applicable. As shown in FIG. 5, the MIMO communication system is constructed to include a transmitter 100 that transmits signals using a plurality of antennas and a receiver 200 that receives signals using a plurality of antennas.

The receiver 200 estimates the channel using the received signal from the transmitter 100. [ Specifically, the receiver 200 estimates the channel by calculating H (m) described above in "2. System and channel model ".

Channels are modeled according to the Kronecker-model. Meanwhile, in order to save memory due to limited memory space, the receiver 200 decomposes the channel matrix using a Toeplitz channel matrix technique based on vectorization in a large scale MIMO communication system. Thereby, the radio channel covariance matrix is modeled as a Kronecker product of small-dimensional matrices.

In the embodiment of the present invention, vectorization based on the Toplitz channel matrix decomposition using Kronecker product is presented, and the Toeplitz Jacket matrix is decomposed into a Cooley-Tukey sparse matrix like the Fourier transform. In order to maximize throughput of the system, we proposed the possibility of switching between antennas.

It is common to adjust between the coding and modulation modes according to the SNR. Embodiments of the present invention have suggested the possibility of proper switching between SM under conditions with high correlation and SM under conditions with low correlation. This switching method increases the performance and requires minimum feedback information because it depends only on the spatial correlation with the two channel characteristics, that is, the average SNR.

5. Toeplitz inverse matrix

The Toeplitz inverse of Eq. (1) can be seen in the following example.

Example 1.

Figure 112015026784704-pat00092
In the case of

Figure 112015026784704-pat00093

Example 2.

Figure 112015026784704-pat00094
In the case of

Figure 112015026784704-pat00095

Example 3.

Figure 112015026784704-pat00096
In the case of

Figure 112015026784704-pat00097

Example 4. For a Circulant Jacket matrix

Figure 112015026784704-pat00098

Figure 112015026784704-pat00099
, Where operator
Figure 112015026784704-pat00100
silver
Figure 112015026784704-pat00101
It means "Sookmyung".

Example 5.

Figure 112015026784704-pat00102
In the case of

Figure 112015026784704-pat00103

Figure 112015026784704-pat00104

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed exemplary embodiments, but, on the contrary, It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the present invention.

100: Transmitter
200: receiver

Claims (8)

The receiver receiving the signal; And
And the receiver estimating the channel using the received signal,
Wherein the estimating step estimates the channel using the following equation,
Figure 112016013303483-pat00105

y (m) is the received signal vector of the m-th channel
x (m) is the transmission signal vector of the m-th channel
z (m) is the noise
Figure 112016013303483-pat00106
≪ RTI ID = 0.0 > m <
Figure 112016013303483-pat00107

I N is a unit matrix
Is the Kronecker product,
The H is decomposed into a product of a Toffler's matrix Tn and diagonal matrices, and is processed by a fast algorithm, as expressed by a Cj (Circulant Jacket) matrix as shown in the following equation.
Figure 112016013303483-pat00122
delete delete delete delete A transmitter for transmitting a signal using a plurality of antennas; And
And a receiver for receiving the signal using a plurality of antennas and estimating the channel using the received signal,
The receiver estimates the channel using the following equation,
Figure 112016013303483-pat00112

y (m) is the received signal vector of the m-th channel
x (m) is the transmission signal vector of the m-th channel
z (m) is the noise
Figure 112016013303483-pat00113
≪ RTI ID = 0.0 > m <
Figure 112016013303483-pat00114

I N is a unit matrix
Is the Kronecker product,
The H is represented by a Cj (Circulant Jacket) matrix and is decomposed into a product of a Toffler's matrix Tn and diagonal matrices and processed by a fast algorithm. Communication system.
Figure 112016013303483-pat00123
delete delete
KR1020150037990A 2015-03-19 2015-03-19 Channel Estimation Method for Large Scale MIMO Systems KR101599189B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020150037990A KR101599189B1 (en) 2015-03-19 2015-03-19 Channel Estimation Method for Large Scale MIMO Systems

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020150037990A KR101599189B1 (en) 2015-03-19 2015-03-19 Channel Estimation Method for Large Scale MIMO Systems

Publications (1)

Publication Number Publication Date
KR101599189B1 true KR101599189B1 (en) 2016-03-04

Family

ID=55536019

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020150037990A KR101599189B1 (en) 2015-03-19 2015-03-19 Channel Estimation Method for Large Scale MIMO Systems

Country Status (1)

Country Link
KR (1) KR101599189B1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101916523B1 (en) 2016-09-28 2018-11-07 전북대학교산학협력단 Signal Processing Method and Apparatus using Genetic RNA Code based Jacket Matrix

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120119935A (en) * 2011-04-22 2012-11-01 한국전자통신연구원 Method and apparatus of detecting signal based on minimum mean square error in multiple input multiple output system
KR101263257B1 (en) 2008-10-24 2013-05-10 퀄컴 인코포레이티드 Method and apparatus for uplink network mimo in a wireless communication system
KR101573001B1 (en) 2009-08-24 2015-11-30 삼성전자주식회사 Receiver and method for using reference singnal thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101263257B1 (en) 2008-10-24 2013-05-10 퀄컴 인코포레이티드 Method and apparatus for uplink network mimo in a wireless communication system
KR101573001B1 (en) 2009-08-24 2015-11-30 삼성전자주식회사 Receiver and method for using reference singnal thereof
KR20120119935A (en) * 2011-04-22 2012-11-01 한국전자통신연구원 Method and apparatus of detecting signal based on minimum mean square error in multiple input multiple output system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
‘Toeplitz channel matrix decomposition with vectorization for very large-scale MIMO’, proc. of ICT Convergence (ICTC), 2013 International Conference on, pp. 974-976. 2013. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101916523B1 (en) 2016-09-28 2018-11-07 전북대학교산학협력단 Signal Processing Method and Apparatus using Genetic RNA Code based Jacket Matrix

Similar Documents

Publication Publication Date Title
Neumann et al. Covariance matrix estimation in massive MIMO
CN107483088B (en) Large-scale MIMO robust precoding transmission method
CN109951214B (en) Signal detection method suitable for large-scale MIMO system
Peken et al. Blind channel estimation for massive MIMO
US20120163488A1 (en) MIMO-OFDM Wireless Communication System
US11646776B2 (en) Massive MIMO beam domain robust precoding transmission method and system
JP5274664B2 (en) Method for precoding a signal for transmission in a wireless MIMO system
CN111162822B (en) Method and system for acquiring large-scale MIMO beam domain statistical channel information
CN108259397B (en) Large-scale MIMO system channel estimation method based on adaptive regularization subspace tracking compressed sensing algorithm
WO2015112883A1 (en) System and method for early termination in iterative null-space directed singular value decomposition for mimo
CN113922848A (en) Signal sending method, channel estimation method, sending end equipment and receiving end equipment
US20220052881A1 (en) Method and apparatus for improving multi-user multiplexing performance, device, and storage medium
Wang et al. Efficient channel statistics estimation for millimeter-wave MIMO systems
CN116192209A (en) Gradient uploading method for air computing federal learning under MIMO channel
CN115053465A (en) Information transmission method and device
CN109802901B (en) 3D MIMO channel estimation method and system based on angle of arrival measurement
Kashyap et al. Frequency-domain interpolation of the zero-forcing matrix in massive MIMO-OFDM
KR101599189B1 (en) Channel Estimation Method for Large Scale MIMO Systems
CN107370533B (en) Method, device and system for carrying out analog CSI feedback
Al-Salihi et al. Bayesian compressed sensing-based channel estimation for massive MIMO systems
CN102025459A (en) Nonparametric estimation ICA-based MIMO-OFDM system blind deconvolution method
JP2016524394A (en) Multi-input multi-output MIMO detection method, apparatus and system for transmission signal
US20230164002A1 (en) User detection technique, and method and apparatus for channel estimation in wireless communication system supporting massive multiple-input multiple-output
CN114629533A (en) Information geometry method and system for large-scale MIMO channel estimation
Li et al. Joint measure matrix and channel estimation for millimeter-wave massive MIMO with hybrid precoding

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20190122

Year of fee payment: 4