Lieb–Liniger model

2014-01-25T12:41:24Z

85.181.53.212: /* Definition and solution of the Model */

{{technical|date=November 2010}}

EMSR stands for Expected Marginal Seat Revenue and is a very popular [[heuristic]] in [[Revenue Management]]. There are two versions: EMSRa<ref>Belobaba, P. P., Air Travel Demand and Airline Seat Inventory Management. Flight Transportation Laboratory. Cambridge, MIT. PhD, 1987</ref> and EMSRb,<ref>Belobaba, P. P., Optimal vs. heuristic methods for nested seat allocation. Presentation at ORSA/TIMS Joint National Meeting, 1992</ref> both of which were introduced by Belobaba. Both methods are for ''n''-class, static, single-resource problems. Because the models are static some assumptions apply: classes are indexed in such a way that the fare for the highest class, <math>r_{1}</math>, is higher than the fare for the next highest class, <math>r_{2}</math>, so <math>r_{1}</math> > <math>r_{2}</math> > ... > <math>r_{n}</math>; demand arrives in a strict low to high order in stages that are indexed with ''j'' as well; demand for class ''j'' is distributed with cdf <math>F_j(x)</math>. For simplicity it is also assumed that demand, capacity and the distributions are continuous, although it is not very difficult to drop this assumption.

==EMSRa==
EMSRa is the first version that Belobaba came up with. The idea behind the heuristic is to add the protection limits that are calculated by applying [[Littlewood's rule]] to successive classes. Suppose that we are in stage ''j+1'' and we want to calculate how much capacity we need to protect for stages ''j, j-1,..., 1''. Then we are actually calculating protection limit <math>y</math><sub>j</sub>. To do so we consider every class in ''j, j-1,..., 1'' and compare that class, indexed with ''k'', with ''j+1'' in isolation. For every combination of ''k'' and ''j+1'' we compute the protection level for that class with [[Littlewood's rule]]:

:<math>P ( D_k > y_k^{j+1}) = \frac{r_{j+1}}{r_k}</math>

The idea of EMSRa then is to add all these protection limits to get the protection limit for <math>y_j</math>.

:<math>y_j = \sum_{k=1}^j y_k^{j+1}</math>

However, there is a problem with this method because it does not take the statistical averaging effect into account. Suppose, for example, that classes ''1'' to ''j'' have the same fare ''r'', then EMSRa will calculate the protection limit for <math>y_{j+1}</math> with

:<math> P( D_k > y_k^{j+1}) = \frac{r_{j+1}}{r}</math>

However, because the fare for all these classes is the same they should be aggregated. EMSRa will calculate protection limits that are too conservative. In other words, it will reserve too many seats for the higher fares, thereby rejecting too many low fare bookings. Although having equal fares is not realistic this will also happen if the difference between fares is small. Therefore EMSRb was invented.

==EMSRb==
One of the most widely used RM heuristics is EMSRb. It is simple and produces under certain conditions close to optimal results. Belobaba reports studies in which both EMSRa and EMSRb were compared. He shows that EMSRb is consistently within 0.5 percent of the optimal solution, while EMSRa under certain conditions can deviate more than 1.5 percent from the optimal solution. However, with mixed order of arrival and frequent reoptimization both methods perform well.<ref>Belobaba, P. P., Optimal vs. heuristic methods for nested seat allocation. Presentation at ORSA/TIMS Joint National Meeting, 1992</ref> There is also study by Polt that shows mixed results.<ref>Polt, S., Back to the roots: New results on leg optimization. In 1999 AGIFORS Reservations and Yield Management Study Group Symposium, London, UK, 1999</ref>

EMSRb is also based on an approximation that compares two classes, but it does take the statistical averaging effect into account. Instead of aggregating protection levels, as EMSRa does, it aggregates demand. Suppose we are again in stage <math>j+1</math> and we want to calculate the protection limit <math>y</math><sub>j</sub>. Then first all future demand for classes ''j, j-1,…, 1'' is aggregated:

:<math>S_j = \sum_{k=1}^j D_k</math>

and the weighted revenues are calculated:

:<math> \overline{r}_j = \frac{\sum_{k=1}^j r_k \cdot D_k}{\sum_{k=1}^j D_k}</math>

Then, again with Littlewood’s rule, the protection limit for classes ''j'' and higher is calculated such that:

:<math> P( S_j > y_j ) = \frac{r_{j+1}}{\overline{r}_j}</math>

Usually demand is considered to be independent and distributed normally with a mean and a variance . Using that the protection limits can be calculated as:

:<math> y_j = \mu_j + z_\alpha \cdot \sigma_j</math>

with the mean and variance of the demand to come as <math>\mu_j = \sum_{k=1}^j \mu_k</math> and <math>\sigma_j^2 = \sum_{k=1}^j \sigma_k^2</math> respectively. <math>z_\alpha</math> is calculated with the inverse of the normal distribution <math>z_\alpha = \phi^{-1}(1-\frac{r_{j+1}}{\overline{r}_j})</math>. This is done for each j, giving the protection limit for every class.

==References==

{{Reflist}}

==See also==
* [[Yield management]]
* [[Littlewood's rule]]

{{DEFAULTSORT:Expected Marginal Seat Revenue}}
[[Category:Pricing]]

formulasearchengine - User contributions [en]

Lieb–Liniger model