Farkasの補題の証明をわかりやすく解説

Farkasの補題を証明します。これは線形計画問題や最適化問題で重要な補題で、凸解析の分離定理を用いて証明します。連立方程式の解の存在と、連立不等式の解の存在のあいだの関係に示唆を与える補題でもあります。

Farkasの補題の証明をわかりやすく解説

Farkasの補題を証明するために、有限次元の凸解析における基本的な命題を復習しておきましょう。
コーシー列・収束列・下限・凸集合・閉包などの実数論あるいは集合位相の標準的な内容を分かっていれば理解することができます。

予備知識1:有限次元射影定理

定理:有限次元射影定理(原点版)

KRd を閉凸集合とする。
xRd
||x||infxK||x||
を満たすものが存在して、一意である。つまり、
infxK||x||
はユニークなミニマイザーをもつ。

証明:inf の定義からxnK
infxK||x||||xn||infxK||x||+1n
を満たすものがとれる。

12xn+12xmK なので、
infxK||x||||12xn+12xm||
であることから、
4(infxK||x||)2||xn+xm||2
である。

中線定理
||xnxm||2+||xn+xm||2=2||xn||2+2||xm||2
を思い出すと、
||xnxm||2=2||xn||2+2||xm||2||xn+xm||22||xn||2+2||xm||24(infxK||x||)2
であるので、極限をとるとlimn||xn||=infxK||x|| であったことを思い出すと、
||xnxm||20
となるのでコーシー列である。

従って、Rn が完備なのでxn は収束列であるので、収束先x が存在する。
K が閉集合であることから、xK である。このx が求めているミニマイザーである。

一意性を示しておく。もし異なるミニマイザーが存在すると、
||12(x+y))||2=12||x||2+12||y||2||xy||2
となり、||x||=||y||,||xy||>0 であることから、
||12(x+y))||<||x||
となってしまい、x がミニマイザーでなくなってしまうので、ミニマイザーは一意である。

補足:原点0p にくるように図形の位置関係を平行移動することで、以下の有名な定理が直ちに従う。

定理:有限次元射影定理

KRd を閉凸集合とする。
任意のpRd に対して、xRd
||px||infxK||px||
を満たすものが存在して、一意である。

予備知識2:有限次元の分離超平面定理の一種

分離超平面の存在に関する命題はさまざまなバージョンが存在するが、次のものを用いることにする。

定理:有限次元の分離超平面定理の一種

pRn とし、K を閉包がpを含まない凸集合とする。
このとき、ξRn で任意のkK に対して
(ξ,p)<(ξ,k)
を満たすものが存在する。

K=Kp={kpkK} と表すことにする。
K の閉包は閉凸集合であるので、
ξKs.t.||ξ||=infkK||k||
が、有限次元射影定理から存在することがわかる。任意にv=kpK をとる。
ξ,vの凸結合(1t)ξ+tvK を考える0t1
(1t)ξ+tv=ξ+t(vξ) であるので、
||ξ||2||ξ+t(vξ)||2=||ξ||2+2t(ξ,vξ)+t2||vξ||2=||ξ||2+2t(ξ,v)2t||ξ||2+t2||vξ||2
が得られる。従って、両辺から||v||2を引くと
02t(ξ,v)2t||ξ||2+t2||vξ||2
が任意の0t1に対して成り立つ。tで割ると(t=0の時は次の式は自明)、
02(ξ,v)2||ξ||2+t||vξ||2
が得られる。t0と極限をとることで、
02(ξ,v)2||ξ||2
が得られる。v=kp であったので、
(ξ,p)+||ξ||2(ξ,k)
が成り立つ。

ここで、p が$K$の閉包に含まれないことから、原点がK=Kpの閉包に含まれないので、
||ξ||>0 が成り立っている。従って、
(ξ,p)<(ξ,p)+||ξ||2(ξ,k)
となり、証明が終了する。

Farkasの補題の証明

補題:Farkasの補題

An×m 行列とし、b を$n$ 次列ベクトルとする。このとき、
(1)非負m次列ベクトルx0Ax=b を満たすものが存在する。
(2)任意のn次列ベクトルyに対して、ytA0 ならばytb0
は同値である。

証明:

(1)ならば(2)を示す。非負$m$次列ベクトルx0Ax=b を満たすものをとる。
ytb=ytAx0
より主張が従う。ただし最後の不等号は、ytA,x が共に非負であることから成り立つ。

(2)ならば(1)を示す。Aの列ベクトルをa1,,am と表記することにする。
(背理法)(1)が成り立たないとすると、ba1,,amの凸錐結合で表せない。即ち、
b{i=1mxiaix1,,xm0}
である。

a1,,amの凸錐結合全体(つまり閉凸錐)はb を含まないので、
有限次元の超平面分離定理から適当なyで任意のa{i=1mxiaix1,,xm0} に対して、
(y,b)<(y,a1),(y,am)
を満たすものが存在する。

特にa としてa=0,a1,a2,,am を考えることで
(y,b)<(y,0)=0
が得られる。また、もしa1,,amの凸錐結合a
(y,a)<0
であるものが存在してしまうと、aR>0倍して、
(y,Ra)<(y,b)
となるように十分大きなR がとれてしまうが、Raも凸錐結合で表されるのでこれはあり得ないので、任意の凸錐結合aに対して
0(y,a)
である。

このことはすなわち、ytA=((y,a1),,(y,am))であることに注意すると、
(y,b)<0,0ytA
を満たすy が存在することを意味している。これは、
(2)任意のn次列ベクトルyに対して、ytA0 ならばytb0
と矛盾する。

以上で証明を終了する。

補足:Farkasの補題は択一定理(Alternative)的な書き方を以下のようにすることができる。

補題:Farkasの補題の択一定理的な書き方

An×m 行列とし、b を$n$ 次列ベクトルとする。このとき、
(1)非負m次列ベクトルx0Ax=b を満たすものが存在する。
(2)n次列ベクトルyで、(y,b)かつ0ytAであるものが存在する。
のいずれか一方のみが成り立つ。

あわせて読みたい記事

記事をシェアして話のネタにする

コメント

コメントする