Methods for accurate use of finite dynamic range measurement data in radio path loss propagation model calibration

Description of equivalent WO 2006099144 A2

METHODS FOR ACCURATE USE OF FINITE DYNAMIC RANGE

MEASUREMENT DATA IN RADIO PATH LOSS PROPAGATION MODEL CALIBRATION

0003Claim of Priority

0004U.S. Provisional Patent Application No. 60/661,128 entitled Methods for Accurate Use of Finite Dynamic Range Measurement Data in Radio Path Loss Propagation Model Calibration, by Nathan Hunt et al, filed March 11, 2005 (Attorney Docket No. NEWW- 01000US0); and

0005U.S. Patent Application No. 11/ , entitled Methods for Accurate Use of Finite

0006Dynamic Range Measurement Data in Radio Path Loss Propagation Model Calibration, by Nathan Hunt et al, filed March 9, 2006 (Attorney Docket No. NEWW-01000US1).

0007Copyright Notice

0008A portion of the disclosure of this patent document contains material that is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the U.S. Patent and

0009Trademark Office patent file or records, but otherwise reserves all copyright rights whatsoever.

0010Field of the Invention

0011This invention relates to the field of calibrating data parameter sets for application to a signal propagation model describing data of limited dynamic range.

0012Background of the Invention

0013The effectiveness of predictive design and optimization of cell phone networks depends on the accuracy of predicted radio signal strength over the geographic area where service is to be provided. The prediction is made with a mathematical model called a propagation model. Propagation models typically predict radio signal strength given the transmitting antenna pattern and digital terrain maps and clutter maps, although they may also use additional descriptions of the physical environment such as vector data and building databases. Propagation models have a number of parameters whose values are determined by calibrating the model with measurement data. Once a model is calibrated for a particular geographic area, the required signal strength and interference predictions can be carried out. Model calibration is required to give specific loss characteristics to the physical environment description. For example, the average properties of clutter categories may vary between markets or be unknown altogether. Model calibration assigns optimal values to each clutter category. The signal strength measurement data used for model calibration is usually either

0014Continuous Wave (CW) data or scanner data. To collect CW data, a clear frequency is selected and a continuous signal is transmitted at that frequency. Data is collected over a drive route. A receiver set to the transmitting frequency collects the data in conjunction with position information from a GPS receiver. Scanner data is also collected on a drive with simultaneous collection of GPS data.

0015Scanner data is distinct from CW data in that the actual pilot signals from the network are detected, rather than the separate monotone signal used for CW. Signals can be collected from many sites simultaneously, depending on the scanner configuration.

0016Model calibration is conventionally carried out by finding the parameter set that gives the minimum least-squares (or OLS, for Ordinary Least Squares) fit to the measurement data. Under certain simple assumptions (i.e. that the errors in the data are normally distributed, and that data is an unbiased sample of the underlying population), the OLS fit is the best prediction that can be made.

0017AU measurement data has a finite dynamic range. This generally introduces bias into the measurement data set. Models created by OLS fit to the data may therefore be biased and will not give the best prediction possible.

0018The effects of the bias due to finite dynamic range are shown in Figures 1 and 2. In Figure 1, received signal strength is shown as a function of distance, along with a simple linear OLS fit to the data. In Figure 2, the same data set is shown, this time with signals below -100 dBm excluded. The resulting OLS fit predicts less attenuation with distance than does the fit to the full data set shown in Figure 1. While the OLS fit in Figure 2 fully reflects the bias in the measurement data set, the Maximum Likelihood (ML) fit (described below) corrects for this bias and gives results nearly equal to those obtained with the full data set.

0019The example data in Figures 1 and 2 is artificial in that the truncation level was applied to the data set, rather than being an inherent feature of it. In real applications, the truncation level or levels are inherent in the data and must be explicitly included in a ML model, as described below. Additionally, a propagation model is usually more complicated than the simple linear fits shown in Figures 1 and 2, and will incorporate diffraction and clutter effects among others, as well as antenna discrimination where appropriate. But the figures serve to illustrate the basic issue of bias introduced by limited dynamic range.

0020Summary of the Invention The invention consists of approaches to calibrating propagation models that removes the bias inherent in using data that is of limited dynamic range. It also covers means of implementing these methods in software.

0021The invention describes treating two types of measurement data, truncated data and censored data. Truncated data is data where observations outside of given limits are lost. A simple example of truncated data is shown in Figure 2, where measurements with values below -100 dBm are not present in the data set. Data may be truncated from below, above, or both. The truncation level or levels may also vary across the data set.

0022Censored data is similar to truncated data, the difference being that the location of the points falling outside the receiver dynamic range is given. For example, in the case of the data in Figure 2, if additionally the distance to each of the observations falling below -100 dBm were known, these would be called censored observations (or left-censored, since they are below a known value). As is the case with truncated data, data may be censored from above (right censored), from below (left censored), or both, and the censoring levels need not be constant across the data set. The invention describes how to obtain truncated data and truncation levels from measurement data. It also describes how to obtain censored data, censoring levels, and censored observations from measurement data. It describes how to obtain optimal models for truncated and censored data sets using Maximum Likelihood (ML).

0023Brief Description of the Drawings

0024Figure 1 is a chart measuring signal vs. log of distance with an OLS linear fit. Figure 2 is a chart measuring signal vs. log of distance, signal level truncated at -100 dBm, with the OLS and truncated ML linear fits to the truncated data. The ML fit corrects for the bias introduced by the truncation, and gives a result quite close to the OLS fit to the full data set shown in Figure 1.

0025Figure 3 is a flow diagram for truncated CW data. Figure 4 is a flow diagram for truncated scanner data. Figure 5 is a flow diagram for censored CW data. Figure 6 is a flow diagram for censored scanner data.

0026Detailed Description of the Invention

0027CW measurements can be treated as truncated data by using the signal sensitivity threshold as the lower truncation level, and the top end of the receiver linear range as the upper truncation level. Any observations outside these limits are excluded from the data set. Figure 3 shows a flow diagram of this process.

0028Scanner data can be treated as truncated data by excluding measurements above the linear range or below the receiver sensitivity. The sensitivity limit for scanner data differs from the limit for CW data in that, for scanner data, the sensitivity is often limited by interference. Therefore the lower truncation levels are determined from the interference level measured at each data point. The lower truncation level at a point in the case of CDMA or

0029UMTS is the measured value of I₀ (in dBm) at a point plus the minimum Ec / 1₀ (in dB) that the scanner can reliably measure. For the case of GSM data, the lower truncation level is the measured value of I (in dBm) at a point plus the minimum value of C / I (in dB) that the scanner can reliably measure. Figure 4 shows a flow diagram of this process.

0030Given a truncated data set, the optimum propagation model is obtained as follows. To help define notation, we start by giving the formulation for the OLS solution. Vector quantities are shown in bold. The OLS solution is given by finding the vector β minimizing the quantity

0031<img file="WO2006099144A2_D0001.tif" />

0032where the sum is taken over all observations i, Vj is the i^th measured value, Xi is a vector of the variables used in the calibration (the dependent variables, such as log distance, diffraction, effective antenna height, clutter terms, etc.), and β is a vector of the parameters to be optimized (constant term or intercept, slope, diffraction multiplier, clutter parameters, etc).

0033The OLS problem can be reformulated as a maximum likelihood (ML) problem, where the log of the likelihood to be maximized is

0034(2) L = Σ_i log (2 π σ) -^1/2 - Σ_i (y_i - Xiβ) ² / (2 σ ²) where y_l5 x, and β are as in equation (1) and σ, an estimate of the dispersion, is an additional model parameter to be optimized along with the components of β. The maximum likelihood solution of (2) will yield the same solution β as (1), with the optimum value of σ equal to the standard deviation of the measurements, y,. The optimum parameter set for the truncated problem is obtained by maximizing the log of the likelihood (Tobin, Breen) L,

0035(3) L = Σ₁ log (2 π σ) -^1/2 - Σ_! (y_I - x_iβ) ² / (2 σ ²)

0036- Σ, log Φ[(c - Xi β) / σ] - Z₁ log Φ[( Xi β - dO / σ]

0037where y_]5 X₁, β and σ are as in (2), the Ci are the upper truncation values, the di are the lower truncation values, and Φ is the standard normal distribution function (Breen).

0038Censored data can be generated from CW data by using the positions of points where the measurement falls below the sensitivity limit of the receiver (left censored observations) or above the receiver's linear range (right censored observations). Censored data is generated from scanner data in the same way, but the lower limit is usually determined by the measured interference at each point. The lower censoring level at a point in the case of CDMA or

0039UMTS is the measured value of Io (in dBm) at a point plus the minimum Ec / 1₀ (in dB) that the scanner can reliably measure. For the case of GSM data, the lower censoring level is the measured value of I (in dBm) at a point plus the minimum value of C / 1 (in dB) that the scanner can reliably measure.

0040Figure 5 shows a flow diagram of the process for censored CW data, while Figure 6 shows a flow diagram of the process for censored scanner data.

0041The optimum parameter set for the censored problem is obtained by maximizing the log of the likelihood, L

0042(4) L = Z₁ log (2 π σ) ^An - X₁ (y, - X₁P) ² / (2 σ ²)

0043- Σ, log Φ[(x, β - c,) / σ] - Σ_k log Φ[( d_{k -} x_k β ) / σ]

0044where the summations over i are sums over the uncensored observations, the summations over j are over the right-censored observations, the summations over k are over the left- censored observations, V₁, Xi, β and σ are as in (3), C₁ and d_k are the upper and lower censoring limits, respectively, and X_k is the vector of dependent variables at the k^th censored observation.

0045References

0046Breen, R. (1996) Regression Models: Censored, Sample Selected, or Truncated Data. Thousand Oaks: Sage Publications Tobin, J. (1958) "Estimation of relationships for limited dependent variables." Econometrica 26: 24-36