OAM multiplex communication system and method for compensating for inter-mode interference

Description

TECHNICAL FIELD

[0001] The present invention relates to an OAM multiplexing communication system and an inter-mode interference compensation method that subject a radio signal to spatial multiplex transmission by using orbital angular momentum (OAM) of an electromagnetic wave.BACKGROUND ART

[0002] There is a technique that subjects a radio signal to spatial multiplex transmission by using OAM, as a technique for improving a transmission capacity in wireless communication. An electric wave having an OAM mode has a feature that an equiphase surface thereof is distributed in a spiral form along a rotational direction centering on a propagation axis of a beam. A mode in which a cycle of the spiral formed by the equiphase surface is 2π×a is referred to as an OAM mode a. Since different OAM modes have orthogonality in a rotational direction, it is possible to subject signals of a plurality of OAM modes to spatial multiplex transmission. For example, since a signal of an OAM mode 1 and a signal of an OAM mode 2 are orthogonal to each other on a space, it is possible to separate the signal of the OAM mode 1 and the signal of the OAM mode 2 from each other by an OAM mode separation circuit on a receiving side, even though the signals of those modes are simultaneously transmitted from a transmitting antenna. As a method for subjecting signals of OAM modes to multiplex transmission, Non-Patent Document 1 reports a method that uses a uniform circular array (UCA) in which a plurality of antenna elements are disposed circularly at an equal interval and a Butler matrix circuit.PRIOR ART DOCUMENT

[0003] Non-Patent Document 1: E. Sasaki, M. Hirabe, T. Maru, N. Zein, "Pragmatic OAM with polarization multiplexing transmission for future 5G ultra-high capacity radio", in proc. of EuMA2016, Oct. 2016.

[0004] WO 2017 / 125969 A1 describes an antenna array with antenna elements arranged in a lattice and with signals multiplexed over a plurality of orthogonal channels by means of wireless signals that have a plurality of orbital angular momentum propagation modes. Multiplexed signals with a number of signals that corresponds to a number of OAM propagation modes are transmitted.

[0005] XIONG WU ET AL: "Inter-mode crosstalk compensation for radio orbital angular momentum multiplexing systems under misaligned condition using multiple-input multiple-output techniques", 2016 15TH INTERNATIONAL CONFERENCE ON OPTICAL COMMUNICATIONS AND NETWORKS (ICOCN), IEEE, 24 September 2016 (2016-09-24), pages 1-3, XP033075089, DOI: 10.1109 / ICOCN.201 6.7875646 describes inter-mode in radio orbital angular momentum mode systems with a misalignment between a transmitter and a receiver. A multiple input multiple output analytical model is used. The misaligned conditions include off axis and non-parallel cases.

[0006] WANGJOO-LEE ET AL: "Experimental results of triply multiplexed microwave orbital angular momentum mode transmission", 2016 INTERNATIONAL CONFERENCE ON INFORMATION AND COMMUNICATION TECHNOLOGY CONVERGENCE (ICTC), IEEE, 19 October 201 6 (2016-10-19), pages 765-767, XP033015558, DOI: 10.1109 / ICTC.201 6.7763290, Describes an OAM system with a 3x4 Butler matrix. Intermodal effects are measured at different distances.

[0007] US 2018 / 069631 A1, Bribes transmitting and orthogonal function, so signal over a communication link including generating a mode crosstalk matrix illustrating mode crosstalk between transmitted modes and adjacent modes within a fiber, selecting adjacent modes to multiplex together based on entries within the generated mode crosstalk matrix being less than or equal to a predetermined value; and multiplexing the transmitted modes and the selected adjacent modes together into the orthogonal function processed signal for transmission on the communications link on the fiber.SUMMARY OF INVENTION PROBLEMS TO BE SOLVED BY THE INVENTION

[0008] In the case of the configuration disclosed in Non-Patent Document 1, signals of respective OAM modes can be separated only in a line-of-sight situation where a transmitting UCA and a receiving UCA are disposed at positions facing each other and there is no reflected wave. However, in an actual operation, it is necessary to consider a case where the transmitting UCA and the receiving UCA are provided at fixed positions deviated from the front facing disposition, or an influence of the reflected wave or the like due to a surrounding environment. When a channel matrix between the transmitting UCA and the receiving UCA is deviated from the ideal front facing disposition due to the influence, an interference component between OAM modes remains in an output stage of the OAM mode separation circuit of the receiving side.

[0009] For example, when the receiving side tries to separate a signal of an OAM mode 1, it is possible to extract only a signal component of the OAM mode 1 when the transmitting UCA and the receiving UCA are in a state of the ideal front facing disposition. On the other hand, when the transmitting UCA and the receiving UCA are provided at fixed positions deviated from the ideal front facing disposition, or when the channel matrix is deviated from the ideal one due to the influence of the reflected wave, a propagation path, an RF circuit and the like, an interference component of another mode such as the adjacent OAM mode 2 or OAM mode 0 is superimposed on the signal component of the OAM mode 1 at a circuit output stage and thus communication quality is deteriorated, even though separation of the signal of the OAM mode 1 is tried.

[0010] An object of the present invention is to provide an OAM multiplexing communication system and an inter-mode interference compensation method that can compensate, with a small computation amount, for inter-mode interference caused by a deviation of an optical axis between a transmitting UCA and a receiving UCA, a tilt, a reflected wave, and the like.MEANS FOR SOLVING THE PROBLEMS

[0011] The object of the invention this solved by the subject matter of the independent claims. Advantageous embodiments are disclosed in the dependent claims.

[0012] A first example is directed to an OAM multiplexing communication system that uses one or more OAM modes and multiplexes signals of one or more sequences for each OAM mode. The OAM multiplexing communication system includes a transmitting station and a receiving station. The transmitting station includes: a transmitting antenna that uses an M-UCA in which a plurality of UCAs having different diameters are concentrically disposed; and an OAM mode generation unit that simultaneously generates the one or more OAM modes from each UCA of the transmitting antenna. The receiving station includes: a receiving antenna that has a configuration equivalent to the M-UCA; an OAM mode separation unit that separates signals received by each UCA of the receiving antenna for each OAM mode; and a channel estimation / interference compensation unit that estimates a channel matrix for each OAM mode separated from the signals received by the each UCA of the receiving antenna, and compensates for inter-mode interference between the OAM modes by using a weight calculated from the channel matrix. The channel estimation / interference compensation unit is configured to select, for each OAM mode, signals of a subject mode and an adjacent mode from among the signals of the OAM modes separated by the OAM mode separation unit, and to compensate for the inter-mode interference by multiplying an approximate weight calculated by using channel matrixes of the subject mode and the adjacent mode.

[0013] In the OAM multiplexing communication system of the first example, the channel estimation / interference compensation unit is configured to determine, as a range of the adjacent mode, a range in which a performance of an equalization process using the approximate weight calculated based on a known signal sequence transmitted from the transmitting station satisfies a predetermined criterion.

[0014] In the OAM multiplexing communication system of the first example, the channel estimation / interference compensation unit is configured to use an approximate channel matrix obtained by approximating the channel matrix as a band matrix to compensate for the inter-mode interference based on Cholesky decomposition of a matrix including the approximate channel matrix.

[0015] A second example is directed to an inter-mode interference compensation method for an OAM multiplexing communication system that uses one or more OAM modes and multiplexes signals of one or more sequences for each OAM mode. The OAM multiplexing communication system includes a transmitting station and a receiving station. The transmitting station includes: a transmitting antenna that uses an M-UCA in which a plurality of UCAs having different diameters are concentrically disposed; and an OAM mode generation unit that simultaneously generates the one or more OAM modes from each UCA of the transmitting antenna. The receiving station includes: a receiving antenna that has a configuration equivalent to the M-UCA; an OAM mode separation unit that separates signals received by each UCA of the receiving antenna for each OAM mode; and a channel estimation / interference compensation unit that estimates a channel matrix for each OAM mode separated from the signals received by the each UCA of the receiving antenna, and compensates for inter-mode interference between the OAM modes by using a weight calculated from the channel matrix. The inter-mode interference compensation method includes: selecting, by the channel estimation / interference compensation unit, for each OAM mode, signals of a subject mode and an adjacent mode from among the signals of the OAM modes separated by the OAM mode separation unit, and compensating for the inter-mode interference by multiplying an approximate weight calculated by using channel matrixes of the subject mode and the adjacent mode.

[0016] In the inter-mode interference compensation method of the second example, the channel estimation / interference compensation unit determines, as a range of the adjacent mode, a range in which a performance of an equalization process using the approximate weight calculated based on a known signal sequence transmitted from the transmitting station satisfies a predetermined criterion.

[0017] In the inter-mode interference compensation method of the second exampleA, the channel estimation / interference compensation unit uses an approximate channel matrix obtained by approximating the channel matrix as a band matrix to compensate for the inter-mode interference based on Cholesky decomposition of a matrix including the approximate channel matrix.EFFECT OF THE INVENTION

[0018] According to the present invention, regarding the inter-mode interference caused by the deviation of the optical axis between the transmitting UCA and the receiving UCA, the tilt, the reflected wave, and the like, it is possible to compensate for the inter-mode interference with a small computation amount by performing the equalization process with restricting a range of interference of the adjacent mode for each OAM mode.BRIEF DESCRIPTION OF THE DRAWINGS

[0019] FIG. 1 is a diagram showing a configuration example of an OAM multiplexing communication system according to the present invention. FIG. 2 is a diagram showing a configuration example of an M-UCA of the OAM multiplexing communication system according to the present invention. FIG. 3 is a diagram showing a configuration example of a channel estimation / interference compensation unit 23. FIG. 4 is a diagram showing an example of a channel matrix. FIG. 5 is a diagram showing an example of processing by the channel estimation / interference compensation unit 23. FIG. 6 is a diagram showing performance characteristics according to adjacent mode considered. FIG. 7 is a diagram showing an example of an approximate channel matrix. EMBODIMENTS

[0020] FIG. 1 is a diagram showing a configuration example of an OAM multiplexing communication system according to the present invention.

[0021] In FIG. 1, a transmitting station 10 includes a plurality of transmitting UCAs 11-1 to 11-N TX as transmitting antennae, and OAM mode generation units 12-1 to 12-N TX each of which receives modulation signals of a plurality of sequences respectively transmitted in OAM modes 1 to L, adjusts phases of the signals to be transmitted as signals of the OAM modes 1 to L, and outputs the signals to antenna elements of each transmitting UCA.

[0022] A receiving station 20 includes a plurality of receiving UCAs 21-1 to 21-N RX as receiving antennae, OAM mode separation units 22-1 to 22-N RX each of which separates the signals of the OAM modes 1 to L from signals received by each receiving UCA, and a channel estimation / interference compensation unit 23 that receives signals r 1 to r L collected for each of the OAM modes 1 to L, and performs channel estimation and interference compensation to output signals r 1 to r L .

[0023] Here, L is the number of OAM modes in use, and 1 to L each indicates a number of the OAM mode. If L = 5, for example, signals of the OAM modes -2, -1, 0, 1, and 2 are multiplexed and transmitted. Any one the OAM modes is hereinafter referred to as an OAM mode k.

[0024] The transmitting UCAs 11-1 to 11-N TX and the receiving UCAs 21-1 to 21-N RX constitute a multi-uniform circular array (Multi-UCA) in which a plurality of UCAs are concentrically disposed as shown in FIG. 2. Herein, a configuration in which four UCAs having different radii with each other are disposed is shown. In other words, N TX = N RX = 4, and the UCAs are referred to as a first UCA, a second UCA, a third UCA, and a fourth UCA in this order from the inner UCAs. As an example, each UCA includes sixteen antenna elements indicated by • in FIG.2, but the number of antenna elements of each UCA is not necessarily required to be the same.

[0025] FIG. 3 shows a configuration example of the channel estimation / interference compensation unit 23. Here, L is assumed to be 5, and the OAM modes 1 to 5 are represented.

[0026] In FIG. 3, the channel estimation / interference compensation unit 23 receives the signals r 1 to r 5 obtained by the OAM mode separation units 22-1 to 22-N RX that separate the received signals for each of the OAM modes 1 to 5 and collect the signals for each of the OAM modes 1 to 5. In other words, the signal r k of the OAM mode k is a set of signals of the OAM mode k received by the receiving UCAs 21-1 to 21-N RX . The signal r 1 of the OAM mode 1 to the signal r 5 of the OAM mode 5 are input not only to a channel estimation / weight calculation unit 31 but also to weight multiplication units 33-1 to 33-5 through adjacent mode selection units 32-1 to 32-5 such that the signals are subjected to interference compensation.

[0027] The channel estimation / weight calculation unit 31 estimates a channel matrix H between the OAM modes combined via the OAM mode generation units 12 of the transmitting station 10, propagation channels, and the OAM mode separation units 22 of the receiving station 20 when the signals of the OAM modes 1 to 5 are multiplexed and transmitted. The channel estimation / weight calculation unit 31 further calculates weights U k,k-p to U k,k+P used for compensating for interference of the signals of each OAM mode. A description of p indicating a range of an adjacent mode with respect to a subject mode will be made later

[0028] The channel matrix H between the transmitting OAM mode and the receiving OAM mode is shown in FIG. 4. Here, an example in a case of the OAM modes -2, -1, 0, 1, and 2 is shown. In an ideal situation where the transmitting UCA and the receiving UCA are disposed to completely face each other, signals of the respective OAM modes are orthogonal to each other, and thus only a diagonal block component H k,k is present and remaining components are 0. However, in a situation where inter-mode interference occurs due to displacement, incompleteness of an analog circuit, and the like, the equalization process with respect to the signals of the OAM mode k is performed by using channel matrixes H k,k-p to H k,k+p in which non-diagonal block components that are p adjacent modes are taken into consideration and the weights U k,k-p to U k,k+p calculated by using inverse matrixes thereof. As a method of calculating the weights, a zero-forcing (ZF) method or a minimum mean-square-error (MMSE) method may be used.

[0029] FIG. 5 shows an example of processing by the channel estimation / interference compensation unit 23.

[0030] (1) of FIG. 5 shows a processing example in which the signal r k of the OAM mode k, which is separated by the OAM mode separation units 22-1 to 22-N RX for each of the OAM modes 1 to 5, is input to the channel estimation / interference compensation unit 23, the signal r k is multiplied by the weights U 1,k to U 5,k associated with interference from a subject OAM mode and the other OAM modes, and r _ k = ∑ U k ,m r m where m = 1 to 5 is output as the signal r k of the OAM mode k subjected to inter-mode interference compensation. In other words, the following results are obtained. r _ 1 = U 1 ,1r1 + U1 ,2r2 + U1 ,3r3 + U1 ,4r4 + U1 ,5r5 r _ 2 = U 2 ,1r1 + U2 ,2r2 + U2 ,3r3 + U2 ,4r4 + U2 ,5r5 r _ 3 = U 3 ,1r1 + U3 ,2r2 + U3 ,3r3 + U3 ,4r4 + U3 ,4r4 r _ 4 = U 4 ,1r1 + U4 ,2r2 + U4 ,3r3 + U4 ,4r4 + U4 ,5r5 r _ 5 = U 5 ,1r1 + U5 ,2r2 + U5 ,3r3 + U5 ,4r4 + U5 ,5r5

[0031] In this case, since the computation process is executed on the squared order of an OAM mode multiplex number and the computation amount is exponentially increased according to the multiplex number, the above processing may be applied only to a case where there is no restriction in a circuit scale.

[0032] Here, since the inter-mode interference dominantly depends on interference from an adjacent mode, only a band-shaped channel matrix indicated by a dashed line that is adj acent to a channel matrix of the subject mode indicated by a thick frame in FIG. 4 is taken into consideration. Thus, an approximate process is possible in which the band-shaped channel matrix and a band-shaped approximate weight calculated by the use of an inverse matrix thereof are used for compensating for the interference.

[0033] (2) of FIG. 5 shows a processing example in which the signal r k of the OAM mode k input to the channel estimation / interference compensation unit 23 is multiplied by the weights U k-p,k to U k+p,k associated with interference from the subject OAM mode and adjacent OAM modes adjacent to the subject OAM mode by p, and r _ k = ∑ U k ,m r m where m = k − p to k + p is output as the signals r k of the OAM mode k subjected to inter-mode interference compensation. In the present example, the number p of the adjacent modes considered is 1.

[0034] The adjacent mode selection unit 32-k shown in FIG. 3 selects the signals r k-1 to r k+1 of OAM modes k-1 to k+1 from among the signals r 1 to rs of the OAM modes 1 to 5. The weight multiplication unit 33-k multiplies the signals r k-1 to r k+1 by the weights U k,k-1 to U k,k+1 correspondingly, and outputs r _ k = U k ,k − 1rk − 1 + Uk ,krk + Uk ,k + 1 rk + 1 as the signals rk of the OAM mode k subj ected to inter-mode interference compensation. In other words, when k = 1 to 5, the following results are obtained. r _ 1 = U 1 ,1r1 + U1 ,2r2 r _ 2 = U 2 ,1r1 + U2 ,2r2 + U2 ,3r3 r _ 3 = U 3 ,2r2 + U3 ,3r3 + U3 ,4r4 r _ 4 = U 4 ,3r3 + U4 ,4r4 + U4 ,5r5 r _ 5 = U 5 ,4r4 + U5 ,5r5

[0035] Such the approximate equalization process of the weights can reduce the computation process as shown in (2) of FIG. 5.

[0036] FIG. 6 shows performance characteristics according to the adjacent mode considered.

[0037] In FIG. 6, a transverse axis represents a tilt angle of an antenna, and a longitudinal axis represents a communication path capacity ratio in which a capacity in a case of MMSE equalization with no approximation is set to 100%.

[0038] The present simulation is performed based on the following conditions. Frequency: 28.5 GHz Numbers N TX and N RX of UCAs: four (each ring includes sixteen elements) Multiplex mode: multiplexing of -2 to +2 for each UCA, and a total of twenty modes Transmission distance: 5 m Equalization weight: MMSE

[0039] In this example, if capacity deterioration is allowed up to 20%, a deviation up to 0.1 degrees can be compensated for in a case of p = 0 where the equalization is performed by using only the subject mode, a deviation up to 0.45 degrees can be compensated for in a case of p = 1 where the equalization is performed by using one adjacent mode, and a deviation up to 1 degree can be compensated for in a case of p = 2 where the equalization is performed by using two adjacent modes. In a case of p = 3 where the the equalization is performed by using three adjacent modes with respect to the number of multiplex modes of five, performance close to the full MIMO equalization process shown in (1) of FIG. 5 can be obtained.

[0040] A method of determining the adjacent mode p to be considered is as follows.

[0041] First, at p = 0, a known signal sequence transmitted from the transmitting station, for example, a preamble for channel estimation is multiplied by the approximate weight for the subject mode only. A cost is calculated based on a Euclid distance between the weight multiplication result and a signal point of the known signal sequence, an EVM, an SINR, an absolute value of a difference between an equalizer output when the equalizer is not approximated and an approximate equalizer output, and the like. If the cost calculated here exceeds a threshold value that is preset by a manager or the like, the weight is determined as the weight for the parameter p. If the cost does not exceed the threshold value, p is incremented by +1, an approximate weight including the subject mode and the adjacent mode is determined, and it is determined whether or not the cost calculated similarly exceeds the threshold value. When the performance of the equalization process exceeds a certain predetermined criterion, a value of the previous p is determined as a parameter for approximation of an equalizer, and weight approximation is performed. Then, the received signal is multiplied by the determined weight to compensate for the the inter-mode interference.

[0042] It should be noted that the value of p may be a fixed value that is preset by a system manager in accordance with accuracy of an antenna position calibrator. For example, when the accuracy of the calibrator is ±0.4 degrees and there is a margin of 20% in a system capacity, p is always set to 1 because capacity deterioration of 20% can be permitted.(Another method)

[0043] In the present method, the interference component is approximated by the band matrix, and the equalization is performed by using the weight.

[0044] In the present method, the equalization process can be performed with the computation amount lower than in a case of full MIMO where the equalization process is performed over the all signal sequences. Calculation required for deriving the inverse matrix required for calculating the weight is not necessary as will be described below. Thus, even when an antenna fluctuates due to vibration caused by wind or the like and thus a channel varies, it is not necessary to calculate the inverse matrix every time for calculating the weight.

[0045] As in the characteristics of the channel matrix in the method described with reference to FIGS. 3 to 5, when the antennae are not in the ideal facing disposition, the channel matrix has a band shape and the interference from the adjacent mode is strong. In a method described below, the channel matrix is approximated as a band matrix instead of approximating the weight, and a process of multiplying the approximated channel matrix by the weight is performed.

[0046] First, a band width p is introduced. For example, if p is 1, the channels are approximated to be a band shape as shown in FIG. 7. When the channel matrix H after approximation is indicated by H, an approximate expression of a received signal R with respect to a transmission signal S is R ∼ H _ S .

[0047] A result of multiplying the received signal R by an approximate MMSE weight is expressed as follows by using the signal to noise power ratio SNR and the unit matrix I. R ′ = H _ H H _ + SNR ⋅ I − 1 H HR

[0048] According to the present method, the above can be calculated without calculating the inverse matrix (H H< H+SNR · I) -1< .

[0049] First, the following calculation is performed. H ∗ = H _ H H _ + SNR ⋅ I

[0050] H* is subjected to Cholesky decomposition. H ∗ = QQ H

[0051] Since H* is the band matrix, Q is also a band matrix having the same band width as that of H* and is a lower triangular matrix. Next, the following calculation is performed. z = H _ H R

[0052] In this calculation, it is enough to perform the calculation process with respect to only the band matrix portion. Next, the following equation is solved with respect to f. Qf = z

[0053] Here, since Q is the lower triangular matrix and the band matrix, the equation can be sequentially solved without deriving an inverse matrix of Q. Any algorithm can be used for solving the equation.

[0054] Next, the following equation is solved with respect to R'. Q H R ′ = f

[0055] Here, snice QH is an upper triangular matrix and the band matrix, the equation can be sequentially solved without deriving an inverse matrix of Q. Any algorithm can be used for solving the equation.

[0056] The computation result Q obtained through the Cholesky decomposition may be applied to equalization of subsequent signal sequences. In other words, calculation of H* and the Cholesky decomposition may be performed only once in the interference elimination process, and then only z and f may be processed in computation of a data portion that subsequently arrives.

[0057] In the present method, the interference compensation process is entirely performed on the receiving side The same interference compensation effect can be achieved by using an apparatus performing the same process on the transmitting side.REFERENCE SIGNS LIST

[0058] 10 TRANSMITTING STATION 11 TRANSMITTING UCA 12 OAM MODE GENERATION UNIT 20 RECEIVING STATION 21 RECEIVING UCA 22 OAM MODE SEPARATION UNIT 23 CHANNEL ESTIMATION / INTERFERENCE COMPENSATION UNIT 31 CHANNEL ESTIMATION / WEIGHT CALCULATION UNIT 32 ADJACENT MODE SELECTION UNIT 33 WEIGHT MULTIPLICATION UNIT

Claims

1. An OAM multiplexing communication system configured to use a plurality of orbital angular momentum, OAM, modes and multiplex signals of the plurality of OAM modes, characterized in that the OAM multiplexing communication system comprises: a transmitting station (10) including: a transmitting antenna (11) including a multi-uniform circular array, M-UCA, in which a plurality of uniform circular arrays, UCAs, having different diameters are concentrically disposed, each UCA having a plurality of antenna elements disposed circularly at an equal interval; and an OAM mode generation unit (12) that is configured to simultaneously generate the plurality of OAM modes from each UCA of the transmitting antenna; and a receiving station (20) including: a receiving antenna (21) including a configuration equivalent to the M-UCA; an OAM mode separation unit (22) that is configured to separate signals received by each UCA of the receiving antenna for each of the plurality of OAM modes; and a channel estimation and interference compensation unit (23) that is configured to estimate a channel matrix for the plurality of OAM modes, each OAM mode separated from the signals received by the each UCA of the receiving antenna, and compensate for inter-mode interference between the plurality of OAM modes by using weights calculated from the channel matrix, wherein the OAM modes are consecutively numbered and the channel estimation and interference compensation unit (23) is configured to select, for each OAM mode, signals of a subject mode and a mode with a number adjacent to the number of the subject mode from among the signals of the OAM modes separated by the OAM mode separation unit (22), and to compensate for the inter-mode interference by multiplying an approximated weight calculated by using channel matrixes of the subject mode and the mode with the number adjacent to the number of the subject mode.

2. The OAM multiplexing communication system according to claim 1, wherein the channel estimation and interference compensation unit (23) is configured to determine a range consisting of several modes which are adjacent to the subject mode and which are used in the inter-mode interference compensation, in which a performance of an equalization process using the approximated weight calculated based on a known signal sequence transmitted from the transmitting station (10) satisfies a predetermined criterion.

3. The OAM multiplexing communication system according to claim 1, wherein the channel estimation and interference compensation unit (23) is configured to use an approximated channel matrix obtained by approximating the channel matrix as a band matrix to compensate for the inter-mode interference based on Cholesky decomposition of a matrix including the approximated channel matrix.

4. An inter-mode interference compensation method for an OAM multiplexing communication system that uses a plurality of angular momentum, OAM, modes and multiplexes signals of the plurality of OAM modes, wherein the OAM multiplexing communication system comprises: a transmitting station (10) including: a transmitting antenna (11) that uses a multi-uniform circular array, M-UCA, in which a plurality of uniform circular arrays, UCAs, having different diameters are concentrically disposed, each UCA having a plurality of antenna elements disposed circularly at an equal interval; and an OAM mode generation unit (12) that simultaneously generates the one or more OAM modes from each UCA of the transmitting antenna; and a receiving station (20) including: a receiving antenna (21) that has a configuration equivalent to the M-UCA; an OAM mode separation unit (22) that separates signals received by each UCA of the receiving antenna for each of the plurality of OAM modes; and a channel estimation and interference compensation unit (23) that estimates a channel matrix for each OAM mode separated from the signals received by the each UCA of the receiving antenna, and compensates for inter-mode interference between the OAM modes by using a weight calculated from the channel matrix, wherein the OAM modes are consecutively numbered and the inter-mode interference compensation method comprises: selecting, by the channel estimation and interference compensation unit (23), for each OAM mode, signals of a subject mode and a mode with a number adjacent to the number of the subject mode from among the signals of the OAM modes separated by the OAM mode separation unit (22), and compensating for the inter-mode interference by multiplying an approximated weight calculated by using channel matrixes of the subject mode and the mode with the number adj acent to the number of the subject mode.

5. The inter-mode interference compensation method according to claim 4, wherein the channel estimation and interference compensation unit (23) determines a range consisting of several modes which are adjacent to the subject mode and which are used in the inter-mode interference compensation, in which a performance of an equalization process using the approximated weight calculated based on a known signal sequence transmitted from the transmitting station (10) satisfies a predetermined criterion.

6. The inter-mode interference compensation method according to claim 4, wherein the channel estimation and interference compensation unit (23) uses an approximate channel matrix obtained by approximating the channel matrix as a band matrix to compensate for the inter-mode interference based on Cholesky decomposition of a matrix including the approximated channel matrix.