A method for calibrating smart antenna array systems in real time

<u style="single">Field of the Technology</u>

0001The present invention generally relates to a smart antenna technology of the wireless communication system, more specifically, to a method for calibrating a smart antenna array system in real time.

<u style="single">Background of the Invention</u>

0002In modem wireless communication systems, especially in a Code Division Multiple Access (CDMA) wireless communication system, smart antenna has been one of the most attractive technologies. By means of smart antenna array and the technology based on digital signal processing, wireless base stations can achieve self-adaptive beam forming of both transmitting and receiving signals, therefore greatly reduce the system interference, increase the system capacity, decrease the transmitting power and improve the receiving sensitivity.

0003In the Chinese patent for invention titled "A Time-Division Duplex Synchronous Code Division Multiple Access Wireless Communication System with Smart Antenna" (<patcit id="pcit0001" dnum="CN97104039"><text>CN 97 104039.7</text></patcit>), the structure of a wireless communication system base station with smart antenna is disclosed. The base station comprises an antenna array with one or more than one antenna elements, radio frequency (RF) cables and RF transceivers correspondingly connected. Based on the signals received from the user terminal in each antenna element of the antenna array, a baseband signal processor can obtain the space vector characteristics and direction of arrival (DOA) of the uplink signals, from which the weights of every link obtained are employed for downlink transmitting beam forming. In this way, all the functionality of a smart antenna is achieved under the circumstances of symmetrical radio wave propagation characterized as the result of time-division duplex communication.

0004In order to transmit and receive signals accurately with the smart antenna, it must be guaranteed that every antenna element, RF cable and transceiver, which the smart antenna array comprises, operate without difference, i.e. every transmitting and receiving link should have the same amplitude and phase response. The procedure and method for amplitude and phase compensation of each transmitting and receiving link are the smart antenna calibration relating to the present invention.

0005As the characteristics of electronic elements, especially active elements, differ from each other, the sensitivities thereof to operation frequency and ambient temperature are different, and the changes in the characteristics of every transmitting and receiving link due to the reasons above are different, the smart antenna calibration should be carried out periodically while the base station is in operation.

0006In the published document of the Chinese patent titled "A method and apparatus for calibrating a smart antenna array" (<patcit id="pcit0002" dnum="CN99111350"><text>CN 99 111350.0</text></patcit>) (See <figref idref="f0001">Figure 1</figref>), a calibrating link is set by an antenna element 201, a couple structure 205, a RF cable 206 and a pilot transceiver 207 sequentially connected. The couple structure 205 sets up a RF couple connection with all the antenna elements 201-1, 201-2 ... 201-N of the smart antenna, and allocates the RF signals to all antenna elements constituting the array according to need. The pilot transceiver 207 has the same structure as the other transceivers 203-1, 203-2...203-N of the base station and uses a common local oscillator 208. The pilot transceiver 207 works coherently with other transceivers and connects with the baseband signal processor 204 via a digital bus. Each antenna element connects to a RF cable and further to a transceiver, and the connected antenna element, RF cable and transceiver form a transmitting link or a receiving link. Ac, A₁, A₂,...A_N in <figref idref="f0001">Fig. 1</figref> represent the connection points between the antenna elements and the RF cables 201-1, 201-2 ... 201-N, respectively; B_C, B₁, B₂...B_N represent the connection points of the pilot transceiver 207 and the radio transceivers 203-1, 203-2...203-N with the baseband signal processor 204, respectively.

0007When making calibration, calibrate the calibration link first by using a network vector analyzer and record the receiving and transmitting transmission coefficients of the calibration link respectively, then perform the receiving calibration and transmitting calibration respectively. In the receiving calibration, the pilot transceiver transmits a signal at a given working frequency, and all the other links in the base station are set in the receiving state. Measure the outputs of all the receiving links and compute the ratio of the receiving transmission-coefficient (vector) of each link to the transmission-coefficient (vector) of a reference link. When the ratio of the amplitudes of the transmission-coefficients equals to 1, record the phase difference of each receiving link from the reference link. In the transmitting calibration, set one link after another of the base station in the transmitting state with all the other links close at the same time and the pilot transceiver will receive the signal of each transmitting link at a given working frequency, respectively; compute the ratio of the transmission-coefficient (vector) of each link during transmission to the transmission-coefficient (vector) of the reference link, and when the ratio of the amplitudes of the transmission-coefficients equals to 1, record the phase difference of every receiving link from the reference link.

0008The patent mentioned above (of which <patcit id="pcit0003" dnum="EP1204161A"><text>EP 1 204 161</text></patcit> is a family member) only relates to the general scheme of the method and apparatus for real-time calibration without a specific engineering implementation thereof, including the calibration sequence used in the transmitting and receiving calibration and the computation by the baseband signal processor, and how to perform the real-time calibration when the smart antenna is in operation. In addition, the transmitting calibration as described above is carried out with one link in the transmitting state at a time while all other links are in the receiving state, which is unfavorable for fast real-time calibration.

0009<patcit id="pcit0004" dnum="EP0805510A"><text>EP 0 805 510</text></patcit> discloses a method of self calibration of an active array system, a module under test being phase-modulated using a special command to increment the phase from pulse to pulse.

<u style="single">Summary of the Invention</u>

0010It is the object of the present invention to provide a real-time calibration method for a smart antenna array so as to periodically calibrate the smart antenna array when a base station is in operation and compute the compensation coefficients for the transmitting and the receiving links of the smart antenna array to be calibrated.

0011The technical solution to achieve the above object is as follows: pre-calibrating each antenna element of a smart antenna array (201-1 to 201-N) to obtain the transmitting compensation coefficient <maths id="math0001"><math display="inline"><msubsup><mi>c</mi><mi>k</mi><mi mathvariant="italic">TX</mi></msubsup></math><img file="EP1585231B1_D0001.tif" /></maths> and the receiving compensation coefficient <maths id="math0002"><math display="inline"><msubsup><mi>C</mi><mi>k</mi><mi mathvariant="italic">RX</mi></msubsup></math><img file="EP1585231B1_D0002.tif" /></maths> of each antenna element (201-1 to 201-N) relative to a calibration antenna element (201), <b>characterized in that:</b><ul id="ul0001" list-style="none"><li>in a transmitting calibration procedure, a plurality of transmitting links transmit each calibration signal simultaneously, and a calibration link receives a combined signal thereof, a baseband signal processor (204) processes the combined signal received by the calibration link to obtain the amplitude and phase response of each transmitting link and computes the compensation coefficient of each transmitting link on the basis of the amplitude and phase response of each transmitting link and the transmitting compensation coefficients <maths id="math0003"><math display="inline"><msubsup><mi>c</mi><mi>k</mi><mi mathvariant="italic">TX</mi></msubsup></math><img file="EP1585231B1_D0003.tif" /></maths> obtained in the pre-calibration for compensating all the downlink data of the base station;</li><li>in the receiving calibration procedure, the calibration link transmits a calibration signal, and the receiving links receive the calibration signal simultaneously, then the baseband signal processor (204) processes the calibration signals received by the receiving links to obtain the amplitude and phase response of each receiving link then computes the compensation coefficient of each receiving link on the basis of the amplitude and phase response of each receiving link and the receiving compensation coefficients <maths id="math0004"><math display="inline"><msubsup><mi>C</mi><mi>k</mi><mi mathvariant="italic">RX</mi></msubsup></math><img file="EP1585231B1_D0004.tif" /></maths> obtained in the pre-calibration for compensating all the uplink data of the base station;</li><li>the calibration signal for each antenna element is generated by a periodic cycling shift of a basic calibration sequence and the calibration signal is a calibration sequence with good anti-white-noise characteristics.</li></ul>

0012The pre-calibrating is carried out after the production of the smart antenna array. The transmitting and receiving compensation coefficients obtained will be stored. After the smart antenna array is installed on the site of the base station, the stored pre-calibration transmitting and receiving compensation coefficients obtained in the pre-calibrating are inputted into the baseband signal processor of the base station.

0013Generating a calibration signal by a periodic cycling shift of a basic calibration sequence includes: taking a binary sequence m_p as the basic calibration sequence with a length P; performing a phase equalization to the sequence m_p to generate <u style="single">m</u>_p, a complex vector for the calibration sequence; expanding the <u style="single">m</u>_p periodically to obtain a new periodical complex vector; obtaining a calibration vector for each antenna element from the <u style="single">m;</u> generating a calibration signal for each antenna element from the calibration vector for each antenna element.

0014The length of said basic calibration sequence is W×N and the length of said calibration sequence is W×N+W-1, where N is the number of antenna elements in the antenna array, and W is the window length in channel estimation for each transmitting or receiving link.

0015Said transmitting calibration and receiving calibration are periodically performed in the idle gap of the mobile communication system.

0016In a TD-SCDMA system, said transmitting calibration and receiving calibration are performed in the protective gap (GP) between the uplink pilot time-slot and the downlink pilot time-slot in a frame.

0017Said computing the compensation coefficient of each transmitting link in the transmitting calibration comprises: first, obtaining the channel impulse response of each transmitting link; second, computing the amplitude and phase response of the path between each transmitting link, including the transceiver, and the calibration link antenna element; third, multiplying the amplitude and phase response with the transmitting compensation coefficient of the corresponding link obtained in pre-calibrating, and then obtaining the transmitting compensation coefficient of each link.

0018Said computing the compensation coefficient of each receiving link in the receiving calibration comprises: first, obtaining the channel impulse response of each receiving link; second, computing the amplitude and phase response of the path between each receiving link, including the calibration link antenna element, and the transceiver; third, multiplying the amplitude and phase response with the receiving compensation coefficient of the corresponding link obtained in the pre-calibrating, and then obtain the receiving compensation coefficient of each link.

0019The real-time calibration method of this invention, for which it is necessary to set a calibration link specially for realizing the calibration function (as described in the background of the invention) which consists of an antenna element, a feeder cable and a pilot transceiver: before the delivery of the antenna array, first pre-calibrating the antenna array to obtain compensation coefficients of each antenna element relative to a calibration antenna element; then storing the compensation coefficients in the network operation and maintenance equipment; after the on-site installation of the smart antenna array, loading the compensation coefficient into the base station.

0020Said calibration performed periodically while the base station is in operation further comprises: in the transmitting calibration, the transmitting links simultaneously transmitting a fixed level calibration sequence, which is received by the calibration link as a combined signal thereof; in the receiving calibration, the calibration link transmitting a fixed level calibration sequence, which is received simultaneously by the receiving links.

0021By computing the received signals according to the computation method provided by this invention, the compensation coefficients of the transmitting and receiving links of the smart antenna array can be obtained so that the real-time calibration can be accomplished.

0022The fixed level calibration sequence employed is generated by a periodic cycling shift to a basic calibration sequence.

0023The method of the present invention has the advantages of short computation time and simple controls. It is especially suitable for a smart antenna array in the third generation mobile communication system with high chip rate.

0024Although the above solution is proposed mainly for CDMA wireless communication systems, the method of the present invention is fully applicable, after simple modification, to frequency division multiple access and time division multiple access wireless communication systems. It can not only be used to calibrate a smart antenna operating in TDD mode, but also a smart antenna operating in FDD mode.

0025By calibrating a smart antenna array with the method of the present invention, beams of the smart antenna array are greatly improved and the transmitting powers are decreased, resulting in a good implementation of the smart antenna technology.

<u style="single">Brief Description of the Drawings</u>

0026<ul id="ul0002" list-style="none" compact="compact"><li><figref idref="f0001">Figure 1</figref> is a schematic diagram for the configuration of a base station with a smart antenna array comprising a calibration link.</li><li><figref idref="f0002">Figure 2</figref> is a schematic diagram for the layout of pre-calibrating of a smart antenna array.</li></ul>

<u style="single">Detailed Description of the Invention</u>

0027The method of the present invention is proposed on the basis of an antenna array which is a passive microwave (radio frequency) network. Characteristics of mutual coupling between each antenna element of this antenna array and the calibration antenna element remain unchanged at a given working frequency provided that the design of the antenna array product has been finalized and the structure thereof is fixed. Therefore, before delivery of the antenna array, each antenna element of the antenna array can be tested relative to the calibration antenna element at a given working frequency or can be pre-calibrated to obtain the compensation coefficient of each antenna element relative to the calibration antenna element which is then stored in the network management database as the pre-calibration data. After the antenna array has been installed on site, the antenna array pre-calibration data are loaded to the base station by the network operation and maintenance equipment, such as OMC_R or LMT. In this way, the antenna array can be calibrated with the real-time calibration method according to the present invention when the antenna array starts into operation.

0028The method of the invention can be employed in a typical time-division duplex (TDD) CDMA base station equipped with a smart antenna. The configuration of the base station is shown in <figref idref="f0001">Fig. 1</figref>. The base station comprises N identical antenna elements 201-1, 201-2 ... 201-N; N identical feeder cables 202-1, 202-2 ... 202-N; N RF transceivers 203-1, 203-2 ... 203-N that work coherently; and an appropriate baseband signal processor 204. In order to implement the calibration, a calibration (reference) link is also established, which consists of a RF coupling structure 205, a calibration antenna element 201, a feeder cable 206 and a pilot transceiver 207, wherein the pilot transceiver 207 works coherently with N transceivers 203-1, 203-2 ... 203-N, and uses a common local oscillator 208. N transceivers 203-1, 203-2 ... 203-N and the pilot transceiver 207 connect with the baseband signal processor 204 via a data bus.

0029The real-time calibration method of the present invention comprises the following key steps: <ul id="ul0003" list-style="none" compact="compact"><li>The First step: before delivery of the smart antenna array, pre-calibrating each antenna element using a radio frequency (microwave) network vector analyzer to obtain the compensation coefficient for each antenna element relative to the calibration antenna element.</li></ul>

0030Refer to <figref idref="f0002">Fig. 2</figref>, in pre-calibration, with one end of the radio frequency network vector analyzer 21 connected to the antenna element of the calibration link 201 through the point Ac and the other end connected to the antenna elements 201-1, 201-2 ... 201-N of the smart antenna array 208 in proper order through the points A₁, A₂...A_N, the radio frequency network vector analyzer 21 performs the transmitting and receiving pre-calibration respectively.

0031If the structure of the smart antenna is rather firm by design, it can be recognized that the channel characteristics between each antenna element 201-1, 201-2 ... 201-N and the calibration antenna unit 201 remain unchanged on the whole with the environmental conditions at a fixed working frequency, given that there is no disruption of the relative location. So it is possible to perform the pre-calibration measurement using a radio frequency network vector analyzer.

0032In the transmitting pre-calibration, a fixed level digital signal is transmitted by each antenna element of 201-1, 201-2 ... 201-N, respectively, and the signal is received by the antenna element of the calibration link 201; the radio frequency network vector analyzer 21 measures and computes the transmitting compensation coefficient <maths id="math0005"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">TX</mi></msubsup></math><img file="EP1585231B1_D0005.tif" /></maths> (TX represents transmitting) between <i>A</i>_i ( i=1 , ...N ) and <i>A_c</i>, thus the transmitting compensation coefficient <maths id="math0006"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">TX</mi></msubsup></math><img file="EP1585231B1_D0006.tif" /></maths> between each antenna element of 201-1, 201-2 ... 201-N and the calibration antenna element 201 is obtained. In the receiving pre-calibration, the antenna element of the calibration link 201 transmits a fixed level digital signal, which is received by each antenna element 201-1, 201-2 ... 201-N; the radio frequency network vector analyzer 21 measures and computes the receiving compensation coefficient <maths id="math0007"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">RX</mi></msubsup></math><img file="EP1585231B1_D0007.tif" /></maths> (RX represents receiving) between <i>A</i>_c and <i>A</i>_i (i=1 , ...N), thus the receiving compensation coefficient <maths id="math0008"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">RX</mi></msubsup></math><img file="EP1585231B1_D0008.tif" /></maths> between the calibration antenna element 201 and each antenna element 201-1, 201-2 ... 201-N is obtained.

0033In general, in a base station of TDD CDMA system, since each transceiver is connected to the same antenna element (i.e. the transmitting and receiving links have a common antenna element), the measured transmitting compensation coefficient of each antenna element equals to the receiving compensation coefficient thereof, i.e. <maths id="math0009"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">TX</mi></msubsup><mo>=</mo></math><img file="EP1585231B1_D0009.tif" /></maths>

0034<maths id="math0010"><math display="inline"><msubsup><mi>c</mi><mi>i</mi><mi mathvariant="italic">RX</mi></msubsup><mn>.</mn></math><img file="EP1585231B1_D0010.tif" /></maths>

0035In an FDD CDMA system, however, when the smart antenna technology is employed, different antenna arrays are usually used for transmitting and receiving, respectively, in order to isolate the transmitting link from the receiving link. Therefore, each antenna element of the two antenna arrays should be measured and pre-calibrated respectively.

0036The second step: inputting the above pre-calibration result (the transmitting compensation coefficient and the receiving compensation coefficient) in the network operation and maintenance equipment. After the antenna array has been installed on site, the antenna array compensation coefficients are loaded to the baseband signal processor of the base station to which the antenna array is connected by the network operation and maintenance equipment, such as OMC_R or LMT.

0037The third step is carried out while the base station starts operation or is in operation. This step comprises: generating a calibration sequence; performing the transmitting calibration; performing the receiving calibration; and computing the transmitting and receiving compensation coefficients.

0038The calibration sequence is generated by a periodic cycling shift of a basic calibration sequence selected with good anti-white-noise characteristics. The length of the basic calibration sequence P is W×N, where N is the number of operating antenna elements of the antenna array, and W is the window length in the channel estimation of each link. The length of the calibration sequence when performing the transmitting and receiving calibration is W×N + W - 1, that is, P + W - 1.

0039Since W relates only to the inconsistency of the hardware time-delay of each antenna element (usually very small), the calibration duration is very short when carrying out calibration by using the method of the present invention. In some systems where the calibration duration is limited, N can take a larger value in order to have a larger antenna gain for the system.

0040The procedure for generating the calibration sequence from the basic calibration sequence is as follows: <ol id="ol0001" compact="compact"><li>(1) Take a binary sequence <b>m_P</b> with the length of P as the basic calibration sequence, where <b>m</b>_P<i>=</i> (<i>m</i>₁, <i>m</i>₂,...,<i>m_P</i>). P = W×N (select a power of 2 as P to simplify the computation);</li><li>(2) In order to avoid an abrupt phase change and increase the calibration accuracy, make a phase equalization of the basic calibration sequence <b>m_P</b> to generate a complex vector of the calibration sequence <b><u style="single">m</u>_P; <u style="single">m</u>_P</b> = (<i><u style="single">m</u></i>₁,<i><u style="single">m</u></i>₂,...,<i><u style="single">m</u>_P</i>), where the element <i><u style="single">m</u>_i</i> is derived from the corresponding element <b><i>m_i</i></b> of the sequence <b>m_P</b>, <b><i><u style="single">m</u>_i</i></b> = (j)^<i>i</i>-1· <i>m_i</i>, (i = 1 ...P), and j is the square root of (-1), in this way the phase equalization is implemented;</li><li>(3) In order to generate the calibration sequence for each antenna element, the basic calibration sequence is periodically expanded to obtain a new complex vector <b><u style="single">m</u></b>, <b><u style="single">m</u></b> = (<i><u style="single">m</u></i>₁,<i><u style="single">m</u></i>₂,...,<u style="single">m</u>_{<i>i</i><sub2>max</sub2>})=(<i><u style="single">m</u></i>₂,<i><u style="single">m</u></i>₃,...,<i><u style="single">m</u>_p</i> ,<i><u style="single">m</u></i>₁,<i><u style="single">m</u></i>₂,...,<i><u style="single">m</u>_p</i>);</li><li>(4) The calibration sequence vector of each antenna element can be obtained from this periodical complex vector, and a fixed level calibration signal is thereby generated: the sequence vector <maths id="math0011"><math display="inline"><msup><munder><mi mathvariant="bold">m</mi><mo>̲</mo></munder><mfenced><mi>k</mi></mfenced></msup><mo>=</mo><mfenced><msubsup><munder><mi>m</mi><mo>̲</mo></munder><mn>1</mn><mfenced><mi>k</mi></mfenced></msubsup><msubsup><munder><mi>m</mi><mo>̲</mo></munder><mn>2</mn><mfenced><mi>k</mi></mfenced></msubsup><mo>…</mo><msubsup><munder><mi>m</mi><mo>̲</mo></munder><msub><mi>L</mi><mi>m</mi></msub><mfenced><mi>k</mi></mfenced></msubsup></mfenced><mo>,</mo></math><img file="EP1585231B1_D0011.tif" /></maths> k = 1 ... N, Lm = P+W-1 (Lm is the window length, i.e. the length of the transmitting calibration sequence).</li></ol>

0041An element of the sequence vector <maths id="math0012"><math display="inline"><msubsup><munder><mi>m</mi><mo>̲</mo></munder><mi>i</mi><mfenced><mi>k</mi></mfenced></msubsup><mo>=</mo><msub><munder><mi>m</mi><mo>̲</mo></munder><mrow><mi>i</mi><mo>+</mo><mfenced><mi>N</mi><mo>-</mo><mi>k</mi></mfenced><mo>⁢</mo><mi>W</mi></mrow></msub><mo>,</mo></math><img file="EP1585231B1_D0012.tif" /></maths><i>i</i> = 1,...,<i>L_m</i> and <i>k =</i> 1,...,<i>N</i> (k represents any antenna element of an operating antenna array).

0042At the same time, a vector S relating to the basic calibration sequence should be computed for this invention as well, which is stored in the baseband signal processor as a constant vector for computing the compensation coefficients when performing the transmitting and receiving calibrations: <maths id="math0013"><math display="block"><mi>S</mi><mo>=</mo><mn>1.</mn><mo>/</mo><mi mathvariant="italic">fft</mi><mfenced><msub><munder><mi mathvariant="bold">m</mi><mo>̲</mo></munder><mi mathvariant="normal">P</mi></msub></mfenced><mo>=</mo><mn>1.</mn><mo>/</mo><mi mathvariant="italic">fft</mi><mfenced><msub><munder><mi>m</mi><mo>̲</mo></munder><mn>1</mn></msub><msub><munder><mi>m</mi><mo>̲</mo></munder><mn>2</mn></msub><mo>…</mo><msub><munder><mi>m</mi><mo>̲</mo></munder><mi>P</mi></msub></mfenced></math><img file="EP1585231B1_D0013.tif" /></maths> where ./ represents a point-division and fft represents the Fast Fourier Transform Algorithm.

0043The selecting the basic calibration sequence step refers to selecting a binary sequence which makes S have a minimum norm and has a length P. The transmitting calibration comprises the following steps: each antenna element transmitting a fixed level calibration sequence simultaneously, and the calibration link receiving the combined signal thereof. With the algorithm provided by the present invention, the baseband signal processor processing the signal data received by the calibration link, computing the amplitude and phase response of each transmitting link, and then computing the compensation coefficient (including the amplitude and the phase compensation) for each transmitting link according to the compensation coefficient (transmitting compensation coefficient) thereof obtained during pre-calibration, by which all the downlink data of the base station are compensated at the baseband signal processor.

0044Taking the base station in <figref idref="f0001">Fig. 1</figref> as an example, N transmitting links transmit calibration sequence vectors <b><u style="single">m</u></b>^(<i>k</i>) with a certain power level at point <i>B_k</i> (k=1, ...N). Through the transceivers 203-1 ... 203-N, the feeder cables 202-1 ... 202-N, the antenna elements 201-1 ... 201-N and the antenna array couple structure 205, the fixed level signals are received by the calibration link antenna element 201. The baseband signal processor computes the received data from the calibration link (201, 206 and 207) to obtain the amplitude and phase response of each transmitting link <i>B</i>_k → <i>A</i>_<i>k</i>· In fact, the amplitude and phase response of the link B_k→ A_k→ A_C→ B_C is needed. Since the amplitude and phase response of the path A_k→ A_C has been obtained in the pre-calibration, only the amplitude and phase response of the path <i>B</i>_k → <i>A_k</i> needs to be computed.

0045Suppose R is the complex vector received from the baseband signal processor after each transmitting calibration sequence signal of the antenna elements 201-1 ... 201-N is accumulated in the calibration link antenna unit 201, the receiving sequence: <maths id="math0014"><math display="block"><mi>R</mi><mo>=</mo><mfenced><msub><mi>r</mi><mn>1</mn></msub><msub><mi>r</mi><mn>2</mn></msub><mo>…</mo><msub><mi>r</mi><mi>l</mi></msub></mfenced><mspace width="1em" /><mi>l</mi><mo>=</mo><mi mathvariant="normal">p</mi><mo mathvariant="normal">+</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">×</mo><mfenced><mi mathvariant="normal">w</mi><mo mathvariant="normal">-</mo><mn mathvariant="normal">1</mn></mfenced><mo>,</mo></math><img file="EP1585231B1_D0014.tif" /></maths>

0046from which a section can be intercepted with a length equaling to that of the basic sequence P, that is, <i><u style="single">R</u>_{<u style="single">p</u>}</i> = (<i><u style="single">r</u></i>₁,<i><u style="single">r</u></i>₂,...,<i><u style="single">r</u>_{<u style="single">p</u>}</i>). Suppose the interception is made in the middle of the sequence (there could be various ways of interception) as represented in the following formula, <i><u style="single">R</u>_{<u style="single">p</u>}</i> = (<i>r</i>_<i>w</i>-1, <i>r_w</i>,..., <i>r</i>_{<i>w</i>+<i>p</i>-2}).

0047A Channel Impulse Response (CIR) sequence with a length of P can be obtained by operation of the following formula: <i>CIR</i> = (<i>c</i>₁,<i>c</i>₂,...,<i>c_p</i>)= <i>ifft</i>(<i>fft(<u style="single">R</u></i>_P.<i>S</i>)), where represents a point multiplication, fft represents the Fast Fourier Transform Algorithm, ifft represents the Inverse Fast Fourier Transform Algorithm, S is the constant vector obtained as described above. <maths id="math0015"><math display="block"><mi>Compute</mi><mspace width="1em" /><msub><mi mathvariant="italic">CIR</mi><mi>k</mi></msub><mo>=</mo><mi>f</mi><mspace width="1em" /><mi>max</mi><mo>⁢</mo><mfenced><msub><mi>c</mi><mrow><mi>w</mi><mo>×</mo><mfenced><mi>k</mi><mo>-</mo><mn>1</mn></mfenced><mo>+</mo><mn>1</mn></mrow></msub><mo>…</mo><msub><mi>c</mi><mrow><mi>w</mi><mo>×</mo><mi>k</mi></mrow></msub></mfenced><mo>,</mo><mi mathvariant="normal">k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>…</mo><mi mathvariant="normal">N</mi><mo>,</mo></math><img file="EP1585231B1_D0015.tif" /></maths><i>f</i> max is an interpolation function to evaluate the peak between the channel estimation results <i>c</i>_{w×(<i>k</i>-1)+1}∼<i>c</i>_{<i>w</i>×<i>k</i>} of the k^th transmitting link (the specific value depends on the required computation accuracy), <i>CIR_k</i> is a complex number comprising the amplitude and phase response of the path <i>B</i>_k → <i>A_c</i> of the k^th link.

0048Multiply the <i>CIR_k</i> with the transmitting compensation coefficient <maths id="math0016"><math display="inline"><msubsup><mi>c</mi><mi>k</mi><mi mathvariant="italic">TX</mi></msubsup></math><img file="EP1585231B1_D0016.tif" /></maths> of the path A_k→ A_C of the k^th link obtained in the pre-calibration, then obtain: <maths id="math0017"><math display="block"><msubsup><mi mathvariant="italic">CIR</mi><mi>k</mi><mi>ʹ</mi></msubsup><mo>=</mo><msub><mi mathvariant="italic">CIR</mi><mi>k</mi></msub><mo>×</mo><msubsup><mi>c</mi><mi>k</mi><mi mathvariant="italic">TX</mi></msubsup><mo>,</mo><mi mathvariant="normal">k</mi><mo>=</mo><mn>1</mn><mspace width="1em" /><mo>…</mo><mspace width="1em" /><mi>N</mi><mo>,</mo></math><img file="EP1585231B1_D0017.tif" /></maths> where, <maths id="math0018"><math display="inline"><msubsup><mi mathvariant="italic">CIR</mi><mi>k</mi><mi>ʹ</mi></msubsup></math><img file="EP1585231B1_D0018.tif" /></maths> is also a complex number, which contains the amplitude and phase response of the path <i>B</i>_k → <i>A_k</i> of the k^th link. By using the amplitude and phase response above, the transmitting compensation coefficient of the k^th link can be obtained.

0049The receiving calibration comprises the following steps: the calibration link transmitting a fixed level calibration sequence signal, which is received by each receiving link simultaneously. The baseband signal processor computing the amplitude and phase response of each receiving link on the basis of the received data at each receiving link, by which and the receiving compensation coefficient obtained in the pre-calibration the compensation coefficient (including the amplitude and phase compensation) of each receiving link is computed and obtained. With the compensation coefficients, all downlink data of the base station can be compensated in the baseband signal processor.

0050Also taking the base station in <figref idref="f0001">Fig. 1</figref> as an example, the calibration link (201, 206 and 207) transmits a calibration vector signal with a certain power level <b><u style="single">m</u></b>^(<i>k</i>) (k = 1, ... N) at point <i>B</i>_c. The signal is received by each receiving link through the couple structure 205, each antenna element of the antenna array 201-1 ... 201-N each feeder cable 202-1, ...202-N, each transceiver 203-1, ...203-N.The baseband signal processor 204 computes the data received from each receiving link to obtain the amplitude and phase response of each receiving link (<i>A_k</i> → <i>B_k</i>). In fact, the amplitude and phase response of the path <i>B_c</i> → <i>A_c</i> → <i>A_k</i> → <i>B_k</i> is needed, since the amplitude and phase response of the path <i>A_c</i> → <i>A_k</i> have been obtained in pre-calibration, only the amplitude and phase response of the path A<i>_k</i> → <i>B_k</i> needs to be computed.

0051Suppose <i>R^k</i> is the complex vector of each link received in the baseband signal processor 204, the receiving sequence: <maths id="math0019"><math display="block"><msup><mi>R</mi><mi>k</mi></msup><mo>=</mo><mfenced><msubsup><mi>r</mi><mn>1</mn><mi>k</mi></msubsup><msubsup><mi>r</mi><mn>2</mn><mi>k</mi></msubsup><mo>…</mo><msubsup><mi>r</mi><mi>l</mi><mi>k</mi></msubsup></mfenced><mo>,</mo><mn>1</mn><mo>=</mo><mi mathvariant="normal">p</mi><mo mathvariant="normal">+</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">×</mo><mfenced><mi mathvariant="normal">w</mi><mo mathvariant="normal">-</mo><mn mathvariant="normal">1</mn></mfenced><mo mathvariant="normal">,</mo><mi mathvariant="normal">k</mi><mo mathvariant="normal">=</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">,</mo><mn mathvariant="normal">..</mn><mi mathvariant="normal">N</mi><mn>.</mn></math><img file="EP1585231B1_D0019.tif" /></maths>

0052Intercept a section from the receiving sequence with a length that equals to the length of the basic calibration sequence P, <maths id="math0020"><math display="inline"><msubsup><munder><mi>R</mi><mo>̲</mo></munder><mi>P</mi><mi>k</mi></msubsup><mo>=</mo><mfenced><msubsup><munder><mi>r</mi><mo>̲</mo></munder><mn>1</mn><mi>k</mi></msubsup><msubsup><munder><mi>r</mi><mo>̲</mo></munder><mn>2</mn><mi>k</mi></msubsup><mo>…</mo><msubsup><munder><mi>r</mi><mo>̲</mo></munder><mi>p</mi><mi>k</mi></msubsup></mfenced><mn>.</mn></math><img file="EP1585231B1_D0020.tif" /></maths> Suppose the interception is made in the middle of the sequence (there could be various ways of interception) as represented in the following formula, <maths id="math0021"><math display="inline"><msub><msup><munder><mi>R</mi><mo>̲</mo></munder><mi>k</mi></msup><mi>p</mi></msub><mo>=</mo><mfenced><msubsup><mi>r</mi><mrow><mi>w</mi><mo>-</mo><mn>1</mn></mrow><mi>k</mi></msubsup><msubsup><mi>r</mi><mi>w</mi><mi>k</mi></msubsup><mo>…</mo><msubsup><mi>r</mi><mrow><mi>w</mi><mo>+</mo><mi>p</mi><mo>-</mo><mn>2</mn></mrow><mi>k</mi></msubsup></mfenced></math><img file="EP1585231B1_D0021.tif" /></maths> k = 1, ..., N.

0053A Channel Impulse Response (CIR) sequence with the length of P can be obtained by operation of the following formula: <maths id="math0022"><math display="inline"><msup><mi mathvariant="italic">CIR</mi><mi>k</mi></msup><mo>=</mo><mfenced><msubsup><mi>c</mi><mn>1</mn><mi>k</mi></msubsup><msubsup><mi>c</mi><mn>2</mn><mi>k</mi></msubsup><mo>…</mo><msubsup><mi>c</mi><mi>p</mi><mi>k</mi></msubsup></mfenced><mo>=</mo><mi mathvariant="italic">ifft</mi><mfenced><mi mathvariant="italic">fft</mi><mfenced><msub><msup><munder><mi>R</mi><mo>̲</mo></munder><mi>k</mi></msup><mi>p</mi></msub><mn>.</mn><mi>S</mi></mfenced></mfenced><mo>,</mo></math><img file="EP1585231B1_D0022.tif" /></maths> k = 1, ...N, where represents a point multiplication, fft represents the Fast Fourier Transform Algorithm, ifft represents the Inverse Fast Fourier Transform Algorithm, and S is the constant vector obtained as described above.

0054Compute <maths id="math0023"><math display="inline"><msub><mi mathvariant="italic">CIR</mi><mi>k</mi></msub><mo>=</mo><mi>f</mi><mspace width="1em" /><mi>max</mi><mfenced><msubsup><mi>c</mi><mn>1</mn><mi>k</mi></msubsup><msubsup><mi>c</mi><mn>2</mn><mi>k</mi></msubsup><mo>…</mo><msubsup><mi>c</mi><mi>p</mi><mi>k</mi></msubsup></mfenced><mo>,</mo></math><img file="EP1585231B1_D0023.tif" /></maths> k = 1, ... N, <i>f</i> max is an interpolation function to evaluate the peak between the channel estimation results <maths id="math0024"><math display="inline"><msubsup><mi>c</mi><mn>1</mn><mi>k</mi></msubsup><mo>∼</mo><msubsup><mi>c</mi><mrow><mi>w</mi><mo>×</mo><mi>k</mi></mrow><mi>k</mi></msubsup></math><img file="EP1585231B1_D0024.tif" /></maths> of the k^th receiving link (the specific value depends on the required computation accuracy), <i>CIR_k</i> is a complex number containing the amplitude and phase response of the path <i>B</i>_k → <i>A_c</i> of the k^th link.

0055Multiply the <i>CIR_k</i> with the receiving compensation coefficient <maths id="math0025"><math display="inline"><msubsup><mi>C</mi><mi>k</mi><mi mathvariant="italic">RX</mi></msubsup></math><img file="EP1585231B1_D0025.tif" /></maths> of the path <i>Ac</i> → <i>A_k</i> of the k^th link obtained in the pre-calibration, and obtain: <maths id="math0026"><math display="block"><msubsup><mi mathvariant="italic">CIR</mi><mi>k</mi><mi>ʹ</mi></msubsup><mo>=</mo><msub><mi mathvariant="italic">CIR</mi><mi>k</mi></msub><mo>×</mo><msubsup><mi>c</mi><mi>k</mi><mi mathvariant="italic">RX</mi></msubsup><mo>,</mo><mi mathvariant="normal">k</mi><mo>=</mo><mn>1</mn><mspace width="1em" /><mo>…</mo><mspace width="1em" /><mi>N</mi><mo>,</mo></math><img file="EP1585231B1_D0026.tif" /></maths> where <maths id="math0027"><math display="inline"><msubsup><mi mathvariant="italic">CIR</mi><mi>k</mi><mi>ʹ</mi></msubsup></math><img file="EP1585231B1_D0027.tif" /></maths> is also a complex number, containing the amplitude and phase response of the path <i>A_k</i> → <i>B_k</i> of the k^th link, by means of which the receiving compensation coefficient of the k^th link can be obtained.

0056The formulas adopted in this invention for computing the compensation coefficient are as follows: <ul id="ul0004" list-style="none" compact="compact"><li>First, compute the mean power of each link, which is carried out for the transmitting link and receiving link respectively, that is, <i>CIR_k^'</i> in the following formula is the result of the transmitting calibration and the receiving calibration, respectively.</li></ul><maths id="math0028"><math display="block"><mi>Mean_power</mi><mo>=</mo><mfenced><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mfenced><mi>abs</mi><mfenced><msubsup><mi>CIR</mi><mi mathvariant="normal">k</mi><mi mathvariant="normal">ʹ</mi></msubsup></mfenced></mfenced><mn>2</mn></msup></mfenced><mo>/</mo><mi mathvariant="normal">N</mi><mo>,</mo><mfenced><mi>abs is an amplitude function</mi></mfenced><mn>.</mn></math><img file="EP1585231B1_D0028.tif" /></maths>

0057In this way, the transmitting compensation coefficient and the receiving compensation coefficient can be computed respectively by the following formula. While computing the compensation power of the receiving link, use the mean power of the receiving link and the <i>CIR_k^'</i> obtained in the receiving calibration, and while computing the compensation power of transmitting link, use the mean power of the transmitting link and the <i>CIR_k^'</i> obtained in the transmitting calibration: <maths id="math0029"><math display="block"><msup><mi>Corr_factor</mi><mi mathvariant="normal">k</mi></msup><mo mathvariant="normal">=</mo><mi>sqrt</mi><mfenced><mi>Mean_power</mi></mfenced><mo mathvariant="normal">/</mo><msubsup><mi>CIR</mi><mi mathvariant="normal">k</mi><mi mathvariant="normal">ʹ</mi></msubsup><mo mathvariant="normal">,</mo><mi mathvariant="normal">k</mi><mo mathvariant="normal">=</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">,</mo><mo mathvariant="normal">…</mo><mi mathvariant="normal">N</mi><mn mathvariant="normal">.</mn></math><img file="EP1585231B1_D0029.tif" /></maths>

0058Although the descriptions above are given with reference to the structure of a TDD CDMA base station, the calibration method is independent of the transmitting calibration and the receiving calibration, so the method of the present invention can also be implemented in an FDD CDMA base station which uses different smart antenna arrays to transmit and receive signals.

0059In the present invention, data transmitted and received are compensated respectively by the baseband signal processor using the transmitting compensation coefficient and the receiving compensation coefficient computed. Thus a software implementation of real-time calibration of a smart antenna array is achieved.

0060In practice, it is not likely for a mobile communication system to run in full load all the time. There would always be some idle gaps which can be used for real-time calibration. For TD-SCDMA system, a third generation mobile communication system, the Gad Period (GP) between the Uplink Pilot Time-Slot (UpPTS) and the Downlink Pilot Time-Slot (DwPTS) in a frame can be used for the real-time calibration.

0061The calibration by this method can be periodically performed while the base station is in operation.

0062Any person skilled in the art of smart antenna calibration with an understanding of the basic principles thereof can easily implement the real-time calibration of a smart antenna array by referring to the method of this invention.