Generalized mean

From Wikipedia, the free encyclopedia

Template:Short description

Page Module:Message box/ambox.css has no content.

File:Generalized means of 1, x.svg
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

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)=(1ni=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=1nwixipi=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

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

File:Generalized Means.svg
A visual depiction of some of the specified cases for n=2. Page Template:Legend/styles.css has no content.
  Harmonic mean: M1(a,b).
Page Template:Legend/styles.css has no content.
  Geometric mean: M0(a,b).
Page Template:Legend/styles.css has no content.
  Arithmetic mean: M1(a,b).
Page Template:Legend/styles.css has no content.
  Quadratic mean: M2(a,b).
Name Exponent Value
Minimum p= min{x1,,xn}
Harmonic mean p=1 n1x1++1xn
Geometric mean p=0 x1xnn
Arithmetic mean p=1 x1++xnn
Root mean square p=2 x12++xn2n
Cubic mean p=3 x13++xn3n3
Maximum p=+ max{x1,,xn}


Page Template:Math proof/styles.css has no content.

Proof of limp0Mp=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 limp0ln(i=1nwixip)p=limp0i=1nwixiplnxij=1nwjxjp1=limp0i=1nwixiplnxij=1nwjxjp=i=1nwilnxij=1nwj=i=1nwilnxi=ln(i=1nxiwi)

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

Template:Proof

Properties

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).Page Template:Block indent/styles.css has no content.
    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.Page Template:Block indent/styles.css has no content.
    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)=bMp(x1,,xn).Page Template:Block indent/styles.css has no content.
    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 bx1,,bxn is equal to b times the generalized mean of the numbers x1, ..., xn.
  4. Mp(x1,,xnk)=Mp[Mp(x1,,xk),Mp(xk+1,,x2k),,Mp(x(n1)k+1,,xnk)].Page Template:Block indent/styles.css has no content.
    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

Template:QM AM GM HM inequality visual proof.svg 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

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

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=1nwixip)1/p=(1i=1nwi1xip)1/p(1i=1nwi1xiq)1/q=(i=1nwixiq)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

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

The proof follows from Jensen's inequality, making use of the fact the logarithm is concave: logi=1nxiwi=i=1nwilogxilogi=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=1nxiwii=1nwixi.

Taking q-th powers of the xi yields i=1nxiqwii=1nwixiqi=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=1nxiqwii=1nwixiq.

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

i=1nxiwi(i=1nwixiq)1/q.

Inequality between any two power means

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/pi=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)(qp1)xqp2 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/pi=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

Script error: No such module "Labelled list hatnote".

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

Mf(x1,,xn)=f1(1ni=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

Signal processing

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

Page Template:Div col/styles.css has no content.

Notes

Page Template:Reflist/styles.css has no content.

Page Template:Reflist/styles.css has no content.

References

Page Template:Reflist/styles.css has no content.

  1. ^ a b Page Module:Citation/CS1/styles.css has no content.Sýkora, Stanislav (2009). "Mathematical means and averages: basic properties". Stan's Library. III. Castano Primo, Italy. doi:10.3247/SL3Math09.001.
  2. ^ a b c P. S. Bullen: Handbook of Means and Their Inequalities. Dordrecht, Netherlands: Kluwer, 2003, pp. 175-177
  3. ^ a b Page Module:Citation/CS1/styles.css has no content.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. ^ Page Module:Citation/CS1/styles.css has no content.Handbook of Means and Their Inequalities (Mathematics and Its Applications).

Further reading

  • Page Module:Citation/CS1/styles.css has no content.Bullen, P. S. (2003). "Chapter III - The Power Means". Handbook of Means and Their Inequalities. Dordrecht, Netherlands: Kluwer. pp. 175–265.