Blind method of equalizing signals in filter bank multi-carrier communications
11483183 · 2022-10-25
Assignee
Inventors
- Muhammad Moinuddin (Jeddah, SA)
- Ubaid M. Al-Saggaf (Jeddah, SA)
- Jawwad Ahmad (Karachi, PK)
- Asmaa Ubaid Al-Saggaf (Jeddah, SA)
- Mohammad S. Alhazmi (Makkah, SA)
Cpc classification
H04L27/2698
ELECTRICITY
International classification
H04L25/03
ELECTRICITY
Abstract
A device, method, and non-transitory computer readable medium to perform a method of blind equalization implemented by a filter bank multi-carrier with offset quadrature amplitude modulation (FBMC-OQAM) transmission system. A first matrix W is obtained by dividing a second matrix V by a receiver waveform matrix G. The second matrix V is obtained by calculating a total objective function J until the total objective function J is either constant or below a threshold error margin. The calculation of the total objective function J includes iterating calculations with a plurality of combination of a frequency bin p and a time slot q. Weights of the obtained first matrix W are applied to an equalizer. A received signal in the equalizer is processed using the applied weights. The weights of the obtained first matrix W are configured to minimize a total outage probability P.sub.out, TOTAL.
Claims
1. A method of blind equalization implemented by processing circuitry of a filter bank multi-carrier with offset quadrature amplitude modulation (FBMC-OQAM) transmission system that includes a transmitter side and a receiver side disposed on opposite sides of a wireless communication channel, comprising: obtaining a first matrix W, which is a matrix of weights to be applied to an equalizer provided in the receiver side of the FBMC-OQAM transmission system, by dividing a second matrix V by a receiver waveform matrix G, wherein the second matrix V is obtained by calculating a total objective function J until the total objective function J is either constant or below a threshold error margin, the calculation of the total objective function J including iterating calculations with a plurality of combinations of a frequency bin p and a time slot q; applying the weights of the obtained first matrix W to the equalizer; and processing a received signal in the equalizer using the applied weights, the received signal being a signal that is received at the receiver side of the FBMC-OQAM transmission system based on a signal that is transmitted over the wireless communication channel from the transmitter side, wherein the weights of the obtained first matrix W are configured to minimize a total outage probability P.sub.out,TOTAL with regard to received signals received after the weights are applied to the equalizer, which the total outage probability P.sub.out,TOTAL is a probabality that data will not be received over the wireless communication channel for a time-frequency batch of symbols.
2. The method of claim 1, wherein calculating the total objective function includes a process comprising: (i) initiating the first matrix W with a plurality of random values; (ii) calculating the second matrix V, the second matrix V being the first matrix W multiplying by the receiver waveform matrix G; (iii) calculating a third matrix A by a second function and a fourth matrix B by a third function, the calculation including the frequency bin p and the time slot q; (iv) performing an eigenvalue decomposition for the fourth matrix B by a fourth function; (v) calculating a fifth matrix C by a fifth function; (vi) performing another eigenvalue decomposition for the fifth matrix C by a sixth function; (vii) calculating a maximum eigenvector u.sub.c,max for the fifth matrix C; and (viii) calculating an optimal solution v.sub.pq,opt by a seventh function.
3. The method of claim 2, wherein the process includes step ix, in which steps (iii)-(viii) are repeated for all p and q.
4. The method of claim 3, wherein the total outage probability P.sub.out,TOTAL is calculated using a first function
5. The method of claim 4, wherein steps i-ix are repeated until the total objective function J is either constant or below a threshold error margin.
6. The method of claim 5, wherein the threshold error margin is |J(i+1)−J(i)|≤ε.
7. The method of claim 2, wherein the second function is A(I.sub.N.Math.g.sub.pq.sup.T)R.sub.H(I.sub.N.Math.g.sub.pq.sup.T).sup.H, wherein I.sub.N is a N-dimensional identity matrix, R.sub.h is a correlational matrix, H is a matrix, and g.sub.pq.sup.T is a pulse vector.
8. The method of claim 2, wherein the third function is B (I.sub.N.Math.g.sub.α.sup.T) R.sub.h(I.sub.N.Math.g.sub.α.sup.T).sup.H+(I.sub.N.Math.g.sub.β.sup.T)R.sub.h(I.sub.N.Math.g.sub.62 .sup.T).sup.H, wherein I.sub.N is a N-dimensional identity matrix, R.sub.h is a correlational matrix, g.sub.α.sup.T and g.sub.β.sup.T are pulse vectors, H is a matrix.
9. The method of claim 2, wherein the fourth function is B=UΛU.sup.H, wherein U is a matrix having eigenvectors and Λ is a diagonal matrix containing eigenvalues of the fourth matrix B.
10. The method of claim 2, wherein the fifth function is C=Λ.sup.−1/2U.sup.HAUΛ.sup.−1/2.
11. The method of claim 2, wherein the sixth function is C=U.sub.cΛ.sub.C U.sub.C.sup.H.
12. The method of claim 2, wherein the seventh function is
13. A filter bank multi-carrier with offset quadrature amplitude modulation (FBMC-OQAM) transmission system that includes a transmitter side and a receiver side disposed on opposite sides of a wireless communication channel and that performs blind equalization, comprising: processing circuitry configured to obtain a first matrix W, which is a matrix of weights to be applied to an equalizer provided in the receiver side of the FBMC-OQAM transmission system, by dividing a second matrix V by a receiver waveform matrix G, wherein the second matrix V is obtained by calculating a total objective function J until the total objective function J is either constant or below a threshold error margin, the calculation of the total objective function J including iterating calculations with a plurality of combinations of a frequency bin p and a time slot q; apply the weights of the obtained first matrix W to the equalizer; and process a received signal in the equalizer using the applied weights, the received signal being a signal that is received at the receiver side of the FBMC-OQAM transmission system based on a signal that is transmitted over the wireless communication channel from the transmitter side, wherein the weights of the obtained first matrix W are configured to minimize a total outage probability P.sub.out,TOTAL with regard to received signals received after the weights are applied to the equalizer, which the total outage probability P.sub.out,TOTAL is a probabality that data will not be received over the wireless communication channel for a time-frequency batch of symbols.
14. The FBMC-OQAM transmission system of claim 13, wherein the processing circuitry calculates the total objective function by performing a process that includes: (i) initiating the first matrix W with a plurality of random values; (ii) calculating the second matrix V, the second matrix V being the first matrix W multiplying by the receiver waveform matrix G; (iii) calculating a third matrix A by a second function and a fourth matrix B by a third function, the calculation including the frequency bin p and the time slot q; (iv) performing an eigenvalue decomposition for the fourth matrix B by a fourth function; (v) calculating a fifth matrix C by a fifth function; (vi) performing another eigenvalue decomposition for the fifth matrix C by a sixth function; (vii) calculating a maximum eigenvector u.sub.c,max for the fifth matrix C; and (viii) calculating an optimal solution v.sub.pq,opt by a seventh function.
15. The FBMC-OQAM transmission system of claim 14, wherein the process includes step ix, in which steps (iii)-(viii) are repeated for all p and q.
16. The FBMC-OQAM transmission system of claim 15, wherein the total outage probability P.sub.out,TOTAL is calculated using a first function
17. The FBMC-OQAM transmission system of claim 16, wherein steps i-ix are repeated until the total objective function J is either constant or below a threshold error margin.
18. The FBMC-OQAM transmission system of claim 17, wherein the threshold error margin is |J(i+1)−J(i)|≤ε.
19. A non-transitory computer readable medium that stores a program that when executed by processing circuitry of a filter bank multi-carrier with offset quadrature amplitude modulation (FBMC-OQAM) transmission system, that includes a transmitter side and a receiver side disposed on opposite sides of a wireless communication channel and that performs blind equalization, causes the processing circuitry to perform a method including: obtaining a first matrix W, which is a matrix of weights to be applied to an equalizer provided in the receiver side of the FBMC-OQAM transmission system, by dividing a second matrix V by a receiver waveform matrix G, wherein the second matrix V is obtained by calculating a total objective function J until the total objective function J is either constant or below a threshold error margin, the calculation of the total objective function J including iterating calculations with a plurality of combinations of a frequency bin p and a time slot q; applying the weights of the obtained first matrix W to the equalizer; and processing a received signal in the equalizer using the applied weights, the received signal being a signal that is received at the receiver side of the FBMC-OQAM transmission system based on a signal that is transmitted over the wireless communication channel from the transmitter side, wherein the weights of the obtained first matrix W are configured to minimize a total outage probability P.sub.out,TOTAL with regard to received signals received after the weights are applied to the equalizer, which the total outage probability P.sub.out,TOTAL is a probabality that data will not be received over the wireless communication channel for a time-frequency batch of symbols.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) A more complete appreciation of this disclosure and many of the attendant advantages thereof will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
DETAILED DESCRIPTION
(12) In the drawings, like reference numerals designate identical or corresponding parts throughout the several views. Further, as used herein, the words “a,” “an” and the like generally carry a meaning of “one or more,” unless stated otherwise.
(13) Furthermore, the terms “approximately,” “approximate,” “about,” and similar terms generally refer to ranges that include the identified value within a margin of 20%, 10%, or preferably 5%, and any values therebetween.
(14) Aspects of this disclosure are directed to a system, device, and method for blind equalization implemented by a filter bank multi-carrier with offset quadrature amplitude modulation (FBMC-OQAM) transmission system (also generally referred to as FBMC-OQAM system). The blind equalization method for the FBMC-OQAM system is based on statistical information of communication channel, i.e., that is, correlation matrix knowledge of the communication channel. According to the aspects of the present disclosure, statistical SINR of the communication channel is maximized to obtain equalizer weights. The blind equalization method is based on minimizing outage probability of the FBMC-OQAM system using maximization of the statistical SINR.
(15)
(16) The FBMC-OQAM system 100 has most of the characteristics of symbol density as that of Orthogonal Frequency Division Multiplexing (OFDM) without cyclic prefix. The complex orthogonality condition gl1,k1(t), gl2,k2(t)
=δ(l2−l1),(k2−k1) of OFDM is substituted by the less strict real orthogonality condition R{
gl1,k1(t), gl2,k2(t)
}=δ(l2−l1),(k2−k1) in FBMC-OQAM system 100.
(17) A transmitted basis pulse for l.sup.th frequency and k.sup.th time symbol denoted by g.sub.l,(t) can be defined as:
(18)
(19) Assuming g.sub.l, represents sampled version of the basis pulse g.sub.l,k(t), then g.sub.l,k is used to denote N samples basis pulse vector for l.sub.th frequency and k.sub.th time symbol. Considering that there are total L frequency sub-carriers and K time symbols. By stacking all the basis pulse vectors in a large transmit matrix G ∈C.sup.N×LK as:
G=[g.sub.1,1. . . g.sub.L,1 g.sub.1,2. . . g.sub.L,K] (2)
(20) and all data symbols in a large transmit symbol vector x ∈C.sup.LK×1 as:
x=[x.sub.1,1. . . x.sub.1,k x.sub.L,1. . . . . . x.sub.L,K].sup.T (3)
(21) the sampled transmit signal s ∈C.sup.N×1 is expressed as
s=Gx (4)
(22) Multipath propagation over time-variant channels is modeled by a time-variant impulse response denoted as h[m.sub.τ, n], where m.sub.τrepresents the delay and n the time position. By writing the impulse response in a time-variant convolution matrix H ∈C.sup.N×N, defined as,
[H].sub.i,j=h[i−j,i] (5)
(23) the received signal can be expressed as:
r=HGx+ñ (6)
(24) where ñ, is a zero mean complex white Gaussian noise vector with correlation matrix P.sub.nI where P.sub.n is noise power. Because of linearity, matrix G can be found even though underlying modulation format is not known in detail. Similarly, sampled receive basis pulses q.sub.l,k∈C.sup.N×1 can also be stacked in a matrix as:
Q=[q.sub.1,1. . . q.sub.1,K q.sub.L,1. . . q.sub.L,K] (7)
(25) For Filter Bank Multi-Carrier (FBMC) system, the receiver uses matched filter, i.e., Q=G. Thus, the received symbols in (6) after pulse de-shaping by Q=G can be written as:
y=G.sup.Hr=G.sup.HHGx+n (8)
(26) where n˜CN (0, P.sub.nGHG).
(27) According to aspects of the present disclosure, the method of blind equalization implemented by the FBMC-OQAM system 100 processes the received symbols y to obtain an estimate of the transmitted symbols.
(28) In conventional techniques of equalization in the FBMC system, the channel induced interference is assumed to be abandoned compared to the noise. Thus, the off-diagonal elements of G.sup.HHG are so minor that they are dominated by noise. Therefore, only the diagonal elements of G.sup.HHG remain and allow to examine the channel according to:
y≈diag{h.sub.1}G.sup.HGx+n (8′)
(29) where h.sub.1∈C.sup.LK×1 describes the one-tap channel. The operator diag{.Math.} creates a diagonal matrix out of a vector. Thus, the conventional one-tap equalizer utilizes knowledge of the communication channel 112 to find the estimate of the transmitted symbols as follows:
(30)
(31) where Re{} represents the real part and Qc(.) is the hard decision operator. Phase equalization of equation (8) in conventional methods is followed by taking the real part, which cancels the imaginary interference in the FBMC-OQAM system 100. Removing the imaginary interference does not remove any important information in an additive white Gaussian noise (AWGN) channel. However, this is not true in fading channels. Thus, the imaginary interference due to loss of the complex orthogonality in the FBMC-OQAM system 100 has a severe effect in fading channels, specifically in selective channels. Also, impact of the off-diagonal elements of G.sup.HHG cannot be ignored as these contribute to interference from other frequency-time symbols. Thus, equalizer design for the FBMC-OQAM system 100 requires dealing with the off-diagonal elements also. Hence, performance of the conventional one-tap channel equalizer and the minimum mean-squared error (MMSE) equalizer is degraded and fails to provide satisfactory results.
(32) According to aspects of the present disclosure, an expression for outage probability for a specific frequency bin (p) and time slot (q) is derived in the presence of an arbitrary equalizer weight matrix denoted by W The equalizer matrix W is multiplied by receiver waveform matrix G The combined receiver matrix is denoted by V defined as:
V=GW (10)
As a result, the received processed signal is given by:
y=V.sup.HHGx+V.sup.Hn (11)
(33) Here, the LK X 1 vector y contains the processed symbols for all time slots and frequency bins. Thus, the pq.sup.th processed symbol takes the following shape:
(34)
(35) Further, the energies of all the terms appearing in equation (12) is evaluated. Considering the energies of the transmitted signal to be unity, i.e., E|x.sub.lk|.sup.2=1, for all l and k, instantaneous SINR for the pq.sup.th processed symbol is given by:
(36)
(37) where Tr(.)represents the trace operator and the terms g.sub.α and g.sub.β represent:
(38)
(39) Further, using the whitening transformation for the channel matrix via H=R.sub.R.sub.
(40)
(41) the SINR expression in equation (13) is reformulated as:
(42)
(43) where
A.sub.1=(ν.sub.pq.sub.
A.sub.2=(ν.sub.pq.sub.
(44) The outage probability of the pq.sup.th processed symbol can now be evaluated as:
P.sub.out.sub.
(45) which is equivalent to write as:
(46)
(47) The above outage probability when evaluated using indefinite quadratic forms approach results in the following expression:
(48)
(49) where λ.sub.t is the t.sup.th eigenvalue of the matrix A.sub.1−γ.sub.thA.sub.2. Now, the total outage probability for a time-frequency batch of LK symbols is computed as follows:
(50)
Design Methodology
(51) According to aspects of the present disclosure, the total outage probability for the whole batch of time-frequency symbols is minimized to design equalizer weights for the FBMC-OQAM system 100.
(52) The expression for the total outage probability derived in equation (23) is a non-convex function of the beamforming weights. Thus, direct minimization of the total outage probability is difficult. Aspects of the present disclosure propose an alternate way to minimize the total outage probability. The outage probability is a function of SINR as shown in equation (20). Thus, the problem of minimizing the total outage probability can be overcome by maximization of SINR. Aspects of the present disclosure propose a technique to maximize the statistical SINR for all the time-frequency symbols of the FBMC-OQAM system 100 in batch processing. To maximize the statistical SINR, expression for the statistical SINR for the time-frequency symbols in the FBMC-OQAM system 100 is derived.
(53) The statistical SINR for the pq.sup.th symbol (denoted as γ.sub.pq) is found by taking expectation on equation (12) and is given by:
(54)
(55) Further, using the linear algebra property: Hg.sub.pq=(I.sub.N.Math.g.sub.pq.sup.τ)h, where hνec(H) is the vectorized version of matrix H, the statistical SINR defined in equation (24) can be expressed as:
(56)
(57) where R.sub.hE(hh.sup.H) is the correlation matrix of channel vector h and I.sub.N is the N-dimensional Identity matrix. Further, by defining the following matrices:
(58)
(59) the statistical SINR can be rewritten as:
(60)
(61) Next, eigenvalue decomposition of the matrix B is applied via B=UΛU.sup.H, where U is the matrix having eigenvectors and Λ is the diagonal matrix containing the eigenvalues of the matrix B. By applying the following variable transformation:
u.sub.pq=Λ.sup.1/2U.sup.Hν.sub.pq (27)
(62) the statistical SINR expressed by equation (26) can be reformulated as:
(63)
(64) where C=Λ.sup.−1/2U.sup.HAUΛ.sup.−1/2.
(65) According to aspects of the present disclosure, the optimum equalizer weight matrix W or equivalently the matrix V whose pq.sup.th column vector is ν.sub.pq are designed. The optimization task is to design the transformed vector u.sub.pq that maximizes the statistical SINR γ.sub.pq. Thus, the following optimization task is defined for the equalizer design in the FBMC-OQAM system 100:
(66)
(67) The optimum value of u.sub.pq is the maximum eigenvector of the matrix C. If C=U.sub.cΛ.sub.CU.sub.C.sup.H is the eigenvalue decomposition of matrix C, the optimum solution of equation (29) is given by:
(68)
(69) Thus, the optimum solution for ν.sub.pq is given by:
(70)
(71) The above method is repeated for all the time-frequency (pq) symbols. The whole process is repeated for certain number of iterations till the change in objective function J is almost zero (or to a very small value, ϵ). The optimum weights of the equalizer design for the FBMC-OQAM system 100 are found through the optimum value of the combined receiver matrix V.
(72) According to aspects of the present disclosure, the equalizer design performs the blind method of equalization for the FBMC-OQAM system 100. The equalizer weights are designed by minimizing the total outage probability which is a function of the SINR. Accordingly, the method of blind equalization implemented by the FBMC-OQAM system 100 is based on maximizing the statistical SINR using the maximum eigenvector technique.
(73) According to aspects of the present disclosure, steps for performing the method of blind equalization implemented by the FBMC-OQAM system 100 includes:
(74) (i) Initializing the equalizer weight matrix W with random values.
(75) (ii) Calculating the combined receiver matrix V using the relation in equation (10).
(76) (iii) Calculating the matrices A and B using the expressions in equation (25) for frequency bin p=1 and time-slot q=1.
(77) (iv) Performing the eigenvalue decomposition of the matrix B via B=UΛU.sup.H and then computing the matrix C using the relation C=Λ.sup.−1/2U.sup.HAUΛ.sup.−1/2.
(78) (v) Performing the eigenvalue decomposition for the matrix C using C=U.sub.cΛ.sub.CU.sub.C.sup.H.
(79) (vi) Computing the maximum eigenvector for the matrix C denoted by u.sub.c,ax.
(80) (vii) Applying the relation in equation (30) to obtain the solution for ν.sub.pq,pt.
(81) (viii) Repeating the steps (iii)-(vii) for all the time-frequency symbols, i.e., frequency bins, p=1, 2 . . . , L and time-slots, q=1, 2, . . . K.
(82) (ix) Repeating the complete process for few iterations till the total objective function J given by
(83)
(84) is almost constant or below a threshold (or acceptable) error margin, c, i.e.,
|J(i+1)−J(i)|≤ε (32)
(85) and storing the optimum solution for the combined receiver matrix V.
(86) (x) Obtaining the equalizer weight matrix W from the solution of the matrix V using equation (10).
(87) (xi) After the optimum solution for the equalizer weight matrix W is obtained, computing the total outage probability P.sub.out,OTAL using the expression in equation (23).
(88) Preliminary Results
(89) Actual implementation of the FBMC-OQAM system 100 follows the steps represented by blocks shown in
(90) The zero mean complex circular Gaussian channel vector h is generated using Jakes model for 2.5 GHz carrier frequency and 500 km/h vehicular speed which corresponds to a maximum Doppler shift of 1.16 kHz. In addition, the channel elements are generated with an exponential model for correlation matrix R.sub.h in which the elements of the correlation matrix are considered a function of correlation coefficient ρ.sub.c;k such that the (i, j).sup.th entry of the correlation matrix is given by:
R.sub.h{i,j}=ρ.sub.h.sup.|i−j|with 0<ρ.sub.h<1
(91) where the parameter ρ.sub.h shows the correlation dependence. For example, ρ.sub.h=0 shows no correlation (or uncorrelated scenario) and ρ.sub.h close to 1 shows that the channel is highly correlated. In the simulation experiments, the value of ρ.sub.h is set to 0.5. The channel matrix H is then constructed from the channel vector h using the relation hνe(H).
(92)
(93) In a first experiment, the total outage probability, P.sub.out,TOTAL of all the time-frequency symbols for L=24 and K=30 obtained via the Monte Carlo simulation is compared with the total outage probability, P.sub.out,TOTAL computed using the analytical expression derived in equation (23). The results of the comparison are shown via the graph in
(94)
(95) In a second experiment, the total outage probability, P.sub.out,OTAL of the FBMC-OQAM system 100 for different equalizers at SNR equal to 20 dB is compared. For comparison, the CMA equalizer and the one-tap equalizer is chosen in addition to blind equalizer of the present disclosure that performs the method of blind equalization according to aspects of the present disclosure. The CMA equalizer is implemented for every sub-carrier using corresponding time sequence. The one-tap equalizer is implemented using the estimation method given in equation (9). The blind equalizer is implemented using the steps (i-xi) for performing the method of blind equalization implemented by the FBMC-OQAM system 100. Initialization of the equalizer weight matrix W is chosen randomly. The total outage probability, P.sub.out,OTAL for all the compared equalizers is computed using the relation in equation (23). The results of the comparison are shown via the graph in
(96)
(97) In a third experiment, the equalization methods of all the equalizers compared in the second experiment are investigated for a fixed threshold value of γ.sup.th=0.0369 and by varying the SNR in the range [1, 45] dB. The results of the comparison are shown via the graph in
(98) According to aspects of the present disclosure, the method of blind equalization implemented by the FBMC-OQAM system 100 is bandwidth efficient as it does not require to send pilot signals for the estimation of instantaneous channel state information (CSI). Instead, the method of blind equalization of the present disclosure relies only on statistical CSI, which can be estimated using received samples based computation.
(99) According to aspects of the present disclosure, the method of blind equalization implemented by the FBMC-OQAM system 100 provides unbiased quality of service (QoS) to all users by imposing appropriate constraints in the optimization task.
(100) According to aspects of the present disclosure, the outage probability is characterized in closed form by employing the approach of characterizing indefinite quadratic forms without imposing approximations.
(101) According to aspects of the present disclosure, the method of blind equalization implemented by the FBMC-OQAM system 100 improves outage performance by maximizing the statistical SINR for each time-frequency symbol.
(102) According to aspects of the present disclosure, the method of blind equalization implemented by the FBMC-OQAM system 100 improves the coverage of a network. Indirect minimization of the total outage probability without sending pilot signals results in better spectral efficiency, which in turn allows accommodating more users.
(103) The method of blind equalization of the present disclosure is based on maximizing the statistical SINR using the maximum eigenvector technique. In contrast to other supervised learning methods of equalization, the proposed method of blind equalization does not require pilot symbols transmission and hence reduces the overhead of transmission. This eventually increases the spectral efficiency, which allows to accommodate more users in the system, thus increasing the network coverage of the system. Moreover, unlike the conventional equalization methods such as Constant Modulus Algorithm (CMA), the proposed method of blind equalization utilizes the channel statistics, i.e., the channel's correlation matrix knowledge. Thus, the proposed method of blind equalization provides better performance than the conventional equalization method.
(104)
(105) At step 502, the method includes obtaining the equalizer weight matrix W by dividing the combined receiver matrix V by the receiver waveform matrix G The combined receiver matrix V is obtained by calculating the total objective function J until the total objective function J is either constant or below the threshold error margin, ε. The calculation of the total objective function J includes iterating calculations with a plurality of combinations of the frequency bin p and the time-slot q.
(106) At step 504, the method includes applying weights of the equalizer weight matrix W to the blind equalizer to perform the method of blind equalization.
(107) At step 506, the method includes processing the received signal in the blind equalizer using the applied weights. The weights of the equalizer weight matrix W are configured to minimize the total outage probability P.sub.out,TOTAL.
(108)
(109) At step 602, the method includes initiating the equalizer weight matrix W with a plurality of random values.
(110) At step 604, the method includes calculating the combined receiver matrix V. The combined receiver matrix V is the equalizer weight matrix W multiplying by the receiver waveform matrix G.
(111) At step 606, the method includes calculating the matrix A and the matrix B by using the expressions in equation (25). The calculation including the frequency bin p and the time slot q.
(112) At step 608, the method includes performing eigenvalue decomposition for the matrix B using the relation B=UΛU.sup.H.
(113) At step 610, the method includes calculating the matrix C using the relation C=Λ.sup.−1/2U.sup.HAUΛ.sup.−1/2.
(114) At step 612, the method includes performing eigenvalue decomposition for the matrix C using C=U.sub.cΛ.sub.CU.sub.C.sup.H.
(115) At step 614, the method includes calculating the maximum eigenvector u.sub.c,ax for the matrix C.
(116) At step 616, the method includes calculating an optimal solution for ν.sub.pq,pt using the relation in equation (30).
(117) The steps 602-616 are repeated for all the time-frequency symbols, i.e., frequency bins, p=1, 2, . . . , L and time-slots, q=1, 2, . . . K.
(118) The first embodiment is illustrated with respect to
(119) The processing circuitry of the FBMC-OQAM system 100 calculates the total objective function by performing a process that includes: (i) initiating the first matrix W with a plurality of random values; (ii) calculating the second matrix V, the second matrix V being the first matrix W multiplying by the receiver waveform matrix G; (iii) calculating a third matrix A by a second function and a fourth matrix B by a third function, the calculation including the frequency bin p and the time slot q; (iv) performing an eigenvalue decomposition for the fourth matrix B by a fourth function; (v) calculating a fifth matrix C by a fifth function; (vi) performing another eigenvalue decomposition for the fifth matrix C by a sixth function; (vii) calculating a maximum eigenvector u.sub.c,max for the fifth matrix C; and (viii) calculating an optimal solution ν.sub.pq,opt by a seventh function.
(120) The process includes step ix, in which steps (iii)-(viii) are repeated for all p and q.
(121) The total outage probability P.sub.out, TOTAL is calculated using a first function
(122)
the p including integers from 1 to L, and the Q including integers from 1 to K.
(123) Steps i-ix are repeated until the total objective function J is either constant or below a threshold error margin.
(124) The threshold error margin is |J(i+1)−J(i)|≤ε.
(125) The second embodiment is illustrated with respect to
(126) Calculating the total objective function includes a process comprising: (i) initiating the first matrix W with a plurality of random values; (ii) calculating the second matrix V, the second matrix V being the first matrix W multiplying by the receiver waveform matrix G; (iii) calculating a third matrix A by a second function and a fourth matrix B by a third function, the calculation including the frequency bin p and the time slot q; (iv) performing an eigenvalue decomposition for the fourth matrix B by a fourth function; (v) calculating a fifth matrix C by a fifth function; (vi) performing another eigenvalue decomposition for the fifth matrix C by a sixth function; (vii) calculating a maximum eigenvector u.sub.c,max for the fifth matrix C; and (viii) calculating an optimal solution v.sub.pq,opt by a seventh function.
(127) The process includes step ix, in which steps (iii)-(viii) are repeated for all p and q. The total outage probability P.sub.out,TOTAL is calculated using a first function
(128)
the p including integers from l to L, and the Q including integers from l to K.
(129) Steps i-ix are repeated until the total objective function J is either constant or below a threshold error margin.
(130) The threshold error margin is |J(i+1)−J(i)|≤ε.
(131) The total outage probability P.sub.out,TOTAL is for a time-frequency batch of LK symbols.
(132) The second function is A (I.sub.N.Math.g.sub.pq.sup.T)R.sub.h(I.sub.N.Math.g.sub.pq.sup.T).sup.H, wherein I.sub.N is a N-dimensional identity matrix, R.sub.h is a correlational matrix, H is a matrix, and g.sub.pq.sup.T is a pulse vector.
(133) The third function is (I.sub.N.Math.g.sub.α.sup.T)R.sub.h(I.sub.N.Math.g.sub.α.sup.T).sup.H+(I.sub.N.Math.g.sub.β.sup.T)R.sub.h(I.sub.N.Math.g.sub.β.sup.T).sup.H+σ.sub.n.sup.2I.sub.N, wherein I.sub.N is a N-dimensional identity matrix, R.sub.h is a correlational matrix, g.sub.α.sup.T, and g.sub.β.sup.T are pulse vectors, and H is a matrix.
(134) The fourth function is B=UΛU.sup.H, wherein U is a matrix having eigenvectors and Λ is a diagonal matrix containing eigenvalues of the fourth matrix B.
(135) The fifth function is C=Λ.sup.−1/2AUΛ.sup.−1/2.
(136) The sixth function is C=U.sub.cΛ.sub.CU.sub.C.sup.H.
(137) The seventh function is
(138)
(139) The third embodiment is illustrated with respect to
(140)
(141) Further, the claims are not limited by the form of the computer-readable media on which the instructions of the inventive process are stored. For example, the instructions may be stored on CDs, DVDs, in FLASH memory, RAM, ROM, PROM, EPROM, EEPROM, hard disk or any other information processing device with which the computing device communicates, such as a server or computer.
(142) Further, the claims may be provided as a utility application, background daemon, or component of an operating system, or combination thereof, executing in conjunction with CPU 701, 703 and an operating system such as Microsoft Windows 7, Microsoft Windows 10, UNIX, Solaris, LINUX, Apple MAC-OS, and other systems known to those skilled in the art.
(143) The hardware elements in order to achieve the computing device may be realized by various circuitry elements, known to those skilled in the art. For example, CPU 701 or CPU 703 may be a Xenon or Core processor from Intel of America or an Opteron processor from AMD of America, or may be other processor types that would be recognized by one of ordinary skill in the art. Alternatively, the CPU 701, 703 may be implemented on an FPGA, ASIC, PLD or using discrete logic circuits, as one of ordinary skill in the art would recognize. Further, CPU 701, 703 may be implemented as multiple processors cooperatively working in parallel to perform the instructions of the inventive processes described above.
(144) The computing device in
(145) The computing device further includes a display controller 708, such as a NVIDIA GeForce GTX or Quadro graphics adaptor from NVIDIA Corporation of America for interfacing with display 710, such as a Hewlett Packard HPL2445w LCD monitor. A general purpose I/O interface 712 interfaces with a keyboard and/or mouse 714 as well as a touch screen panel 716 on or separate from display 710. General purpose I/O interface also connects to a variety of peripherals 718 including printers and scanners, such as an OfficeJet or DeskJet from Hewlett Packard.
(146) A sound controller 720 is also provided in the computing device such as Sound Blaster X-Fi Titanium from Creative, to interface with speakers/microphone 722 thereby providing sounds and/or music.
(147) The general purpose storage controller 724 connects the storage medium disk 704 with communication bus 726, which may be an ISA, EISA, VESA, PCI, or similar, for interconnecting all of the components of the computing device. A description of the general features and functionality of the display 710, keyboard and/or mouse 714, as well as the display controller 708, storage controller 724, network controller 706, sound controller 720, and general purpose I/O interface 712 is omitted herein for brevity as these features are known.
(148) The exemplary circuit elements described in the context of the present disclosure may be replaced with other elements and structured differently than the examples provided herein. Moreover, circuitry configured to perform features described herein may be implemented in multiple circuit units (e.g., chips), or the features may be combined in circuitry on a single chipset, as shown on
(149)
(150) In
(151) For example,
(152) Referring again to
(153) The PCI devices may include, for example, Ethernet adapters, add-in cards, and PC cards for notebook computers. The Hard disk drive 860 and CD-ROM 866 can use, for example, an integrated drive electronics (IDE) or serial advanced technology attachment (SATA) interface. In one implementation the I/O bus can include a super I/O (SIO) device.
(154) Further, the hard disk drive (HDD) 860 and optical drive 866 can also be coupled to the SB/ICH 820 through a system bus. In one implementation, a keyboard 870, a mouse 872, a parallel port 878, and a serial port 876 can be connected to the system bus through the I/O bus. Other peripherals and devices that can be connected to the SB/ICH 820 using a mass storage controller such as SATA or PATA , an Ethernet port, an ISA bus, a LPC bridge, SMBus, a DMA controller, and an Audio Codec.
(155) Moreover, the present disclosure is not limited to the specific circuit elements described herein, nor is the present disclosure limited to the specific sizing and classification of these elements. For example, the skilled artisan will appreciate that the circuitry described herein may be adapted based on changes on battery sizing and chemistry, or based on the requirements of the intended back-up load to be powered.
(156) The functions and features described herein may also be executed by various distributed components of a system. For example, one or more processors may execute these system functions, wherein the processors are distributed across multiple components communicating in a network. The distributed components may include one or more client and server machines, which may share processing, as shown by
(157) The above-described hardware description is a non-limiting example of corresponding structure for performing the functionality described herein.
(158) Obviously, numerous modifications and variations of the present disclosure are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims, the invention may be practiced otherwise than as specifically described herein.