一般化平均

いくつかの一般化平均のプロットMp1×{\displaystyle M_{p}(1,x)}

数学において、一般化平均(またはオットー・ヘルダーに由来する冪平均ヘルダー平均[ 1 ]は、数値集合を集計するための関数群である。これには、ピタゴラス平均算術平均幾何平均調和平均)が特別な場合として含まれる。

意味

pが非ゼロの実数で、正の実数である場合、これらの正の実数の指数pの一般化平均またはべき乗平均は[ 2 ] [ 3 ]である。×1×n{\displaystyle x_{1},\dots ,x_{n}}

Mp×1×n1n1n×p1/p{\displaystyle M_{p}(x_{1},\dots ,x_{n})=\left({\frac {1}{n}}\sum _{i=1}^{n}x_{i}^{p}\right)^{{1}/{p}}.}

pノルムを参照)。p = 0の場合、 pを幾何平均(以下で証明されるように、指数が0に近づく平均の極限)と等しく設定します。

M0×1×n1n×1/n{\displaystyle M_{0}(x_{1},\dots ,x_{n})=\left(\prod _{i=1}^{n}x_{i}\right)^{1/n}.}

さらに、正の重みw iのシーケンスに対して、重み付き累乗平均を[ 2 ]と定義し、 p = 0の ときは重み付き幾何平均に等しくなります。 Mp(x1,,xn)=(i=1nwixipi=1nwi)1/p{\displaystyle M_{p}(x_{1},\dots ,x_{n})=\left({\frac {\sum _{i=1}^{n}w_{i}x_{i}^{p}}{\sum _{i=1}^{n}w_{i}}}\right)^{{1}/{p}}}

M0(x1,,xn)=(i=1nxiwi)1/i=1nwi.{\displaystyle M_{0}(x_{1},\dots ,x_{n})=\left(\prod _{i=1}^{n}x_{i}^{w_{i}}\right)^{1/\sum _{i=1}^{n}w_{i}}.}

重み付けされていない平均は、すべてのw i = 1 の設定に相当します。

特殊なケース

のいくつかの値では、平均はよく知られた平均に対応します。 p{\displaystyle p}Mp(x1,,xn){\displaystyle M_{p}(x_{1},\dots ,x_{n})}

指定されたケースの一部を視覚的に表したものです。n=2{\displaystyle n=2}
  調和平均: .M1(a,b){\displaystyle M_{-1}(a,b)}
  幾何平均: .M0(a,b){\displaystyle M_{0}(a,b)}
  算術平均: .M1(a,b){\displaystyle M_{1}(a,b)}
  二次平均: .M2(a,b){\displaystyle M_{2}(a,b)}
名前 指数 価値
最小p={\displaystyle p=-\infty }min{x1,,xn}{\displaystyle \min\{x_{1},\dots ,x_{n}\}}
調和平均p=1{\displaystyle p=-1}n1x1++1xn{\displaystyle {\frac {n}{{\frac {1}{x_{1}}}+\dots +{\frac {1}{x_{n}}}}}}
幾何平均p=0{\displaystyle p=0}x1xnn{\displaystyle {\sqrt[{n}]{x_{1}\dots x_{n}}}}
算術平均p=1{\displaystyle p=1}x1++xnn{\displaystyle {\frac {x_{1}+\dots +x_{n}}{n}}}
二乗平均平方根p=2{\displaystyle p=2}x12++xn2n{\displaystyle {\sqrt {\frac {x_{1}^{2}+\dots +x_{n}^{2}}{n}}}}
立方平均p=3{\displaystyle p=3}x13++xn3n3{\displaystyle {\sqrt[{3}]{\frac {x_{1}^{3}+\dots +x_{n}^{3}}{n}}}}
最大p=+{\displaystyle p=+\infty }max{x1,,xn}{\displaystyle \max\{x_{1},\dots ,x_{n}\}}

(幾何平均)の証明limp0Mp=M0{\textstyle \lim _{p\to 0}M_{p}=M_{0}}

証明のために、一般性を失うことなく、次 のことを仮定する 。wi[0,1]{\displaystyle w_{i}\in [0,1]}i=1nwi=1.{\displaystyle \sum _{i=1}^{n}w_{i}=1.}

指数関数 の定義を次のように書き直すことができる。Mp{\displaystyle M_{p}}

Mp(x1,,xn)=exp(ln[(i=1nwixip)1/p])=exp(ln(i=1nwixip)p){\displaystyle M_{p}(x_{1},\dots ,x_{n})=\exp {\left(\ln {\left[\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{1/p}\right]}\right)}=\exp {\left({\frac {\ln {\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)}}{p}}\right)}}

p → 0の極限において、指数関数の偏角にロピタルの定理を適用することができる。ただしp ≠ 0 であり、 w iの和は1に等しい(一般性を損なうことなく)と仮定する。 [ 4 ]分子と分母をpについて微分すると、 pR{\displaystyle p\in \mathbb {R} }limp0ln(i=1nwixip)p=limp0i=1nwixiplnxij=1nwjxjp1=limp0i=1nwixiplnxij=1nwjxjp=i=1nwilnxij=1nwj=i=1nwilnxi=ln(i=1nxiwi){\displaystyle {\begin{aligned}\lim _{p\to 0}{\frac {\ln {\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)}}{p}}&=\lim _{p\to 0}{\frac {\frac {\sum _{i=1}^{n}w_{i}x_{i}^{p}\ln {x_{i}}}{\sum _{j=1}^{n}w_{j}x_{j}^{p}}}{1}}\\&=\lim _{p\to 0}{\frac {\sum _{i=1}^{n}w_{i}x_{i}^{p}\ln {x_{i}}}{\sum _{j=1}^{n}w_{j}x_{j}^{p}}}\\&={\frac {\sum _{i=1}^{n}w_{i}\ln {x_{i}}}{\sum _{j=1}^{n}w_{j}}}\\&=\sum _{i=1}^{n}w_{i}\ln {x_{i}}\\&=\ln {\left(\prod _{i=1}^{n}x_{i}^{w_{i}}\right)}\end{aligned}}}

指数関数の連続性により、上記の関係式に代入して 希望する結果を得ることができます。[ 2 ]limp0Mp(x1,,xn)=exp(ln(i=1nxiwi))=i=1nxiwi=M0(x1,,xn){\displaystyle \lim _{p\to 0}M_{p}(x_{1},\dots ,x_{n})=\exp {\left(\ln {\left(\prod _{i=1}^{n}x_{i}^{w_{i}}\right)}\right)}=\prod _{i=1}^{n}x_{i}^{w_{i}}=M_{0}(x_{1},\dots ,x_{n})}

証明とlimpMp=M{\textstyle \lim _{p\to \infty }M_{p}=M_{\infty }}limpMp=M{\textstyle \lim _{p\to -\infty }M_{p}=M_{-\infty }}

(おそらくラベルを付け直して用語を結合した後で) と仮定する。すると x1xn{\displaystyle x_{1}\geq \dots \geq x_{n}}

limpMp(x1,,xn)=limp(i=1nwixip)1/p=x1limp(i=1nwi(xix1)p)1/p=x1=M(x1,,xn).{\displaystyle {\begin{aligned}\lim _{p\to \infty }M_{p}(x_{1},\dots ,x_{n})&=\lim _{p\to \infty }\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{1/p}\\&=x_{1}\lim _{p\to \infty }\left(\sum _{i=1}^{n}w_{i}\left({\frac {x_{i}}{x_{1}}}\right)^{p}\right)^{1/p}\\&=x_{1}=M_{\infty }(x_{1},\dots ,x_{n}).\end{aligned}}}

の式は次のようになります。 M{\displaystyle M_{-\infty }}M(x1,,xn)=1M(1/x1,,1/xn)=xn.{\displaystyle M_{-\infty }(x_{1},\dots ,x_{n})={\frac {1}{M_{\infty }(1/x_{1},\dots ,1/x_{n})}}=x_{n}.}

プロパティ

を正の実数列とすると、以下の性質が成り立つ:[ 1 ]x1,,xn{\displaystyle x_{1},\dots ,x_{n}}

  1. min(x1,,xn)Mp(x1,,xn)max(x1,,xn){\displaystyle \min(x_{1},\dots ,x_{n})\leq M_{p}(x_{1},\dots ,x_{n})\leq \max(x_{1},\dots ,x_{n})}
    各一般化平均は常にx値の最小値と最大値の間にあります。
  2. Mp(x1,,xn)=Mp(P(x1,,xn)){\displaystyle M_{p}(x_{1},\dots ,x_{n})=M_{p}(P(x_{1},\dots ,x_{n}))}、ここでは順列演算子です。P{\displaystyle P}
    それぞれの一般化平均はその引数の対称関数です。一般化平均の引数を並べ替えても、その値は変わりません。
  3. Mp(bx1,,bxn)=bMp(x1,,xn){\displaystyle M_{p}(bx_{1},\dots ,bx_{n})=b\cdot M_{p}(x_{1},\dots ,x_{n})}
    ほとんどの平均と同様に、一般化平均は引数x 1 , ..., x nの同次関数です。つまり、b が正の実数である場合、指数pを持つ一般化平均は、引数x 1 , ..., x nの一般化平均のb倍に等しくなります。bx1,,bxn{\displaystyle b\cdot x_{1},\dots ,b\cdot x_{n}}
  4. Mp(x1,,xnk)=Mp[Mp(x1,,xk),Mp(xk+1,,x2k),,Mp(x(n1)k+1,,xnk)]{\displaystyle M_{p}(x_{1},\dots ,x_{n\cdot k})=M_{p}\left[M_{p}(x_{1},\dots ,x_{k}),M_{p}(x_{k+1},\dots ,x_{2\cdot k}),\dots ,M_{p}(x_{(n-1)\cdot k+1},\dots ,x_{n\cdot k})\right]}
    準算術平均と同様に、平均の計算は等しいサイズのサブブロックに分割することができます。これにより、必要に応じて分割統治アルゴリズムを用いて平均を計算することができます。

一般化平均不平等

2つの異なる正の数aとbの最大値 ( ab ) >二乗平均平方根( RMS )または二次平均( QM ) >算術平均( AM ) >幾何平均( GM ) >調和平均( HM ) >最小値 ( ab )を言葉なしで証明する[1 ]

一般に、p < qの場合、 x 1 = x 2 = ... = x n の場合にのみ、 2 つの平均は等しくなります。 Mp(x1,,xn)Mq(x1,,xn){\displaystyle M_{p}(x_{1},\dots ,x_{n})\leq M_{q}(x_{1},\dots ,x_{n})}

この不等式は、 pqの実数値、および正と負の無限大値 に対しても当てはまります。

これは、すべての実数pに対して であるという事実から成り、 これはJensen の不等式 を使用して証明できます。 pMp(x1,,xn)0{\displaystyle {\frac {\partial }{\partial p}}M_{p}(x_{1},\dots ,x_{n})\geq 0}

特に、pが{−1, 0, 1}の範囲内にある場合、一般化平均不等式は、ピタゴラス平均不等式だけでなく、算術平均と幾何平均の不等式も意味します。

重み付き不等式の証明

重み付きべき乗平均不等式を証明します。証明のために、一般性を損なうことなく、以下の仮定を置きます。 wi[0,1]i=1nwi=1{\displaystyle {\begin{aligned}w_{i}\in [0,1]\\\sum _{i=1}^{n}w_{i}=1\end{aligned}}}

重み付けされていない累乗平均の証明は、w i = 1/ nを代入することによって簡単に得られます。

反対符号の平均間の不等式の同値性

指数pqのべき乗平均間の平均が成り立つと仮定します。 これを適用すると、次のようになります。 (i=1nwixip)1/p(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{1/p}\geq \left(\sum _{i=1}^{n}w_{i}x_{i}^{q}\right)^{1/q}}(i=1nwixip)1/p(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}{\frac {w_{i}}{x_{i}^{p}}}\right)^{1/p}\geq \left(\sum _{i=1}^{n}{\frac {w_{i}}{x_{i}^{q}}}\right)^{1/q}}

両辺を−1乗します(正の実数では厳密に減少する関数です)。 (i=1nwixip)1/p=(1i=1nwi1xip)1/p(1i=1nwi1xiq)1/q=(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}w_{i}x_{i}^{-p}\right)^{-1/p}=\left({\frac {1}{\sum _{i=1}^{n}w_{i}{\frac {1}{x_{i}^{p}}}}}\right)^{1/p}\leq \left({\frac {1}{\sum _{i=1}^{n}w_{i}{\frac {1}{x_{i}^{q}}}}}\right)^{1/q}=\left(\sum _{i=1}^{n}w_{i}x_{i}^{-q}\right)^{-1/q}}

指数が− pqである平均値の不等式が得られ、同じ推論を逆に適用することで不等式​​が同等であることが証明され、これは後の証明のいくつかで使用されます。

幾何平均

任意のq > 0および非負の重みの合計が1の 場合、次の不等式が成り立ちます。(i=1nwixiq)1/qi=1nxiwi(i=1nwixiq)1/q.{\displaystyle \left(\sum _{i=1}^{n}w_{i}x_{i}^{-q}\right)^{-1/q}\leq \prod _{i=1}^{n}x_{i}^{w_{i}}\leq \left(\sum _{i=1}^{n}w_{i}x_{i}^{q}\right)^{1/q}.}

証明は、対数が凹であるという事実を利用して、ジェンセンの不等式から導かれます。 logi=1nxiwi=i=1nwilogxilogi=1nwixi.{\displaystyle \log \prod _{i=1}^{n}x_{i}^{w_{i}}=\sum _{i=1}^{n}w_{i}\log x_{i}\leq \log \sum _{i=1}^{n}w_{i}x_{i}.}

指数関数を両辺に 適用し、厳密に増加関数として不等式の符号が保存されることを観察すると、次の式が得られます。i=1nxiwii=1nwixi.{\displaystyle \prod _{i=1}^{n}x_{i}^{w_{i}}\leq \sum _{i=1}^{n}w_{i}x_{i}.}

x iq乗をとると、 i=1nxiqwii=1nwixiqi=1nxiwi(i=1nwixiq)1/q.{\displaystyle {\begin{aligned}&\prod _{i=1}^{n}x_{i}^{q{\cdot }w_{i}}\leq \sum _{i=1}^{n}w_{i}x_{i}^{q}\\&\prod _{i=1}^{n}x_{i}^{w_{i}}\leq \left(\sum _{i=1}^{n}w_{i}x_{i}^{q}\right)^{1/q}.\end{aligned}}}

したがって、 qが正の不等式についてはこれで完了です。 q が負の場合も、最後のステップで符号が入れ替わる点を除いて同じです。

i=1nxiqwii=1nwixiq.{\displaystyle \prod _{i=1}^{n}x_{i}^{-q{\cdot }w_{i}}\leq \sum _{i=1}^{n}w_{i}x_{i}^{-q}.}

もちろん、各辺を負の数-1/ qで乗じると、不等式の方向が入れ替わります。

i=1nxiwi(i=1nwixiq)1/q.{\displaystyle \prod _{i=1}^{n}x_{i}^{w_{i}}\geq \left(\sum _{i=1}^{n}w_{i}x_{i}^{-q}\right)^{-1/q}.}

任意の2つの累乗平均間の不平等

任意のp < qに対して次の不等式が成り立つ ことを証明します。p負でqが正の 場合、この不等式は上で証明したものと等しくなります。 (i=1nwixip)1/p(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{1/p}\leq \left(\sum _{i=1}^{n}w_{i}x_{i}^{q}\right)^{1/q}}(i=1nwixip)1/pi=1nxiwi(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}w_{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}}

pq が正であることの証明は次のとおりです。次の関数を定義します: f  : R +R +fはべき関数なので、2 次導関数があります。 これは、 q > pであるため、 f のドメイン内で厳密に正であり、fが凸であることがわかります。 f(x)=xqp{\displaystyle f(x)=x^{\frac {q}{p}}}f(x)=(qp)(qp1)xqp2{\displaystyle f''(x)=\left({\frac {q}{p}}\right)\left({\frac {q}{p}}-1\right)x^{{\frac {q}{p}}-2}}

これとジェンセンの不等式を使うと、次の式が得られます。 両辺を1/ q 乗すると( 1/ qは正なので増加関数です)、証明すべき不等式が得られます。 f(i=1nwixip)i=1nwif(xip)(i=1nwixip)q/pi=1nwixiq{\displaystyle {\begin{aligned}f\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)&\leq \sum _{i=1}^{n}w_{i}f(x_{i}^{p})\\[3pt]\left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{q/p}&\leq \sum _{i=1}^{n}w_{i}x_{i}^{q}\end{aligned}}}

(i=1nwixip)1/p(i=1nwixiq)1/q{\displaystyle \left(\sum _{i=1}^{n}w_{i}x_{i}^{p}\right)^{1/p}\leq \left(\sum _{i=1}^{n}w_{i}x_{i}^{q}\right)^{1/q}}

前に示した同値性を使用して、負のpqをそれぞれ−q−pに置き換えることで、不等式を証明できます。

一般化f平均

べき乗平均はさらに一般化f平均に一般化できます。

Mf(x1,,xn)=f1(1ni=1nf(xi)){\displaystyle M_{f}(x_{1},\dots ,x_{n})=f^{-1}\left({{\frac {1}{n}}\cdot \sum _{i=1}^{n}{f(x_{i})}}\right)}

これはf ( x ) = log( x )の極限を用いずに幾何平均をカバーします。べき乗平均はf ( x ) = x pで得られます。これらの平均の性質は de Carvalho (2016) で研究されています。[ 3 ]

アプリケーション

信号処理

べき乗平均は非線形移動平均として機能し、 pが小さい場合には小さな信号値にシフトし、 pが大きい場合には大きな信号値を強調します。移動算術平均の効率的な実装が と呼ばれる場合、移動べき乗平均は次のHaskellsmoothコードに従って実装できます。

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

参照

注記

  1. ^ NM = a、 PM = bとすると、 AM = abAM、半径r = AQ = AG となります。ピタゴラスの定理を用いると、 QM² = AQ² + AM² ∴ QM = √ AQ² + AM² = QMとなります。ピタゴラスの定理を用いると、 AM² = AG² + GM² ∴ GM = √ AM² − AG² = GMとなります。相似三角形を用いると、HM/GM = GM/午前 ∴ HM = GM²/午前 = HM

参考文献

  1. ^ a b Sýkora, Stanislav (2009). 「数学的平均と平均値:基本的な性質」. 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 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. ^平均とその不等式ハンドブック (数学とその応用)

さらに読む

  • Bullen, PS (2003). 「第3章 パワー・ミーンズ」. 『ミーンズとその不等式ハンドブック』 . ドルドレヒト、オランダ: Kluwer. pp.  175– 265.