首页 > 其他分享 >康托洛维奇不等式

康托洛维奇不等式

时间:2023-03-24 16:00:14浏览次数:41  
标签:dots xi frac limits 不等式 洛维奇 mathop 康托 lambda

康托洛维奇不等式是数值优化中收敛性分析的一个常用工具:

康托洛维奇不等式:设\(Q\)为正定对称阵,\(x \in \mathbb{R}^n\),则有

\[\frac{(x^Tx)^2}{(x^TQx)(x^TQ^{-1}x)}\geq\frac{4aA}{(a+A)^2} \]

其中\(a,A\)分别为\(Q\)的最小和最大特征值。

证明:设\(Q\)的特征值位\(\lambda_1,\lambda_2,\dots,\lambda_n\),且$$0<a\leq \lambda_1\leq \lambda_2 \leq \dots \leq \lambda_n=A$$

设\(D^TQD=diag\{\lambda_1,\lambda_2,\dots,\lambda_n\},D^TD=I_n,x=Dy,y=(y_1,y_2,\dots,y_n)^T\),
于是:

\[\begin{align*} \frac{(x^Tx)^2}{(x^TQx)(x^TQ^{-1}x)} & =\frac{(y^TD^TDy)^2}{(y^TD^TQDy)((y^TD^TQ^{-1}Dy)}\\ & =\frac{(\sum_{i=1}^n y_i^2)^2}{(\sum_{i=1}^n (\lambda_iy_i^2))(\sum_{i=1}^n (y_i^2/\lambda_i))}\\ & = \frac{1/(\sum_{i=1}^n \xi_i \lambda_i)}{\sum_{i=1}^n (\xi_i/\lambda_i)} \equiv \frac{\phi(\xi)}{\psi(\xi)} \end{align*}\]

其中\(\xi_i=y_i^2 / \sum_{i=1}^n y_i^2\)
已知\(0<a\leq \lambda_1\leq \lambda_2 \leq \dots \leq \lambda_n=A\),
令\(\xi=(\xi_1,\xi_2,\dots,\xi_n)^T\),\(\Lambda=(\lambda_1,\lambda_2,\dots,\lambda_n)^T,\mathop{\Lambda}\limits^{-}=(\frac{1}{\lambda_1},\frac{1}{\lambda_2},\dots,\frac{1}{\lambda_n})^T\)
于是问题就变成了

\[\begin{align*} \mathop{min}\limits_{\xi} & \quad \frac{1/(\Lambda^T\xi)}{(\mathop{\Lambda}\limits^{-})^T\xi}\\ s.t. & \quad \xi_1+\xi_2+\dots+\xi_n=1 \end{align*}\]

再设\(\forall \lambda \in [\lambda_1,\lambda_n],\Gamma(\lambda)=\{\xi : \Lambda^T\xi=\lambda,\xi_1+\xi_2+\dots+\xi_n=1\}\)
于是问题可以表示为:

\[\mathop{min}\limits_{\lambda \in [\lambda_1,\lambda_n]} \mathop{min}\limits_{\xi \in \Gamma(\lambda)} \frac{1/(\Lambda^T\xi)}{(\mathop{\Lambda}\limits^{-})^T\xi} \]

如果先固定\(\lambda \in [\lambda_1,\lambda_n]\),则
\(\mathop{min}\limits_{\xi \in \Gamma(\lambda)} \frac{1/(\Lambda^T\xi)}{(\mathop{\Lambda}\limits^{-})^T\xi} = \frac{1}{(\mathop{\Lambda}\limits^{-})^T\xi}\iff \mathop{max}\limits_{\xi \in \Gamma(\lambda)}(\mathop{\Lambda}\limits^{-})^T\xi=\xi_1\frac{1}{\lambda_2}+\xi_2\frac{1}{\lambda_2}+\dots+\xi_n\frac{1}{\lambda_n}\)
进一步发现:\(\xi_1\frac{1}{\lambda_2}+\xi_2\frac{1}{\lambda_2}+\dots+\xi_n\frac{1}{\lambda_n} \leq \frac{\xi_1+\xi_2+\dots+\xi_{n-1}}{\lambda_1}+\frac{\xi_n}{\lambda_n}=\frac{1-\xi_n}{\lambda_1}+\frac{\xi_n}{\lambda_n}\)
而上述不等式的等号可以取到,事实上,只需要使得:"\(\xi_1=\xi_n=1,\xi_2=\xi_3=\dots=\xi_{n-1}=0,\xi_1\lambda_1+\xi_n\lambda_n=\lambda\)"即可.
于是原问题相当于

\[\mathop{min}\limits_{\lambda \in [\lambda_1,\lambda_n],\xi_1+\xi_n=1} \quad \frac{1/\lambda}{\xi_1/\lambda_1+\xi_n/\lambda_n} \]

\[\frac{1/\lambda}{\xi_1/\lambda_1+\xi_n/\lambda_n}=\frac{1/\lambda}{(\lambda_1+\lambda_n-\lambda)/(\lambda_1\lambda_n)}=\frac{\lambda_1\lambda_n}{\lambda(\lambda_1+\lambda_n-\lambda) }\geq \frac{4\lambda_1\lambda_n}{(\lambda_1+\lambda_n)^2} \]

上述等号在\(\lambda=\frac{1}{2}(\lambda_1+\lambda_n)\)时成立.
由于\(\lambda_1=a,\lambda_n=A\),所以得到康托洛维奇不等式

\[\frac{(x^Tx)^2}{(x^TQx)(x^TQ^{-1}x)}\geq\frac{4aA}{(a+A)^2} \]

同时本文也是对叶荫宇老师的《线性与非线性规划》一书8.2节对康托洛维奇不等式证明的补充:
image
image

标签:dots,xi,frac,limits,不等式,洛维奇,mathop,康托,lambda
From: https://www.cnblogs.com/wjma2719/p/17251945.html

相关文章

  • 1499. 满足不等式的最大值
    题目描述数组points在x轴上是严格单调增,需要求一个不等式x1+y2+x2-x1的最大值?要求是x2-x1不能超过kf1-分析不等式+单调队列基本分析怎么能让值最大?对当前x2和y2......
  • 康托展开式与逆康托展开式
    康托展开官方简介:康托展开是一个全排列到一个自然数的双射,常用于构建哈希表时的空间压缩。康托展开的实质是计算当前排列在所有由小到大全排列中的顺序,因此是可逆的。......
  • 均值不等式学习笔记
    从平均数说起我们都知道\(n\)个数的平均数表示为:\[\frac{a_1+a_2+a_3+\cdotsa_n}{n}\]这种最常见的平均数被称为“算术平均数”(ArithmeticMean)。还有一种常用的平均......
  • 1235. 付账问题(均值不等式贪心)
    https://www.acwing.com/problem/content/1237/均值不等式贪心#include<iostream>#include<algorithm>#include<cmath>usingnamespacestd;constintN=5e5+10;lo......
  • 决策单调性与四边形不等式
    参考自:《决策单调性与四边形不等式》,彭思进,感谢Itst的耐心答疑/bx《决策单调性与四边形不等式-学习笔记》,p_b_p_b的博客一、决策单调性与四边形不等式定义1.......
  • 康托展开
    康托展开名词解释:康托展开是一个全排列到一个自然数的双射,常用于构建哈希表时的空间压缩。---康托展开的实质是计算当前排列在所有由小到大全排列中的顺序,因此是可逆的......
  • 柯西不等式也是重要的消元方式(涉及2次放缩)
    已知\(a\),\(b\in\textbf{R}\),函数\(f(x)=\text{e}^x-a\sinx\),\(g(x)=b\sqrtx\).若\(y=f(x)\)和\(y=g(x)\)有公共点.求证:\(a^2+b^2>\text{e}\).分析:\(\text{e}^x-a\sinx=......
  • 康托展开和康托逆展开
    #include<bits/stdc++.h>usingnamespacestd;intf[20];intmain(){ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);f[0]=1;for(inti=1;......
  • 算法学习笔记(54)——绝对值不等式
    绝对值不等式题目链接:AcWing104.货仓选址\[\begin{align*}f(x)&=\lvertx_1-x\rvert+\lvertx_2-x\rvert+\cdots+\lvertx_n-x\rvert\\&=(\lve......
  • 算法学习笔记(53)——排序不等式
    排序不等式题目链接:AcWing913.排队打水让最磨叽的人最后打水。如图所示,第一个同学被等了6次,第二个同学被等了5次,以此类推...\[总时间=t_1\times(n-1)+t_2\t......