An Improved Anti-Interference Precoding of Large-Scale Fading System Based on Channel Inversion

: Recently, a large-scale fading precoding (LSFP) for the wireless massive multiple-input, multiple-output (MIMO) systems has been proposed. In this precoding, the channel information of all the cells using re-use pilot sequences is processed jointly, and pilot contamination and interference due to a certain number of antennas are e ﬀ ectively eliminated. Additionally, recent studies have found that research in the asymptotic ﬁeld can be applied to the wireless large-scale MIMO systems. In the LSFP, pilot contamination and signal interference will be completely eliminated when a number of antennas at a base station tend to be unlimited. In this research found that the LSFP method can eliminate most pilot contamination and interference in practical applications only when the number of antennas of the base station reaches hundreds of orders, which greatly increases the equipment construction cost. On the other hand, channel inversion denotes a multi-user channel modulation technology, where a vector signal generated between a user and a base station is used to form an inverse channel matrix so that the channels of each user are balanced during the transmission. In this paper, the channel inversion technology is used in the LSFP. The improved LSFP can e ﬀ ectively reduce the number of antennas required by the base station without a ﬀ ecting the performance of eliminating the pilot contamination and interference. It is shown that when the number of antennas of a base station tends to be unlimited, the improved LSFP can eliminate pilot contamination and signal interference. The simulation results show that in the same practical scenario, when the base station is equipped with the same number of antennas, the improved method can more e ﬀ ectively improve the anti-contamination and anti-interference performance over conventional LSFP. the MMSE pre-equalization algorithm is applied to the optimized precoding algorithm, and the advanced optimal LSFP based on the MMSE is obtained. Furthermore, the asymptotic regime is analyzed when the number of antennas of a base station approaches inﬁnity, and it is shown that the advanced zero-forcing precoding algorithm can completely eliminate all sources of interference and additive noise. Finally, the SINR expressions for anti-interference precoding are derived and the above-improved algorithm is simulated. The simulation results show that the anti-interference precoding that use the channel inversion technology achieves the same communication transmission rate as conventional LSFP when the number of base antennas reducing. That means the improved precoding scheme have an improvement on anti-interference and improving channel transmission rate without increasing the number of base antennas in the practical cellular communication network.


Introduction
A massive MIMO (multiple-input multiple-output) system entails a large number (tens or hundreds) of base station antennas serving a much smaller number of terminals. These systems demonstrate large gains in spectral and energy efficiency compared with the conventional MIMO technology [1]. Massive MIMO (multiple-input multiple-output) is no longer a "wild" or "promising" concept for future cellular networks. In 2018, the key ingredients of Massive MIMO have made it into the 5G standard. The development of fully digital Massive MIMO arrays for mmWave frequencies is well underway [2]. However, its research is prohibitively complicated and costly. In a few years, maybe Massive MIMO with fully digital transceivers will be a mainstream feature at both sub-6 GHz and mmWave frequencies.
regularization with the pseudo-noise in the MU-MIMO transmission. The performance evaluation can be carried out through computer simulation. In [10], researchers applied the regularized channel inversion (RCI) combined with one of the 5G waveform candidates, namely generalized frequency division multiplexing (GFDM), in multiuser MIMO downlink systems and presents how the RCI in the multiuser MIMO-GFDM systems performs such good results as the RCI implemented in the multiuser MIMO-OFDM (Orthogonal Frequency Division Multiplexing) systems based on bit error rate (BER) parameter without causing a significant increase in system complexity.
The design idea of the precoding LSFP based on the zero-forcing algorithm is relatively simple, but it can only effectively eliminate interference such as pilot pollution when the number of base station antennas is very large. Optimizing the design of the precoding LSFP can ensure that interference is effectively eliminated. However, it needs to increase the power of the transmitting device to achieve efficient communication.
Aiming at the problem that the number of base station antennas of existing LSFP precoding is too large, the construction cost is too high. Based on the zero-forcing precoding algorithm and the optimized precoding algorithm, this paper applies channel inversion to improve and optimize the precoding. In the case of fewer base station antennas, interference such as pilot pollution and noise is suppressed, and the channel transmission rate is not affected at the same time, thereby achieving the purpose of further increasing the system rate and increasing the SINR.
In this work, channel inversion technology is used to improve the standard large scale fading precoding (LSFP), obtaining an improved anti-interference precoding, which can effectively eliminate pilot contamination, noise and interference using a small number of antennas. First, we applied the channel inversion technique to derive the theoretical formula for the zero-forcing large scale fading precoding algorithm to obtain its SINR expression. Under ideal circumstances, we analyzed and proved that the improved zero-forcing precoding can achieve the ideal anti-interference performance. The communication transmission rate performance was verified in the simulation. Secondly, the channel inversion technique was used to derive the theoretical formula of the optimized precoding algorithm. The SINR expression was obtained. The performance of the communication transmission rate was verified in the simulation. Finally, the regularized channel inversion technique was used to derive the theoretical formula for the optimized precoding algorithm. Its SINR expression was obtained. The performance of the communication transmission rate was verified in simulation. The SINR simulation expression of LSAS precoding of massive MIMO system was verified by simulation. The simulation compares the precoding communication rate performance before and after improvement. The simulation results show that reducing the number of antennas has no effect on the communication transmission rate. After applying the improved method to downlink transmission precoding, the communication transmission rate can be effectively improved.
The remaining of this article is organized as follows. Section 2 describes the system model used in this paper, including the channel model and orthogonal frequency division multiplexing scheme. Section 3 employs the channel inversion to improve the LSFP zero-forcing precoding, optimizes the precoding algorithm and performs the theoretical analysis based on the improved anti-interference precoding matrix, and lastly, obtains the SINR and channel transmission rate. Section 4 sets the infinite number of antennas in the asymptotic field, and performs the progressive analysis of the improved anti-interference precoding to determine the channel anti-interference transcoding technology combined with the anti-jamming precoding that can produce benefits. Section 5 presents the simulations of the channel transmission rate in the limited antenna field and compares the performance of the improved LSFP precoding and conventional LSFP precoding. Section 6 concludes the paper.
Throughout the paper, we use A T , A * , Tr[A] and A (m,n) for matrix transposition, Hermitian transposition, the trace of matrix and the (m, n)-th entry of matrix A, respectively.

System Model
In this paper, the LSAS is set as an OFDM based cellular mobile communication network. In this kind of cellular network, the cells within a group use different OFDM orthogonal sub-channels, and the number of cells in the group is set to L [7]. The cell represents a regular hexagonal cell network. Each cell is equipped with a multi-antenna base station, and the base station can provide services to multiple users in the cell. The number of antennas of the base station is set to M, and each antenna is an omnidirectional antenna. The number of users in the cell is set to K, and each user terminal is equipped with an omnidirectional antenna. Accordingly, a cell constitutes a multi-user MIMO system, where K independent users can form a virtual antenna array with number of antennas K. The configuration of the base station and the user end constitutes a MIMO downlink system with the size of KM, or a MIMO uplink system with the size of MK. In the traditional single-user MIMO systems, the total rate capacity of the uplink and downlink channels of a multi-user MIMO system is proportional to min(M, K).
In the OFDM-based cellular communication network, due to frequency reuse, the channel spectrum utilization is very high and the system communication performance is greatly improved [11]. We further assume that Orthogonal Frequency-Division Multiplexing (OFDM) is used and that in a given subcarrier we have a flat-fading channel [7]. In order to simplify the presentation of the multi-user MIMO system and reduce complex settings, we set that all the subcarriers in the communication frequency band can be used for communication between base stations and users.

Channel Model Settings
In a cellular network OFDM-based system, one block group allocates all orthogonal subcarrier channels, and the other block group reuses subcarrier channels. For a cell in a cellular network, which uses an OFDM subcarrier channel, the channel factor is expressed as a correlation operation between a large-scale attenuation factor and a small-scale attenuation factor, which is expressed as in Figure 1.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 16 In this paper, the LSAS is set as an OFDM based cellular mobile communication network. In this kind of cellular network, the cells within a group use different OFDM orthogonal sub-channels, and the number of cells in the group is set to L [7]. The cell represents a regular hexagonal cell network. Each cell is equipped with a multi-antenna base station, and the base station can provide services to multiple users in the cell. The number of antennas of the base station is set to M, and each antenna is an omnidirectional antenna. The number of users in the cell is set to K, and each user terminal is equipped with an omnidirectional antenna. Accordingly, a cell constitutes a multi-user MIMO system, where K independent users can form a virtual antenna array with number of antennas K. The configuration of the base station and the user end constitutes a MIMO downlink system with the size of KM, or a MIMO uplink system with the size of MK. In the traditional single-user MIMO systems, the total rate capacity of the uplink and downlink channels of a multi-user MIMO system is proportional to min(M, K).
In the OFDM-based cellular communication network, due to frequency reuse, the channel spectrum utilization is very high and the system communication performance is greatly improved [11]. We further assume that Orthogonal Frequency-Division Multiplexing (OFDM) is used and that in a given subcarrier we have a flat-fading channel [7]. In order to simplify the presentation of the multi-user MIMO system and reduce complex settings, we set that all the subcarriers in the communication frequency band can be used for communication between base stations and users.

Channel Model Settings
In a cellular network OFDM-based system, one block group allocates all orthogonal subcarrier channels, and the other block group reuses subcarrier channels. For a cell in a cellular network, which uses an OFDM subcarrier channel, the channel factor is expressed as a correlation operation between a large-scale attenuation factor and a small-scale attenuation factor, which is expressed as in Figure  1. Cell j indicates the j-th cell in the area group, and the number of antennas of a cell base station is M; cell l indicates the l-th cell in the area group, and the number of cell users is K. The channel factor between the m-th antenna of a base station in the j-th cell and the k-th user in the l-th cell can be expressed as: where ℎ denotes a fast-fading factor which is a small-scale fading factor. It is mainly caused by a multipath propagation and modeled as . . . (0,1) random variables. Cell j indicates the j-th cell in the area group, and the number of antennas of a cell base station is M; cell l indicates the l-th cell in the area group, and the number of cell users is K. The channel factor between the m-th antenna of a base station in the j-th cell and the k-th user in the l-th cell can be expressed as: where h kl mj denotes a fast-fading factor which is a small-scale fading factor. It is mainly caused by a multipath propagation and modeled as i.i.d. CN(0, 1) random variables. Furthermore, β kl j denotes a slow fading factor which is a large-scale fading factor, mainly caused by the building blocking. We model the large-scale fading coefficient β kl j as a log-normal distribution based on the Urban Macro model. The specific value of β kl j is a positive real number, which is determined by shadowing and distance between the corresponding user and base station and it is given by: where ψ denotes the shadow fading coefficient with i.i.d. N 0, σ 2 ψ distribution and d kl j denotes the distance, measured in kilometer, between the j-th BS and the k-th user in the l-th cell.

OFDM Duplex Setting
In this article, a multi-user MIMO system that used the OFDMA (Orthogonal Frequency Division Multiple Access) technology is studied. In the OFDMA, each user is assigned to a group of subcarriers, and the number of subcarriers for each user in each cell is defined by the system scheduling. For the sake of simplicity, it is assumed that the system model uniformly sets K users in each cell, that is, K subcarriers. The communication link between a base station and a user is divided into downlink and uplink duplex modes. We consider LSAS (large scale antenna system) with the TDD (Time Division Duplexing) protocol. The TDD protocol was set up with three steps: Tu channel was used for uplink training; Tpchannel was used for processing precoding; Td channel was used for downlink training.
In TDD, the uplink and downlink transmissions were allocated to different communication time slots and could not be performed simultaneously; thus, the transmission and reception of a user signal was completed through the assigned subcarriers. However, they were performed at different times, and the corresponding time slots could be dynamically and flexibly allocated. Due to the reciprocity of uplink transmission and downlink transmission, the MIMO technology and antenna technology could be combined with time division duplexing.

LSFP
In the LSFP (large scale fading precoding) network, several assumptions were made, and they are as follows.
(1) Any base station in a cell can accurately estimate and obtain a large-scale fading factor between itself and any user in the group. (2) Any base station in a cell can send a large-scale fading factor matrix to any user in the cell.
(3) All base stations in a cellular network are connected to the same network hub controller.
The large-scale fading factor of the entire cellular network is visible to this network hub controller. (4) Any downlink signal in a a cellular network is visible to this network hub controller.
By adjusting the hardware of the base station equipment and using the reciprocity of uplink and downlink [12,13], the network hub controller in a cluster generates a LSFP matrix denoted as W k based on the large-scale fading factor β kl j [14]: The was is specifically divided into the precoding of downlink information from a base station to the user, and the decoding of uplink information from a user to the base station.

Downlink LSFP
In Step 1, base station of a cell performs channel estimation according to the uplink signal and uplink pilot training sequence sent by a user to obtain a large-scale attenuation factor and sends it to the network hub controller of the cellular network.
In Step 2, the network hub controller receives the large-scale fading factors sent by all the base stations in the area group and forms a large-scale fading factor precoding matrix.
In Step 3, the network hub controller sends the large-scale fading factor precoding matrix to each base station.
In Step 4, before the base station sends the information to a user, it performs a precoding operation on the signal according to the large-scale fading factor precoding matrix.
In the Downlink PCP (pilot contamination precoding) proposed for the downlink pilot contamination precoding [5], the PCP matrix is expressed as: where B k denotes a matrix composed of a large-scale fading factor β kl j , that is, the zero-forcing pilot contamination precoding (ZF-PCP).
In the LSFP proposed in [8], the LSFP matrix is expressed as: where γ A denotes a normalization factor that defines the power limit and B k denotes a matrix composed of the large-scale fading factor. β kl j , that is, the zero-forcing large-scale fading precoding (ZF-LSFP). By applying the asymptotic theory, it has been proven that systems using zero-forcing precoding can achieve communication performance. At a limited number of antennas, the ZF-PCP and ZF-LSFP with zero-forcing precoding achieve average performance, and there are also pilot contamination, noise and interference in the LSAS. In [4,7], an optimized precoding scheme named the Optimized LSFP (O-LSFP) was proposed; this scheme aims to optimize the precoding matrix in the field of limited antennas. Lemma 1.
W denotes an O-LSFP matrix, B denotes a large-scale fading factor matrix and P kl denotes a matrix obtained by a large-scale fading factor β kl j operation.
The SINR expression, after applying the O-LSFP matrix is as follows: The simulation results presented in [4,7] show that the number of antennas at a base station must be at the hundred-level to eliminate most of the pilot contamination, noise and interference. However, such a configuration of antennas will inevitably increase the cost of the base station equipment. In [7], the simulation assumed that each base station is equipped with M omnidirectional antennas and each user has a single omnidirectional antenna. In practice, there will be unnecessary power consumption, failing to utilize antenna resources fully and increasing the operating costs. Therefore, in this paper, the LSFP is improved to reduce the pilot contamination, noise and interference performance without reducing the number of antennas at a base station. At the same time, if the number of antennas at a base station antennas stays unchanged, the improved LSFP precoding scheme can improve the performance. Communication transmission rate to improve LSAS communication performance. It means that the improved LSFP precoding scheme can more effectively improve the anti-contamination and anti-interference performance over conventional LSFP as the number of antennas at a base station antennas stays unchanged.

Uplink LSFP
Similarly, the uplink transmission process consists of the following steps. In Step 1, network hub controller of a group receives the large-scale fading factors sent by all base stations in the group and forms a large-scale fading factor decoding matrix.
In Step 2, network hub controller sends the large-scale fading factor decoding matrix to each base station. In Step 3, when a base station receives the information sent by the user, the signal is decoded according to the large-scale fading factor decoding matrix.

Advanced LSFP
In this paper, channel inversion technology was applied to the ZF-LSFP and O-LSFP to obtain an improved anti-interference precoding named the Advanced LSFP (A-LSFP). Channel inversion belongs to precoding technology. In the downlink communication, channel inversion can be used to eliminate data interference at a base station. In this work, channel inversion and regularized channel inversion techniques are implemented in the multi-cell LSAS. Channel inversion is similar to the zero-forcing pre-equalization technique in the single-user MIMO systems, while regularized channel inversion is similar to the MMSE pre-equalization technique in the single-user MIMO systems.
In the pre-equalization technology, weighting matrix W represents the pre-equalizing weighting matrix. When the zero-forcing pre-equalization is used, the corresponding weighting matrix W is expressed as: where β denotes a constant that keeps the total power constant after pre-equalization.
Appl. Sci. 2020, 10, 867 8 of 16 By applying the channel inversion technology to the LSAS and substituting β into the LSFP matrix, the Advanced ZF-LSFP (A-ZF-LSFP) is obtained, which can be expressed as follows: where W k denotes a large-scale fading factor matrix for channel inversion processing and γ A denotes a normalization factor that guarantees power limitation. β represents a constant that keeps the total power constant after pre-equalization. Substitute the ZF-LSFP (A-ZF-LSFP) to which channel inversion technology is applied: The SINR of the A-ZF-LSFP can be expressed as: In the LSAS, the channel inversion technology is applied to optimize precoding. The Optimized LSFP (O-LSFP), and the Advanced O-LSFP (AO-LSFP) can be obtained: W denotes the A-O-LSFP precoding matrix for channel inversion processing. Applying W to Lemmas 1 and 2, it is easy to get the SINR expression of A-O-LSFP. The SINR of the A-O-LSFP can be expressed as: where B represents a large-scale fading factor matrix and P kl represents a matrix obtained by using a large-scale fading factor β kl j operation.
In the MMSE pre-equalization, the corresponding weighting matrix W can be expressed as: where β denotes a constant that keeps the total power constant after pre-equalization. By applying regular channel inversion technology to the LSAS, substituting β into the LSFP matrix, and using the MMSE pre-equalization weighting matrix W to represent the O-LSFP matrix, the MMSE Advanced O-LSFP (MMSE-AO-LSFP) can be obtained as follows: where W k denotes a matrix obtained by performing regular channel inversion processing on a large-scale fading factor matrix. W denotes the MMSE-A-O-LSFP precoding matrix for regularized channel inversion processing.

Lemma 3.
In this work, the theorem presented in [10] was used, and it is as follows: The large-scale fading factor β kl j can be applied to obtain the data transmission rate, which is expressed as follows: At the receiving end, in order to compensate for the amplification effect of a transmitter, the receiver uses the automatic gain control technology to adjust the received signal so that the received signal is multiplied by 1 β , and then the received signal of the same strength is restored.

Performance Analysis
When the number of antennas of a base station approaches infinity, the pre-coded A-ZF-LSFP and AO-LSFP can achieve excellent communication performance after channel inversion, and in the same network group, all users will reach the theoretically maximal network transmission rate.
First, before the asymptotic analysis, the most suitable LSFP matrix is determined. According to the previous analysis, the channel inversion technology is applied to LSFP to obtain a suitable precoding matrix. The zero-forcing A-ZF-LSFP matrix for channel inversion processing is obtained as: where γ A denotes the normalization factor that guaranteed the power limit and ϕ A = γ A * β denotes the normalization factor for the zero-forcing precoding A-ZF-LSFP.
Among them: Therefore, when M → ∞ , J kl 1 = 0, J kl 1 = Mρ f ρ r τ we could get the idealized SINR of the A-ZF-LSFP: In the OFDM communication system, it means that noise and interference can be removed when the number of antennas tends to be infinite.
In practical applications, the number of antennas in the LSAS is often tens or hundreds. In the following section, the channel transmission rate of the optimized precoding before and after channel inversion is compared in detail to prove that the optimized LSFP after channel inversion processing can achieve the same anti-contamination and anti-interference performance as equipped with a smaller number of base antennas than conventional LSFP.

Simulation
The performance of anti-interference precoding was analyzed in terms of the communication transmission rate. The simulation graph represents the cumulative distribution function (CDF), which can reflect the received signal-to-interference and noise ratio and the information transmission throughput rate [15]. In the simulation, the performance of communication transmission rate of improved and conventional precoding schemes, which includes ZF-LSFP, AZF-LSFP, O-LSFP, AO-LSFP and MMSE-AO-LSFP, is compared under different base antenna numbers.
The simulation model was based on Section 2 channel model setting. We model the large-scale fading coefficient β kl j as log-normal distribution based on the Urban Macro model that followed the 3GPP (3rd Generation Partnership Project) standard. Cell base stations and user settings in Figure 2: the number of cells contained in a single cell group was L = 19, the number of users contained in a single cell was K = 10, and the number of cell base station antennas was M. The cell was randomly distributed to ten mobile users. The cell radius was r = 1 km. The cell base station was set in the center. The mobile user was randomly located in a circular area having the base station as a center, a radius of 62.5 m. The cell numbers in the area group were set to L , l = 1, ⋯ , 19. Each cell center point was set with two-dimensional coordinates and the L coordinate is L (0,0). According to the cell radius, the coordinates of all cell center points L (m , n ) can be generated.
The distribution map of users of a single cell is presented in Figure 3 left, where each dot represents a user. The user distribution range was set in the simulation, and the generated points were randomly distributed. Randomly selecting a hexagonal cell, and the center point is the base station position. The positions of 10 users in this cell are randomly generated. The user positions are evenly distributed in the cell. We eliminate the extreme case that all users are distributed at the cell edge.
A random distribution diagram of a total of 190 users grouped into 19 cells within the cellular network to which they belonged is presented in Figure 3 right. The relative distance between a user and 19 base stations was obtained using the same principle as in spatial channel model for multiple input multiple output (MIMO) simulations of Technical specification group radio access network. There are 19 communities in the cellular network, with a total of 19 × 10 users. The location of 190 users is randomly generated in a hexagonal cell with the center point coordinate of (x, y), and a group of 10 users Substitute the coordinates of 19 center points L (m , n ), l = 1, ⋯ , 19 in turn, and obtain the coordinate of all users in each cell. Then, we can calculate the distance according to coordinates of the j-th BS and the k-th user in the l-th cell. The cell numbers in the area group were set to L 1 , l = 1, · · · , 19. Each cell center point was set with two-dimensional coordinates and the L 1 coordinate is L 1 (0,0). According to the cell radius, the coordinates of all cell center points L 1 (m l , n l ) can be generated.
The distribution map of users of a single cell is presented in Figure 3 left, where each dot represents a user. The user distribution range was set in the simulation, and the generated points were randomly distributed. Randomly selecting a hexagonal cell, and the center point is the base station position. The positions of 10 users in this cell are randomly generated. The user positions are evenly distributed in the cell. We eliminate the extreme case that all users are distributed at the cell edge. Considering the movement of end users, it was assumed that users move freely in the cell to which they belong, but irrespective of the extreme case where all users move to the edge of the cell. The distribution of users in the cell is still roughly evenly distributed. The positions of 190 users were randomly generated in the hexagonal cell with the center point coordinate of (x, y) for many times. The user position distribution changes as shown in the point in Figure 4, which simulates the free movement of users in each cell. A random distribution diagram of a total of 190 users grouped into 19 cells within the cellular network to which they belonged is presented in Figure 3 right. The relative distance between a user and 19 base stations was obtained using the same principle as in spatial channel model for multiple input multiple output (MIMO) simulations of Technical specification group radio access network.
There are 19 communities in the cellular network, with a total of 19 × 10 users. The location of 190 users is randomly generated in a hexagonal cell with the center point coordinate of (x,y), and a group of 10 users Substitute the coordinates of 19 center points L 1 (m l , n l ), l = 1, · · · , 19 in turn, and obtain the coordinate of all users in each cell. Then, we can calculate the distance d kl j according to coordinates of the j-th BS and the k-th user in the l-th cell.
Considering the movement of end users, it was assumed that users move freely in the cell to which they belong, but irrespective of the extreme case where all users move to the edge of the cell. The distribution of users in the cell is still roughly evenly distributed. The positions of 190 users were randomly generated in the hexagonal cell with the center point coordinate of (x,y) for many times. The user position distribution changes as shown in the point in Figure 4, which simulates the free movement of users in each cell. Considering the movement of end users, it was assumed that users move freely in the cell to which they belong, but irrespective of the extreme case where all users move to the edge of the cell. The distribution of users in the cell is still roughly evenly distributed. The positions of 190 users were randomly generated in the hexagonal cell with the center point coordinate of (x, y) for many times. The user position distribution changes as shown in the point in Figure 4, which simulates the free movement of users in each cell. According to the set system model, a large-scale fading factor β was obtained using the 3GPP urban macro model, which is given by (2). Other simulation settings: the transmit powers of the BS According to the set system model, a large-scale fading factor β kl j was obtained using the 3GPP urban macro model, which is given by (2). Other simulation settings: the transmit powers of the BS and users were ρ r = 23 dBm and ρ f = 23 dBm, respectively. The noise figures of the BS and user were 5 and 10, respectively; the bandwidth was 20 MHZ.
In the simulation, different numbers of base station antennas were used, M was set to 100, 1000, 10,000, and 1,000,000. Channel transmission rate can be calculated by (36). The simulation results show the lowest achievable rate in a cellular network group.
In Figure 5, the simulation results are compared with those of the ZF-LSFP that are provided in [4]. The comparison results show that when M reached the value of 10 6 , the achievable channel transmission rate of the A-ZF-LSFP reached the value of 10 0 , and the performance of channel transmission rate increased by 60%. When M was between 10 2 and 10 4 , the channel transmission rate of the A-ZF-LSFP increased by 90%.
In Figure 5, the simulation results are compared with those of the ZF-LSFP that are provided in [4]. The comparison results show that when M reached the value of 10 , the achievable channel transmission rate of the A-ZF-LSFP reached the value of 10 , and the performance of channel transmission rate increased by 60%. When M was between 10 and 10 , the channel transmission rate of the A-ZF-LSFP increased by 90%. As Figure 6 shows, when the number of antennas was 100, the probability that the unimproved optimized precoding could achieve a channel transmission rate greater than 10 was about 0.8, and the probability of reaching the value of 10 was zero; when the number of antennas was 1000, the probability that the O-LSFP could achieve a channel transmission rate greater than 10 was about 0.9, and the probability of reaching the value of 10 was close to zero. As Figure 6 shows, when the number of antennas was 100, the probability that the unimproved optimized precoding could achieve a channel transmission rate greater than 10 0 was about 0.8, and the probability of reaching the value of 10 1 was zero; when the number of antennas was 1000, the probability that the O-LSFP could achieve a channel transmission rate greater than 10 0 was about 0.9, and the probability of reaching the value of 10 1 was close to zero. As Figure 7 shows, when the number of antennas was 100, the probability that the optimized precoding with channel inversion processing could achieve a channel transmission rate greater than 10 was about 0.9, and the probability of reaching the value of 10 was 0.05; when the number of antennas was 1000, the probability that the A-O-LSFP could achieve a channel transmission rate greater than 10 was about 0.96, and the probability of reaching the value of 10 was about 0.5. In the simulations, when the number of antennas was a few hundreds, the transmission rate of the A-O-LSFP was increased by at least 10%. As Figure 7 shows, when the number of antennas was 100, the probability that the optimized precoding with channel inversion processing could achieve a channel transmission rate greater than 10 0 was about 0.9, and the probability of reaching the value of 10 1 was 0.05; when the number of antennas was 1000, the probability that the A-O-LSFP could achieve a channel transmission rate greater than 10 0 was about 0.96, and the probability of reaching the value of 10 1 was about 0.5. In the simulations, when the number of antennas was a few hundreds, the transmission rate of the A-O-LSFP was increased by at least 10%.
precoding with channel inversion processing could achieve a channel transmission rate greater than 10 was about 0.9, and the probability of reaching the value of 10 was 0.05; when the number of antennas was 1000, the probability that the A-O-LSFP could achieve a channel transmission rate greater than 10 was about 0.96, and the probability of reaching the value of 10 was about 0.5. In the simulations, when the number of antennas was a few hundreds, the transmission rate of the A-O-LSFP was increased by at least 10%.  Figure 8 shows that when the number of antennas was 100, the probability that the optimized precoding with regular channel inversion processing could achieve a channel transmission rate greater than 10 was about 0.7, and the probability that it could reach the value of 10 was 0.15; when the number of antennas was 1000, the probability that the MMSE-A-O-LSFP could achieve a channel transmission rate greater than 10 was about 0.9, and the probability of reaching the value of 10 was about 0.25. The simulation results showed that when the number of antennas was about 100, the MMSE-A-O-LSFP had a probability of 15% to reach 10 times the transmission rate gain, but there was still a 10% attenuation probability. The MMSE-A-O-LSFP had a probability of reaching 10  Figure 8 shows that when the number of antennas was 100, the probability that the optimized precoding with regular channel inversion processing could achieve a channel transmission rate greater than 10 0 was about 0.7, and the probability that it could reach the value of 10 1 was 0.15; when the number of antennas was 1000, the probability that the MMSE-A-O-LSFP could achieve a channel transmission rate greater than 10 0 was about 0.9, and the probability of reaching the value of 10 1 was about 0.25. The simulation results showed that when the number of antennas was about 100, the MMSE-A-O-LSFP had a probability of 15% to reach 10 times the transmission rate gain, but there was still a 10% attenuation probability. The MMSE-A-O-LSFP had a probability of reaching 10 times the transmission rate gain of 25%, and there was no attenuation of the transmission rate gain, and the improved performance is completely improved.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 15 of 16 times the transmission rate gain of 25%, and there was no attenuation of the transmission rate gain, and the improved performance is completely improved.

Conclusions
This paper proposed and studied the anti-interference precoding algorithms that apply the channel inversion technology and use the large-scale fading factor matrix. First, channel-inversion processing was applied to the zero-forcing precoding matrix, and the advanced zero-forcing precoding algorithm is developed. Next, the advanced optimized LSFP is presented, which is based on the optimized precoding algorithm and uses the regularized pre-equalization technology. In addition, the MMSE pre-equalization algorithm is applied to the optimized precoding algorithm, and the advanced optimal LSFP based on the MMSE is obtained. Furthermore, the asymptotic regime is analyzed when the number of antennas of a base station approaches infinity, and it is shown that the advanced zero-forcing precoding algorithm can completely eliminate all sources of interference and additive noise. Finally, the SINR expressions for anti-interference precoding are derived and the

Conclusions
This paper proposed and studied the anti-interference precoding algorithms that apply the channel inversion technology and use the large-scale fading factor matrix. First, channel-inversion processing was applied to the zero-forcing precoding matrix, and the advanced zero-forcing precoding algorithm is developed. Next, the advanced optimized LSFP is presented, which is based on the optimized precoding algorithm and uses the regularized pre-equalization technology. In addition, the MMSE pre-equalization algorithm is applied to the optimized precoding algorithm, and the advanced optimal LSFP based on the MMSE is obtained. Furthermore, the asymptotic regime is analyzed when the number of antennas of a base station approaches infinity, and it is shown that the advanced zero-forcing precoding algorithm can completely eliminate all sources of interference and additive noise. Finally, the SINR expressions for anti-interference precoding are derived and the above-improved algorithm is simulated. The simulation results show that the anti-interference precoding that use the channel inversion technology achieves the same communication transmission rate as conventional LSFP when the number of base antennas reducing. That means the improved precoding scheme have an improvement on anti-interference and improving channel transmission rate without increasing the number of base antennas in the practical cellular communication network.