Jump to content

Generalized mean

From Wikipedia, the free encyclopedia
Plot of several generalized means Mp(1,x)

In mathematics, generalized means (or power mean or Hölder mean from Otto Hölder)[1] are a family of functions for aggregating sets of numbers. These include as special cases the Pythagorean means (arithmetic, geometric, and harmonic means).

Definition

[edit]

If p is a non-zero real number, and x1,…,xn are positive real numbers, then the generalized mean or power mean with exponent p of these positive real numbers is[2][3]

Mp(x1,…,xn)=(1n∑i=1nxip)1/p.

(See p-norm). For p = 0 we set it equal to the geometric mean (which is the limit of means with exponents approaching zero, as proved below):

M0(x1,…,xn)=(∏i=1nxi)1/n.

Furthermore, for a sequence of positive weights wi we define the weighted power mean as[2] Mp(x1,…,xn)=(∑i=1nwixip∑i=1nwi)1/p and when p = 0, it is equal to the weighted geometric mean:

M0(x1,…,xn)=(∏i=1nxiwi)1/∑i=1nwi.

The unweighted means correspond to setting all wi = 1.

Special cases

[edit]

For some values of p, the mean Mp(x1,…,xn) corresponds to a well known mean.

A visual depiction of some of the specified cases for n=2.
  Harmonic mean: M−1(a,b).
  Geometric mean: M0(a,b).
  Arithmetic mean: M1(a,b).
  Quadratic mean: M2(a,b).
Name Exponent Value
Minimum p=−∞ min⁡{x1,…,xn}
Harmonic mean p=−1 n1x1+⋯+1xn
Geometric mean p=0 x1…xnn
Arithmetic mean p=1 x1+⋯+xnn
Root mean square p=2 x12+⋯+xn2n
Cubic mean p=3 x13+⋯+xn3n3
Maximum p=+∞ max⁡{x1,…,xn}


Proof of limp→0Mp=M0 (geometric mean)

For the purpose of the proof, we will assume without loss of generality that wi∈[0,1] and ∑i=1nwi=1.

We can rewrite the definition of Mp using the exponential function as

Mp(x1,…,xn)=exp⁡(ln⁡[(∑i=1nwixip)1/p])=exp⁡(ln⁡(∑i=1nwixip)p)

In the limit p → 0, we can apply L'Hôpital's rule to the argument of the exponential function. We assume that p∈ℝ but p ≠ 0, and that the sum of wi is equal to 1 (without loss in generality);[4] differentiating the numerator and denominator with respect to p, we have limp→0ln⁡(∑i=1nwixip)p=limp→0∑i=1nwixipln⁡xi∑j=1nwjxjp1=limp→0∑i=1nwixipln⁡xi∑j=1nwjxjp=∑i=1nwiln⁡xi∑j=1nwj=∑i=1nwiln⁡xi=ln⁡(∏i=1nxiwi)

By the continuity of the exponential function, we can substitute back into the above relation to obtain limp→0Mp(x1,…,xn)=exp⁡(ln⁡(∏i=1nxiwi))=∏i=1nxiwi=M0(x1,…,xn) as desired.[2]

Proof of limp→∞Mp=M∞ and limp→−∞Mp=M−∞

Assume (possibly after relabeling and combining terms together) that x1≥⋯≥xn. Then

limp→∞Mp(x1,…,xn)=limp→∞(∑i=1nwixip)1/p=x1limp→∞(∑i=1nwi(xix1)p)1/p=x1=M∞(x1,…,xn).

The formula for M−∞ follows from M−∞(x1,…,xn)=1M∞(1/x1,…,1/xn)=xn.

Properties

[edit]

Let x1,…,xn be a sequence of positive real numbers, then the following properties hold:[1]

  1. min⁡(x1,…,xn)≤Mp(x1,…,xn)≤max⁡(x1,…,xn).
    Each generalized mean always lies between the smallest and largest of the x values.
  2. Mp(x1,…,xn)=Mp(P(x1,…,xn)), where P is a permutation operator.
    Each generalized mean is a symmetric function of its arguments; permuting the arguments of a generalized mean does not change its value.
  3. Mp(bx1,…,bxn)=b⋅Mp(x1,…,xn).
    Like most means, the generalized mean is a homogeneous function of its arguments x1, ..., xn. That is, if b is a positive real number, then the generalized mean with exponent p of the numbers b⋅x1,…,b⋅xn is equal to b times the generalized mean of the numbers x1, ..., xn.
  4. Mp(x1,…,xn⋅k)=Mp[Mp(x1,…,xk),Mp(xk+1,…,x2⋅k),…,Mp(x(n−1)⋅k+1,…,xn⋅k)].
    Like the quasi-arithmetic means, the computation of the mean can be split into computations of equal sized sub-blocks. This enables use of a divide and conquer algorithm to calculate the means, when desirable.

Generalized mean inequality

[edit]
Geometric proof without words that max (a,b) > root mean square (RMS) or quadratic mean (QM) > arithmetic mean (AM) > geometric mean (GM) > harmonic mean (HM) > min (a,b) of two distinct positive numbers a and b[note 1]

In general, if p < q, then Mp(x1,…,xn)≤Mq(x1,…,xn) and the two means are equal if and only if x1 = x2 = ... = xn.

The inequality is true for real values of p and q, as well as positive and negative infinity values.

It follows from the fact that, for all real p, ∂∂pMp(x1,…,xn)≥0 which can be proved using Jensen's inequality.

In particular, for p in {−1, 0, 1}, the generalized mean inequality implies the Pythagorean means inequality as well as the inequality of arithmetic and geometric means.

Proof of the weighted inequality

[edit]

We will prove the weighted power mean inequality. For the purpose of the proof we will assume the following without loss of generality: wi∈[0,1]∑i=1nwi=1

The proof for unweighted power means can be easily obtained by substituting wi = 1/n.

Equivalence of inequalities between means of opposite signs

[edit]

Suppose an average between power means with exponents p and q holds: (∑i=1nwixip)1/p≥(∑i=1nwixiq)1/q applying this, then: (∑i=1nwixip)1/p≥(∑i=1nwixiq)1/q

We raise both sides to the power of −1 (strictly decreasing function in positive reals): (∑i=1nwixi−p)−1/p=(1∑i=1nwi1xip)1/p≤(1∑i=1nwi1xiq)1/q=(∑i=1nwixi−q)−1/q

We get the inequality for means with exponents −p and −q, and we can use the same reasoning backwards, thus proving the inequalities to be equivalent, which will be used in some of the later proofs.

Geometric mean

[edit]

For any q > 0 and non-negative weights summing to 1, the following inequality holds: (∑i=1nwixi−q)−1/q≤∏i=1nxiwi≤(∑i=1nwixiq)1/q.

The proof follows from Jensen's inequality, making use of the fact the logarithm is concave: log⁡∏i=1nxiwi=∑i=1nwilog⁡xi≤log⁡∑i=1nwixi.

By applying the exponential function to both sides and observing that as a strictly increasing function it preserves the sign of the inequality, we get ∏i=1nxiwi≤∑i=1nwixi.

Taking q-th powers of the xi yields ∏i=1nxiq⋅wi≤∑i=1nwixiq∏i=1nxiwi≤(∑i=1nwixiq)1/q.

Thus, we are done for the inequality with positive q; the case for negatives is identical but for the swapped signs in the last step:

∏i=1nxi−q⋅wi≤∑i=1nwixi−q.

Of course, taking each side to the power of a negative number -1/q swaps the direction of the inequality.

∏i=1nxiwi≥(∑i=1nwixi−q)−1/q.

Inequality between any two power means

[edit]

We are to prove that for any p < q the following inequality holds: (∑i=1nwixip)1/p≤(∑i=1nwixiq)1/q if p is negative, and q is positive, the inequality is equivalent to the one proved above: (∑i=1nwixip)1/p≤∏i=1nxiwi≤(∑i=1nwixiq)1/q

The proof for positive p and q is as follows: Define the following function: f : R+ → R+ f(x)=xqp. f is a power function, so it does have a second derivative: f″(x)=(qp)(qp−1)xqp−2 which is strictly positive within the domain of f, since q > p, so we know f is convex.

Using this, and the Jensen's inequality we get: f(∑i=1nwixip)≤∑i=1nwif(xip)(∑i=1nwixip)q/p≤∑i=1nwixiq after raising both side to the power of 1/q (an increasing function, since 1/q is positive) we get the inequality which was to be proven:

(∑i=1nwixip)1/p≤(∑i=1nwixiq)1/q

Using the previously shown equivalence we can prove the inequality for negative p and q by replacing them with −q and −p, respectively.

Generalized f-mean

[edit]

The power mean could be generalized further to the generalized f-mean:

Mf(x1,…,xn)=f−1(1n⋅∑i=1nf(xi))

This covers the geometric mean without using a limit with f(x) = log(x). The power mean is obtained for f(x) = xp. Properties of these means are studied in de Carvalho (2016).[3]

Applications

[edit]

Signal processing

[edit]

A power mean serves a non-linear moving average which is shifted towards small signal values for small p and emphasizes big signal values for big p. Given an efficient implementation of a moving arithmetic mean called smooth one can implement a moving power mean according to the following Haskell code.

powerSmooth :: Floating a => ([a] -> [a]) -> a -> [a] -> [a]
powerSmooth smooth p = map (** recip p) . smooth . map (**p)

See also

[edit]

Notes

[edit]
  1. ↑ If NM = a and PM = b. AM = AM of a and b, and radius r = AQ = AG.
    Using Pythagoras' theorem, QM2 = AQ2 + AM2 ∴ QM = AQ2+AM2 = QM.
    Using Pythagoras' theorem, AM2 = AG2 + GM2 ∴ GM = AM2−AG2 = GM.
    Using similar triangles, ⁠HM/GM⁠ = ⁠GM/AM⁠ ∴ HM = ⁠GM2/AM⁠ = HM.

References

[edit]
  1. 1 2 Sýkora, Stanislav (2009). "Mathematical means and averages: basic properties". Stan's Library. III. Castano Primo, Italy. doi:10.3247/SL3Math09.001.
  2. 1 2 3 P. S. Bullen: Handbook of Means and Their Inequalities. Dordrecht, Netherlands: Kluwer, 2003, pp. 175-177
  3. 1 2 de Carvalho, Miguel (2016). "Mean, what do you Mean?". The American Statistician. 70 (3): 764‒776. doi:10.1080/00031305.2016.1148632. hdl:20.500.11820/fd7a8991-69a4-4fe5-876f-abcd2957a88c.
  4. ↑ Handbook of Means and Their Inequalities (Mathematics and Its Applications).

Further reading

[edit]
  • Bullen, P. S. (2003). "Chapter III - The Power Means". Handbook of Means and Their Inequalities. Dordrecht, Netherlands: Kluwer. pp. 175–265.
[edit]