Generalized mean: Difference between revisions

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search
imported>SimLibrarian
m https, image caption style (MOS:CAPFRAG)
imported>Stqckfish
m standardized spelling
 
Line 4: Line 4:
[[File:Generalized means of 1, x.svg|400px|thumb|right|Plot of several generalized means <math>M_p(1, x)</math>]]
[[File:Generalized means of 1, x.svg|400px|thumb|right|Plot of several generalized means <math>M_p(1, x)</math>]]


In [[mathematics]], '''generalised means''' (or '''power mean''' or '''Hölder mean''' from [[Otto Hölder]])<ref name=sykora/> are a family of functions for aggregating sets of numbers. These include as special cases the [[Pythagorean means]] ([[arithmetic mean|arithmetic]], [[geometric mean|geometric]], and [[harmonic mean|harmonic]] [[mean]]s).
In [[mathematics]], '''generalized means''' (or '''power mean''' or '''Hölder mean''' from [[Otto Hölder]])<ref name=sykora/> are a family of functions for aggregating sets of numbers. These include as special cases the [[Pythagorean means]] ([[arithmetic mean|arithmetic]], [[geometric mean|geometric]], and [[harmonic mean|harmonic]] [[mean]]s).


==Definition==
==Definition==
If {{mvar|p}} is a non-zero [[real number]], and <math>x_1, \dots, x_n</math> are positive real numbers, then the '''generalized mean''' or '''power mean''' with exponent {{mvar|p}} of these positive real numbers is<ref name="Bullen1"/><ref name = "dC2016">{{cite journal|last=de Carvalho|first=Miguel|title=Mean, what do you Mean?|journal=[[The American Statistician]]|year=2016|volume=70|issue=3|pages=764‒776|doi=10.1080/00031305.2016.1148632|url=https://zenodo.org/record/895400|hdl=20.500.11820/fd7a8991-69a4-4fe5-876f-abcd2957a88c|hdl-access=free}}</ref>
If {{mvar|p}} is a non-zero [[real number]], and <math>x_1, \dots, x_n</math> are [[positive real numbers]], then the '''generalized mean''' or '''power mean''' with exponent {{mvar|p}} of these positive real numbers is<ref name="Bullen1"/><ref name = "dC2016">{{cite journal|last=de Carvalho|first=Miguel|title=Mean, what do you Mean?|journal=[[The American Statistician]]|year=2016|volume=70|issue=3|pages=764‒776|doi=10.1080/00031305.2016.1148632|url=https://zenodo.org/record/895400|hdl=20.500.11820/fd7a8991-69a4-4fe5-876f-abcd2957a88c|hdl-access=free}}</ref>


<math display=block>M_p(x_1,\dots,x_n) = \left( \frac{1}{n} \sum_{i=1}^n x_i^p \right)^{{1}/{p}} .</math>
<math display=block>M_p(x_1,\dots,x_n) = \left( \frac{1}{n} \sum_{i=1}^n x_i^p \right)^{{1}/{p}} .</math>
Line 24: Line 24:


== Special cases ==
== Special cases ==
A few particular values of {{mvar|p}} yield special cases with their own names:<ref name="mw">{{MathWorld|title=Power Mean|urlname=PowerMean}} (retrieved 2019-08-17)</ref>
 
;[[minimum]] :<math>M_{-\infty}(x_1,\dots,x_n) = \lim_{p\to-\infty} M_p(x_1,\dots,x_n) = \min \{x_1,\dots,x_n\}</math>
For some values of <math>p</math>, the mean <math>M_p(x_1, \dots, x_n)</math> corresponds to a well known mean.
;[[Image:MathematicalMeans.svg|thumb|A visual depiction of some of the specified cases for {{math|1=''n'' = 2}} with {{math|1=''a'' = ''x''{{sub|1}} = ''M''{{sub|∞}}}} and {{math|1=''b'' = ''x''{{sub|2}} = ''M''{{sub|−∞}}}}: {{legend|magenta|harmonic mean, {{math|''H'' {{=}} ''M''{{sub|−1}}(''a'', ''b'')}},}} {{legend|blue|geometric mean, {{math|''G'' {{=}} ''M''{{sub|0}}(''a'', ''b'')}}}} {{legend|red|arithmetic mean, {{math|''A'' {{=}} ''M''{{sub|1}}(''a'', ''b'')}}}} {{legend|lime|quadratic mean, {{math|''Q'' {{=}} ''M''{{sub|2}}(''a'', ''b'')}}}}]][[harmonic mean]] :<math>M_{-1}(x_1,\dots,x_n) = \frac{n}{\frac{1}{x_1}+\dots+\frac{1}{x_n}}</math>
 
;[[geometric mean]] <math>M_0(x_1,\dots,x_n) = \lim_{p\to0} M_p(x_1,\dots,x_n) = \sqrt[n]{x_1\cdot\dots\cdot x_n}</math>
[[File:Generalized Means.svg|thumb|A visual depiction of some of the specified cases for <math>n = 2</math>. {{legend|magenta|Harmonic mean: <math>M_{-1}(a, b)</math>.}} {{legend|blue|Geometric mean: <math>M_0(a, b)</math>.}} {{legend|red|Arithmetic mean: <math>M_1(a, b)</math>.}} {{legend|lime|Quadratic mean: <math>M_2(a, b)</math>.}}]]
;[[arithmetic mean]] :<math>M_1(x_1,\dots,x_n) = \frac{x_1 + \dots + x_n}{n}</math>
 
;[[root mean square]]{{anchor|Quadratic}}<br/>or quadratic mean<ref>{{cite book |last1=Thompson |first1=Sylvanus P. |title=Calculus Made Easy |date=1965 |publisher=Macmillan International Higher Education |isbn=9781349004874 |page=185 |url=https://books.google.com/books?id=6VJdDwAAQBAJ&pg=PA185 |access-date=5 July 2020 }}{{Dead link|date=May 2024 |bot=InternetArchiveBot |fix-attempted=yes }}</ref><ref>{{cite book |last1=Jones |first1=Alan R. |title=Probability, Statistics and Other Frightening Stuff |date=2018 |publisher=Routledge |isbn=9781351661386 |page=48 |url=https://books.google.com/books?id=OvtsDwAAQBAJ&pg=PA48 |access-date=5 July 2020}}</ref> :<math>M_2(x_1,\dots,x_n) = \sqrt{\frac{x_1^2 + \dots + x_n^2}{n}}</math>
{| class="wikitable"
;[[cubic mean]] :<math>M_3(x_1,\dots,x_n) = \sqrt[3]{\frac{x_1^3 + \dots + x_n^3}{n}}</math>
|-
;[[maximum]] :<math>M_{+\infty}(x_1,\dots,x_n) = \lim_{p\to\infty} M_p(x_1,\dots,x_n) = \max \{x_1,\dots,x_n\}</math>
! Name
! Exponent
! Value
|-
| [[Minimum]]
| <math>p = -\infty</math>
| <math>\min \{x_1, \dots, x_n\}</math>
|-
| [[Harmonic mean]]
| <math>p = -1</math>
| <math>\frac{n}{\frac{1}{x_1}+\dots+\frac{1}{x_n}}</math>
|-
| [[Geometric mean]]
| <math>p = 0</math>
| <math>\sqrt[n]{x_1\dots x_n}</math>
|-
| [[Arithmetic mean]]
| <math>p = 1</math>
| <math>\frac{x_1 + \dots + x_n}{n}</math>
|-
| [[Root mean square]]
| <math>p = 2</math>
| <math>\sqrt{\frac{x_1^2 + \dots + x_n^2}{n}}</math>
|-
| [[Cubic mean]]
| <math>p = 3</math>
| <math>\sqrt[3]{\frac{x_1^3 + \dots + x_n^3}{n}}</math>
|-
| [[Maximum]]
| <math>p = +\infty</math>
| <math>\max\{x_1, \dots, x_n\}</math>
|}
 
 


{{Math proof|title=Proof of <math display="inline"> \lim_{p \to 0} M_p = M_0 </math> (geometric mean)|proof=For the purpose of the proof, we will assume without loss of generality that  
{{Math proof|title=Proof of <math display="inline"> \lim_{p \to 0} M_p = M_0 </math> (geometric mean)|proof=For the purpose of the proof, we will assume without loss of generality that  
Line 96: Line 129:


==Proof of the weighted inequality==
==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:
We will prove the weighted power mean inequality. For the purpose of the proof we will assume the following [[without loss of generality]]:
<math display=block>\begin{align}
<math display="block">\begin{align}
   w_i \in [0, 1] \\
   w_i \in [0, 1] \\
   \sum_{i=1}^nw_i = 1
   \sum_{i=1}^nw_i = 1
Line 147: Line 180:
<math display="block">\left(\sum_{i=1}^nw_i x_i^p\right)^{1/p} \leq \prod_{i=1}^n x_i^{w_i} \leq \left(\sum_{i=1}^n w_i x_i^q\right)^{1/q}</math>
<math display="block">\left(\sum_{i=1}^nw_i x_i^p\right)^{1/p} \leq \prod_{i=1}^n x_i^{w_i} \leq \left(\sum_{i=1}^n w_i x_i^q\right)^{1/q}</math>


The proof for positive {{mvar|p}} and {{mvar|q}} is as follows: Define the following function: {{math|''f'' : '''R'''<sub>+</sub> → '''R'''<sub>+</sub>}} <math>f(x)=x^{\frac{q}{p}}</math>. {{mvar|f}} is a power function, so it does have a second derivative:
The proof for positive {{mvar|p}} and {{mvar|q}} is as follows: Define the following function: {{math|''f'' : '''R'''<sub>+</sub> → '''R'''<sub>+</sub>}} <math>f(x)=x^{\frac{q}{p}}</math>. {{mvar|f}} is a power function, so it does have a [[second derivative]]:
<math display="block">f''(x) = \left(\frac{q}{p} \right) \left( \frac{q}{p}-1 \right)x^{\frac{q}{p}-2}</math>
<math display="block">f''(x) = \left(\frac{q}{p} \right) \left( \frac{q}{p}-1 \right)x^{\frac{q}{p}-2}</math>
which is strictly positive within the domain of {{mvar|f}}, since {{math|''q'' > ''p''}}, so we know {{mvar|f}} is convex.
which is strictly positive within the domain of {{mvar|f}}, since {{math|''q'' > ''p''}}, so we know {{mvar|f}} is convex.

Latest revision as of 16:44, 4 August 2025

Template:Short description Template:More citations needed

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 Template:Mvar is a non-zero real number, and x1,,xn are positive real numbers, then the generalized mean or power mean with exponent Template:Mvar of these positive real numbers is[2][3]

Mp(x1,,xn)=(1ni=1nxip)1/p.

(See [[Norm (mathematics)#p-norm|Template:Mvar-norm]]). For Template:Math 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 Template:Mvar we define the weighted power mean as[2] Mp(x1,,xn)=(i=1nwixipi=1nwi)1/p and when Template:Math, it is equal to the weighted geometric mean:

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

The unweighted means correspond to setting all Template:Math.

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. <templatestyles src="Legend/styles.css" />
  Harmonic mean: M1(a,b).
<templatestyles src="Legend/styles.css" />
  Geometric mean: M0(a,b).
<templatestyles src="Legend/styles.css" />
  Arithmetic mean: M1(a,b).
<templatestyles src="Legend/styles.css" />
  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}


Template:Math proof

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).Template:Block indent
  2. Mp(x1,,xn)=Mp(P(x1,,xn)), where P is a permutation operator.Template:Block indent
  3. Mp(bx1,,bxn)=bMp(x1,,xn).Template:Block indent
  4. Mp(x1,,xnk)=Mp[Mp(x1,,xk),Mp(xk+1,,x2k),,Mp(x(n1)k+1,,xnk)].Template:Block indent

Generalized mean inequality

Template:QM AM GM HM inequality visual proof.svg In general, if Template:Math, then Mp(x1,,xn)Mq(x1,,xn) and the two means are equal if and only if Template:Math.

The inequality is true for real values of Template:Mvar and Template:Mvar, as well as positive and negative infinity values.

It follows from the fact that, for all real Template:Mvar, pMp(x1,,xn)0 which can be proved using Jensen's inequality.

In particular, for Template:Mvar in Template:Math, 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 Template:Math.

Equivalence of inequalities between means of opposite signs

Suppose an average between power means with exponents Template:Mvar and Template:Mvar 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 Template:Math and Template:Math, 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 Template:Math 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 Template:Mvar-th powers of the Template:Mvar yields i=1nxiqwii=1nwixiqi=1nxiwi(i=1nwixiq)1/q.

Thus, we are done for the inequality with positive Template:Mvar; 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 Template:Math 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 Template:Math the following inequality holds: (i=1nwixip)1/p(i=1nwixiq)1/q if Template:Mvar is negative, and Template:Mvar is positive, the inequality is equivalent to the one proved above: (i=1nwixip)1/pi=1nxiwi(i=1nwixiq)1/q

The proof for positive Template:Mvar and Template:Mvar is as follows: Define the following function: Template:Math f(x)=xqp. Template:Mvar is a power function, so it does have a second derivative: f(x)=(qp)(qp1)xqp2 which is strictly positive within the domain of Template:Mvar, since Template:Math, so we know Template:Mvar 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 Template:Math (an increasing function, since Template:Math 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 Template:Mvar and Template:Mvar by replacing them with Template:Mvar and Template:Mvar, respectively.

Generalized f-mean

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

The power mean could be generalized further to the [[generalized f-mean|generalized Template:Mvar-mean]]:

Mf(x1,,xn)=f1(1ni=1nf(xi))

This covers the geometric mean without using a limit with Template:Math. The power mean is obtained for Template:Mvar. 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 Template:Mvar and emphasizes big signal values for big Template:Mvar. 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

Template:Cols

Template:Colend

Notes

Template:Notelist Template:Reflist

References

Template:Reflist

Further reading

  • Script error: No such module "citation/CS1".

External links

  1. a b Script error: No such module "Citation/CS1".
  2. a b Cite error: Invalid <ref> tag; no text was provided for refs named Bullen1
  3. a b Script error: No such module "Citation/CS1".