Determining distance to source of passive intermodulation product (PIM) in a distributed base station

Abstract

A method for determining the distance to source of a passive intermodulation (PIM) product for a distributed base station is provided. The method includes generating a test signal; transmitting the test signal on a front haul communication link between a baseband unit and a remote radio head of the distributed base station; receiving the PIM product as a reflected signal at the remote radio head; multiplying the reflected signal by a conjugate of the test signal to create a sinewave; determining a frequency of the sinewave; and determining the distance to the source of the PIM based on the determined frequency of the sinewave.

Claims

1. A method for determining the distance to a source of a passive intermodulation (PIM) product for a distributed base station, the method comprising: generating a test signal wherein the test signal comprises a chirp defined by: $\mathrm{chirp}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{t}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}$ with an initial frequency (f.sub.t) that is set to a negative value and a frequency change rate (m) is positive, A is a constant, and n is a discrete time index; transmitting the test signal on a front haul communication link between a baseband unit and a remote radio head of the distributed base station; receiving the PIM product as a reflected signal at the remote radio head; multiplying the reflected signal by a conjugate of the test signal to create a sinewave; determining a frequency of the sinewave; and determining the distance to the source of the PIM based on the determined frequency of the sinewave.

2. The method of claim 1, wherein generating the test signal comprises generating a chirp signal and a continuous wave (CW) signal.

3. The method of claim 1, wherein generating the test signal includes generating a PN sequence, a Zadoff-Chu sequence or any sequence that shows good correlation properties.

4. The method of claim 1, wherein transmitting the test signal comprises transmitting the test signal at a high enough power to generated an intermodulation (IM) product in a component of the remote radio head in order to estimate the RRH round trip delay and the delay of the front haul link.

5. The method of claim 1, wherein determining the distance to the source comprises determining round trip travel time between the remote radio head and the PIM source from the frequency of the reflected signal.

6. A method for determining the distance to a source of a passive intermodulation (PIM) product for a distributed base station, the method comprising: generating a test signal wherein the test signal comprises a chirp defined by: $\mathrm{chirp}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{t}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}$ with an initial frequency (f.sub.t) that is set to a negative value and a frequency change rate (m) is positive, A is a constant, and n is a discrete time index; transmitting the test signal on a front haul communication link between a baseband unit and a remote radio head of the distributed base station; receiving the PIM product as a reflected signal at the remote radio head; analyzing an Eigen component of the reflected signal; and determining the distance to the source of the PIM based on the analysis of the Eigen component.

7. The method of claim 6, wherein analyzing an Eigen component comprises: multiplying the reflected signal by a conjugate of the test signal to create a sinewave; and determining a frequency of the sinewave, wherein the distance to the source of the PIM is based on the frequency of the sinewave.

8. A tester for determining the distance from a distributed base station to a source of passive intermodulation (PIM), the tester comprising: a signal generator configured to be coupled to a front haul communication link between a baseband unit and a remote radio head of the distributed base station, wherein the signal generator is configured to generate a test signal, wherein the test signal comprises a chirp defined by: $\mathrm{chirp}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{t}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}$ with an initial frequency (f.sub.t) that is set to a negative value and a frequency change rate (m) is positive, A is a constant, and n is a discrete time index, and to transmit the test signal on the front haul communication link; a digital signal processor, configured to be coupled to the front haul communication link, wherein the digital signal processor is operable to cause the tester to do the following: receiving the PIM product as a reflected signal at the remote radio head; multiplying the reflected signal by a conjugate of the test signal to create a sinewave; determining a frequency of the sinewave; and determining the distance to the source of the PIM based on the determined frequency of the sinewave.

9. The tester of claim 8, wherein the digital signal processor includes a chirp product module that is configured to multiply the reflected signal by a conjugate of a reference chirp signal.

10. The tester of claim 9, wherein the digital signal processor includes a cross correlation module, coupled to the chirp product module, that is configured to apply cross correlation to the sinewave from the chirp product module.

11. The tester of claim 10, wherein the digital signal processor includes a frequency estimation circuit, coupled to the cross correlation module, that is configured to determine the frequency of the sinewave.

12. The tester of claim 11, wherein the frequency estimation circuit includes: an Eigen decomposition module, configured to separate the output of the cross correlation module into a noise subspace and a signal subspace; a vector extractor, coupled to the Eigen decomposition module, configured to extract a first vector that corresponds to the signal subspace and at least one vector corresponding to the noise subspace; and a least square estimator, coupled to the vector extractor, configured to extract the frequency of the first vector.

13. The tester of claim 11, wherein the digital signal processor further comprises a summation module that is adapted to sum at least two sample bunches of signals received at the digital signal processor.

Description

DRAWINGS

(1) FIG. 1 is a block diagram of one exemplary embodiment of a tester for a distributed base station system within which the techniques for determining the distance to a PIM source described here can be used.

(2) FIG. 2 is block diagram of an embodiment of a tester coupled to a remote radio head that illustrates a general approach to determining a distance to a PIM source.

(3) FIG. 3 is a timing diagram illustrating the relative timing of signals transmitted in the tester of FIG. 2.

(4) FIG. 4 is a block diagram of one embodiment of a DSP chain in a tester of FIG. 2.

(5) FIG. 5 is an alternate embodiment of a portion of the DSP chain in the tester of FIG. 2.

(6) FIG. 6 is another alternate embodiment of a portion of the DSP chain in the tester of FIG. 2.

(7) FIG. 7 is another alternate embodiment of a portion of the DSP chain in the tester of FIG. 2.

(8) FIG. 8 is another alternate embodiment of a portion of the DSP chain in the tester of FIG. 2.

(9) FIG. 9 is another alternate embodiment of a portion of the DSP chain in the tester of FIG. 2.

(10) FIGS. 10 and 11 are illustrations of chirp signals.

DETAILED DESCRIPTION

(11) Embodiments of the present invention provide a technique to measure a distance from a remote radio head (RRH) to a source of a passive intermodulation (PIM) product. The distance to PIM consists in measuring the distance between the antenna connector and the source of the PIM product. Advantageously, embodiments of the present invention enable measurement of the distance to the PIM source over the front-haul communication link, e.g., a CPRI/OBSAI/ORI interface, between the RRH and a baseband unit (BBU). As a result, RF maintenance and troubleshooting does not require a technician to climb up a cell tower to connect the test equipment to the antenna port.

(12) Embodiments of the present invention measure distance to a PIM source by transmitting a test signal or predefined pattern (e.g., a chirp) at the transmit antenna port and analyzing a reflecting signal, e.g., the PIM product. The interest of the chirp is its good correlation property, product property and low sensitivity to self-interference in the radio channel between reflections of the test signal (chirp).

(13) The test equipment used to determine the distance to the PIM source is placed between the RRH and the BBU and is connected to the physical optical interface ports of the RRH and the BBU. The test equipment transmits the test signal on the front-haul communication link (e.g., CPRI, OBSAI, ORI) to be transmitted by the RRH. The RRH receives the reflected signal (PIM product) in the uplink band and passes the received signal to the BBU over the front-haul communication link. The capture of the reflected signal is done at the CPRI/OBSAI interface carrying I/Q samples received from the RRH module and transmitted to the BBU module. The test equipment operates during site maintenance. It can transfer any type of I/Q samples to the antenna and receive any type of I/Q samples captured from the antenna. The test equipment can detect external or internal PIM sources. External PIM sources are issued from corroded objects, such as fences, metal roofs or tower rusty bolts in sight from the antenna. Internal PIM sources are issued from imperfect RF feed line and connectors (e.g., metal chips floating inside connectors).

(14) This disclosure provides multiple solutions to the problem of determining the distance to the PIM source using the transmitted test signal that is launched from the test equipment inserted between the BBU and the RRH. In one embodiment, an Eigen component of the reflected test signal is analyzed to determine the distance to the PIM source. In an alternative embodiment, a cross correlation of the reflected signal with its corresponding transmitted test signal (reference pattern) is analyzed to produce the distance measurement. In this embodiment, the test signal or reference pattern is a chirp signal, a PN sequence, a Zadoff-Chu sequence or any other appropriate sequence that shows good correlation properties. In another embodiment, the distance to the PIM source is determined using an enhanced version of reflectometry for PIM source default.

(15) Advantageously, the solutions described herein do not require knowledge of the proprietary operation and maintenance messages at the digital interface to the front-haul communication link to obtain various delay of transceiver loop (RRH, fiber, antenna cable). Further, the embodiments described herein enable determination of the distance to a single PIM source with a good accuracy. Embodiments of the present invention also extract information on the distance to more than one PIM source if there is sufficient separation between the PIM sources.

(16) Distributed Base Station

(17) FIG. 1 is a block diagram of one exemplary embodiment of a tester 100 for a distributed base station system, indicated generally at 102. In the exemplary embodiment shown in FIG. 1, the system 102 comprises a plurality of baseband units (BBU) 104-1 to 104-N and a plurality of remote radio heads (RRH) 106-1 to 106-N that communicate over a plurality of wireless radio frequency (RF) channels with one or more wireless units 108 (such as mobile telephones, smartphones, tablets, wireless modems for laptops or other computers or for other devices such as wireless sensors or other Internet of Things (IOT) or machine-to-machine (M2M) devices) using one or more standard wireless air interfaces. The exemplary embodiment of system 102 shown in FIG. 1 may support several air interfaces, e.g., three air interfaces including, but not limited to, Long-Term Evolution (LTE) 4G air interface described in the Third Generation Partnership Project (3GPP) Technical Specification (TS) 36.211 Evolved Universal Terrestrial Radio Access (E-UTRA); Physical Channels and Modulation specification produced by the 3GPP, Advanced Wireless Services (AWS-1), Personal Communications Services (PCS), CLR, GSM, WiMax, and others. It is to be understood that other air interfaces can be used.

(18) Each BBU 104 is communicatively coupled to the core network 110 of a wireless service provider using a suitable bi-directional backhaul communication link 111 and interface (for example, using a wireless or wired ETHERNET connection and using the LTE S1 interface). The backhaul communication link 111 can also be used for base station-to-base station communications using the LTE X2 interface.

(19) Each BBU 104 is communicatively coupled to a corresponding RRH 106 using a bi-directional front-haul communication link 112. In the exemplary embodiment shown in FIG. 1, the bi-directional front-haul communication link 112 is implemented using a plurality of pairs of optical fibers, where, in each pair, one optical fiber is used for downlink communications from the BBU 104 to the RRH 106 and the other optical fiber is used for uplink communications from the RRH 106 to the BBU 104. Further, as shown in FIG. 1, the bi-directional communication link 112 are split into two parts to allow a tester 100 (described in more detail below) to be inserted between the BBUs 104 and the RRHs 106. Namely, the bi-directional communication link 112 include a first part 112a connecting BBU 104 to tester 100 and a second part 112b connecting tester 100 to a respective RRH 106. It is to be understood that the front-haul communication link 112 can be implemented in other ways. The exemplary embodiment shown in FIG. 1 is described here as using a CPRI interface for communications between each BBU 104 and the corresponding RRH 106 over the front-haul communication link 112. It is to be understood, however, that a different front-haul interface could be used (for example, the OBSAI or ORI interface).

(20) As noted above, each BBU 104 performs baseband processing for the particular air interface that is being used to wirelessly communicate over its assigned RF channel, and the RRH 106 performs radio frequency processing to convert baseband data output from the BBU 104 to radio frequency signals for radiating from one or more antennas 114 that are connected to the RRH 106 at antenna port 113 via coaxial cable 115 and to produce baseband data for the associated BBU 104 from radio frequency signals that are received at the RRH 106 via one or more antennas 114. Advantageously, tester 100 enables determination of the distance between antenna port 113 of RRH 106 and PIM source 103. Although PIM source 103 is illustrated as being external to system 102, as noted above, PIM source 103 may be internal to the system 102, e.g., in antenna connector 113.

(21) During normal operation of the system 102, in the downlink direction, the BBUs 104 generate downlink baseband IQ data to encode frames of downlink user and control information received from the core network for communication to the wireless units 108 over the appropriate wireless RF channels. The downlink baseband IQ data is communicated from the BBUs 104 to the RRHs 106 over the respective front-haul communication link 112. The RRHs 106 receive the downlink baseband IQ data and generate one or more downlink analog radio frequency signals that are radiated from the one or more antennas 114 for reception by the wireless units 108. The wireless units 108 perform baseband processing, in accordance with the air interface, on the received downlink analog RF downlink signals in order to recover the frames of downlink user and control information.

(22) During normal operation of the system 102, in the uplink direction, the wireless units 108 generate, in accordance with the air interface, uplink analog radio frequency signals that encode uplink user and control information that is to be communicated to the core network 110 and transmits the generated uplink analog RF signals over the wireless RF channel. The uplink analog RF signals are received by one or more antennas 114 connected to the RRHs 106. The RRH 106 that receives the uplink analog RF signal produces uplink baseband IQ data from the received uplink analog RF signals. The uplink baseband IQ data is communicated from the RRH 106 to the associated BBU 104 over the front-haul communication link 112. The BBU 104 receives the uplink baseband IQ data and performs baseband processing, in accordance with the air interface, on the uplink baseband IQ data in order to recover the uplink user and control information transmitted from the wireless units 108. The BBU 104 communicates the recovered uplink user and control information to the core network 110 over the backhaul communication link 111 using the backhaul interface.

(23) The RRHs 106 are typically installed remotely from its corresponding BBU 104, near the antennas 114 and is mounted to a structure 116 (such as a tower, pole, building, tree, or other structure). For example, the RRH 104 can be mounted near the top of the structure 116 and the BBU 104 can be located on the ground, where the optical fibers used to implement the front-haul communication link 112 run up the structure 116 to couple the BBU 104 to the RRU 106. Although FIG. 1 shows the RRH 106 mounted near the top of structure 116, the RRH 106 can be mounted at other positions relative to the structure 116, for example, approximately midway between the bottom and top of the structure 116.

(24) During normal operation, signals transmitted from RRH 106 may interact with PIM source 103 and produce PIM products. PIM products 103 are received by RRH 106 and can interfere with the normal operation of system 102. As such, during maintenance operations, tester 100 is used to determine the distance from the RRH 106 to help identify the PIM sources so that corrective action can be taken.

(25) Tester

(26) As shown in FIG. 1, tester 100 can be coupled into the font-haul communication link 112 between the BBUs 104 and the RRHs 106. In this manner, tester 100 is enabled to insert test signals in the downlink frames of baseband IQ data transmitted to the RRH. Also, tester 100 is enabled to receive signals in the uplink frames of baseband IQ data communicated from the RRHs 106. This capturing of baseband data enables tester 100 to determine the distance from an antenna port 113 to a PIM source 103 using the techniques described in more detail below.

(27) A user can interact with the software 130 executing on the tester 100 using a user device 136, e.g., smartphone, tablet, or computer. The user device 136 is communicatively coupled to the tester 100. In the exemplary embodiment shown in FIG. 1, the tester 100 includes one or more wired interfaces 138 (for example, an ETHERNET interface and/or a USB interface) and wireless interfaces 140 (for example, a Wi-Fi wireless interface) to communicatively couple the tester 100 to a local area network or directly to the user device 136. Moreover, a remotely located user device 136 can access the tester 100 via a connection established over the local area network and/or a public network such as the Internet. In one embodiment, the software 130 implements a webserver that is operable to present a browser-based user interface that enables a user to use a general-purpose Internet browser installed on the user device 136 to interact with the software 130 on the tester 100.

(28) Also, although the embodiments described above are described as using antenna carriers in downlink CPRI frames, it is to be understood that the techniques described here can be used with other streams of baseband IQ data (for example, streams of baseband IQ data communicated over an OBSAI or ORI interface).

(29) Common Functionality

(30) FIG. 2 is a block diagram that illustrates a general approach to determining the distance to a PIM source 103 from an antenna connector 113 of remote radio head 106 according to the various embodiments of the present invention. Under this approach, signal generator 201 in tester 100 generates a test signal in appropriate IQ containers of the front-haul communication link (CPRI/OBSAI/ORI interface) 112. In one embodiment, the test signal includes a chirp signal (discussed in more detail below) and a continuous wave (CW) signal. The chirp signal and the continuous wave signal are transmitted as baseband data over the front-haul communication link 112 to RRH 106.

(31) RRH 106 includes downlink components 214 and uplink components 216 that are each coupled to diplexer 205. Downlink components 214 include the series connection of digital to analog converter (DAC) 202, mixer 203 and power amplifier (PA) 204. Power amplifier 204 is coupled to diplexer 205. Antenna 114 is coupled to antenna port 113 at diplexer 205. Diplexer 205 is also coupled to the uplink components 216. The uplink components 216 include the series connection of low noise amplifier (LNA) 206, mixer 207, and analog to digital converter (ADC) 208. ADC 208 is coupled to a digital signal processing (DSP) chain 209 in tester 100 via front-haul communication link 112 (e.g., CPRI/OBSAI/ORI connection).

(32) In operation, tester 100 determines the distance to a PIM source 103. To do this, tester 100 launches the test signal toward the RRH 106 on the front-haul communication link 112. These test signals, and their reflections, traverse two different paths before arriving back at the tester 100. The test signals and reflections on these two paths enable determination of the distance to PIM source 103. On a first path 210, the test signal travels to the PIM source 103 and produces a PIM product that is reflected and returns to the tester 100. On a second path 212, the test signal produces an intermodulation (IM) product that is reflected back to the tester 100 within the RRH 106.

(33) The manner in which the distance to the PIM source 103 is determined depends on the sampling rate used in RRH 106. In one embodiment, a low sampling rate is used in the RRH 106. In other embodiments, a high sampling rate is used in the RRH 106.

(34) In the case of the low sampling rate, a three-step process is used to determine the distance to the PIM source 103. First, a test signal is launched to determine a coarse timing estimation of the chirp arrival. The purpose of the coarse timing is described in more detail below. Second, further test signals are launched to generate an estimation of the PIM relative timing. At this second step, the signal level for the test signal is set so that no saturation occurs in RRH 106 (neither in PA 204 nor in DAC 202) such that only the PIM reflection will be reflected back to the RRH 106 and the delay of the PIM reflection can be estimated with a dedicated algorithm. Finally, additional test signals are launched that saturate circuitry in the RRH 106 (e.g., either PA 204 or DAC 202). These test signals are used to detect the relative timing of the IM reflection signal (reference timing). At this stage both IM and PIM reflections arrive at the RRH 106 and optical PIM tester 100. An algorithm is employed in the DSP chain 209 of optical PIM tester 100, like a canceler or DC blocker, to eliminate the PIM reflection during IM reflection extraction. Once the PIM reflection is eliminated, the IM relative timing can be estimated with the dedicated algorithm. Finally the difference between both IM and PIM relative timing is used in the distance to PIM calculation.

(35) In the case of a high sampling rate in the RRH 106, a two-step process is used to determine the distance to the PIM source 103. First, a test signal is launched to determine a coarse timing estimation of the chirp arrival. Secondly, further test signals are launched to generate a fine estimation of the PIM and IM relative timing using dedicated a pseudo spectrum analysis algorithm.

(36) Returning to the low sampling rate embodiment, in the first step, the transmitted chirp sequence is repeated such that the tester can use the multiple sequences for signal to noise ratio improvement. The coarse timing determination is implemented with a classical correlator with an expected transmitted sequence. To have a better chance to find the coarse timing, the PA 204 of RRH 106 is driven to saturation. The accuracy of the coarse timing is expected at worse to be around +/1 sampling timing (e.g. 1/7.56 MHz). Once the coarse timing of the repeated chirp is identified, the second and third steps can proceed.

(37) In the second and third steps, the duration of the chirp sequence (Tchirp) and the repeating interval is known. The coarse timing of the received chirp is known thanks to the first step (coarse timing estimation). As a result, at this stage it is possible to sum at least two Tchirp sample bunches to improve signal to noise ratio (because those chirp sequences are repeated on a regular basis). The number of samples collected within a duration of Tchirp is N=Tchirp*fs where fs is the receiver sampling rate.

(38) Further, in the second and third steps, test signals (a number of chirp signals) are transmitted sequentially on first path 210 and second path 212 to enable comparison of the timing on the respective paths. Advantageously, the difference between the first path 210 and the second path 212 corresponds to the distance between the antenna port 113 of RRH 106 and PIM source 103.

(39) This is shown graphically in FIG. 3. The test signal, e.g., chirp 302, is transmitted from tester 100 at time t=0. At time t1, the IM signal 304 generated on path 212 is received at tester 100. Finally, at time t2, the PIM product 306 produced by the test signal 302 traversing first path 210 is received at tester 100. The time difference of arrival (t2t1) is the time it takes for the test signal 302 to travel to the PIM source 103 from the antenna port 113 and for the reflection (PIM product) 306 to return to the antenna port 113. The distance to the PIM source from the RRH 106 is calculated as follows:

(40) $\begin{array}{cc}D=v\frac{\left(t2-t1\right)}{2}& \left(1\right)\end{array}$

(41) where v wave speed of the media between the antenna port 113 of RRH 106 and PIM source 103. If the media is composed from different materials, the formula is a composite formula taking into account the length of different materials and associated propagation speed.

(42) Returning to FIG. 2, signal generator 201 launches a test signal, e.g., the chirp and CW signals, in the front-haul communication link 112 toward DAC 202. DAC 202 converts the received test signal to an analog signal. The analog signal is mixed with a carrier frequency at mixer 203 and amplified at power amplifier 204. The amplified analog signal is then passed to the antenna 114 by diplexer 205. In the presence of a non-linearity (PIM source 103) a PIM product is reflected back to RRH 106 and received at antenna 114. The reflected signal (PIM product) is passed by the diplexer to the uplink components 216 where it is amplified (LNA 206), down converted (mixer 207) and converted to a digital signal (ADC 208). The digital signal is passed over front-haul communication link 112 to DSP chain 209 of tester 100 (reflection 306 received at time t2 in FIG. 3) for processing to determine the distance to the PIM source 103.

(43) As discussed above, a third test (third step) or reference signal is used in determining the distance to the PIM source 103. This third test signal is also generated by signal generator 201 and launched onto the front-haul communication link 112 (CPRI/OBSAI interface). This third test signal is designed to generate intermodulation within the RRH 106. The third test signal can be a chirp with amplitude that saturates power amplifier 204 or DAC 202. A nonlinear signal can also be convenient to generate IM in the transceiver if it is not smoothed by a digital frontend filter of RRH 106. As a result of intermodulation between tones, the PA 204 generates an IM signal that is also composed of a chirp waveform at frequency falling in the uplink band of the RRH 106. This last signal crosses the diplexer 205 to the LNA 206. The chirp wave can therefore be analyzed by the DSP chain 209 to capture a reference timing (e.g., time t1 in FIG. 3). This IM reference timing is used to compute the PIM source distance as discussed above.

(44) The foregoing description provides the general framework for determining the distance from the RRH 106 to the PIM source 103. As discussed with respect to FIG. 2, DSP chain 209 processes the reflected signals to determine this distance based on the reflected signals. Various exemplary embodiments are described below for this DSP chain 209 that can be used to extract out the distance to the PIM source.

(45) Distance Determination Using Frequency Estimation

(46) FIG. 4 is a block diagram of one embodiment of a DSP chain 209 for determining the distance to a PIM source 103. In this embodiment, the distance is determined using a frequency estimation. At a high level, this embodiment converts the received reflections to a sinewave with a frequency that is proportional to the delay of the reflected signal. The difference between the time delays of these two signals is the same as t2t1 in FIG. 3 above, and can thus be used to determine the distance to the PIM source according to equation (1) above.

(47) As discussed above with respect to FIG. 2, signal generator 201 generates a test signal to be used to determine the distance to the PIM source. This test signal, in one embodiment, includes a chirp signal and a continuous wave (CW) signal. The expression of the transmitted chirp signal is:

(48) $\begin{array}{cc}\mathrm{chirp}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{t}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}& \left(2\right)\end{array}$

(49) The expression of the transmitted CW signal is:
CW(n)=e.sup.j2f.sup.cw.sup.n(3)

(50) It is noted that f.sub.CW, f.sub.r, f.sub.t are normalized frequencies. They are normalized with the frequency sampling fs. In one embodiment, both chirp signal and CW signal are transmitted at the same time to produce intermodulation at the PIM source 103 and the IM signal at the analog nonlinearity (e.g., PA 204).

(51) At the DSP chain 209, the received chirp signal is:

(52) $\begin{array}{cc}s\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{r}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}+w\left(n\right)& \left(4\right)\end{array}$

(53) In equation (4), s(n) represents the received noisy chirp with the following parameters: m: the chirp frequency change rate f.sub.r: the chirp initial frequency n: the discrete time index A: constant complex component representing the chirp amplitude and phase rotation due to PIM/IM reflection. w(n): the received noise

(54) In the DSP chain 209, the reflection (noisy chirp, s(n)) is received at chirp product module 402. The chirp product module 402 converts the chirp signal into a sinewave by multiplying the reflection (noisy chirp) with the conjugate of a reference chirp expression. The chirp product formulation is as follows:

(55) $\begin{array}{cc}s\left(n\right){\mathrm{chirp}}^{*}\left(n\right)=\left({\mathrm{Ae}}^{\left\{j2\left(f_r\left(n-\right)+\frac{1}{2}{m\left(n-\right)}^2\right)\right\}}+w\left(n\right)\right)e^{{}^{\left\{-j2\left(f_{r}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}}& \left(5\right)\end{array}$

(56) The exact position in time of the received chirp is not known. Thus, the product will apply with an unknown time delay . This uncertainty is modeled by including the parameter in equation (5). After some simplification, the formula after chirp product can be written as:

(57) $\begin{array}{cc}s\left(n\right){\mathrm{chirp}}^{*}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(-mn+\frac{1}{2}m{}^2-f_r\right)\right\}}+w^{}\left(n\right)& \left(6\right)\end{array}$

(58) Equation (6) represents a sinusoidal signal with frequency proportional to the time delay uncertainty.

(59) In the context of distance to PIM measurement, the ADC 208 captures two chirps with two different time uncertainties: 1 and 2 (a first chirp from IM reflection and a second from PIM reflection). After chirp product module 402, the signal is as provided in equation (7):

(60) $\begin{array}{cc}x\left(n\right)=s\left(n\right){\mathrm{chirp}}^{*}\left(n\right)=A_{\mathrm{IM}}{e}^{\left\{j2\left(-m{}_{1}n+\frac{1}{2}m{}_1^2-f_r{}_1\right)\right\}}+A_{\mathrm{PIM}}{e}^{\left\{j2\left(-m{}_{2}n+\frac{1}{2}m{}_2^2-f_r{}_2\right)\right\}}+w^{}\left(n\right)& \left(7\right)\end{array}$

(61) This signal is composed of two sinewaves whose frequencies are proportional to the time delay uncertainty. Note that the phase rotation at the PIM/IM source (phase of complex component A.sub.IM and A.sub.PIM) has no influence on the resulting frequency values m2 and m1. The round trip time duration (antenna port 113 to PIM source 103 to antenna port 113) can be expressed as: 21. As a result the estimation of the round trip time duration consists in estimating the frequencies of the noisy composite sinewave signal and then extracting the distance to the PIM source 103 based on the estimated frequencies.

(62) The embodiment of FIG. 4 includes an optional summation component 401 that can be used to improve the signal-to-noise ratio. When the summation component 401 is included, it is possible to sum at least two Tchirp sample bunches to improve signal to noise ratio (because those chirp sequence are repeated on a regular basis). For the sake of processing load, summation component 401 precedes the chirp product module 402 in FIG. 4. The number of samples collected within a duration of Tchirp is N=Tchirp*fs where fs is the receiver sampling rate.

(63) To further improve signal-to-noise ratio a cross correlation applies on the composite sinewave signal by cross correlation module 404. Two vectors of signal are extracted at two different time instants sufficiently spaced such that noise is assumed uncorrelated. The subsequent signal vectors are x1(n)=x(n) and x2(n)=x(n+D). The cross correlation expression is as follows:
r.sub.x.sub.1.sub.x.sub.2(k)=E[x.sub.1*(n)x.sub.2(n+k)]=A.sub.IM.sup.2e.sup.{j2(m.sup.1.sup.(k+D))}+A.sub.PIM.sup.2e.sup.{j2(m.sup.2.sup.(k+D))}+w(k)(8)

(64) Parameter D is a constant parameter. w(k) is the noise component whose value is equal zero thanks to expectation operation. In the implementation, the noise component will not equal zero as the expectation operation will not apply on an infinite sample. Thus, to eliminate the residual noise, Eigen decomposition is applied to the correlation matrix. Note that the cross correlation formula represents the summation of two sinusoids whose frequency are proportional to the respective 1 and 2 delays.

(65) R.sub.x1x2 represents all cross correlations for different time lag k. Column vector and row vector of the correlation matrix contains the noisy sinewaves whose frequencies needs to be determined. R.sub.x1x2 is structured as a Toeplitz matrix:

(66) $\begin{array}{cc}{R}_{x_1{x}_2=}\left[\begin{array}{cccc}{r}_{x_1{x}_2}\left(0\right)& r_{x_1{x}_2}\left(-1\right)& \mathrm{.Math.}& r_{x_1{x}_2}\left(-p+1\right)\\ r_{x_1{x}_2}\left(1\right)& r_{x_1{x}_2}\left(0\right)& \mathrm{.Math.}& r_{x_1{x}_2}\left(-p+2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ r_{x_1{x}_2}\left(p-1\right)& r_{x_1{x}_2}\left(p-2\right)& \mathrm{.Math.}& r_{x_1{x}_2}\left(0\right)\end{array}\right]& \left(9\right)\end{array}$

(67) The Rx1x2 is not necessarily symmetric. The next step consists in extracting precise values of the tones.

(68) The process for determining distance to a PIM source using frequency estimation differs based on the sampling rate used by the distributed base station under test. Those I/Q samples can be different due to the LTE bandwidth used at the air interface. There are 4 sampling rates in LTE (7.68 MHz, 15.36 MHz, 23.04 MHz, 30.72 MHz). In one embodiment, the tester 100 is able to estimate the distance to the PIM source 103 at +/2 meters. As a result of the speed of light, this is corresponding to 1/2.25*Ts of possible induced delay accuracy requirement (Ts is the sampling period with 30.72 MHz sampling frequency). It is difficult to obtain this accuracy with a sampling at I/Q CPRI interface of 7.68 MHz, 15.36 MHz, and 23.04 MHz. As a result two different algorithms are proposed for the different CPRI interface sampling configuration: configuration with low sampling rate (7.68 MHz, 15.36 MHz and 23.04 MHz) and one configuration with high sampling rate (30.72 MHz).

(69) Frequency Estimation with Low Sampling Rate

(70) If the I/Q sampling period is 7.68 MHz, 15.36 MHz or 23.04 MHz (sampling period greater than delay estimation accuracy requirement) and if the PIM and IM source are closely located, the delay estimation process is performed in three steps. These are the same three steps as discussed above with respect to FIG. 2.

(71) Step 1

(72) Step 1 is described above and is used to extract the coarse timing estimate.

(73) Step 2

(74) The block diagram of FIG. 4 includes frequency estimator 406 that is part of the signal processing chain 209 and is dedicated to the extraction of the frequency tones from each intermodulation reflection.

(75) Step 2 consists in examining only one PIM reflection (no PA saturation forced at the RRH 106). First assume only one PIM source or a grouped of PIM sources located close to each other. The instant the chirp product applies on the received signal is arbitrary. Nevertheless, if the product timing applies aligned with the received chirp reflection, a DC component will result making impossible the Eigen decomposition. Therefore, to avoid DC component as a result of the chirp product, a rough chirp timing estimation is performed in step 1.

(76) Eigen decomposition (Singular Value Decomposition can also be used) is applied on the cross correlation matrix at Eigen decomposition module 410 to separate the noisy signal into a noise subspace and a signal subspace. The signal subspace is the sinusoidal or group of close sinusoidal information. And the noise subspace contains information on the noise w(k).
R.sub.x.sub.1.sub.x.sub.2=VV.sup.1(10)

(77) V is Eigen vector matrix. is the diagonal matrix of Eigen values. The first Eigen vector is extracted by vector extractor 412 and corresponds to the signal subspace and the last Eigen vectors are the noise subspace. The Eigen vector does not need to be orthogonal between each other to proceed with frequency estimation. The Eigen value can be positive, negative or complex.

(78) For high signal-to-noise ratio (SNR), the computed first eigenvector is a pure sinewave. For low SNR, the extracted eigenvector is further cleaned from high frequency noise. An optional low pass filter 414 is used in that case. As discussed above, the sinewave frequency is proportional to the delay. A linear least square method is implemented by least square estimator 416 to extract frequency as described below. The signal Eigen vector comprises the complex signal of a sinewave. To obtain the frequency the trigonometric problem is converted into a linear problem. The sinewave trigonometric function is:
v=Acos(wt+)(11)

(79) where v is the signal Eigen vector, t is a time vector, w is the angular frequency and is the initial phase of the resulting sinewave.

(80) The linearization of the sinewave trigonometric function is given in the equation (11). The angular frequency w can be found as follows:

(81) $\begin{array}{cc}\left[\begin{array}{c}w\\ \end{array}\right]={\left[\begin{array}{cc}{t}_0& 1\\ \mathrm{.Math.}& \mathrm{.Math.}\\ t_N& 1\end{array}\right]}^{-1}\mathrm{angle}\left(v\right)& \left(12\right)\end{array}$

(82) Note that the inverse matrix

(83) ${\left[\begin{array}{cc}{t}_0& 1\\ \mathrm{.Math.}& \mathrm{.Math.}\\ t_N& 1\end{array}\right]}^{-1}$
can be predetermined. As a result, no matrix inversion is needed in the implementation in that case.

(84) To improve reliability of the delay some statistic on the estimated frequency may be done. Step 2 can be repeated several times and the corresponding result can be analyzed using the histogram of the measured frequencies. The histogram is centered to the true frequency. A weighted mean is then applied to find the estimation of the true frequency. The PIM relative delay is then obtained.

(85) Step 3

(86) Once the estimation of the PIM relative delay is found, the third step can take place. The third step consists in sending the same pair of signals (CW and Chirp) with power increased to saturate the PA 204. The PIM signal may be stronger than the IM resulting in blinding the IM frequency estimation. To avoid this, the timing chirp product may be adjusted to mitigate the PIM effect over the measured IM. To do so the timing chirp can be set aligned to the instant corresponding to the PIM time arrival instant (as previously computed in step 2). The resulting signal after chirp product will be a DC component that can be removed with a DC blocker (notch filter at 0 Hz) or DC subtractor (the estimation of the DC is done from samples captured in step 2). Once the chirp product timing adjusted and the DC component removed, the IM reflection delay estimation described in FIG. 4 can proceed. Step 3 outcome is the IM chirp timing.

(87) Alternative Step 3

(88) An alternative of the above step 3 algorithm (adjust timing and DC blocker) is described hereafter. The DC blocker can be replaced with a canceler as shown in FIG. 5. After DAC/PA saturation, the first eigenvector is composed of two sinewaves resulting from both PIM and IM reflections. The PIM reflection timing is known perfectly, therefore it is possible to eliminate its contribution from the received data. A canceller 500 can be employed for this purpose.

(89) Canceler 500 includes a signal synthesizer 502 that synthetizes an estimation of PIM correlation samples from Eigen component found during step 2. The synthesized formula is:
{circumflex over (R)}.sub.PIM=V.sub.1.sub.1V.sub.1.sup.1(13)
Where: V.sub.1: is the first eigenvector from step 1. .sub.1: is the first eigenvalue from step 1. {circumflex over (R)}.sub.PIM: is the cross correlation matrix estimation carrying PIM signal subspace information only.

(90) A gain is applied on the synthesized signal decomposition to equalize with the step 3 power increased (power increase to saturate the PA/DAC). The estimated cross-correlation PIM matrix is subtracted from step 3 data capture such that only IM and noise signal remain. Then, the eigenvector decomposition can apply to estimate the frequency and IM reflection delay.

(91) Second Alternative to Step 3

(92) Another alternative of the above step 3 algorithm is described hereafter. In this embodiment, the signal synthesis block of canceler 500 is replaced with a sample mean estimator (SME) 602 as a way to eliminate the 0 mean centered noise signal and to keep only the PIM information. Step 3 runs without timing adjustment matching the PIM delay but step 3 runs with saturation in the RRH (e.g., PA 204 or DAC 202). The cross correlation vector estimated from step 2 is:

(93) 0${\hat{R}}_{\mathrm{PIM}}\left(k\right)=\frac{1}{N}\underset{i}{\overset{N}{.Math.}}{r}_{x1x2,i}\left(k\right)$
Where: r.sub.x1x2,i: is the cross correlation vector from step 2 data capture i. N: is the number of captures performed from step 2 algorithm {circumflex over (R)}.sub.PIM(k): is the estimation of the PIM cross correlation signal at lag k cleaned from noise component.
Alternative Embodiment to Eigen Decomposition for Low Sampling Rate

(94) The Eigen decomposition technique described above is an effective solution but comes with significant processing cost. The processing complexity depends on the matrix structure and is rated as a O(n.sup.3) or O(n.sup.2 log(n)) algorithm. For embedded solution and large matrix decomposition, the Eigen decomposition algorithm may be replaced by another method to separate noise from the user signal data. A simple method already described above is the Sample Mean Estimation algorithm. Compared to Eigen decomposition this method requires a significant number of received signal occurrences such that the SME algorithm mitigates noise by summation.

(95) Step 2 uses the sample mean estimator (SME) plus frequency estimation based Delay estimator block diagram of FIG. 7 to estimate PIM cross correlation vector and PIM delay.

(96) Step 3 aims to eliminate the contribution of PIM signal from PA/DAC saturation data captures (step 3 data capture). The PIM signal is eliminated by subtracting the step 2 PIM cross correlation vector from step 3 cross correlation vector as described in Step 3 block diagram of FIG. 8. Once the PIM signal is eliminated, the resulting data samples comprise the IM signal plus noise. To remove noise, these data are forwarded to SME (noise cleaning algorithm) plus frequency estimator resulting in the IM delay associated frequency. The delay computation block diagram converts the computed frequency into the IM delay timing.

(97) Statistic may apply on the estimated frequency values to improve estimation reliability.

(98) The step 2 and step 3 can be summarized in FIG. 8.

(99) Because the SME algorithm is based on long term integration principle, it is sensitive to ADDA sampling jitter and local oscillator drift.

(100) Frequency Estimation with High Sampling Rate (30.72 MHz)

(101) Under this condition several frequency tones (PIM sources) can be estimated in one step. The frequencies estimation process is detailed in FIG. 9.

(102) In this embodiment, DSP chain 209 includes a frequency estimation circuit 900. Frequency estimation circuit 900 includes a matrix converter 902 that converts the cross correlation matrix to a Hermitian matrix formulation according to the following equation:
R=(R.sub.x.sub.1.sub.x.sub.2+R.sub.x.sub.1.sub.x.sub.2.sup.H)/2(14)

(103) It is assumed that the signal is interfered only with white Gaussian noise. As a result the Eigen vectors of signal sub-space are orthogonal to Eigen vectors of noise sub-space. The different sinewaves resulting from the chirp product are most of the time mixed together in one or more Eigen vector. Therefore, it is not possible to apply the previous method on each Eigen vector. A pseudo spectrum method can be applied here. The intended pseudo spectrum can be the method called MUSIC, or modified pseudo spectrum, Eigen Vector method, Eigen AIC, Eigen MDL, minimum variance spectral estimator.

(104) Some tests show negative Eigen value, as a result it is recommended this vector shall have a constructive contribution into the spectrum estimation than a destructive behavior. As result in the formula of the EigenVector (EV) noise subspace frequency estimator an absolute operator applies on the Eigen values. The enhanced pseudo spectrum EV formula is therefore:

(105) $\begin{array}{cc}{P}_{\mathrm{EEV}}\left(e^{\mathrm{jw}}\right)=\frac{1}{\underset{i=p+1}{\overset{M}{.Math.}}\frac{1}{.Math.{}_i.Math.}{.Math.e^H{v}_i.Math.}^2}& \left(15\right)\end{array}$
P.sub.EEV(e.sup.jw): is the enhance Eigen vector pseudo density spectrum.
M: the matrix correlation size.
p: is the number of vector in the signal sub-space.
v.sub.i: is the Eigen vector i (noise sub-space).
e.sup.H: is the Hermitian reference complex exponential vector of angular frequency w.

(106) The resulting spectrum can be sum-up after each iteration to improve the estimation. As an alternative of the summation and if the number of PIM/IM is estimated, it is possible to build a histogram of the detected frequency.

(107) The first significant peak in the histogram gives the relative delay frequency for the IM reflection (earliest arrived chirp). The next significant peaks give relative delay frequencies of the PIM reflections.

(108) Distance Computation

(109) Once the frequencies are extracted for the PIM and IM, the distance to the PIM source 103 can be determined. In all instances, the distance to the PIM source is computed from the relative frequencies as follows:

(110) Relative delay of IM source: .sub.IM=f.sub.IM/m

(111) Relative delay of PIM source: .sub.PIM=f.sub.PIM/m

(112) Distance to PIM according: (.sub.PIM.sub.PIM)/2*v

(113) Chirp Structure at Transmitter

(114) The chirp signal samples in the tester is built with sampling period that is higher than the time resolution necessary to recover accurate delay estimation. A commonly used chirp structure configuration corresponds to the signal shown in FIG. 10.

(115) Assume one transmitted chirp that is captured at the tester with a sampling position that is not corresponding to the sampling position at chirp creation. A chirp product applied on it will not provide a constant frequency sinewave. In other words, the chirp product property is lost if the receive part does not sample it at the exact position the chirp was created.

(116) In the system described above, if two reflected chirps are received with overlapping timing less than the baseband sampling period, the chirp product property will be lost. This deficiency will cause an unreliable estimation of delay spacing (21).

(117) Instead of the commonly used configuration of FIG. 10, in one embodiment, the chirp wave in FIG. 11 is used. The chirp structure formula is unchanged:

(118) $\begin{array}{cc}\mathrm{chirp}\left(n\right)={\mathrm{Ae}}^{\left\{j2\left(f_{t}n+\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}& \left(16\right)\end{array}$
but the initial frequency shall be set to a negative value and the frequency change rate shall be positive. With a proper number of samples, the chirp product property is not altered from baseband sampling as the number of high frequency components are reduced. This chirp structure is used as the transmitted signal in the DL CPRI link. To keep the chirp product property the configuration parameters are:

(119) The initial frequency is computed as follows:

(120) $\begin{array}{cc}{f}_t=-m\frac{\mathrm{chirpSize}-1}{f_s}& \left(17\right)\end{array}$
Where:
chirp size is chirpSize=correlaticnMatrixSize*2(18)
the change rate is

(121) $m=\frac{1}{2\mathrm{chirpSize}}$
to guaranty chirp product property integrity all along the chirp signal.
f.sub.t: the initial frequency of the transmitted chirp signal
m: the frequency change rate of the chirp signal
f.sub.s: the I/Q frequency sampling at the CPRI/OBSAI interface

(122) The correlation squared matrix size is recommended to be around 6464.

(123) Chirp Structure at Receiver

(124) The chirp product operation on a received chirp signal will produce a sinewave. A reference chirp is used by the chirp product as indicated in (5). The reference chirp expression employed at the transmit chain differs from the one at the receive chain. Due to transmission in the RF followed by non-linear reflection, the chirp configuration is altered. Here only the third order reflection is considered. However, higher order reflections could exist with similar properties.

(125) The modified reference chirp expression to apply at the receiver:

(126) $\begin{array}{cc}\mathrm{referenceChirp}\left(n\right)=e^{\left\{j2\left(f_{r}n-\frac{1}{2}{\mathrm{mn}}^2\right)\right\}}& \left(19\right)\end{array}$
With the initial frequency equals:
f.sub.r=2f.sub.CWf.sub.t+2f.sub.dl.sub._.sub.cwf.sub.dl.sub._.sub.chirpf.sub.ul(20)
Where: f.sub.CW: sinusoidal CW frequency f.sub.t: initial frequency of the transmitted chirp f.sub.dl.sub._.sub.chirp: LTE downlink carrier carrying the chirp signal f.sub.dl.sub._.sub.cw: LTE downlink carrier carrying the CW signal (if different from this used to carry the chirp) f.sub.ul: LTE uplink carrier carrying the reflected chirp
Distance Determination Using Matched Filter

(127) Another category of algorithm that may address the delay estimation problem is the one based on matched filter. A predetermine sequence is employed that constitutes the matched filer coefficient. The targeted predetermined pattern can be a chirp signal, a PN sequence or a Zadoff Chu sequence. A correlation process is applied with the matched filter along the received samples. The correlation output results into a window of correlation peaks. The windows length corresponds to the pattern sequence length. Throughout the received process the correlation result may be sum-up to improve the detection performance.

(128) Distance Determination Using Adapted Reflectometry

(129) One issue with reflectometry method is its sensibility to phase rotation at the PIM reflection point. Prior to applying reflectometry, a correlation method can be used to estimate the phase rotation of the pattern sequence or modulated signal reflected at the PIM point. Once known, the received signal can be phase equalized to recover from the defective phase rotation. Then, the reflectometry method can apply to estimate the delay.

(130) The methods and techniques described here may be implemented in digital electronic circuitry, or with a programmable processor (for example, a special-purpose processor or a general-purpose processor such as a computer) firmware, software, or in combinations of them. Apparatus embodying these techniques may include appropriate input and output devices, a programmable processor, and a storage medium tangibly embodying program instructions for execution by the programmable processor. A process embodying these techniques may be performed by a programmable processor executing a program of instructions to perform desired functions by operating on input data and generating appropriate output. The techniques may advantageously be implemented in one or more programs that are executable on a programmable system including at least one programmable processor coupled to receive data and instructions from, and to transmit data and instructions to, a data storage system, at least one input device, and at least one output device. Generally, a processor will receive instructions and data from a read-only memory and/or a random access memory. Storage devices suitable for tangibly embodying computer program instructions and data include all forms of non-volatile memory, including by way of example semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks such as internal hard disks and removable disks; magneto-optical disks; and DVD disks. Any of the foregoing may be supplemented by, or incorporated in, specially-designed application-specific integrated circuits (ASICs) or Field Programmable Gate Arrays (FGPAs).

(131) A number of embodiments of the invention defined by the following claims have been described. Nevertheless, it will be understood that various modifications to the described embodiments may be made without departing from the spirit and scope of the claimed invention. Accordingly, other embodiments are within the scope of the following claims.