Method and apparatus for modeling the propagation of wireless signals in buildings
Abstract
Techniques for predicting RF signal propagation in indoor environments are disclosed. A typical technique comprises six distinct phases. In one phase, the mean wall separation, d¯, of one floor of a building is estimated. In a second phase, a reflection coefficient, s, is estimated for the floor in general. In a third phase, a number of trial RF propagation measurements are made to gather empirical data about the RF propagation characteristics of the floor. In the fourth phase, the mean wall separation, d¯, the reflection coefficient, s, and the trial RF propagation measurements are fit, using well-known techniques, into a wireless propagation model, such asP(r) = P02πDπξ2re-r/ξ, where P(r) is the measured or predicted power at a distance, r, from the transmitter, and ξ is the parameter that is fit to the empirical data. In the fifth phase the wireless propagation model is used to predict the RF signal strength throughout the floor from a base station at a given location, and in the six phase one or more base stations are installed in the building based on the results predicted by the wireless propagation model. The fifth and sixth phases can be thereafter repeated for other building with sufficiently similar characteristics.

Term
Term ended
Projected expiry passed 15 September 2018, 8 years ago.
- Priority
- Filed
- Published
- Projected expiry
- Today
29 claims: 4 independent, 25 dependent
- 1A method for placing a wireless base station inside a building, said method comprising the steps of:estimating a mean wall separation, d ¯ , of a floor of said building;creating a wireless propagation model based on said mean wall separation, d ¯ ;and placing said wireless base station in said building at a location based on said wireless propagation model.
- 8A building comprising one or more floors, with a wireless base station or stations being disposed on a given floor, with the number and location of said wireless base stations selected such that the signal from said wireless base station or stations has a power level essentially everywhere on said given floor above a threshold;CHARACTERIZED IN THAT the number and location of said wireless base station or stations is based on a wireless propagation model that involves estimation of a mean wall separation, d ¯ , of said given floor.
- 15A method for placing a wireless base station inside a building, said method comprising the steps of estimating a mean wall separation, d ¯ x , in a first direction of a floor of said building;estimating a mean wall separation, d ¯ y , in a second direction of said floor, wherein said second direction is orthogonal to said first direction;creating a wireless propagation model based on said mean wall separation, d ¯ x , in said first direction of said floor and said mean wall separation, d ¯ y , in said second direction of said floor;and placing said wireless base station in said building at a location based on said wireless propagation model.
- 22A building comprising one or more floors, with a wireless base station or stations being disposed on a given floor, with the number and location of said wireless base stations selected such that the signal from said wireless base station or stations has a power level essentially everywhere on said given floor above a threshold;CHARACTERIZED IN THAT the number and location of said wireless base station or stations is based on a wireless propagation model that involves estimation of a mean wall separation, d ¯ x , in a first direction of said given floor and estimation of a mean wall separation, d ¯ y , in a first direction of said given floor.
Independent claims4
57 paragraphs, as filed
<u>Field of the Invention</u>
0001The present invention relates to wireless telecommunications in general, and, more particularly, to a technique for modeling the propagation of wireless signals in buildings.
<u>Background of the Invention</u>
0002Indoor wireless telecommunications has been the subject of intense investigation in recent years for both voice and data communication. One particular area of investigation is how to ensure the adequate propagation of a wireless signal from a base station to a wireless terminal through a typical indoor environment. The walls, furniture and other objects in atypical indoor environment scatter the wireless signal and thus produce a complex multipath channel in which the signal propagation characteristics are substantially more difficult to predict than those in outdoor contexts.
0003Typically, the signal propagation characteristics are considered when designing and installing an indoor wireless telecommunications system. In particular, the signal propagation characteristics are advantageously considered when determining how many base stations are needed to provide coverage for a building and where in the building those base stations should be located. Because base stations are typically expensive to install and operate, it is advantageous to be able to determine how to provide the necessary coverage for the building with the fewest number of base stations. To do his, several techniques have been developed for modeling the propagation of wireless signals indoors.
0004One technique in the prior art for measuring and modeling indoor signal propagation is adapted from the power-law decay model used in modeling outdoor environments. The power-law decay model assumes that the base station's antenna is high above the ground and that there is line-of-sight propagation to the wireless terminal. In this case, the mean power, <i>P</i>, received at the wireless terminal decays in inverse proportion to the square of the distance from the transmitter,<maths id="math0001" num="(Eq. 1)"><math display="block"><mrow><mtext mathvariant="italic">P∝</mtext><mtext></mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><msup><mrow><mtext mathvariant="italic">r</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext>,</mtext></mrow></math><img file="EP0903880A2_D0001.tif" /></maths> up to some break-point. Beyond that breakpoint, the mean power at the wireless terminal decays in inverse proportion to the fourth power of the distance from the transmitter:<maths id="math0002" num="(Eq. 2)"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>∝ </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><msup><mrow><mtext mathvariant="italic">r</mtext></mrow><mrow><mtext>4</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP0903880A2_D0002.tif" /></maths> The location of the break-point is determined by the location at which the ground bounce signal interferes with the line-of-sight signal. For indoor environments, Equation 1 has been adapted to Equation 3, where γ is fit to empirical trial measurements of the building of interest, in well-known fashion:<maths id="math0003" num="(Eq. 3)"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>∝ </mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><msup><mrow><mtext mathvariant="italic">r</mtext></mrow><mrow><mtext>γ</mtext></mrow></msup></mrow></mfrac></mrow></math><img file="EP0903880A2_D0003.tif" /></maths> Unfortunately, Equation 3 typically does not provide a satisfactory model of the building of interest for aiding a designer in designing a wireless system for the building. The use of Equation 3 does, however, have advantages. First, the empirical data needed to fit γ is easily obtainable, and second the fitting of γ and subsequent use of Equation 3 as a model is computationally straightforward.
0005H.L. Bertoni, W. Honcharenko, L.R Maciel, and H.H. Xia, "UHF Propagation Prediction for Wireless Personal Communications," Proc. of the IEEE, vol. 82, pp. 1333-1359 (Sep. 1994); and J.B. Anderson, T.S. Rappaport, and S. Yoshida, "Propagation Measurements and Models for Wireless Communications Channels," IEEE Communic. Mag., pp. 42-49 (Jan. 1995) both provide an excellent overview of wireless propagation models in the prior art.
0006Another technique in the prior art for measuring and modeling indoor RF propagation is adapted from ray-tracing techniques, as taught by S.J. Fortune et al., U.S. Patent 5,450,615, issued September 12, 1995. According to this technique the RF propagation within a building is predicted by modeling an RF signal as a plurality of rays that pass through objects in the building or are reflected off of objects in the building or both. An advantage of this technique is that the resulting wireless propagation model can be very effective for modeling the RF propagation characteristics of the building of interest. There are, however, two disadvantages of this technique. The first disadvantage is that the technique requires an extremely detailed floorplan of the building of interest including the major RF obstacles in the building. For a typical office building this can be burdensome. The second disadvantage is that the technique is computationally intense, and, therefore, typically requires fast expensive computers.
0007Therefore, here exists the need for a wireless propagation model that is effective, that does not require a great deal of empirical data about the building to be gathered and that can be reasonably implemented on a typical desktop computer.
<u>Summary of the Invention</u>
0008Some embodiments of the present invention are capable of modeling the propagation of RF signals in an indoor environment without the restrictions and disadvantages of techniques in the prior art. In particular, some embodiments are extremely effective, yet require only a modest amount of architectural information and computational power.
0009These advantages may be found in some embodiments of the present invention that comprise six distinct phases. In one phase, the mean wall separation, <maths id="math0004"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0004.tif" /></maths><i> ,</i> of one floor of a building is estimated. In a second phase, a reflection coefficient, <i>s</i>, is estimated for the floor in general. In a third phase, a number of trial RF propagation measurements are made to gather empirical data about the RF propagation characteristics of the floor. In the fourth phase, the mean wall separation, <maths id="math0005"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0005.tif" /></maths><i> ,</i> the reflection coefficient s, and the trial RF propagation measurements are fit, using well-known techniques, into a wireless propagation model, such as<maths id="math0006"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">r</mtext><mtext>) = </mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext mathvariant="italic">0</mtext></mrow></msub></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">πD</mtext></mrow></mfrac><msqrt><mfrac><mrow><mtext>πξ</mtext></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">r</mtext></mrow></mfrac></msqrt><msup><mrow><mtext mathvariant="italic">e</mtext></mrow><mrow><mtext>-r/ξ</mtext></mrow></msup><mtext>,</mtext></mrow></math><img file="EP0903880A2_D0006.tif" /></maths> where <i>P</i>(<i>r</i>) is the measured or predicted power at a distance, <i>r</i>, from the transmitter,<maths id="math0007"><img file="EP0903880A2_D0007.tif" /></maths> and ξ is the parameter that is fit to the empirical data. In the fifth phase the wireless propagation model is used to predict the RF signal strength throughout the floor from a base station at a given location, and in the six phase one or more base stations are installed in the building based on the results predicted by the wireless propagation model. The fifth and sixth phases can be thereafter repeated for other building with sufficiently similar characteristics.
<u>Brief Description of the Drawings</u>
0010FIG. 1 depicts a flowchart of the illustrative embodiment of the present invention.
0011FIG. 2 depicts a outline of one floor of atypical office building wherein the location of a trial transmitter is represented by a triangle and each location at which an empirical signal strength measurement is made is represented by a square.
0012FIG. 3 depicts a map of the wireless coverage for the floor of FIG. 2 from a single base station according to one wireless propagation model.
0013FIG. 4 depicts a map of the wireless coverage for the floor of FIG. 2 from a single base station according to a second wireless propagation model.
0014FIG. 5 depicts a map of the wireless coverage for the floor of FIG. 2 from two base stations according to the first wireless propagation model.
0015FIG. 6 depicts a map of the wireless coverage for the floor of FIG. 2 from two base stations according to the second wireless propagation model.
0016FIG. 7 depicts an outline of the floor in FIG. 2 and the location of two base stations that have been placed on that floor.
<u>Detailed Description</u>
0017FIG. 1 depicts a flowchart of the illustrative embodiment of the present invention, which can employ one of the wireless propagation models described below to assist in requires the number and placement of wireless base stations in a building. Each of the wireless propagation models described below belong to one of three classes. The three classes are: (1) the Bulk Diffusion Models, (2) the Boundary Diffusion Models, and (3) the Statistical Boltzmann Models.
0018The remainder of the detailed description teaches how each of the wireless propagation models is applied to one floor of an illustrative multi-story office building to predict that the signal strength from one or more base stations on the floor is at least -65 dBm everywhere on the floor. Although it will be clear to those skilled in the art how to apply each wireless propagation model to an entire multi-story building as a whole, the building materials between adjacent floors of a multi-story building typically attenuate RF signals to the extent that there is little coverage provided by a base station on one floor to a wireless terminal on another floor. Therefore, it is advantageous to apply the illustrative embodiment of the present invention to a multi-story building one floor at a time.
0019FIG. 2 depicts an outline of the outer walls of one floor, floor 200, of the illustrative office building. As defined by the outer walls, floor 200 is rectangular and 100 meters long in the x-direction and 50 meters long in the y-direction. Although the outer wall of the building is rectangular, it will be clear to those skilled in the art how to make and use embodiments of the present invention for any size or shape building.
0020Referring to step 101 in FIG. 1, the mean wall separation of the floor of interest is estimated. Depending on which wireless propagation model is to be used, either an isotropic estimate or an anisotropic estimate is made.
0021The Bulk Diffusion Models use an isotropic estimate of the mean wall separation, <maths id="math0008"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0008.tif" /></maths><i> ,</i> and the Boundary Diffusion Models and the Statistical Boltzmann Models advantageously use anisotropic estimates, <maths id="math0009"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0009.tif" /></maths><sub><i>x</i></sub><i>,</i> and <maths id="math0010"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0010.tif" /></maths><sub><i>y</i></sub><i>,</i> each in orthogonal directions.
0022In accordance with the illustrative embodiment, a detailed floorplan of floor 200 is not needed to estimate the mean wall separation. Instead, the mean wall separation can be estimated by measuring the dimensions of a few typical rooms in floor 200, in well-known fashion. Illustratively, the isotropic mean wall separation, <maths id="math0011"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0011.tif" /></maths> , of floor 200 is 4.5 meters, and the anisotropic mean wall separation in the x-direction, <maths id="math0012"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0012.tif" /></maths><sub><i>X</i></sub><i>,</i> of floor 200 is 4 meters and the mean wall separation in the y-direction, <maths id="math0013"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0013.tif" /></maths><sub><i>y</i></sub><i>,</i> is 5 meters.
0023At step 103 in FIG. 1, one or more reflection coefficients are estimated for the floor of interest. Each reflection coefficient indicates what percentage of the energy of an RF signal, at the frequency of interest, is typically reflected off of a wall in comparison to the amount of energy that is transmitted through the wall.
0024The Bulk Diffusion Models require a single isotropic estimate of reflection, <i>s</i>, which is based on the type of building materials used in the construction of both the interior and exterior walls of floor 200. It will be clear to those skilled in the art how to estimate the isotropic reflection coefficient, <i>s</i>, for a given set of building materials. Illustratively, the interior walls of floor 200 are constructed of sheetrock and the exterior walls of floor 200 are glass. Therefore, a reasonable estimate of reflection, <i>s</i>, is 0.4, as is well-known in the art
0025The Boundary Diffusion Models and the Statistical Boltzmann Models advantageously require an anisotropic estimate of reflection, <i>s(θ)</i> , of the interior walls and an anisotropic estimate of reflection, <i>B(θ)</i>, of the exterior walls. The anisotropic nature of these estimates consider that the amount of energy in an RF signal that is reflected from a wall depends on the angle between the direction of propagation and the direction normal to the wall at the point of contact. Separate estimates of reflection are advantageously made for both the interior walls and for the exterior walls so that the typical disparity in their construction can be considered. It will be clear to those skilled in the art how to estimate values for <i>s(θ)</i> and <i>B(θ)</i> for given building materials. When, as in this example, the interior walls are mostly sheetrock and the exterior walls are glass, <i>s(θ)</i> can be reasonably estimated as sin<sup>2</sup> θ and <i>B(θ)</i> reasonably estimated as 0.4 + 0.6sin<sup>2</sup> θ, where θ is the angle between the direction of propagation and the direction normal to the wall at the point of contact.
0026At step 105, empirical measurements of RF propagation throughout floor 200 are made which are advantageously used to help create the wireless propagation model of floor 200. Advantageously, a trial transmitter is placed at a location, <i>x</i><sub><i>0,</i></sub><i> y</i><sub><i>0,</i></sub> on floor 200 that radiates a trial signal at power, <i>P</i><sub><i>0</i></sub><i> ,</i> at the frequency of interest. In the illustrative embodiment, the power, <i>P</i><sub><i>0</i></sub><i> ,</i> of the trial transmitter is +20 dBm. Advantageously, the location, <i>x</i><sub><i>0</i></sub><sub>,</sub><i>y</i><sub><i>0</i></sub><i>,</i> is asymmetric with respect to the exterior walls of the building. In FIG. 2, trial transmitter 201 is shown at location <i>x</i><sub><i>0</i></sub>=30, <i>y</i><sub><i>0</i></sub>=20.
0027Advantageously, the strength of the trial signal, <i>P</i><sub><i>i</i></sub>, is measured with a signal strength detector at <i>n</i> locations, <i>x</i><sub><i>i</i></sub><i>, y</i><sub><i>i</i></sub>, for <i>i</i> = 1 to <i>n</i>. Advantageously, the <i>n</i> locations are asymmetrically located around the trial transmitter and are fairly distributed throughout floor 200.
0028When reasonable care has been taken to estimate the mean wall separation and reflection coefficient(s) for floor 200, the trial signal, <i>P</i><sub><i>i</i></sub><i>,</i> need be only measured at a few locations. In contrast, when reasonable care has not been taken in estimating the mean wall separation and reflection coefficient(s) for floor 200, the trial signal, <i>P</i><sub><i>i</i></sub><i>,</i> should be measured at a larger number of locations to compensate. Furthermore, if no estimate of the mean wall separation or the reflection coefficient(s) is made for the building of interest, and estimates of a typical building are used instead then measuring the trial signal, <i>P</i><sub><i>i</i></sub><i>,</i> at a large number of locations may provide enough empirical data to enable the creation of a satisfactory wireless propagation model for the floor of interest. It will be clear to those skilled in the art how to set up the trial transmitter and how to measure the strength of the trial signal at various locations on floor 200.
0029For the illustrative floor 200, Table 1 lists the <i>n</i>=6 locations at which the trial signal, <i>P</i><sub><i>i</i></sub>, is measured and the resulting measurements. <tables id="tabl0001" num="0001"><table frame="all"><title>Table 1 -</title><tgroup cols="4" colsep="1" rowsep="1"><colspec colnum="1" colname="col1" colwidth="39.37mm" /><colspec colnum="2" colname="col2" colwidth="39.37mm" /><colspec colnum="3" colname="col3" colwidth="39.37mm" /><colspec colnum="4" colname="col4" colwidth="39.37mm" /><thead valign="top"><row rowsep="1"><entry namest="col1" nameend="col4" align="center">Empirical Signal Strength Measurements</entry></row><row><entry namest="col1" nameend="col1" align="center"><i>i</i></entry><entry namest="col2" nameend="col2" align="center"><i>x</i><sub><i>i</i></sub></entry><entry namest="col3" nameend="col3" align="center"><i>y</i><sub><i>i</i></sub></entry><entry namest="col4" nameend="col4" align="center"><i>P</i><sub><i>i</i></sub></entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="center">1</entry><entry namest="col2" nameend="col2" align="center">20</entry><entry namest="col3" nameend="col3" align="center">40</entry><entry namest="col4" nameend="col4" align="center">-50dBm</entry></row><row><entry namest="col1" nameend="col1" align="center">2</entry><entry namest="col2" nameend="col2" align="center">40</entry><entry namest="col3" nameend="col3" align="center">40</entry><entry namest="col4" nameend="col4" align="center">-55 dBm</entry></row><row><entry namest="col1" nameend="col1" align="center">3</entry><entry namest="col2" nameend="col2" align="center">60</entry><entry namest="col3" nameend="col3" align="center">40</entry><entry namest="col4" nameend="col4" align="center">-65 dBm</entry></row><row><entry namest="col1" nameend="col1" align="center">4</entry><entry namest="col2" nameend="col2" align="center">80</entry><entry namest="col3" nameend="col3" align="center">40</entry><entry namest="col4" nameend="col4" align="center">-80 dBm</entry></row><row><entry namest="col1" nameend="col1" align="center">5</entry><entry namest="col2" nameend="col2" align="center">60</entry><entry namest="col3" nameend="col3" align="center">20</entry><entry namest="col4" nameend="col4" align="center">-55dBm</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="center">6</entry><entry namest="col2" nameend="col2" align="center">80</entry><entry namest="col3" nameend="col3" align="center">20</entry><entry namest="col4" nameend="col4" align="center">-65 dBm</entry></row></tbody></tgroup></table></tables>
0030FIG. 2 shows the location of trial transmitter 201 and the respective locations at which the trial signal was measured are depicted by squares in floor 200.
0031At step 109, the mean wall separation, the reflection coefficient(s) and the empirical data from step 107 are fit into one or more of the wireless propagation models.
0032<b><u>Bulk Diffusion Models</u> -</b> Two specific Bulk Diffusion Models are presented, both of which are isotropic. The first Bulk Diffusion Model is based on:<maths id="math0014" num="(Eq. 4)"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">r</mtext><mtext>) = </mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext mathvariant="italic">0</mtext></mrow></msub></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">πD</mtext></mrow></mfrac><msqrt><mfrac><mrow><mtext>πξ</mtext></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">r</mtext></mrow></mfrac></msqrt><msup><mrow><mtext mathvariant="italic">e</mtext></mrow><mrow><mtext>-r/ξ</mtext></mrow></msup></mrow></math><img file="EP0903880A2_D0014.tif" /></maths> wherein <i>P(r)</i> is the measured or predicted power at a distance, <i>r,</i> from the transmitter,<maths id="math0015"><img file="EP0903880A2_D0015.tif" /></maths> and ξ is the parameter that is fit to the data in Table 1, in well-known fashion. In other words, Equation 4 is a least-squares regression equation <i>of P(r)</i> on <i>r</i> and the mean wall separation, <maths id="math0016"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d ,</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0016.tif" /></maths> the reflection coefficient <maths id="math0017"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">s</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0017.tif" /></maths>, and the empirical data in Table 1 are used to find the best fit value for ξ. For the values <i>P</i><sub><i>0</i></sub> = +20 dBm, <maths id="math0018"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0018.tif" /></maths> =4.5, <i>s</i> = 0.4 and the values in Table 1, the best fit value for ξ is 24. Alternatively, when the characteristic decay lenght, ξ is large <i>D</i> = <maths id="math0019"><math display="inline"><mrow><mfrac><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">s</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0019.tif" /></maths> can be used.
0033The second Bulk Diffusion Model is based on:<maths id="math0020"><img file="EP0903880A2_D0020.tif" /></maths> wherein <i>P(r)</i> is the measured or predicted power at a distance, <i>r</i>, from the transmitter,<maths id="math0021"><img file="EP0903880A2_D0021.tif" /></maths><i>K</i><sub><i>0</i></sub> is a modified Bessel function and ξ is the parameter that is fit to the data in Table 1, in well-known fashion. Alternatively, when the characteristic decay length, ξ, is large D = <maths id="math0022"><math display="inline"><mrow><mfrac><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">2s</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0022.tif" /></maths> can be used.
0034<b><u>Boundary Diffusion Models</u></b>-Two specific Boundary Diffusion Models are presented, both of which are anisotropic. In general, the Boundary Diffusion Models are expected to be more effective than the Bulk Diffusion Models. The first Boundary Diffusion Model is based on the solution to the partial differential equation:<maths id="math0023" num="(Eq. 6)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mtext> (θ)</mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">∂</mtext></mrow><mrow><mtext mathvariant="italic">2</mtext></mrow></msup><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,θ</mtext><mtext>)</mtext></mrow><mrow><msup><mrow><mtext mathvariant="italic">∂X</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mtext>(</mtext><mtext mathvariant="italic">θ</mtext><mtext>) </mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">∂</mtext></mrow><mrow><mtext mathvariant="italic">2</mtext></mrow></msup><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,θ</mtext><mtext>)</mtext></mrow><mrow><msup><mrow><mtext mathvariant="italic">∂y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> - </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,θ</mtext><mtext>)</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext mathvariant="italic">δ</mtext><msub><mrow><mtext>(x-x</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><msub><mrow><mtext>)δ(y-y</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>) = 0</mtext></mrow></math><img file="EP0903880A2_D0023.tif" /></maths> wherein <i>x</i><sub><i>0</i></sub><sub>,</sub><i>y</i><sub><i>0</i></sub> is the location of the transmitter, <i>P(x,y,</i>θ) is the measured or predicted power at location <i>x</i>, <i>y,</i> propagating in the direction θ with respect to the x-axis, δ is the Dirac delta function, <chemistry id="chem0001" num="0001"><img file="EP0903880A2_D0024.tif" /></chemistry> and ξ<sub>1</sub> is the parameter to be fit to the data in Table 1, in well-known fashion; with the boundary condition<maths id="math0024"><math display="block"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">n</mtext></mrow><mo>ˆ</mo></mover><mtext mathvariant="italic"> ·</mtext><mtext></mtext><mover accent="true"><mrow><mtext>∇</mtext></mrow><mo>→</mo></mover><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,θ</mtext><mtext> ) </mtext><mtext mathvariant="italic">= -</mtext><mtext></mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">K</mtext></mrow><mrow><mtext mathvariant="italic">ext</mtext></mrow></msub></mrow><mrow><mfenced open="|" close="|"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">n</mtext></mrow><mo>ˆ</mo></mover><mtext> · </mtext><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">l</mtext></mrow><mo>→</mo></mover></mrow><mrow><mtext mathvariant="italic">wall</mtext></mrow></msub><mtext> (</mtext><mtext mathvariant="italic">θ</mtext><mtext>)</mtext></mrow></mfenced></mrow></mfrac><mtext></mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,θ</mtext><mtext> ), wherein </mtext><msub><mrow><mtext mathvariant="italic">K</mtext></mrow><mrow><mtext mathvariant="italic">ext</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>2(1-</mtext><mtext mathvariant="italic">B</mtext><mtext>(θ))</mtext></mrow><mrow><mtext>1+</mtext><mtext mathvariant="italic">B</mtext><mtext>(θ)(1-2</mtext><mtext mathvariant="italic">s</mtext><mtext>(θ))</mtext></mrow></mfrac><mtext>,</mtext></mrow></math><img file="EP0903880A2_D0025.tif" /></maths><maths id="math0025"><img file="EP0903880A2_D0026.tif" /></maths> and n̂ is the unit vector normal to the boundary at the boundary point x, y. For the values, <maths id="math0026"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0027.tif" /></maths><sub><i>X</i></sub>= 4 meters, <maths id="math0027"><math display="inline"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow></math><img file="EP0903880A2_D0028.tif" /></maths><sub><i>y</i></sub><i>=</i> 5 meters, <i>s</i>(θ) = sin<sup>2</sup> θ, <i>B</i>(θ) = 0.4 + 0.6sin<sup>2</sup><i>θ</i> and the values in Table 1, the best fit value for ξ<sub>1</sub> is 16.
0035Alternatively, when the characteristic decay length ξ<sub>1</sub> is large,<maths id="math0028"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mtext>(</mtext><mtext mathvariant="italic">θ</mtext><mtext>) = </mtext><mfrac><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mfenced open="|" close="|"><mrow><mtext mathvariant="italic">cosθ</mtext></mrow></mfenced></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">s</mtext><mtext>(θ)</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0029.tif" /></maths> and<maths id="math0029"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mtext>(</mtext><mtext mathvariant="italic">θ</mtext><mtext>) = </mtext><mfrac><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mfenced open="|" close="|"><mrow><mtext>sin θ</mtext></mrow></mfenced></mrow><mrow><mtext mathvariant="italic">2s</mtext><mtext>(</mtext><mtext mathvariant="italic">π</mtext><mtext>/2 - </mtext><mtext mathvariant="italic">θ</mtext><mtext>)</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0030.tif" /></maths> can be used.
0036The second Boundary Diffusion Model is based on the solution to the partial differential equation:<maths id="math0030" num="(Eq. 7)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mtext>(θ) </mtext><mfrac><mrow><msup><mrow><mtext>∂</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ)</mtext></mrow><mrow><mtext>∂</mtext><msup><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mtext>(θ) </mtext><mfrac><mrow><msup><mrow><mtext>∂</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ)</mtext></mrow><mrow><mtext>∂</mtext><msup><mrow><mtext mathvariant="italic">y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></mfrac><mtext> - </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ)</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><mspace linebreak="newline" /><mtext>- </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ) - 〈</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ)〉</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub></mrow></mfrac><mtext> + </mtext><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>δ(</mtext><mtext mathvariant="italic">x</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>)δ(</mtext><mtext mathvariant="italic">y</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">y</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>) = 0</mtext></mrow></math><img file="EP0903880A2_D0031.tif" /></maths> wherein <i>x</i><sub><i>0</i></sub>, <i>y</i><sub><i>0</i></sub> is the location of the transmitter, <i>P</i>(<i>x, y, θ</i> ) is the measured or predicted power at location x, <i>y,</i> propagating in the direction θ with respect to the x-axis, δ is the Dirac delta function,<maths id="math0031"><img file="EP0903880A2_D0032.tif" /></maths><maths id="math0032"><img file="EP0903880A2_D0033.tif" /></maths> and ξ<sub>1</sub> and ξ<sub>2</sub> are parameters to be fit to the data in Table 1, in well-known fashion; with the boundary condition<maths id="math0033"><math display="block"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">n</mtext></mrow><mo>ˆ</mo></mover><mtext mathvariant="italic"> ·</mtext><mtext></mtext><mover accent="true"><mrow><mtext>∇</mtext></mrow><mo>→</mo></mover><mtext></mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>,</mtext><mtext mathvariant="italic">y</mtext><mtext>,θ) = - </mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">K</mtext></mrow><mrow><mtext mathvariant="italic">ext</mtext></mrow></msub></mrow><mrow><mfenced open="|" close="|"><mrow><mover accent="true"><mrow><mtext mathvariant="italic">n</mtext></mrow><mo>ˆ</mo></mover><mtext> · </mtext><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">l</mtext></mrow><mo>→</mo></mover></mrow><mrow><mtext mathvariant="italic">wall</mtext></mrow></msub><mtext> (θ)</mtext></mrow></mfenced></mrow></mfrac><mtext></mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ), wherein </mtext><msub><mrow><mtext mathvariant="italic">K</mtext></mrow><mrow><mtext mathvariant="italic">ext</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext>2(1-B(θ))</mtext></mrow><mrow><mtext>1+B(θ)(1-2s(θ))</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0034.tif" /></maths> and<maths id="math0034"><img file="EP0903880A2_D0035.tif" /></maths>
0037Alternatively, when the characteristic decay length ξ<sub>1</sub> and ξ<sub>2</sub> are large,<maths id="math0035"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mtext> (θ) = </mtext><mfrac><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub><mfenced open="|" close="|"><mrow><mtext>cos θ</mtext></mrow></mfenced></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">s</mtext><mtext>(θ)</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0036.tif" /></maths> and<maths id="math0036"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">D</mtext></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mtext>(θ) = </mtext><mfrac><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub><mfenced open="|" close="|"><mrow><mtext>sin θ</mtext></mrow></mfenced></mrow><mrow><mtext>2</mtext><mtext mathvariant="italic">s</mtext><mtext>(π/2 - θ)</mtext></mrow></mfrac></mrow></math><img file="EP0903880A2_D0037.tif" /></maths> can be used.
0038The solution to equation 6 or 7 yields the anisotropic function <i>P(x, y, θ)</i> which predicts the power at location <i>x, y,</i> propagating in the direction θ with respect to the x-axis. Each final wireless propagation model, <i>P(x, y),</i> is advantageously the sum of all the power traveling in all directions at location <i>x, y.</i> Therefore, the anisotropic function <i>P(x, y,θ)</i> is related to the total power by:<maths id="math0037"><img file="EP0903880A2_D0038.tif" /></maths>
0039<b><u>Statistical Boltzmann Models</u></b>-Two specific Statistical Boltzmann Models are presented, both of which are anisotropic. In general, the Statistical Boltzmann Models are expected to be more effective than either the Bulk Diffusion Models or the Boundary Diffusion Models, especially for buildings that have an irregular shape. The first Statistical Boltzmann Model is based on the solution to the partial differential equation:<maths id="math0038" num="(Eq. 8)"><math display="block"><mrow><mover accent="true"><mrow><mtext>θ</mtext></mrow><mo>ˆ</mo></mover><mtext> ·</mtext><mover accent="true"><mrow><mtext>∇</mtext></mrow><mo>→</mo></mover><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ) + </mtext><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>δ(</mtext><mtext mathvariant="italic">x</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>)δ(</mtext><mtext mathvariant="italic">y</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">y</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>) + </mtext><mfrac><mrow><mtext mathvariant="italic">s</mtext><mtext>(θ)</mtext><mfenced open="|" close="|"><mrow><mtext>cosθ</mtext></mrow></mfenced></mrow><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">x</mtext></mrow></msub></mrow></mfrac><mtext> [</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>,</mtext><mtext mathvariant="italic">y</mtext><mtext>,θ) - </mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>,</mtext><mtext mathvariant="italic">y</mtext><mtext>,π-θ)] + </mtext><mfrac><mrow><mtext mathvariant="italic">s</mtext><mtext>(π/2-θ)</mtext><mfenced open="|" close="|"><mrow><mtext>sinθ</mtext></mrow></mfenced></mrow><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub></mrow></mfrac><mtext>[</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ) - </mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,-θ)] + </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ)</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><mtext> = 0</mtext></mrow></math><img file="EP0903880A2_D0039.tif" /></maths> wherein <i>x</i><sub><i>0</i></sub><i>, y</i><sub><i>0</i></sub> is the location of the transmitter, <i>P(x, y, θ)</i> is the measured or predicted power at location <i>x, y</i>, propagating in the direction θ with respect to the x-axis, δ is the Dirac delta function, θ̂ is a unit vector in the direction θ and ξ<sub>1</sub> is parameter to be fit to the data in Table 1, in well-known fashion; with the boundary condition: <i>P</i>(<i>x, y</i>, θ) = <i>B</i>(θ)<i>P</i>(<i>x,y</i>,ϕ), where ϕ is angle of the vector θ̂ - 2<i>n̂</i>(<i>n̂</i> · θ̂), wherein <i>n̂</i> is a unit vector normal to the boundary at the boundary point <i>x</i>, <i>y.</i>
0040The second Statistical Boltzmann Model is based on the solution to the partial differential equation:<maths id="math0039" num="(Eq. 9)"><math display="block"><mrow><mtext>θ ·</mtext><mover accent="true"><mrow><mtext>∇</mtext></mrow><mo>→</mo></mover><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ) + </mtext><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>δ(</mtext><mtext mathvariant="italic">x</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>)δ(</mtext><mtext mathvariant="italic">y</mtext><mtext>-</mtext><msub><mrow><mtext mathvariant="italic">y</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>) + </mtext><mfrac><mrow><mtext mathvariant="italic">s</mtext><mtext>(θ)</mtext><mfenced open="|" close="|"><mrow><mtext>cosθ</mtext></mrow></mfenced></mrow><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext>x</mtext></mrow></msub></mrow></mfrac><mtext> [</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ) - </mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>π-θ)]</mtext><mspace linebreak="newline" /><mtext>+ </mtext><mfrac><mrow><mtext mathvariant="italic">s</mtext><mtext>(π/2-θ)</mtext><mfenced open="|" close="|"><mrow><mtext>sinθ</mtext></mrow></mfenced></mrow><mrow><msub><mrow><mover accent="true"><mrow><mtext mathvariant="italic">d</mtext></mrow><mo>¯</mo></mover></mrow><mrow><mtext mathvariant="italic">y</mtext></mrow></msub></mrow></mfrac><mtext>[</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ) - </mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>-θ)] + </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ)</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>1</mtext></mrow></msub></mrow></mfrac><mtext> + </mtext><mfrac><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">x,y,</mtext><mtext>θ) - 〈P(</mtext><mtext mathvariant="italic">x,y</mtext><mtext>,θ)〉</mtext></mrow><mrow><msub><mrow><mtext>ξ</mtext></mrow><mrow><mtext>2</mtext></mrow></msub></mrow></mfrac><mtext>= 0</mtext></mrow></math><img file="EP0903880A2_D0040.tif" /></maths> wherein <i>x</i><sub><i>0</i></sub>, <i>y</i><sub><i>0</i></sub> is the location of the transmitter, <i>P(x, y,</i> θ) is the measured or predicted power at location x, <i>y,</i> propagating in the direction θ with respect to the x-axis, δ is the Dirac delta function,θ̂ is a unit vector in the direction θ,<maths id="math0040"><img file="EP0903880A2_D0041.tif" /></maths> and ξ<sub>1</sub>, and ξ<sub>2</sub> are parameters to be fit to the data in Table 1, in well-known fashion; with the boundary condition: <i>P(x,y,θ)</i> = <i>B(θ)P(x,y,ϕ)</i><sub><i>,</i></sub> where ϕ is angle of the vector <i>θ̂ - 2n̂(n̂·θ̂)</i>, wherein <i>n̂</i> is a unit vector normal to the boundary at the boundary point <i>x,y.</i>
0041Each solution to equation 8 or 9 yields the anisotropic function <i>P(x,y,θ</i> ) which predicts the power at location <i>x</i>, <i>y,</i> propagating in the direction θ with respect to the x-axis. The final wireless propagation model, <i>P(x, y),</i> is advantageously the sum of all the power traveling in all directions at location x, <i>y.</i> Therefore, the anisotropic function <i>P(x, y,</i> θ) is related to the total power by:<maths id="math0041"><img file="EP0903880A2_D0042.tif" /></maths>
0042When step 107 is complete, a wireless propagation model, <i>P(x, y),</i> for floor 200 is complete.
0043At step 109, a guess is made as to the potential location for placing a wireless base station on floor 200. For economic reasons, it is advantageous if all of floor 200 can be covered with a single base station. If all of floor 200 can be covered with a single base station, it is reasonable to guess that the best location for that base station would be the exact center of floor 200 <i>(x</i><sub><i>0</i></sub> = 50, <i>y</i><sub><i>0</i></sub><i> = 25).</i>
0044At step III, the wireless propagation model <i>P(x, y),</i> is used to predict the mean power throughout floor 200 with the base station at location <i>x</i><sub><i>0,</i></sub><i> y</i><sub><i>0</i></sub><i>,</i> in well-known fashion.
0045FIG. 3 depicts a topographical map of floor 200 that predicts the signal strength (in dBm) throughout floor 200 based on a single base station at <i>x</i><sub><i>0</i></sub> = 50, <i>y</i><sub><i>0</i></sub> = 25 and based on the first Bulk Diffusion Model of Equation 4 :<maths id="math0042"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">r</mtext><mtext>) = </mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">P</mtext></mrow><mrow><mtext>0</mtext></mrow></msub></mrow><mrow><mtext>5.76</mtext></mrow></mfrac><msqrt><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext mathvariant="italic">r</mtext></mrow></mfrac></msqrt><msup><mrow><mtext mathvariant="italic">e</mtext></mrow><mrow><mtext>-</mtext><mtext mathvariant="italic">r</mtext><mtext>/24</mtext></mrow></msup><mtext>,</mtext></mrow></math><img file="EP0903880A2_D0043.tif" /></maths> where<maths id="math0043"><math display="block"><mrow><mtext mathvariant="italic">r</mtext><mtext> = </mtext><msqrt><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext> - </mtext><mtext mathvariant="italic">x</mtext><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext> + (</mtext><msub><mrow><mtext mathvariant="italic">y</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext> - </mtext><mtext mathvariant="italic">y</mtext><msup><mrow><mtext>)</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></msqrt><mtext>.</mtext></mrow></math><img file="EP0903880A2_D0044.tif" /></maths> As discussed above, adequate coverage is achieved in the illustrative example when the signal strength in all of floor 200 is above -65 dBm. As can be seen in FIG. 3, the first Bulk Diffusion Model predicts that there will be gaps in coverage at the ends of floor 200 if only one base station is used.
0046For comparison, FIG. 4 depicts a topological map of floor 200 that predicts the signal strength (in dBm) throughout floor 200 based on a single base station at <i>x</i><sub><i>0</i></sub> = 50, <i>y</i><sub><i>0</i></sub> = 25 and based on the first Boundary Diffusion Model of Equation 6. As can be seen in FIG. 4, the first Boundary Diffusion Model predicts substantially different isopower lines than does the first Bulk Diffusion Model, but both predict that there will be gaps in coverage at the ends of floor 200 if only one base station is used.
0047At step 113, the question is posed whether the wireless propagation model predicts regions with insufficient coverage. Although FIG. 3 and FIG. 4 predict different signal strengths, both predict that a single base station is insufficient to cover floor 200, and, therefore, control passes back to step 109.
0048At step 109, another guess is made as to the possible locations for placing base stations on floor 200. For economical reasons, it is advantageously if floor 200 can be covered with two base stations. If two base stations are to cover floor 200, a reasonable guess is to locate each at the center of one half of floor 200 (i.e., one base station at <i>x</i><sub><i>0</i></sub> = 25, <i>y</i><sub><i>0</i></sub> = 25 and the second base station at <i>x</i><sub><i>0</i></sub> = 75, <i>y</i><sub><i>0</i></sub> = 25).
0049At step 111, the wireless propagation model, <i>P(x, y),</i> is used to predict the mean power throughout floor 200 with the base stations at <i>x</i><sub><i>0</i></sub> = 25, <i>y</i><sub><i>0</i></sub> = 25 <i>and x</i><sub><i>0</i></sub><i> =</i> 75, <i>y</i><sub><i>0</i></sub> = 25, where the total power at each location on floor 200 is equal to the sum of the power from each base station. FIG. 5 depicts a map that predicts the signal strength throughout floor 200 based on the first Bulk Diffusion Model and FIG. 6 depicts a similar map based on the first Boundary Diffusion Model.
0050At step 113, the question is posed whether the wireless propagation model predicts regions with insufficient coverage. Although FIG. 5 and FIG. 6 predict substantially different isopower lines, both models predict that two base stations (one at <i>x</i><sub><i>0</i></sub> = 25, <i>y</i><sub><i>0</i></sub> = 25 and one at <i>x</i><sub><i>0</i></sub> = 75, <i>y</i><sub><i>0</i></sub> = 25) provide sufficient coverage for floor 200.
0051At step 115, two wireless base stations are placed at <i>x</i><sub><i>0</i></sub> = 25, <i>y</i><sub><i>0</i></sub> = 25 and <i>x</i><sub><i>0</i></sub> = 75, <i>y</i><sub><i>0</i></sub> = 25 in floor 200. FIG. 7 depicts floor 200 with base station 701 placed at <i>x</i><sub><i>0</i></sub> = 25, <i>y</i><sub><i>0</i></sub> = 25 and base station <i>702 at x</i><sub><i>0</i></sub><i> = 75,y</i><sub><i>0</i></sub><i> = 25.</i>
0052It is to be understood that the above-described embodiments are merely illustrative of the present invention and that many variations may be devised by those skilled in the art without departing from the scope of the invention. It is therefore intended that such variations be included within the scope of the following claims.
111 sheets
Sheet 1 Sheet 2 Sheet 3 Sheet 4 Sheet 5 Sheet 6 Sheet 7 Sheet 8 Sheet 9 Sheet 10 Sheet 11 Sheet 12 Sheet 13 Sheet 14 Sheet 15 Sheet 16 Sheet 17 Sheet 18 Sheet 19 Sheet 20 Sheet 21 Sheet 22 Sheet 23 Sheet 24 Sheet 25 Sheet 26 Sheet 27 Sheet 28 Sheet 29 Sheet 30 Sheet 31 Sheet 32 Sheet 33 Sheet 34 Sheet 35 Sheet 36 Sheet 37 Sheet 38 Sheet 39 Sheet 40 Sheet 41 Sheet 42 Sheet 43 Sheet 44 Sheet 45 Sheet 46 Sheet 47 Sheet 48 Sheet 49 Sheet 50 Sheet 51 Sheet 52 Sheet 53 Sheet 54 Sheet 55 Sheet 56 Sheet 57 Sheet 58 Sheet 59 Sheet 60 Sheet 61 Sheet 62 Sheet 63 Sheet 64 Sheet 65 Sheet 66 Sheet 67 Sheet 68 Sheet 69 Sheet 70 Sheet 71 Sheet 72 Sheet 73 Sheet 74 Sheet 75 Sheet 76 Sheet 77 Sheet 78 Sheet 79 Sheet 80 Sheet 81 Sheet 82 Sheet 83 Sheet 84 Sheet 85 Sheet 86 Sheet 87 Sheet 88 Sheet 89 Sheet 90 Sheet 91 Sheet 92 Sheet 93 Sheet 94 Sheet 95 Sheet 96 Sheet 97 Sheet 98 Sheet 99 Sheet 100 Sheet 101 Sheet 102 Sheet 103 Sheet 104 Sheet 105 Sheet 106 Sheet 107 Sheet 108 Sheet 109 Sheet 110 Sheet 111
Every citation, both ways
| Document | Relation | Office | Cited during |
|---|---|---|---|
| US8233906B2 | Cited by | United States of America | Applicant |
| EP1181574A1 | Cited by | European Patent Office (EPO) | Search report |
| WO03019826A1 | Cited by | World Intellectual Property Organization (WIPO) | International search |
| WO03019827A1 | Cited by | World Intellectual Property Organization (WIPO) | International search |
| EP1869578A2 | Cited by | European Patent Office (EPO) | Search report |
| US7796983B2 | Cited by | United States of America | Applicant |
| US7583961B2 | Cited by | United States of America | Applicant |
| WO2007023074A1 | Cited by | World Intellectual Property Organization (WIPO) | International search |
| EP1869578A4 | Cited by | European Patent Office (EPO) | Search report |
| EP1592273A3 | Cited by | European Patent Office (EPO) | Search report |
| EP1592273A2 | Cited by | European Patent Office (EPO) | Search report |
| WO2006116766A1 | Cited by | World Intellectual Property Organization (WIPO) | International search |
| EP1181574A4 | Cited by | European Patent Office (EPO) | Search report |
| US5450615A | Cites | United States of America | Search report |
| US5491644A | Cites | United States of America | Search report |
5 members in 3 offices; this record represents the family
Priority claims2
| Document | Office | Kind | Date |
|---|---|---|---|
| 933493 | United States of America | – | |
| 93349397 | United States of America | A |
Members5
| Document | Office | Kind | |
|---|---|---|---|
| EP0903880A2This record | European Patent Office (EPO) | A2 | |
| JPH11164354A | Japan | A | |
| US6119009A | United States of America | A | |
| EP0903880A3 | European Patent Office (EPO) | A3 | |
| US6678525B1 | United States of America | B1 |
12 legal events, as the office reported them to INPADOC
Over the term
Point at a mark for the eventEvents
| Event | Code | |
|---|---|---|
| Application withdrawnWithdrawn18W | 18W | |
| Information on the status of an ep patent application or granted ep patentGrantedSTATUS: THE APPLICATION HAS BEEN WITHDRAWNSTAA | STAA | |
| Designation fees paidDE FR GBAKX | AKX | |
| Request for examination filed17P | 17P | |
| Designated contracting statesAK | AK | |
| Request for extension of the european patentAL;LT;LV;MK;RO;SIAX | AX | |
| Information provided on ipc code assigned before grant7H 04B 17/00 A, 7H 04Q 7/36 BRIC1 | RIC1 | |
| Search report despatchedORIGINAL CODE: 0009013PUAL | PUAL | |
| Designated contracting statesAK | AK | |
| Designated contracting statesAK | AK | |
| Request for extension of the european patentAL;LT;LV;MK;RO;SIAX | AX | |
| Public reference made under article 153(3) epc to a published international application that has entered the european phaseORIGINAL CODE: 0009012PUAI | PUAI |
Numbers
- Publication
- 0903880
- Application
- 983074774
Titles3
- German
- Verfahren und Vorrichtung zur Modellierung der drahtlosen Signalenslaufzeit in Gebäuden
- English
- Method and apparatus for modeling the propagation of wireless signals in buildings
- French
- Procédé et dispositif pour la modelisation de la propagation des signaux sans fil dans des bâtiments
Classification
- CPC, 4
- H04W16/20
- H04B17/318
- H04B17/391
- H04B17/3913
- IPC, 3
- H04B7 26
- H04B17 00
- H04W16 20
Designated states2
- Contracting states, 1
- Sweden
- Extension states, 1
- Slovenia