幂平均

来自testwiki
跳转到导航 跳转到搜索

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

數個不同冪平均 Mp(1,x) 的圖表。

定义

p 是一非零实数,可定义非負实数 x1,,xn p次幂平均

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

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

M0(x1,,xn)=(i=1nxi)1n

性质

  • 和所有平均一样,幂平均是各参数 x1,,xn 的一次齐次函数。即若 b 是一个正实数,则 bx1,,bxn 指数为 p 的幂平均等于 bx1,,xn 的幂平均。
  • 几何算术平均一样,这种平均的计算可以分解成同样大小的子块来计算。
Mp(x1,,xnk)=Mp(Mp(x1,,xk),Mp(xk+1,,x2k),,Mp(x(n1)k+1,,xnk))

幂平均不等式

一般地,如果 p<q,则 Mp(x1,,xn)Mq(x1,,xn) 且这两个平均相等当且仅当 x1=x2==xn。这由事实

p Mp(x1,,xn)p0,

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

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

特例

特例n=2时的图形描述。
M(x1,,xn)=limpMp(x1,,xn)=min{x1,,xn} 最小值
M1(x1,,xn)=n1x1++1xn 调和平均
M0(x1,,xn)=limp0Mp(x1,,xn)=x1xnn 几何平均
M1(x1,,xn)=x1++xnn 算术平均
M2(x1,,xn)=x12++xn2n 平方平均Template:Anchor
M3(x1,,xn)=x13++xn3n3 立方平均
M+(x1,,xn)=limpMp(x1,,xn)=max{x1,,xn} 最大值

幂平均不等式的证明

不同符号的不等式之等价

假设指数 pq 的幂平均间有不等式:

i=1nwixippi=1nwixiqq

i=1nwixippi=1nwixiqq.

我们在两边取倒数(正实数上的严格递减函数,不等号反向):

i=1nwixipp=1i=1nwi1xipp1i=1nwi1xiqq=i=1nwixiqq,

我们得到了关于 pq 的幂平均不等式,同样的推理可以倒推,从而证明了两个不等式等价,这在后面的证明中将用到。

几何平均

对任何 q,指数为 q 的幂平均与几何平均之间的不等式为:

i=1nxiwii=1nwixiqq
i=1nwixiqqi=1nxiwi

(第一个不等式对正数 q,第二个对负数)

我们在两边取 q 次幂:

i=1nxiwiqi=1nwixiq

两种情形我们都得到关于 xiq 的加权算术几何平均不等式,这可以用琴生不等式证明,利用对数函数是凸函数的事实:

i=1nwilog(xi)log(i=1nwixi)
log(i=1nxiwi)log(i=1nwixi)

两边取指数函数(严格递增),我们得到了不等式:

i=1nxiwii=1nwixi.

从而对任何正数 q,下式成立:

i=1nwixiqqi=1nxiwii=1nwixiqq.

因为此不等式对任何 q 成立,足够小同样成立,可以将证明(利用洛必达法则),当 q 趋于 0 时,左右两边趋于几何平均,q 趋于 0 时的幂平均是几何平均:

limq0i=1nwixiqq=i=1nxiwi

幂平均不等式

我们将证明对任何 p<q 如下不等式成立:

i=1nwixippi=1nwixiqq.

如果 p 是负数且 q 是正数,不等式等价于上面已证过的

i=1nwixippi=1nxiwii=1nwixiqq

对正数 pq 的证明如下:定义函数 f:++f(x)=xqpf是一个幂函数,所以有二阶导数:f(x)=(qp)(qp1)xqp2,在 f 的定义域内严格正,因为 q>p,从而我们知道 f 是凸的。

利用这一点以及琴生不等式,我们得到:

f(i=1nwixip)i=1nwif(xip)
i=1nwixippqi=1nwixiq

两边取 1q 次幂(递增函数,因 1q 为正数)我们得到了欲证之不等式:

i=1nwixippi=1nwixiqq

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

最小值与最大值

此段最后将证明当指数 p 趋于 +,其幂平均的幂平均分别趋于最小值与最大值。定义指数为 + 的幂平均为最大值与最小值。从而应该有:

min(x1,x2,,xn)i=1nwixiqqmax(x1,x2,,xn)

对最大值证明如下:不失一般性假设序列 xi 非减且全不为零。则不等式等价于:

i=1nwixiqqx1

两边取 q 次幂,我们得到不等式(取决于 q 的符号):

i=1nwixiqx1q

q>0 为 ≤,q<0 为 ≥

两边同时减去 w1x1q 我们得到:

i=2nwixiq(1w1)x1q

除以 (1w1)

i=2nwi(1w1)xiqx1q

1w1 不为零,从而:

i=2nwi(1w1)=1

减去 x1q 剩下:

i=2nwi(1w1)(xiqx1q)0

这是显然的,因为 x1 大于或等于任何 xi,从而

xiqx1q0

对最小值证明几乎相同,只不过将 x1w1 换作 xnwn,证毕。

另一方面,当 q 大于零时,由简单的推理以及上面的不等式有

w1x1qq<i=1nwixiqqx1,

q 趋于 + 时,左边同样趋于 x1,由夹逼定理知中间项幂平均趋于 x1。最小值的证明完全类似。

广义 f-平均

Template:Main

幂平均可以推广到更一般的Template:Link-en

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

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

应用

信号处理

幂平均作为一个非线性移动平均。对於小 p 值,幂平均比较偏重小信号值,对於大 p 值,幂平均则会强调大信号值。给予一个高效率移动算术平均的实施函数,称为 smooth ,工程师可以按照下述 Haskell 代码,设计一个移动幂平均实施函数:

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

参见条目

外部链接