PHASE NOISE SUPPRESSION METHOD FOR A MULTIPLE-INPUT MULTIPLE-OUTPUT (MIMO) SYSTEM WITH A PLURALITY OF CO-REFERENCE CHANNELS

Abstract

A phase noise suppression method for a multiple-input multiple-output (MIMO) system with a plurality of co-reference channels includes: dividing the phase noise of each channel in the MIMO system into common phase noise and independent phase noise, and constructing a certain number of joint phase states for the independent phase noise; inserting a pilot sequence into the sent signal based on a preset cycle, obtaining the common phase noise based on the pilot at receiver, and performing compensation; and performing signal demodulation on each joint state of the independent phase noise, and comparing the posterior log likelihood values to select the optimal result to output. The above method can significantly improve the phase noise suppression performance of the MIMO system with a plurality of co-reference channels, thereby providing support for improving the system capacity by using MIMO technology.

Claims

1. A phase noise suppression method for a multiple-input multiple-output (MIMO) system with a plurality of co-reference channels, comprising the following steps: S11: dividing phase noise of each channel in the MIMO system into common phase noise and independent phase noise, constructing a certain number of joint phase states for the independent phase noise, and calculating a log likelihood value and a multi-channel decorrelation coefficient matrix that correspond to each joint phase state, wherein step S11 comprises: S111: dividing the phase noise of each channel in the MIMO system into the common phase noise and the independent phase noise, and establishing a signal model under impact of the phase noise in the MIMO system with the plurality of co-reference channels by considering the number of antennas at transmitter as N.sub.t and the number of antennas at receiver as N.sub.r:
r(k)=e.sup.jθ(k)[ΔΘ.sub.r(k)HΔΘ.sub.t(k)Px(k)]+n(k)
θ.sub.t,j(k)=θ.sub.t(k)+Δθ.sub.t,j(k) (1≤j≤N.sub.t)
θ.sub.r,i(k)=θ.sub.r(k)+Δθ.sub.rj(k)(1≤i≤N.sub.r)
θ(k)=θ.sub.t(k)+θ.sub.r(k) wherein θ.sub.t,j (k) represents phase noise of the j.sup.th channel at transmitter at time point k, and comprises common phase noise θ.sub.t(k) and independent phase noise Δθ.sub.t,j(k) of the j.sup.th channel; θ.sub.r,i(k) represents phase noise of the i.sup.th channel at receiver at time point k, and comprises common phase noise θ.sub.r(k) and independent phase noise Δθ.sub.r,i(k) of the i.sup.th channel; θ(k) represents the sum of the common phase noise θ.sub.t(k) at transmitter and the common phase noise θ.sub.r(k) receiver at time point k; x.sub.(k) represents a sent signal of each channel at transmitter at time point k, r(k) represents a received signal of each channel at receiver at time point k, and n (k) represents additive white Gaussian noise at time point k; P represents an N.sub.t×N.sub.t precoding matrix, satisfying P.sup.HP=I.sub.N.sub.t.sub.×N.sub.t, wherein a superscript H represents that conjugate transposition is performed, and I.sub.N.sub.t.sub.×N.sub.t represents an N.sub.t×N.sub.t identity matrix; H represents an N.sub.r×N.sub.t channel response matrix; and ΔΘ.sub.t(k) represents a diagonal matrix formed by an exponent e.sup.jΔθ.sup.i,j.sup.(k) of independent phase noise of each channel at transmitter at time point k, and ΔΘ.sub.r(k) represents a diagonal matrix formed by the exponent e.sup.jΔθ.sup.rj.sup.(k) of independent phase noise of each channel at receiver at time point k; S112: forming a vector Δθ(k)=[Δθ.sub.f,1(k) Δθ.sub.f,2(k) L Δθ.sub.t,N.sub.t(k) Δθ.sub.r,1(k) Δθ.sub.r,2(k) L Δθ.sub.r,N.sub.r(k)].sup.T=[Δθ.sub.1(k) Δθ.sub.2(k) L Δθ.sub.N.sub.t.sub.+N.sub.r(k)].sup.T by the independent phase noise of each channel at transmitter and the independent phase noise of each channel at receiver, and quantizing each element Δθ.sub.m(k) (1≤m≤N.sub.t+N.sub.r) in the vector into Z.sub.m phase states that are denoted as a phase state φ.sub.m=[φ.sub.m,1 φ.sub.m,2 L φ.sub.m,Z.sub.m].sup.T, wherein φ.sub.m,1<φ.sub.m,2<L<φ.sub.m,Z.sub.m; S113: selecting a state from each element based on the phase state of each element of the independent phase noise to form a joint state ψ.sub.z (1≤z≤Z), wherein $Z=\underset{m=1}{\overset{N_t+N_r}{.Math.}}{Z}_m$ represents a total number of joint states; S114: calculating the log likelihood value corresponding to each joint state ψ.sub.z according to the following formula: ${\lambda }_z=\log f\left({\psi }_z\right)=\stackrel{N_t+N_r}{\underset{m=1}{.Math.}}\log f\left({\psi }_{z,m}\right)=-\underset{m=1}{\overset{N_t+N_r}{.Math.}}\frac{{\psi }_{z,m}^2}{2{\sigma }_m^2}+\upsilon$ wherein υ represents a common coefficient, λ.sub.z represents the log likelihood value corresponding to ψ.sub.z, f(.Math.) represents a likelihood function, ψ.sub.z,m represents the m.sup.th element in the joint state ψ.sub.z, and σ.sub.m.sup.2 represents a variance; and S115: calculating the multi-channel decorrelation coefficient matrix corresponding to each joint state ψ.sub.z according to the following formula:
W.sub.z=(H.sub.z.sup.NH.sub.z+σ.sub.n.sup.2I.sub.N.sub.t.sub.×N.sub.t).sup.−1H.sub.z.sup.N wherein W.sub.z represents the multi-channel decorrelation coefficient matrix corresponding to ψ.sub.z, H.sub.z represents an equivalent response matrix corresponding to ψ.sub.z, and σ.sub.n.sup.2 represents power of the additive white Gaussian noise; S12: inserting a pilot sequence into the sent signal based on a preset cycle, estimating common phase noise based on the pilot at receiver, obtaining common phase noise of effective data other than the pilot through interpolation fitting, and performing compensation, wherein S12 comprises: S121: inserting the pilot sequence into the sent signal based on the preset cycle, and. establishing a signal model of common phase noise between adjacent pilots:
r(αP+P)=e.sup.j[{circumflex over (θ)}(αP)+Δθ(αP)]H(αP+P)x(αP+P)+n(αP+P)
θ(αP+P)=θ(αP)+Δθ(αP) wherein pilot density is 1/P that is, the pilot is corresponded when k=αP, and the effective data is corresponded at other time points, wherein α is an integer; r(αP+P) represents a received signal at pilot position αP+P, θ(αP) and {circumflex over (θ)}(αP) represent common phase noise at pilot position αP and the corresponding estimated value respectively, and Δθ(αP) represents a change value of common phase noise at pilot position αP+P relative to common phase noise at pilot position αP ; and x(αP+P) represents a known sent signal at the pilot position αP+P, n(αP+P) represents additive white Gaussian noise at pilot position αP+P, and H(αP+P) represents an unknown equivalent response matrix; S122, for each joint state of the independent phase noise, estimating a change value of the common phase noise between the adjacent pilots, and calculating the corresponding log likelihood value according to the following formula:
Δθ.sub.z(αP)=−∠[e.sup.j{circumflex over (θ)}(αP)r.sup.H(αP+P)H.sub.zx(αP+P)], (1≤z≤Z) wherein Δθ.sub.z(αP) represents an initial estimated value of the change value of the z.sup.th piece of common phase noise at pilot position αP+P relative to position αP initial estimated values of change values of Z pieces of common phase noise are obtained in total, and Z represents the total number of joint states; S123: obtaining log likelihood values corresponding to the Z initial estimated values with respect to Δθ(αP) based on the change value of the common phase noise, the log likelihood value of each joint state of the independent phase noise, and an error of the received signal ${\overline{\Lambda }}_z\left(\mathrm{aP}\right)={\lambda }_z-\frac{\Delta {\theta }_z^2\left(\mathrm{aP}\right)}{2{\sigma }_P^2}-\frac{{.Math.r(\left(\mathrm{aP}+P\right)-e^{j\left[\hat{\theta }\left(\mathrm{aP}\right)+\Delta {\theta }_z\left(\mathrm{aP}\right)\right]}{H}_{z}x\left(\mathrm{aP}+P\right).Math.}^2}{{\sigma }_n^2},\left(1\le z\le Z\right)$ wherein |.Math.|.sup.2 represents a square sum of each element module, σ.sub.P.sup.2 represents a variance of Δθ(αP), H.sub.z represents the equivalent response matrix corresponding to ψ.sub.z, σ.sub.n.sup.2 represents the power of the additive white Gaussian noise, λ.sub.z represents the log likelihood value corresponding to each joint phase state, and ∧.sub.z (αP) represents the log likelihood value corresponding to the z.sup.th initial estimated value Δθ.sub.z(αP) with respect to Δθ(αP); comparing the Z log likelihood values, selecting the maximum log likelihood value denoted as z′.sub.αP=arg max ∧.sub.z(αP), and then obtaining the final estimated value of the common phase noise at pilot position αP+P according to the following formula:
{circumflex over (θ)}(αP+P)={circumflex over (θ)}(αP)+Δθ.sub.z′.sub.αP(αP) obtaining the common phase noise of the effective data through interpolation; and S124: suppressing the common phase noise through phase reversal:
r′(k)=e.sup.−jθ(k)r(k) wherein {circumflex over (θ)}(k) represents an estimated value of the common phase noise at time point k after interpolation, r(k) represents the received signal of each channel at receiver at time point k before suppression, and r′(k) represents a result obtained after the common phase noise is suppressed; and S13: for each joint phase state of the independent phase noise, based on the multi-channel decorrelation coefficient matrix, performing signal demodulation on a result Obtained after compensating for the common phase noise in S12, and calculating a posterior log likelihood value; and taking a joint phase state with the maximum posterior log likelihood value as estimation result of the independent phase noise of the plurality of channels, and taking the corresponding demodulation result as a final output; wherein step S13 comprises: S131: calculating, according to the following formula, a decorrelation result based on the multi-channel decorrelation coefficient matrix W.sub.z corresponding to each joint state of the independent phase noise and the result r′(k) obtained after the common phase noise is suppressed:
y.sub.z(k)=W.sub.zr′(k), (1≤z≤Z) wherein y.sub.z (k) represents a decorrelation result corresponding to the z.sup.th joint state; S132: performing decision demodulation on the decorrelation result to obtain the demodulation result %.sub.z, wherein 1≤z≤Z; S133: for the demodulation result %.sub.z (k) corresponding to each joint state, calculating the posterior log likelihood value according to the following formula: ${\Lambda }_z\left(k\right)={\lambda }_z-\frac{{.Math.r^{\prime }\left(k\right)-H_z{x}_z\left(k\right).Math.}^2}{{\sigma }_n^2}$ wherein λ.sub.z represents the log likelihood value corresponding to each joint phase state, H.sub.z represents the equivalent response matrix corresponding to ψ.sub.z, σ.sub.n.sup.2 represents the power of the additive white Gaussian noise, and ∧.sub.z (k) represents the posterior log likelihood value corresponding to each joint phase state; and S134: selecting the maximum posterior log likelihood value denoted as z′.sub.k=arg max ∧.sub.z(k), and then obtaining the estimated value ψ.sub.z′.sub.z of the independent phase noise at time point k, and a final demodulation output structure {circumflex over (x)}(k)=%.sub.z′.sub.k(k).

2. The phase noise suppression method for the MIMO system with the plurality of co-reference channels according to claim 1, wherein a quantization criterion in step S112 is as follows: $\begin{array}{cc}\underset{{\phi }_m}{\min }{\int }_{-\infty }^{+\infty }{\underset{z\left(1\le z\le Z_m\right)}{\min }\left(x-{\phi }_{m,z}\right)}^2{f}_m\left(x\right)\mathrm{dx}=\underset{{\phi }_m}{\min }\underset{z=1}{\stackrel{Z_m}{.Math.}}{\int }_{\frac{{\phi }_{m,z-1}+{\phi }_{m,z}}{2}}^{\frac{{\phi }_{m,z}+{\phi }_{m,z+1}}{2}}{\left(x-{\phi }_{m,z}\right)}^2{f}_m\left(x\right)\mathrm{dx}& \end{array}$ wherein φ.sub.m,0=−∞, φ.sub.m,Z.sub.m.sub.+1=+∞, and f.sub.m(x) represents a distribution function of an element Δθ.sub.m(k); the minimization objective in the quantization criterion is a piecewise integral function, and solving is performed through offline search based on the given variance σ.sub.m.sup.2 to obtain a value of each discrete phase state φ.sub.m,z(1≤z≤Z.sub.m).

3. The phase noise suppression method for the MIMO system with the plurality of co-reference channels according to claim 1, wherein an interpolation formula in step S123 is as follows: $\hat{\theta }\left(\alpha P+p\right)=\frac{P-p}{P}\hat{\theta }\left(aP\right)+\frac{p}{P}\hat{\theta }\left(\mathrm{aP}+P\right),\left(1\le p\le P-1\right)$ wherein αP+p represents the p.sup.th interpolation point between αP and αP+P, and {circumflex over (θ)}(αP+p) represents an interpolation result.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0015] FIG. 1 is a flowchart of a phase noise suppression method for an MIMO system with a plurality of co-reference channels according to the present disclosure;

[0016] FIG. 2 is a flowchart of step 1 according to the present disclosure;

[0017] FIG. 3 is a schematic diagram of a system model under impact of phase noise in the MIMO system with a plurality of co-reference channels;

[0018] FIG. 4 is a flowchart of step 2 according to the present disclosure;

[0019] FIG. 5 is a schematic diagram of an implementation of obtaining common phase noise through pilot insertion and interpolating fitting according to the present disclosure; and.

[0020] FIG. 6 compares BER performance in the technical solutions of the present disclosure and existing technical solutions.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0021] The present disclosure is further described below with reference to the accompanying drawings and embodiments.

[0022] The present disclosure provides a phase noise suppression method for an MIMO system with a plurality of co-reference channels. As shown in PG. 1, the method includes: dividing phase noise of each channel in the MIMO system into common phase noise and independent phase noise, and constructing a certain number of joint phase states for the independent phase noise; inserting a pilot sequence into the sent signal based on a preset cycle, obtaining the common phase noise based on the pilot at receiver, and performing compensation; and performing signal demodulation on each joint state of the independent phase noise, and comparing posterior log likelihood values to select the optimal result to output. The above method can significantly improve phase noise suppression performance of the MIMO system with the plurality of co-reference channels, thereby providing support for improving the system capacity by using MEVIO technology.

[0023] A specific implementation process includes the following steps:

[0024] Step 1: Phase noise of each channel in the MIMO system with a plurality of co-reference channels is divided into common phase noise and independent phase noise, a certain number of joint phase states is constructed for the independent phase noise, and a log likelihood value and a multi-channel decorrelation coefficient matrix that correspond to each joint phase state are calculated, as shown in FIG. 2.

[0025] FIG. 3 shows the system model under impact of the phase noise in the MIMO system with the plurality of co-reference channels. Considering that the number of antennas at transmitter and the number of antennas at receiver are N.sub.t and N.sub.r respectively, it may be assumed that the number of signal streams is equal to N.sub.t and N.sub.t≤N.sub.r, and the received signal under the impact of phase noise of the plurality of channels is expressed as formula (1):

r(k)=Θ.sub.r(k)HΘ.sub.t(k)Px(k)+n(k)   (1)

[0026] Herein, x(k)=[x.sub.1(k) x.sub.2(k) L x.sub.N.sub.i(k)].sup.T represents a sent signal of each channel at time point k, and the superscript T represents that transposition is performed; P represents an N.sub.t×N.sub.t precoding matrix, satisfying P.sup.HP=I.sub.N.sub.t.sub.×N.sub.t, where the superscript H represents that conjugate transposition is performed, and I.sub.N.sub.t.sub.×N, represents an N.sub.t×N.sub.t identity matrix;

[00001]${\Theta }_t\left(k\right)=\mathrm{diag}\left\{{\left[\begin{array}{cccc}{e}^{j{\theta }_{t,1}\left(k\right)}& e^{j{\theta }_{t,2}\left(k\right)}& L& e^{j{\theta }_{t,N_t}\left(k\right)}\end{array}\right]}^T\right\}$

represents a diagonal matrix formed by e.sup.jθ.sup.t.sup.,(k), namely, the exponent of phase noise at transmitter of each channel at time; point k; H represents an N.sub.r×N.sub.t channel response matrix, where the element in the i.sup.th row and the j.sup.th column, namely, h.sub.i,j, represents the channel gain from the j.sup.th transmitting antenna to the i.sup.th receiving antenna;

[00002]${\Theta }_r\left(k\right)=\mathrm{diag}\left\{{\left[\begin{array}{cccc}{e}^{j{\theta }_{r,1}\left(k\right)}& e^{j{\theta }_{r,2}\left(k\right)}& L& e^{j{\theta }_{r,N_r}\left(k\right)}\end{array}\right]}^T\right\}$

represents a diagonal matrix formed by e.sup.jθ.sup.r.sup.,(k), namely, the exponent of phase noise at receiver of each channel at time point k; and

[00003]$r\left(k\right)={\left[\begin{array}{cccc}{r}_1\left(k\right)& r_2\left(k\right)& L& {r_{N_r}\left(k\right)]}^T\mathrm{and}n\left(k\right)\begin{array}{cccc}=[n_1\left(k\right)& n_2\left(k\right)& L& n_{N_r}\left(k\right)\end{array}\end{array}\right]}^T$

represent a received signal and additive white Gaussian noise of each channel at time point k respectively.

[0027] In the case of co-reference, phase noise of each channel at transmitter and receiver includes common phase noise and independent phase noise, where the common phase noise is from a reference signal, and the independent phase noise is additionally generated in the frequency multiplication process. In other words, the phase noise of each channel at transmitter and the phase noise of each channel at receiver can be expressed as θ.sub.t,j(k)=θ.sub.t(k)+Δθ.sub.t,k(k) (1≤j≤N.sub.t) and θ.sub.r,i(k)=θ.sub.r(k)+Δθ.sub.r,i(k) (1≤i≤N.sub.r) respectively, where θ.sub.t,j(k) and θ.sub.r,i(k) represent phase noise of the j.sup.th channel at transmitter and phase noise of the i.sup.th channel at receiver at time point k respectively, θ.sub.t(k) and θ.sub.r(k) represent the common phase noise of each channel at transmitter and the common phase noise of each channel at receiver at time point k respectively, and Δθ.sub.t,j(k) and Δθ.sub.r,i(k) represent the independent phase noise of the j.sup.th channel at transmitter and the independent phase noise of the i.sup.th channel at receiver at time point k respectively. Therefore, formula (1) can be further expressed as formula (2):

r(k)=e.sup.jθ(k)[ΔΘ.sub.r(k)HΔΘ.sub.t(k)Px(k)]+n(k)   (2)

[0028] In the above formula, θ(k)=θ.sub.t(k)+θ.sub.r(k) represents the sum of the common phase noise at transmitter and the common phase noise at receiver at time point k. After loop filtering with a low-pass characteristic in a phase-locked loop, the common phase noise is mainly distributed in a low-frequency region, in other words, the common phase noise between adjacent symbols in the time domain is strongly correlated, After the common phase noise at pilot is estimated by periodically inserting the pilot, the common phase noise of effective data can be obtained through interpolation fitting, and then corresponding compensation can be performed.

[00004]${\mathrm{\Delta \Theta }}_t\left(k\right)=\mathrm{diag}\left\{{\left[\begin{array}{cccc}{e}^{j{\mathrm{\Delta \theta }}_{t,1}\left(k\right)}& e^{j{\mathrm{\Delta \theta }}_{t,2}\left(k\right)}& L& e^{j\Delta {\theta }_{t,N_t}\left(k\right)}\end{array}\right]}^T\right\}$

represents a diagonal matrix formed by e.sup.jΔθ.sup.t,j.sup.(k), namely, the exponent of the independent phase noise of each channel of the transmitting antenna at time point k; and

[00005]${\mathrm{\Delta \Theta }}_r\left(k\right)=\mathrm{diag}\left\{{\left[\begin{array}{cccc}{e}^{j{\mathrm{\Delta \theta }}_{r,1}\left(k\right)}& e^{j\Delta {\theta }_{r,2}\left(k\right)}& L& e^{j{\mathrm{\Delta \theta }}_{r,N_r}\left(k\right)}\end{array}\right]}^T\right\}$

represents a diagonal matrix formed by e.sup.jΔθ.sup.r,j.sup.(k), namely, the exponent of the independent phase noise of each channel of the receiving antenna at time point k. Different from the common phase noise, the independent phase noise of each channel is mainly distributed in a high-frequency region, and fast changes. lindependent phase noise between adjacent symbols are poorly correlated, which are like white Gaussian noise and cannot be accurately estimated through conventional interpolation fitting or by using a prediction algorithm.

[0029] The following formula is obtained:

[00006]$\begin{array}{cccccccc}\mathrm{\Delta \theta }\left(k\right)=[{\mathrm{\Delta \theta }}_{t,1}\left(k\right)& {\mathrm{\Delta \theta }}_{t,2}\left(k\right)& L& {\mathrm{\Delta \theta }}_{t,N_t}\left(k\right)& {\mathrm{\Delta \theta }}_{r,1}\left(k\right)& {\mathrm{\Delta \theta }}_{r,2}\left(k\right)& L& {{\mathrm{\Delta \theta }}_{r,N_r}\left(k\right)]}^T\\ =[{\mathrm{\Delta \theta }}_1\left(k\right)& {\mathrm{\Delta \theta }}_2\left(k\right)& L& {{\mathrm{\Delta \theta }}_{N_t+N_r}\left(k\right)]}^T& & & & \end{array}$

[0030] The present disclosure quantizes each element of the independent phase noise, in other words, it is considered that each element is distributed only in a. few phase states. It is assumed that the distribution function of the m.sup.th element Δθ.sub.m(k) (1≤m≤N.sub.t+N.sub.t) in Δθ(k) is f.sub.m(x), and Δθ.sub.m(k) is quantized into Z.sub.m phase states, which are denoted as

[00007]${\phi }_m=[\begin{array}{cccc}{\phi }_{m,1}& {\phi }_{m,2}& L& {{\phi }_{m,Z_m}]}^T\end{array},$

where φ.sub.m,l<φ.sub.m,2<L<φ.sub.m,Z.sub.m. A quantization criterion is as follows:

[00008]$\begin{array}{cc}\underset{{\phi }_m}{\min }{\int }_{-\infty }^{+\infty }{\underset{z\left(1\le z\le Z_m\right)}{\min }\left(x-{\phi }_{m,z}\right)}^2{f}_m\left(x\right)\mathrm{dx}& \left(3\right)\end{array}$$=\underset{{\phi }_m}{\min }\underset{z=1}{\overset{Z_m}{.Math.}}{\int }_{\frac{{\phi }_{m,z-1}+{\phi }_{m,z}}{2}}^{\frac{{\phi }_{m,z}+{\phi }_{m,z+1}}{2}}{\left(x-{\phi }_{m,z}\right)}^2{f}_m\left(x\right)\mathrm{dx}$

[0031] In the above formula, φ.sub.m,0=−∞, φ.sub.m,Z.sub.m.sub.+1=+∞.The independent phase noise is like the white Gaussian noise. A distribution interval of the white Gaussian noise is (−∞, +∞), and the periodic phase element should be distributed between −π and π. However, the independent phase noise is small and only fluctuates around zero. Therefore, the integral region between −∞ and +∞ is also equivalent to that between −π and π. The minimization objective in formula (3) is a piecewise integral function. To satisfy the minimization, each discrete phase state φ.sub.m,Z (1≤z≤Z.sub.m) needs to be equal to the mean value of the corresponding integral segment. It may be assumed that Δθ.sub.m(k) follows the zero-mean Gaussian distribution with a variance being σ.sub.m.sup.2. In this case, each discrete phase state can be obtained according to the following formula:

[00009]$\begin{array}{cc}{\phi }_{m,z}=\frac{{\int }_{\frac{{\phi }_{m,z-1}+{\phi }_{m,z}}{2}}^{\frac{{\phi }_{m,z}+{\phi }_{m,z+1}}{2}}{\mathrm{xf}}_m\left(x\right)\mathrm{dx}}{{\int }_{\frac{{\phi }_{m,z-1}+{\phi }_{m,z}}{2}}^{\frac{{\phi }_{m,z}+{\phi }_{m,z+1}}{2}}{f}_m\left(x\right)\mathrm{dx}}=\sqrt{\frac{2}{\pi }}{\sigma }_m\frac{\exp \left\{-\frac{{\left({\phi }_{m,z-1}+{\phi }_{m,z}\right)}^2}{8{\sigma }_m^2}\right\}-\exp \left\{-\frac{{\left({\phi }_{m,z}+{\phi }_{m,z+1}\right)}^2}{8{\sigma }_m^2}\right\}}{\mathrm{erf}\left(\frac{{\phi }_{m,z}+{\phi }_{m,z+1}}{2\sqrt{2}{\sigma }_m}\right)-\mathrm{erf}\left(\frac{{\phi }_{m,z-1}+{\phi }_{m,z}}{2\sqrt{2}{\sigma }_m}\right)}& \left(4\right)\end{array}$

[0032] In the above formula,

[00010]$\mathrm{erf}\left(x\right)=\frac{2}{\sqrt{\pi }}{\int }_0^x{e}^{-{\eta }^2}d\eta ,$

and exp(x)=e.sup.x. Formula (4) provides a method for solving based on φ.sub.m,Z+1 based on φ.sub.m,Z−1 and φ.sub.m,Z. Therefore, when a given variance σ.sub.m.sup.2 is provided, the value of each discrete phase state φ.sub.m,z (1≤z≤Z.sub.m) accurately calculated offline by searching φ.sub.m,1. Table 1 gives the values of φ.sub.m,z/σ.sub.m that correspond to different numbers of phase states.

TABLE-US-00001 TABLE 1 Discrete values corresponding to different numbers of phase states Z.sub.m φ.sub.m, z/σ.sub.m (1 ≤ z ≤ Z.sub.m) 2 −0.798, 0.798 3 −1.224, 0, 1.224 4 −1.510, −0.453, 0.453, 1.510 5 −1.724, −0.765, 0, 0.765, 1.724 6 −1.894, −1.000, −0.318, 0.318, 1.000, 1.894 7 −2.033, −1.188, −0.561, 0, 0.561, 1.188, 2.033 8 −2.152, −1.344, −0.756, −0.245, 0.245, 0.756, 1.344, 2.152

[0033] After all N.sub.t+N.sub.r elements of the independent phase noise are quantized, a state is selected from each element to construct a joint state. The total number of joint states is

[00011]$Z=\underset{m=1}{\overset{N_i+N_r}{.Math.}}{Z}_m.$

Herein, the number Z.sub.m of quantized states of each element may be the same or different, depending on the valiance of each element.

[0034] For convenience of description, the element of the z.sup.th selection from the phase state vector φ.sub.m, corresponding to the element Δθ.sub.m(k) is denoted as ψ.sub.z,m (where ψ.sub.z,m is an element in φ.sub.m (1≤m≤N.sub.t+N.sub.r)) to construct the z.sup.th joint state, which is denoted as

[00012]${\psi }_z=[\begin{array}{cccc}{\psi }_{z,1}& {\psi }_{z,2}& L& {{\psi }_{z,N_t,N_r}]}^T.\end{array}$

[0035] A boa likelihood value corresponding to the joint state ψ.sub.z is expressed as follows:

[00013]$\begin{array}{cc}{\lambda }_z=\log f\left({\psi }_z\right)=\underset{m=1}{\overset{N_t+N_r}{.Math.}}\log f\left({\psi }_{z,m}\right)=-\underset{m=1}{\overset{N_t+N_r}{.Math.}}\frac{{\psi }_{z,m}^2}{2{\sigma }_m^2}+\upsilon & \left(5\right)\end{array}$

[0036] in the above formula, υ represents a common coefficient that can be ignored, and f (.) represents a likelihood function.

[0037] When the common phase noise is corrected, according to formula (2), a multi-channel decorrelation coefficient matrix corresponding to the joint state ψ.sub.z of the independent phase noise can be obtained based on the minimum mean square error ° MESE) algorithm, namely

[00014]$\begin{array}{cc}{w}_z={\left(H_z^H{H}_z+{\sigma }_n^2{I}_{N_i\times N_i}\right)}^{-1}{H}_z^H& \left(6\right)\end{array}$

[0038] In the above formula,

[00015]$H_z=\mathrm{diag}\left\{[\begin{array}{cccc}{e}^{j{\psi }_{z,N_l+1}}& e^{j{\psi }_{z,N_l+2}}& L& {e^{j{\psi }_{z,N_l+N_r}}]}^T\end{array}\right\}H\mathrm{diag}\{[\begin{array}{cccc}{e}^{j{\psi }_{z,1}}& e^{j{\psi }_{z,2}}& L& {e^{j{\psi }_{z,N_l}}]}^T\}P\end{array}$

represents an equivalent response matrix corresponding to ψ.sub.z, and σ.sub.n.sup.2 represents power of the additive white Gaussian noise.

[0039] In the present disclosure, the joint state ψ.sub.z of the independent phase noise of the plurality of channels, the log likelihood value λ.sub.z, and the multi-channel decorrelation coefficient matrix

[00016]$W_z\left(1\le z\le Z=\underset{m=1}{\overset{N_t+N_r}{.Math.}}{Z}_m\right)$

are independent of the symbol sequence number k and can be calculated offline. A key to determine the complexity is the selected number of quantized states Z.sub.m on each element of independent phase noise. A larger value of Z, corresponds to a smaller quantization error, namely, a smaller suppression error of the independent phase noise, but results in a larger number of joint states Z, namely, higher complexity. Therefore, it is necessary to select the number of quantized states Z.sub.m on each element of the independent phase noise under a tradeoff between performance and complexity.

[0040] Step 2: A pilot sequence is inserted into the sent signal based on a preset cycle, the common phase noise at the pilot is estimated at receiver, the common phase noise of the effective data is obtained through interpolation fitting, and suppression is performed, as shown in FIG. 4,

[0041] During suppression of the phase noise of the plurality of channels, the common phase noise needs to be suppressed first. In the present disclosure, the common phase noise at the pilot is estimated based on the pilot, then the common phase noise of the effective data is obtained through interpolation, and then compensation is performed, as shown in FIG. 5. The pilot density is 1/P, that is, the pilot is corresponded when k=αP (α is an integer), and the effective data is corresponded at other time points. The change of the common phase noise at the pilot can be regarded as a Wiener process, as shown in the following formula (7), where Δθ(αP) represents a change value of the common phase noise at the position αP+P relative to the position αP, which follows the zero-mean Gaussian distribution, and its variance is denoted as a σ.sub.P.sup.2. An estimation process of the common phase noise can be equivalent to an estimation process of the change value Δθ(αP).

θ(αP+P)=θ(αP)+Δθ(αP)

[0042] Assuming that an estimated value of the common phase noise at the position αP is {circumflex over (θ)}(αP), by denoting H(k)=ΔΘ.sub.r(k)HΔΘ.sub.t(k)P, the received symbol of the pilot at the position αP+P is r(αP+P)=e.sup.j[θ(αP)+Δθ(αP)] H(αP+P)x(αP+P)+n(αP+P), where x(αP +P) is a known symbol and H(αP+P) is unknown. However, corresponding Δθ(αP) can be estimated and solved for each joint state constructed by the independent phase noise. A result shown in the following formula (8) can be obtained based on the least square (LS) algorithm, where the operator ∠ represents that the phase is taken from a complex data:

[00017]$\begin{array}{cc}{\mathrm{\Delta \theta }}_z\left(aP\right)=-\angle \left[e^{j\hat{\theta }\left(\mathrm{aP}\right)}{r}^H\left(\mathrm{aP}+P\right)H_{z}x\left(\mathrm{aP}+P\right)\right],\left(1\le z\le Z\right)& \left(8\right)\end{array}$

[0043] In the above formula, Δθ.sub.z (αP) represents an initial estimated value of the change value of the zth piece of common phase noise at the pilot position αP +P relative to the position αP. After initial estimated values of change values of Z pieces of common phase noise are obtained in total, the optimal estimation result needs to be selected through comparison. The present disclosure obtains log likelihood values corresponding to the Z initial estimated values with respect to Δθ(αP) based on the change value of the common phase noise, the log likelihood value of each joint state of the independent phase noise, and the error of the received signal, as shown in the following formula:

[00018]$\begin{array}{cc}{\overline{\Lambda }}_z\left(\mathrm{aP}\right)={\lambda }_z-\frac{\Delta {\theta }_z^2\left(\mathrm{aP}\right)}{2{\sigma }_P^2}-\frac{{.Math.r(\left(\mathrm{aP}+P\right)-e^{j\left[\hat{\theta }\left(\mathrm{aP}\right)+\Delta {\theta }_z\left(\mathrm{aP}\right)\right]}{H}_{z}x\left(\mathrm{aP}+P\right).Math.}^2}{{\sigma }_n^2},\left(1\le z\le Z\right)& \left(9\right)\end{array}$

[0044] In the above formula, |.Math.|.sup.2 represents the square sum of each element module. After that, the maximum log likelihood value is selected and denoted as z′.sub.αP=arg max ∧.sub.z(αP). Therefore, an estimated value of the common phase noise at the pilot position αP+P is obtained according to the following formula:

.sub.z{circumflex over (θ)}(αP+P)={circumflex over (θ)}(αP)+Δθ.sub.z′.sub.op(αP)   (10)

[0045] Herein, an estimate (corresponding to α=−1) of the common phase noise of the initial pilot with k=0 can be considered as {circumflex over (θ)}(0)=Δθ.sub.z′.sub.op(0).

[0046] After the estimated value of the common phase noise at the pilot position is obtained, the corresponding value at the effective data can be obtained through interpolation. For example, a linear interpolation process can be expressed as follows:

[00019]$\begin{array}{cc}\hat{\theta }\left(\alpha P+p\right)=\frac{P-p}{P}\hat{\theta }\left(aP\right)+\frac{p}{P}\hat{\theta }\left(\mathrm{aP}+P\right),\left(1\le p\le P-1\right)& \left(1\right)\end{array}$

[0047] In the above formula, αP+p represents the p.sup.th interpolation point between αP and αP+P.

[0048] Through interpolation, the common phase noise at time point k can be obtained, and then the common phase noise can be suppressed through phase reversal, as shown in formula (12):

r′(k)=e.sup.−j{circumflex over (θ)}(k)r(k)   (12)

[0049] In the above formula, r′(k) represents a result obtained after the common phase noise is suppressed, and {circumflex over (θ)}(k) represents an estimated value of the common phase noise at time point k.

[0050] Step 3: For each joint phase state of the independent phase noise, signal demodulation is performed, and a posterior log likelihood value is calculated; and the joint phase state with the maximum posterior log likelihood value is taken as an estimation result of the independent phase noise of the plurality of channels, and the corresponding demodulati©n result is taken as the final output.

[0051] After the common phase noise is compensated, according to the following formula, a decorrelation result can be obtained based on the multi-channel decorrelation coefficient matrix W.sub.z corresponding to each joint state of the independent phase noise:

y.sub.z(k)=W.sub.zr′(k), (1≤z≤Z, k≠αP)   (13)

[0052] In the above formula, y.sub.z(k) represents a decorrelation result corresponding to the z.sup.th joint state ψ.sub.z. Then, decision demodulation is performed on y.sub.z (k) to obtain x0.sub.z(k).

[0053] For the demodulation result (k) corresponding to each joint state ((1≤z≤Z, k≠αP)), a posterior log likelihood value is calculated according to the following formula:

[00020]$\begin{array}{cc}{\Lambda }_z\left(k\right)={\lambda }_z-\frac{{.Math.r^{\prime }\left(k\right)-H_z{x}_z\left(k\right).Math.}^2}{{\sigma }_n^2},\left(1\le z\le Z,k\ne \mathrm{aP}\right)& \left(14\right)\end{array}$

[0054] After that, the maximum posterior log likelihood value is selected and denoted as z′.sub.k=arg max ∧.sub.z(k). Therefore, the estimated value of the independent phase noise at time point k is ψ.sub.z′.sub.r, and the demodulation output is {circumflex over (x)}(k)=x0.sub.′.sub.k(k).

[0055] The present disclosure divides the phase noise of each channel in the MIMO system into the common phase noise and the independent phase noise for processing. Compared with the existing methods for the case in which independent local oscillators are used or the case in which a local oscillator is shared, the present disclosure is more suitable for the actual co-reference architecture, and can ensure better phase noise suppression performance, thereby providing support for improving the transmission capacity through spatial multiplexing.

[0056] To further explain a beneficial effect of the present disclosure compared with the prior art, a performance simulation example is provided herein. The number of transmitting antennas and the number of receiving antennas are the same, namely, N.sub.t=N.sub.r=2 the communication frequency is set to 15 GHz, the transmission distance is 1.6 Km, and the spacings between two antennas at receiver and transmitter are both 3 m (equal to 0.75 times of the optimal antenna spacing), the normalized channel response matrix is shown in the following formula (15), a preceding matrix is obtained through singular value decomposition, and the result is shown in the following formula (16). After singular value decomposition, the demodulation performance of each channel in the MLMO system is different. Herein, the modulation modes of two channels are configured as 64 QAM and 16 QAM respectively. A Wiener model is set for the common phase noise at transmitter and receiver, where the power spectral density is −90 dBc/Hz at 100 kHz from the central frequency, the white Gaussian noise mode is set for the independent phase noise, and the standard deviation is 2 degree. The number of quantized states of each element of the independent phase noise is Z.sub.m=4(m=1,2,3,4), the value of each state can be obtained based on Table 1, and the corresponding number of joint states is 256. In addition, the pilot spacing is set to 32.

[00021]$\begin{array}{cc}H=\left[\begin{array}{cc}0.7071& 0.4486+0.5466i\\ 0.4486+0.5466i& 0.7071\end{array}\right]& \left(15\right)\end{array}$$\begin{array}{cc}P=\left[\begin{array}{cc}-0.7071& 0.7071\\ -0.7071& -0.7071\end{array}\right]& \left(16\right)\end{array}$

[0057] FIG. 6 compares the BER performance of the two channels in the phase noise suppression method for the case in which independent local oscillators are assumed, the phase noise suppression method for the case in which a shared local oscillator is assumed, and the technical solutions in the present disclosure, where different SNRs are used, and no error correction code is considered. It can be seen from the figure that the technical solutions of the present disclosure significantly improve the BER performance of the two channels. A main problem of the phase noise suppression method for the case in which independent local oscillators are assumed is that internal characteristics of the co-reference architecture cannot be fully mined. The phase noise of each channel is considered as completely independent, while the phase noise suppression method for the case in which a shared local oscillator is assumed cannot process the independent phase noise. For channel 1, the BER performance in the phase noise suppression method for the case in which independent local oscillators are assumed and the phase noise suppression method for the case in which a shared local oscillator is assumed has a high error floor, which makes it difficult to ensure reliability. In contrast, the technical solutions of the present disclosure can reduce the error floor by more than 10 times and can ensure high reliability after adopting an error correction code. Channel 2 has a lower modulation order and has better BER performance than channel 1. However, the technical solutions of the present disclosure still have great advantages, and at BER=10.sup.−4, the gains are about 4 dB and 1 dB respectively compared with the phase noise suppression method for the case in which independent local oscillators are assumed and the phase noise suppression method for the case in which a shared local oscillator is assumed.

[0058] The examples listed above are merely specific embodiments of the present disclosure. Apparently, the present disclosure is not limited to the above embodiments, and many variations are possible. All deformations that can be directly derived or associated by those skilled in the art from the present disclosure should be deemed as falling within the protection scope of the present disclosure.