幂平均

幂平均(),又稱广义平均()或赫尔德平均(),是一族從數列到實數的函數。幂平均函數的特殊情況包括毕达哥拉斯平均(算术、几何、调和平均),因此可視為毕达哥拉斯平均的一種推廣。

定义
若 p 是一非零实数,可定义非負实数 x_1,\dots,x_n
*p次幂平均*为
:
M_p(x_1,\dots,x_n) = \left( \frac{1}{n} \cdot \sum_{i=1}^n x_{i}^p \right)^{\frac{1}{p}}\,

冪平均在p=0等於幾何平均(冪平均函數p收斂於0時的收斂於幾何平均)

:
M_0(x_1,\dots,x_n) = \left( \prod_{i=1}^n x_{i} \right)^{\frac{1}{n}}\,

性质

  • 和所有平均一样,幂平均是各参数 x_1,\dots,x_n 的一次齐次函数。即若 b 是一个正实数,则 b\cdot x_1,\dots, b\cdot x_n 指数为 p 的幂平均等于 b 倍 x_1,\dots, x_n 的幂平均。
  • 与几何算术平均一样,这种平均的计算可以分解成同样大小的子块来计算。

::
M_p(x_1,\dots,x_{n\cdot k}) =
M_p(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}))

幂平均不等式
一般地,如果 p ,则 M_p(x_1,\dots,x_n) \le M_q(x_1,\dots,x_n) 且这两个平均相等当且仅当 x_1 = x_2 = \cdots = x_n。这由事实

: \forall p\in\mathbb{R}\ \frac{\partial M_p(x_1,\dots,x_n)}{\partial p}\geq 0,

得出,上述不等式可由琴生不等式证明。

特别地,对 p\in\{-1, 0, 1\},幂平均不等式蕴含了毕达哥拉斯平均不等式以及算术几何平均不等式。

特例
幂平均不等式的证明
不同符号的不等式之等价
假设指数 p 与 q 的幂平均间有不等式:
:\sqrt[p]{\sum_{i=1}^nw_ix_i^p}\leq \sqrt[q]{\sum_{i=1}^nw_ix_i^q}

:\sqrt[p]{\sum_{i=1}^n\frac{w_i}{x_i^p}}\leq \sqrt[q]{\sum_{i=1}^n\frac{w_i}{x_i^q}}.
我们在两边取倒数(正实数上的严格递减函数,不等号反向):
:\sqrt[-p]{\sum_{i=1}^nw_ix_i^{-p}}=\sqrt[p]{\frac{1}{\sum_{i=1}^nw_i\frac{1}{x_i^p}}}\geq \sqrt[q]{\frac{1}{\sum_{i=1}^nw_i\frac{1}{x_i^q}}}=\sqrt[-q]{\sum_{i=1}^nw_ix_i^{-q}},
我们得到了关于 -p 与 -q 的幂平均不等式,同样的推理可以倒推,从而证明了两个不等式等价,这在后面的证明中将用到。

几何平均
对任何 q,指数为 q 的幂平均与几何平均之间的不等式为:
:\prod_{i=1}^nx_i^{w_i} \leq \sqrt[q]{\sum_{i=1}^nw_ix_i^q}
:\sqrt[q]{\sum_{i=1}^nw_ix_i^q}\leq \prod_{i=1}^nx_i^{w_i}
(第一个不等式对正数 q,第二个对负数)

我们在两边取 q 次幂:
:\prod_{i=1}^nx_i^{w_i\cdot q} \leq \sum_{i=1}^nw_ix_i^q
两种情形我们都得到关于 x_i^q 的加权算术几何平均不等式,这可以用琴生不等式证明,利用对数函数是凸函数的事实:
:\sum_{i=1}^nw_i\log(x_i) \leq \log\left(\sum_{i=1}^nw_ix_i\right)
:\log\left(\prod_{i=1}^nx_i^{w_i}\right) \leq \log\left(\sum_{i=1}^nw_ix_i\right)
两边取指数函数(严格递增),我们得到了不等式:
:\prod_{i=1}^nx_i^{w_i} \leq \sum_{i=1}^nw_ix_i.

从而对任何正数 q,下式成立:
:\sqrt[-q]{\sum_{i=1}^nw_ix_i^{-q}}\leq \prod_{i=1}^nx_i^{w_i} \leq \sqrt[q]{\sum_{i=1}^nw_ix_i^q}.
因为此不等式对任何 q 成立,足够小同样成立,可以将证明(利用洛必达法则),当 q 趋于 0 时,左右两边趋于几何平均,q 趋于 0 时的幂平均是几何平均:
:\lim_{q\rightarrow 0}\sqrt[q]{\sum_{i=1}^nw_ix_i^{q}}=\prod_{i=1}^nx_i^{w_i}

幂平均不等式
我们将证明对任何 p 如下不等式成立:
:\sqrt[p]{\sum_{i=1}^nw_ix_i^p}\leq \sqrt[q]{\sum_{i=1}^nw_ix_i^q}.
如果 p 是负数且 q 是正数,不等式等价于上面已证过的
:\sqrt[p]{\sum_{i=1}^nw_ix_i^p}\leq \prod_{i=1}^nx_i^{w_i} \leq\sqrt[q]{\sum_{i=1}^nw_ix_i^q}
对正数 pq 的证明如下:定义函数 f:{\mathbb R_+}\rightarrow{\mathbb R_+}, f(x)=x^{\frac{q}{p}}。 f是一个幂函数,所以有二阶导数:f(x)=\left(\frac{q}{p}\right)\left(\frac{q}{p}-1\right)x^{\frac{q}{p}-2},在 f 的定义域内严格正,因为 q>p,从而我们知道 f* 是凸的。

利用这一点以及琴生不等式,我们得到:
:f(\sum_{i=1}^nw_ix_i^p)\leq\sum_{i=1}^nw_if(x_i^p)
:\sqrt[\frac{p}{q}]{\sum_{i=1}^nw_ix_i^p}\leq\sum_{i=1}^nw_ix_i^q
两边取 \frac{1}{q} 次幂(递增函数,因 \frac{1}{q} 为正数)我们得到了欲证之不等式:

:\sqrt[p]{\sum_{i=1}^nw_ix_i^p}\leq\sqrt[q]{\sum_{i=1}^nw_ix_i^q}。

最后使用先前证过的等价性,我们得到了关于负数 pq 的不等式,证毕。

最小值与最大值
此段最后将证明当指数 p 趋于 -\infty 与 +\infty,其幂平均的幂平均分别趋于最小值与最大值。定义指数为 -\infty 与 +\infty 的幂平均为最大值与最小值。从而应该有:
:\min (x_1,x_2,\ldots ,x_n)\leq \sqrt[q]{\sum_{i=1}^nw_ix_i^q}\leq \max (x_1,x_2,\ldots ,x_n)
对最大值证明如下:不失一般性假设序列 x_i 非减且全不为零。则不等式等价于:
:\sqrt[q]{\sum_{i=1}^nw_ix_i^q}\leq x_1
两边取 q 次幂,我们得到不等式(取决于 q 的符号):
:\sum_{i=1}^nw_ix_i^q\leq {\color{red} \geq} x_1^q
若 q>0 为 ≤, 若 q 为 ≥。

两边同时减去 w_1x_1^q 我们得到:
:\sum_{i=2}^nw_ix_i^q\leq {\color{red} \geq} (1-w_1)x_1^q
除以 (1-w_1):
:\sum_{i=2}^n\frac{w_i}{(1-w_1)}x_i^q\leq {\color{red} \geq} x_1^q
1-w_1 不为零,从而:
:\sum_{i=2}^n\frac{w_i}{(1-w_1)}=1
减去 {x_1}^q 剩下:
:\sum_{i=2}^n\frac{w_i}{(1-w_1)}(x_i^q-x_1^q)\leq {\color{red} \geq} 0
这是显然的,因为 x_1 大于或等于任何 x_i,从而
:x_i^q-x_1^q\leq {\color{red} \geq} 0

对最小值证明几乎相同,只不过将 x_1、w_1 换作 x_n、w_n,证毕。

另一方面,当 q 大于零时,由简单的推理以及上面的不等式有
:\sqrt[q]{w_1x_1^q}
令 q 趋于 +\infty 时,左边同样趋于 x_1,由夹逼定理知中间项幂平均趋于 x_1。最小值的证明完全类似。

广义 f-平均
幂平均可以推广到更一般的:
: M_f(x_1,\dots,x_n) = f^{-1}
\left[{\frac{1}{n}\cdot\sum_{i=1}^n{f(x_i)}}\right]

例如这包括了几何平均而勿需使用极限。幂平均是由 f\left(x\right)=x^p 得到的。

应用
信号处理
幂平均作为一个非线性移动平均。对於小 p 值,幂平均比较偏重小信号值,对於大 p 值,幂平均则会强调大信号值。给予一个高效率移动算术平均的实施函数,称为 smooth ,工程师可以按照下述 Haskell 代码,设计一个移动幂平均实施函数:
powerSmooth :: Floating a => ([a] -> [a]) -> a -> [a] -> [a]
powerSmooth smooth p =
map ( recip p) . smooth . map (p)

  • 对於大 p 值,这可作为一个整流信号的包封檢測器()。
  • 对於小 p 值,这可作为一个質量譜的基线侦测器()。

参见条目

  • 平均
  • 平方平均
  • 算术平均
  • 几何平均
  • 调和平均
  • 海伦平均(Heronian mean)
  • 莱默平均(Lehmer mean)——也是一个与幂有关的平均数。
  • 算术几何平均不等式

外部链接
*[http://mathworld.wolfram.com/PowerMean.html Power mean at MathWorld]
*[http://people.revoledu.com/kardi/tutorial/BasicMath/Average/Generalized%20mean.html Examples of Generalized Mean]
*[http://planetmath.org/?op=getobj&from=objects&name=ProofOfGeneralMeansInequality A proof of the Generalized Mean] on PlanetMath
*[http://w3.math.sinica.edu.tw/math_media/d284/28406.pdf 平均論]

评论 (0)

  • 还没有评论,来抢沙发吧。