Method and system for measuring interference signal in carrier signal

Abstract

Disclosed in the present disclosure are a method and system for measuring interference signals in carrier signals. The method comprises: receiving multiple valid samples of the carrier signals; based on the multiple valid samples of the carrier signals, gaining value of constant part and mean square error of zero mean part in the interference signals through iterative calculation; calculating the total power value of the interference signals. Since considering the constant part in the interference signals, the measurement method is suitable in the case that the interference signals in the carrier signals is non-zero mean value, and the obtained measurement result of the interference signals will be more accurate than the result through the existing method.

Claims

1. A method for measuring an interference signal in a carrier signal, wherein: receiving, by a receiver, multiple valid samples of the carrier signal, the carrier signal including a transmitted signal and an interference signal, the interference signal having a constant component and a zero-mean component; based on the multiple valid samples of the carrier signal, obtaining, by the receiver, a quantity value of the constant component and a mean square error of the zero-mean component of the interference signal through an iterative calculation; obtaining, by the receiver, a total power value of the interference signal from the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal, so that the quality of carrier channel can be assessed; and performing an iterative calculation using I.sub.0 and .sub.0.sup.2 to obtain a quantity value I.sub.K+1 of the constant component and a mean square error .sub.k+1.sup.2 of the zero-mean component of the interference signal, wherein I.sub.k+1 is the (K+1)th iteration value of the constant component of the interference signal, .sub.k+1.sup.2 is the (K+1)th iteration value of the mean square error of the zero-mean component of the interference signal; and, wherein in the iterative calculation: I.sub.0 is an iteration initial value of the constant component of the interference signal and is calculated according to the following equation: $I_0=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n,$ wherein N is the number of the valid samples of the carrier signal being received, Xn is a sample value received at a time n; and .sub.0.sup.2 is an iteration initial value of the mean square error of the zero-mean component of the interference signal and is calculated according to the following equation: ${}_0^2=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{\left(x_n-I_0\right)}^2+1,$ wherein N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n.

2. The method according to claim 1, wherein, the transmitted signal is a BPSK-modulated signal, in the iterative calculation, I.sub.k+i and o.sub.k+i.sup.2 are calculated according to the following equations:
Ik+1=I01Nn=1N tan h((xnIk).Math.sigma.Math.k2).Math.sigma.Math.k+12=1+1Nn=1N((xnIk)2tan h((xnIk).Math.sigma.Math.k2)2(xnIk)) Wherein, n and k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n, k is a iteration number; I.sub.k is the K iteration value of the constant component of the interference signal, o.sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

3. The method according to claim 2, wherein further comprising: comparing, by the receiver, |i.sub.k+iI.sub.k| against a first threshold, comparing, by the receiver, |<7.sub.k+i.sup.2Ok.sup.2| against a second threshold; when |I.sub.k+iIk| is less than the first threshold and |ok+i.sup.2Ok.sup.2| is less than the second threshold, terminating, by the receiver, the iterative calculation, wherein |.sub.k+i and Ok+i.sup.2 are the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal obtained in the K+I iteration, respectively.

4. The method according to claim 3, wherein that the total power value of the interference signal is the sum of I.sub.k+1.sup.2 and .sub.k+1.sup.2.

5. The method according to claim 1, wherein, the transmitted signal is a QAM-modulated signal, in the iterative calculation, I.sub.k+i and o.sub.k+i.sup.2 are calculated according to the following equations:
Ik+1=I01Nn=1Nm=1M(2m1)Sm,n(Ik,.Math.sigma.Math.k2)m=1 MCm,n(Ik,.Math.sigma.Math.k2).Math.sigma.Math.k+12=.Math.sigma.Math.021Nn=1N(m=1M(4m(m1)Cm,n(Ik,.Math.sigma.Math.k2)2(2m1)(xnIk)Sm,n(Ik,.Math.sigma.Math.k2))m=1MCm,n(Ik,.Math.sigma.Math.k2)) Wherein, M represents a constellation order of the real and imaginary parts in QAM, Cm, n (I, .Math.sigma.Math.2)=exp(2m(m1).Math.sigma.Math.2) cos h((2m1)(xnI).Math.sigma.Math.2)Sm,n(I,.Math.sigma.Math.2)=exp(2m(m1).Math.sigma.Math.2)sin h((2m1)(xnI).Math.sigma.Math.2) Wherein, n, m, k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n; I.sub.k is the K iteration value of the constant component of the interference signal, o.sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

6. The method according to claim 5, wherein further comprising: comparing, by the receiver, |I.sub.k+1I.sub.k| against a first threshold, comparing, by the receiver, |.sub.k+1.sup.2.sub.k.sup.2| against a second threshold; when |I.sub.k+1I.sub.k| is less than the first threshold and |.sub.k+1.sup.2.sub.k.sup.2| is less than the second threshold, terminating, by the receiver, the iterative calculation, wherein I.sub.k+1 and .sub.k+1.sup.2 are the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal obtained in the K+1 iteration, respectively.

7. The method according to claim 6, wherein that the total power value of the interference signal is the sum of I.sub.k+1.sup.2 and .sub.k+1.sup.2.

8. A system for measuring interference signal in a carrier signal, characterized in comprising at least one processor configured to: receive multiple valid samples of the carrier signal, the carrier signal including a transmitted signal and an interference signal, the interference signal having a constant component and a zero-mean component; based on the multiple valid samples of the carrier signal, obtain a quantity value of the constant component and a mean square error of the zero-mean component of the interference signal through an iterative calculation; obtain a total power value of the interference signal from the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal, so that the quality of carrier channel can be assessed; and performing an iterative calculation using I.sub.0 and .sub.0.sup.2 to obtain a quantity value I.sub.K+1 of the constant component and a mean square error .sub.k+1.sup.2 of the zero-mean component of the interference signal, wherein I.sub.k+1 is the (K+1)th iteration value of the constant component of the interference signal, .sub.k+1.sup.2 is the (K+1)th iteration value of the mean square error of the zero-mean component of the interference signal; and, wherein in the iterative calculation: I.sub.0 is an iteration initial value of the constant component of the interference signal and is calculated according to the following equation: $I_0=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n,$ wherein N is the number of the valid samples of the carrier signal being received, Xn is a sample value received at a time n; and .sub.0.sup.2 is an iteration initial value of the mean square error of the zero-mean component of the interference signal and is calculated according to the following equation: ${}_0^2=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{\left(x_n-I_0\right)}^2+1,$ wherein N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n.

9. The system according to claim 8, wherein, for a transmitted signal which is a BPSK-modulated signal, the at least one processor is further configured to calculate I.sub.k+i and c.sub.k+i.sup.2 according to the following equations respectively:
Ik+1=I01Nn=1N tan h((xnIk).Math.sigma.Math.k2).Math.sigma.Math.k+12=1+1 Nn=1N((xnIk)2tan h((xnIk).Math.sigma.Math.k2)2(xnIk)) wherein, n and k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n; I.sub.k is the K iteration value of the constant component of the interference signal, .Math.sigma.Math..sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal; for a transmitted signal which is a QAM-modulated signal, the iterative calculation unit calculates I.sub.k+1 and .Math.sigma.Math..sub.k+1.sup.2 according to the following equations respectively:
Ik+1=I01Nn=1Nm=1M(2m1)Sm,n(Ik,.Math.sigma.Math.k2)m=1MCm,n(Ik,.Math.sigma.Math.k2).Math.sigma.Math.k+12=.Math.sigma.Math.021Nn=1N(m=1M(4m(m1)Cm,n(Ik,.Math.sigma.Math.k2)2(2m1)(xnIk)Sm,n(Ik,.Math.sigma.Math.k2))m=1MCm,n(Ik,.Math.sigma.Math.k2)) wherein, M represents a constellation order of the real and imaginary parts in QAM, Cm,n(I, .Math.sigma.Math.2)=exp(2m (m1).Math.sigma.Math.2) cos h((2m1)(xnI).Math.sigma.Math.2) Sm,n(I, .Math.sigma.Math.2)=exp (2m (m1).Math.sigma.Math.2) sin h((2m1)(xnI).Math.sigma.Math.2) wherein, n, m, k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n; I.sub.k is the K iteration value of the constant component of the interference signal, or.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

10. The system according to claim 9, wherein that the at least one processor is further configured to: compare |I.sub.k+1I.sub.k| against a first threshold and compare |.sub.k+1.sup.2.sub.k.sup.2| against a second threshold; when |I.sub.k+1I.sub.k| is less than the first threshold and |.sub.k+1.sup.2.sub.k.sup.2| is less than the second threshold, terminate the iterative calculation, wherein I.sub.k+1 and .sub.k+1.sup.2 are the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal obtained in the K+1 iteration, respectively; the at least one processor is further configured to calculate the sum of I.sub.k+1.sup.2 and .sub.k+1.sup.2 as the total power value of the interference signal.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the disclosure and, together with the description, serve to explain the principles of the disclosure.

(2) With reference to the drawings, this disclosure may be understood more clearly according to the following description, wherein:

(3) FIG. 1 shows a flowchart of an embodiment of a method that is provided in this disclosure;

(4) FIG. 2A shows estimation errors of the method provided in this disclosure and a contrast method;

(5) FIG. 2B shows an iteration convergence test of the method provided in this disclosure;

(6) FIG. 3 shows a structural diagram of an embodiment of a system provided in this disclosure.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

(7) Various exemplary embodiments of the present disclosure will now be described in detail with reference to the drawings. It should be noted that the relative arrangement of the components and steps, the numerical expressions, and numerical values set forth in these embodiments do not limit the scope of the present disclosure unless it is specifically stated otherwise.

(8) Meanwhile, it should be appreciated that, for the convenience of description, various parts shown in those drawings are not necessarily drawn on scale.

(9) The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the disclosure, its application, or uses.

(10) Techniques, methods and apparatus as known by one of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate.

(11) In all of the examples illustrated and discussed herein, any specific values should be interpreted to be illustrative only and non-limiting. Thus, other examples of the exemplary embodiments could have different values.

(12) Notice that similar reference numerals and letters refer to similar items in the following figures, and thus once an item is defined in one figure, it is possible that it need not be further discussed for following figures.

(13) FIG. 1 shows a flowchart of an embodiment of a method that is provided in this disclosure.

(14) At step S101, multiple valid samples of a carrier signal are received.

(15) A process may be performed for each carrier wave separately. A carrier signal may be expressed with the following signal model:
x(n)=s(n)+u(n)=s(n)+I+v(n)

(16) Wherein, n is a time index, x(n) represent a received carrier signal, s(n) represents a transmitted signal, u(n) represents an interference signal with a non-zero mean value, wherein the interference signal is an interference signal contained in the carrier signal in addition to an effective signal. The received carrier signal may comprise a transmitted signal s(n) and an interference signal u(n). Wherein the interference signal u(n) may comprise a constant component I and a zero-mean component v(n). I is a constant representing the mean of non-zero mean interference, which is also referred to as deterministic interference.

(17) For each carrier wave, channel gain, i.e., the transmitted signal s(n), may be obtained with any channel estimation algorithm that is well known by those skilled in the art.

(18) Multiple effective samples of a carrier signal are received to measure interference contained in the carrier signal. In an embodiment provided in this disclosure, receiving multiple effective samples of a carrier signal is only based on a set of effective symbols that are received without being aware of the content of the transmitted signal s(n). The number of effective samples may be determined based on actual needs, and does not have any limitation thereof.

(19) Further, in practical applications, according to requirements of practical applications, the effectiveness of samples and conditions under which the samples are obtained may be further considered. Effective samples are sample point signals of the carrier signal that are correctly received and identified. For example, it is necessary to keep statistical interference/noise characteristics unchanged in an interference/noise estimation process, that is, the interference must have constant mean and variance. Thereby, it is required that the samples are obtained from the same carrier wave and the same user, spanning over a time window that is not too large, so as to guarantee unchanged channel interference/noise characteristics of the samples. Also, there is a requirement about the transmitted data, for example, the transmitted data must show a random characteristic.

(20) At step S102, based on the multiple effective samples of the received carrier signal, a quantity value of the constant component and a mean square error of the zero-mean component of the interference signal are obtained through an iterative calculation.

(21) In an embodiment of this disclosure, a quantity value I of the constant component of the interference signal and a variance of the zero-mean component v(n) of the interference signal are estimated respectively, wherein the variance of v(n) is .sup.2, after that a total power value of the interference signal is further calculated.

(22) Different estimation methods may be used for transmitted signals in different modulation manners. Below, an illustration will be given in which Binary Phase Shift Keying (BPSK) and Quadrature Amplitude Modulation (QAM) are taken as examples. Those skilled in the art benefiting from the concept of this disclosure may come up with estimation methods for transmitted signals modulated in other manners to further calculate the total power of the interference signal.

(23) For a transmitted signal that is BPSK modulated, the zero-mean component of the interference signal v(n) is a Gaussian noise. When the transmitted signal is a BPSK signal with a constellation {1,+1}, the probability density function of the received signal may be represented as:

(24) $f\left(x\right)=\frac{1}{2}\frac{1}{\sqrt{2{}^2}}\left[\exp \left(-\frac{{\left(x-1-I\right)}^2}{2{}^2}\right)+\exp \left(-\frac{{\left(x+1-I\right)}^2}{2{}^2}\right)\right]$

(25) A receiver receives N samples x(1), x(2), . . . x(N) of the received signal, and then obtain a log-likelihood function L(x)=ln {f(x.sub.1)f(x.sub.2) . . . f(x.sub.N)} represented as follows:

(26) $L\left(x\right)=\frac{N}{2}\ln \left(\frac{1}{2{}^2}\right)+\underset{n=1}{\overset{N}{.Math.}}\left(-\frac{{\left(x_n-I\right)}^2+1}{2{}^2}+\ln \cosh \left(\frac{2\left(x_n-I\right)}{2{}^2}\right)\right)$

(27) According to the Maximum likelihood method (ML), the above log-likelihood function is maximized to obtain corresponding I and .sup.2 as an estimation result. Based on the ML criterion, a more accurate estimation may be obtained.

(28) Taking partial derivatives of L(x) with respect to I and .sup.2 respectively and setting the partial derivatives to zero, the following linear equations containing I and .sup.2 may be obtained:

(29) $I=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\tanh \left(\frac{\left(x_n-I\right)}{{}^2}\right)$${}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\left({\left(x_n-I\right)}^2-\tanh \left(\frac{\left(x_n-I\right)}{{}^2}\right)2\left(x_n-I\right)\right)$

(30) In order to obtain a closed-form estimation of I and .sup.2, an iterative calculation is performed to solve the equations. The iterative calculation mainly comprises two steps: initialization and iterative updating.

(31) At the initialization step, to simplify the calculation, the last items on the right side of the above two equations are ignored, to obtain an iteration initial value I.sub.0 of the constant component of the interference signal and an iteration initial value .sub.0.sup.2 of the mean square error of the zero-mean component of the interference signal. The calculation equations are as follows:

(32) $I_0=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n$${}_0^2=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{\left(x_n-I_0\right)}^2+1$

(33) Wherein, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n.

(34) Then, a quantity value I.sub.K+1 of the constant component and a mean square error .sub.k+1.sup.2 of the zero-mean component of the interference signal are obtained through the iterative calculation using I.sub.0 and .sub.0.sup.2, wherein I.sub.k+1 is the K+1 iteration value of the constant component of the interference signal, .sub.k+1.sup.2 is the K+1 iteration value of the mean square error of the zero-mean component of the interference signal.

(35) Particularly, an iteration update result is calculated using two equations as follows:

(36) $I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\tanh \left(\frac{\left(x_n-I_k\right)}{{}_k^2}\right)$${}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\left({\left(x_n-I_k\right)}^2-\tanh \left(\frac{\left(x_n-I_k\right)}{{}_k^2}\right)2\left(x_n-I_k\right)\right)$

(37) Wherein, n and k are natural numbers, N is the number of the valid samples of the received carrier signal, n is the serial number of sample point, Xn is a sample value received at a time n, k is a iteration number; I.sub.k is the K iteration value of the constant component of the interference signal, .sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

(38) For a transmitted signal that is QAM modulated, the real part and the imaginary part of the received signal may be processed separately. Both parts may be processed in the same manner. The real part and the imaginary part of the QAM signal have the same constellation {2M+1, . . . 1,+1 . . . 2M1}, wherein M represents the order of the constellation of the QAM real part and the imaginary part. Both the probability density functions of real part and the imaginary part of the received signal may be represented as:

(39) $f\left(x\right)=\frac{1}{2M}\frac{1}{\sqrt{2{}^2}}\left[\underset{m=1}{\overset{M}{.Math.}}\exp \left(-\frac{{\left(x-\left(2m-1\right)-I\right)}^2}{2{}^2}\right)+\exp \left(-\frac{{\left(x+\left(2m-1\right)-I\right)}^2}{2{}^2}\right)\right]$

(40) The obtained log-likelihood function is represented as follows:

(41) $L\left(x\right)=-N\ln M+\frac{N}{2}\ln \left(\frac{1}{2{}^2}\right)+\underset{n=1}{\overset{N}{.Math.}}\left(-\frac{{\left(x_n-I\right)}^2+1}{2{}^2}+\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I,{}^2\right)\right)$

(42) Wherein

(43) $C_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\cosh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$
is an estimation result. n, m are natural numbers.

(44) Similarly, taking partial derivatives of L(x) with respect to I and .sup.2 respectively and setting the partial derivatives to zero, the following linear equations containing I and .sup.2 may be obtained:

(45) $I=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\frac{\underset{m=1}{\overset{M}{.Math.}}\left(2m-1\right)S_{m,n}\left(I,{}^2\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I,{}^2\right)}$${}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}({\left(x_n-I\right)}^2-\frac{\underset{m=1}{\overset{M}{.Math.}}\left(4m\left(m-1\right)C_{m,n}\left(I,{}^2\right)-2\left(2m-1\right)\left(x_n-I\right)S_{m,n}\left(I,{}^2\right)\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I,{}^2\right)})$
wherein,

(46) 0$S_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\sinh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$

(47) Similarly, with an iteration calculation, the iterative calculation mainly comprises two steps: initialization and iterative updating. At the initialization step, to simplify the calculation, the last items on the right side of the above two equations are ignored, to obtain an iteration initial value I.sub.0 of the constant component of the interference signal and an iteration initial value .sub.0.sup.2 of the mean square error of the zero-mean component of the interference signal similar to BPSK modulation. The calculation equations are as follows:

(48) $I_0=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n$${}_0^2=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{\left(x_n-I_0\right)}^2+1$

(49) An iterative calculation is performed using I.sub.0 and .sub.0.sup.2, in which I.sub.0 and .sub.0.sup.2 are calculated according to the following equations respectively:

(50) $I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\frac{\underset{m=1}{\overset{M}{.Math.}}\left(2m-1\right)S_{m,n}\left(I_k,{}_k^2\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I_k,{}_k^2\right)}$${}_{k+1}^2={}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\left(\frac{\underset{m=1}{\overset{M}{.Math.}}\left(4m\left(m-1\right)C_{m,n}\left(I_k,{}_k^2\right)-2\left(2m-1\right)\left(x_n-I_k\right)S_{m,n}\left(I_k,{}_k^2\right)\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I_k,{}_k^2\right)}\right)$

(51) Wherein, M is the constellation order of the QAM real part and the imagery part, n, m, k are natural numbers, and m is an intermediate variable.

(52) $C_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\cosh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$$S_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\sinh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$

(53) N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n. I.sub.k is the K iteration value of the constant component of the interference signal, .sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

(54) Notice that the transmitted signal may be modulated in different manners. The modulation methods include but not limited to BPSK and QAM. Thus, for a transmitted signal that is modulated in a different manner, expressions that are used in the iterative calculation of the K iteration value I.sub.k of the constant component and the K iteration value .sub.k.sup.2 of the mean square error of the zero-mean component may be different.

(55) For the completion of the iterative calculation, an iteration number may be specified. When the iteration number is reached, the iterative calculation terminates. Or a determination may be made based on an iteration result. The iterative calculation terminates if the iteration result satisfies a certain condition. As to the embodiment in which an iteration number is specified, it may be implemented according to methods that are well known by those skilled in the art. Taking the last method, i.e., making a determination based on an iteration result, as an example, a detail description will be given in this embodiment.

(56) In the iteration calculation, a first determination and a second determination may be made on an iteration result. Particularly, a first threshold and a second threshold may be specified to compare |I.sub.k+1I.sub.k| against the first threshold and compare |.sub.k+1.sup.2.sub.k.sup.2| against the second threshold.

(57) When |I.sub.k+1I.sub.k| is less than the first threshold and |.sub.k+1.sup.2.sub.k.sup.2| is less than the second threshold, it may be considered that a condition of convergence has been met and the iterative calculation may terminate. I.sub.k+1 and .sub.k+1.sup.2 are the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal obtained through the iteration calculation respectively.

(58) The thresholds may be reasonably set in the iteration calculation according to computing conditions and error requirements. For example, the thresholds may have values ranging from 10.sup.3 to 10.sup.5. For example, in this embodiment, the threshold used herein is 10.sup.3.

(59) In order to clearly describe the technical solution of this disclosure and its effect, a detail description of the method of this disclosure for measuring interference signal in a carrier signal and its effect will be given below with reference to FIG. 2A to FIG. 2B. FIG. 2A shows estimation errors of a method provided in this disclosure and a contrast method, FIG. 2B shows an iteration convergence test of the method provided in this disclosure.

(60) In this example, a measurement result of computer simulation is given. Simulation parameters are as follows: transmitted symbols are binary random signals, deterministic interference I=0.5, noise variance=0.2512, corresponding SNR=6 dB, simulation number=200. Signals that are transmitted each time and noises are randomly generated, and then an average error of parameter estimations is obtained.

(61) Estimation error shown by the vertical axis of FIG. 2A is measured as mean square error, i.e., estimation error=E{(estimation resultactual parameter).sup.2}. Wherein, E{ } represents the average of multiple repetitive tests.

(62) In the contrast method shown in FIG. 2A, only random noise error is considered. Calculation is performed in an application based on statistics. Particularly, the constant component of the interference signal is not considered, instead of an iterative calculation, the initial value of the mean square error of the zero-mean component is taken as its estimation result.

(63) In FIG. 2A, line A represents estimation errors obtained with an existing method, line B represents estimation errors of the estimation method provided in this disclosure. It can be seen from FIG. 2A, estimation errors obtained with the contrast method is far larger than the estimation errors of measurements performed according to this disclosure.

(64) In FIG. 2B, line C represents estimated quantities of the constant component of the interference signal, and line D represents estimated mean square errors of the zero-mean component of the interference signal. It can be seen from FIG. 2B, the iterative calculation adopted in this disclosure converges upon 10 iterations.

(65) At step S103, a total power value of the interference signal is obtained from the quantity of the constant component of the interference signal and the mean-square error of the zero-mean component of the interference signal.

(66) Particularly, the total power value of the interference signal may be the sum of I.sub.k+1.sup.2 and .sub.k+1.sup.2.

(67) With the carrier interference measurement method of this disclosure, the calculated total power value of the interference signal is the sum of the quantity of the constant component of the interference signal and the mean-square error of the zero-mean component of the interference signal. This measurement method is directed to a situation in which interference signal in the carrier signal has a non-zero mean value to take the constant component of the interference signal into account. Thus, a more accurate interference signal measurement result may be obtained.

(68) A system for measuring interference signal contained in a carrier signal is further disclosed in this disclosure. FIG. 3 shows a structural diagram of an embodiment of a system provided in this disclosure.

(69) The system may comprise a carrier signal receiving unit 11, an iterative calculation unit 12 and a total power calculation unit 13.

(70) The carrier signal receiving unit 11 receives multiple valid samples of the carrier signal. The carrier signal may comprise a transmitted signal and an interference signal.

(71) The interference signal may comprise a constant component and a zero-mean component.

(72) The iterative calculation unit 12 may obtain a quantity value of the constant component and a mean square error of the zero-mean component of the interference signal through an iterative calculation based on the multiple valid samples received by the carrier signal receiving unit 11.

(73) The iterative calculation unit 12 may calculate an iteration initial value I.sub.0 of the constant component of the interference signal according to the following equation:

(74) $I_0=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{x}_n$

(75) Wherein, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n.

(76) The iterative calculation unit 12 may calculate an initial iteration value .sub.0.sup.2 of the mean square error of the zero-mean component of the interference signal according to the following equation:

(77) ${}_0^2=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{\left(x_n-I_0\right)}^2+1$

(78) Wherein, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n.

(79) The iterative calculation unit 12 may perform an iterative calculation using I.sub.0 and .sub.0.sup.2 to obtain a quantity value I.sub.K+1 of the constant component and a mean square error .sub.k+1.sup.2 of the zero-mean component of the interference signal, wherein I.sub.k+1 is the K+1 iteration value of the constant component of the interference signal, .sub.k+1.sup.2 is the K+1 iteration value of the mean square error of the zero-mean component of the interference signal.

(80) For transmitted signals that are modulated in different manners, the K iteration values of the constant components and the K+1 iteration values .sub.k.sup.2 of the mean square errors of the zero-mean components of the interference signals may be different.

(81) Preferably, for a transmitted signal which is a BPSK-modulated signal, in the iterative calculation:

(82) $I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\tanh \left(\frac{\left(x_n-I_k\right)}{{}_k^2}\right)$${}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\left({\left(x_n-I_k\right)}^2-\tanh \left(\frac{\left(x_n-I_k\right)}{{}_k^2}\right)2\left(x_n-I_k\right)\right)$

(83) Wherein, n and k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n; I.sub.k is the K iteration value of the constant component of the interference signal, .sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

(84) For a transmitted signal which is a QAM-modulated signal, in the iterative calculation:

(85) $I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\frac{\underset{m=1}{\overset{M}{.Math.}}\left(2m-1\right)S_{m,n}\left(I_k,{}_k^2\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I_k,{}_k^2\right)}$${}_{k+1}^2={}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}\left(\frac{\underset{m=1}{\overset{M}{.Math.}}\left(4m\left(m-1\right)C_{m,n}\left(I_k,{}_k^2\right)-2\left(2m-1\right)\left(x_n-I_k\right)S_{m,n}\left(I_k,{}_k^2\right)\right)}{\underset{m=1}{\overset{M}{.Math.}}{C}_{m,n}\left(I_k,{}_k^2\right)}\right)$

(86) Wherein, M represents a constellation order of the real and imaginary parts in QAM,

(87) $C_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\cosh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$$S_{m,n}\left(I,{}^2\right)=\exp \left(-\frac{2m\left(m-1\right)}{{}^2}\right)\sinh \left(\frac{\left(2m-1\right)\left(x_n-I\right)}{{}^2}\right)$

(88) Wherein, n, m, k are natural numbers, N is the number of the valid samples of the received carrier signal, Xn is a sample value received at a time n; I.sub.k is the K iteration value of the constant component of the interference signal, .sub.k.sup.2 is the K iteration value of the mean square error of the zero-mean component of the interference signal.

(89) Notice that the modulation methods of the transmitted signal are not limited to the above two manners. Correspondingly, in the iteration calculation, the expression of the quantity value of the constant component and the expression of the mean-square error of the zero-mean component may be adjusted accordingly.

(90) To terminate the iterative calculation, an iteration number may be specified so that the iterative calculation terminates when the iteration number is reached, or a determination may be made on the result of the iterative calculation and the iterative calculation terminates if the result satisfies a certain condition.

(91) For instant, in this embodiment, the iterative calculation unit 12 may compare |I.sub.k+1I.sub.k| against a first threshold and compare |.sub.k+1.sup.2.sub.k.sup.2| against a second threshold. When |I.sub.k+1I.sub.k| is less than the first threshold and |.sub.k+1.sup.2.sub.k.sup.2| is less than the second threshold, the iterative calculation terminates, wherein I.sub.k+1 and .sub.k+1.sup.2 are the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal obtained in the K+1 iteration, respectively.

(92) The total power calculation unit 13 may obtain a total power value of the interference signal from the quantity value of the constant component and the mean square error of the zero-mean component of the interference signal.

(93) Particularly, the total power calculation unit 13 may calculate the sum of I.sub.K+1.sup.2 and .sub.k+1.sup.2 as the total power value of the interference signal.

(94) Various embodiments of this disclosure have been described in a progressive manner, wherein each embodiment has emphases different to other embodiments. A mutual reference can be made to various embodiments for the same or similar portions thereof. With respect to the system embodiment, due to the correspondence to the method embodiments thereof, only a simple description is given, refer to those method embodiments for related portions.

(95) The method and system of this disclosure may be realized in many manners. For example, the method and system of this disclosure may be realized by software, hardware, firmware or any combination thereof. The above sequence of steps in the method is merely for illustration, and thus unless specified otherwise, the steps in the method of this disclosure are not limited to the sequence particularly described above. Further, in some embodiments, this disclosure can be implemented as programs recorded in a record medium, including machine readable instructions for implementing the method of this disclosure. Therefore, those record mediums storing programs for implementing the method of this disclosure are also covered in this disclosure.

(96) Although some particular embodiments of this disclosure have been described in detail with examples, those skilled in the art may understand that the above examples are merely for illustration but not limitation on the scope of this disclosure. Those skilled in the art may understand various modifications can be made to the above embodiments without departing from the scope of this disclosure. The scope of this disclosure is defined by attached claims.