有界區(qū)域上二階線性橢圓方程解的存在性

感覺學(xué)得十分混亂, 所以有必要寫一篇筆記來捋一捋自己的思路, 盡量做到以簡(jiǎn)馭繁. 這篇筆記只考慮實(shí)系數(shù)方程.

我們考慮的區(qū)域\mathbb{R}^d中的一個(gè)有界區(qū)域\Omega, 即連通開集. 現(xiàn)在暫時(shí)不要求邊界光滑, 所以正方體這種也在我們的考慮之列.

我們考慮一個(gè)映射L:{\cal D}'(\Omega)\rightarrow{\cal D}'(\Omega), u\mapsto Lu=\sum_{i,j}a_{ij}u_{x_ix_j}+\sum_ib_iu_{x_i}+cu, 其中實(shí)函數(shù)a_{ij},b_i,c\in C^\infty(\Omega). 系數(shù)的假定可能有些過強(qiáng), 但是為了語(yǔ)言上的方便我愿意做出這種犧牲, 尤其是為了使我們能夠把L理解成這種映射的形式而不必去討論它的定義域. 另外對(duì)我而言, 在應(yīng)用中我也尚未遇到過不光滑的系數(shù), 所以這樣得到的結(jié)論暫時(shí)是完全夠用的.

我們還要假定L是一致橢圓的, 即\exists\theta>0, 使得\forall x\in\Omega, \forall\xi\in\mathbb{R}^n, 有\sum_{i,j}a_{ij}(x)\xi_i\xi_j\ge\theta|\xi|^2.

先說一下什么叫解. 我們先引入一個(gè)最廣義的解的定義, 我們稱之為分布解. 設(shè)f\in{\cal D}'(\Omega), 如果Lu=f(注意這個(gè)等號(hào)是{\cal D}'(\Omega)中的等號(hào)), 那么我們稱u是一個(gè)分布解.

對(duì)方程\Delta u=0而言, 它是泛函J(u)=\int_\Omega|\nabla u|^2dx的Euler-Lagrange方程. 出于這個(gè)物理上的考慮, 我們希望所尋找的解應(yīng)該至少落在H^1(\Omega)里.

注意到如果u\in H^m(\Omega), 由于L是二階的微分算子, 所以有Lu\in H^{m-2}(\Omega). 如果u\in H^1(\Omega), 那么Lu\in H^{-1}(\Omega), 所以如果Lu=f\mathcal{D}'(\Omega)中成立的話, 我們至少應(yīng)該假定f\in H^{-1}(\Omega).

現(xiàn)在我們來考慮邊界條件. 如果邊界條件是u|_{\partial\Omega}=0(即Dirichlet問題), 那么我們簡(jiǎn)單要求u\in H_0^1(\Omega)即可. 其他邊界條件暫時(shí)先不考慮.

把上面的討論總結(jié)起來, 就是

定義1(弱解). 設(shè)a_{ij},b_i,c\in C^\infty(\overline{\Omega}), L:\mathcal{D}'(\Omega)\rightarrow\mathcal{D}'(\Omega)如前定義, f\in H^{-1}(\Omega). 如果u\in H_0^1(\Omega)滿足Lu=f(\mathcal{D}'(\Omega)中的等號(hào)), 那么稱u是問題
\left\{\begin{aligned} &Lu=f\\ &u|_{\partial\Omega}=0 \end{aligned}\right.
的弱解.

注意到等式兩端都是H^{-1}(\Omega)內(nèi)的元素, 故經(jīng)過極限步驟, 我們知道試驗(yàn)函數(shù)可以取在H_0^1(\Omega)內(nèi), 這樣就得到和Evans一模一樣的定義.

顯然L|_{H_0^1}(為方便我們也直接記作L, 下文中的L們到底是L還是L|_{H_0^1}, 讀者可以聯(lián)系上下文選擇合理的那個(gè), 我每寫下一個(gè)L時(shí)總會(huì)斟酌使得總有一個(gè)是合理的)是H_0^1H^{-1}的有界算子. 為了說明對(duì)任何f\in H^{-1}, 都存在u\in H_0^1滿足Lu=f, 我們只需要說明L是滿的.

為了說明L是滿的, 我們需要Lax-Milgram定理.

定理2(實(shí)的Lax-Milgram定理).設(shè)H\mathbb{R}上的Hilbert空間, BH\times H上的雙線性形式, 并且存在\alpha,\beta>0使得\forall u,v\in H, \left|B[u,v]\right| \le\alpha\|u\|\|v\|, |B[u,u]|\ge\beta\|u\|^2, 則對(duì)任何f\in H^*, 存在唯一的u\in H使得\forall v, f(v)=B[u,v].
證明.Riesz表示定理和Lax-Milgram定理. \blacksquare

在我們面對(duì)的情形里, 我們考慮的Hilbert空間為H_0^1, 雙線性形式為B[u,v]=\langle -Lu,v\rangle. 因?yàn)?img class="math-inline" src="https://math.jianshu.com/math?formula=u%5Cin%20H_0%5E1" alt="u\in H_0^1" mathimg="1">時(shí), Lu\in H^{-1}, 所以這個(gè)配對(duì)是沒問題的. 只要我們能驗(yàn)證|\langle Lu,v\rangle|\le\alpha\|u\|_{H_0^1}\|v\|_{H_0^1}, |\langle -Lu,u\rangle|\ge\beta\|u\|_{H_0^1}, 那么由上述定理, 對(duì)任何f\in H^{-1}, 都存在唯一的u\in H_0^1使得Lu=f.

現(xiàn)在我們就來驗(yàn)證這個(gè)雙線性形式的連續(xù)性和強(qiáng)制性是否成立. 連續(xù)性由\|Lu\|_{H^{-1}}\lesssim_{a_{ij},b_i,c}\|u\|_{H_0^1}得到, 關(guān)于強(qiáng)制性, 我們有:

定理3(G?rding不等式). 存在C=C(L)>0使得\forall u\in H_0^1, 有
\langle -Lu,u\rangle\ge\frac{\theta}{2}\|u\|_{H_0^1}^2-C\|u\|_{L^2}^2
這里的\thetaL的一致橢圓性的那個(gè)系數(shù).

這看起來離強(qiáng)制性還差一點(diǎn). 但是我們注意到當(dāng)\lambda>C(L)時(shí), 有\langle(-L+\lambda)u,u\rangle\ge\frac{\theta}{2}\|u\|_{H_0^1}^2, 即-L+\lambda伴隨的雙線性形式是有強(qiáng)制性的, 由此我們立刻得到如下定理:

定理4(第一存在性定理). 當(dāng)\lambda>C(\lambda)時(shí)(這里的C(\lambda)是G?rding不等式中所定義的), 對(duì)任何f\in H^{-1}, 都存在著唯一的u\in H_0^1使得(L-\lambda)u=f.

這個(gè)定理說明L-\lambda:H_0^1\rightarrow H^{-1}是滿的, 另外注意到L-\lambda是有界的, 并且再由\langle(-L+\lambda)u,u\rangle\ge\frac{\theta}{2}\|u\|_{H_0^1}^2得到\|(L-\lambda)u\|_{H_{-1}}\ge\frac{\theta}{2}\|u\|_{H_0^1}, 從而L-\lambda是兩個(gè)Banach空間之間的線性同胚.

第一存在性定理看起來不錯(cuò), 但它沒有回答Lu=f的可解性. 為了回答這個(gè)問題, 我們接下來將會(huì)假設(shè)f\in L^2, 而不再是原先的f\in H^{-1}(因?yàn)闆]有能力做到H^{-1}里的可解性).

在此之前, 我們還要引入L的形式伴隨L^*:\mathcal{D}'(\Omega) \rightarrow\mathcal{D}'(\Omega), 定義為L^*u=\sum_{i,j}(a_{ij}v)_{x_ix_j}-\sum_i(b_iv)_{x_i}+cv. 分部積分告訴我們, 當(dāng)u,v一個(gè)屬于\mathcal{D}'(\Omega)一個(gè)屬于C_c^\infty(\Omega)時(shí)(無所謂順序), \langle Lu,v\rangle=\langle u,L^*v\rangle. 當(dāng)u,v\in H_0^1時(shí), 可以取\{v_n\}\subset H_0^1且在H_0^1v_n\rightarrow v. 來看\langle Lu,v_n\rangle=\langle u,L^*v_n\rangle, 該式左邊由于Lu\in H^{-1}可以知道是趨于\langle Lu,v\rangle的, 右邊由于L^*H_0^1H^{-1}的有界線性算子知也趨于\langle u,L^*v\rangle, 故\langle Lu,v\rangle=\langle u,L^*v\rangle.

我們以下總假設(shè)\lambda>\max(C(L),C(L^*)), 那么此時(shí)L-\lambda,L^*-\lambdaH_0^1H^{-1}的線性同胚, (L-\lambda)^{-1},(L^*-\lambda)^{-1}H^{-1}H_0^1的線性同胚, 把它們限制到L^2\subset H^{-1}上, 并且和i:H_0^1\hookrightarrow L^2復(fù)合起來, 得到i(L-\lambda)^{-1}|_{L^2},i(L^*-\lambda)^{-1}|_{L^2}, 把這兩個(gè)算子分別記作M,M^*: L^2\rightarrow L^2. 由緊嵌入定理, 當(dāng)\partial\Omega\in C^1時(shí), i是緊的, 從而M,M^*是緊的.

現(xiàn)在我們把Lu=f改寫為(L-\lambda)u+\lambda u=f, 兩邊作用(L-\lambda)^{-1}得到u+\lambda Mu=Mf.

引理5. 設(shè)f\in L^2. 以下陳述是等價(jià)的:
(1)存在u\in H_0^1使得Lu=f;
(2)存在u\in L^2使得u+\lambda Mu=Mf.
證明. (1)推(2): 把Lu=f改寫為(L-\lambda)u+\lambda u=f, 注意兩邊都屬于H^{-1}, 故可以兩邊作用(L-\lambda)^{-1}得到u+\lambda Mu=Mf.
(2)推(1): 由條件有u+\lambda(L-\lambda)^{-1}u=(L-\lambda)^{-1}f, 我們很想兩邊作用L-\lambda, 為此得先檢查是不是每個(gè)元素都屬于H_0^1. 顯然\lambda Mu,Mf\in H_0^1, 故u\in H_0^1, 這樣我們可以放心作用L-\lambda得到Lu=f. \blacksquare

更一般地, 我們有

引理6. 設(shè)f\in L^2, 則:
(1)u\in H_0^1,Lu=f\Leftrightarrow u\in L^2, u+\lambda Mu=Mf;
(2)u\in H_0^1,Lu=f\Leftrightarrow 0\in L^2, u+\lambda Mu=0;
(3)u\in H_0^1,L^*u=f\Leftrightarrow u\in L^2, u+\lambda M^*u=M^*f;
(4)u\in H_0^1,L^*u=0\Leftrightarrow u\in L^2, u+\lambda M^*u=0;
證明. 和上一個(gè)引理的證明完全類似. \blacksquare

由于M,M^*是分別通過L,L^*來定義的, 現(xiàn)在我們要說明, 作為L^2L^2的有界線性算子, M^*恰恰是M的伴隨.

引理7. M^*M的伴隨, 即對(duì)任何f,g\in L^2, 有(Mf,g)_{L^2}=(f,M^*g)_{L^2}.
證明. (Mf,g)_{L^2}=((L-\lambda)^{-1}f,(L^*-\lambda)(L^*-\lambda)^{-1}g)_{L^2} =((L^*-\lambda)(L-\lambda)^{-1}f,(L^*-\lambda)^{-1}g)_{L^2}=(f,M^*g)_{L^2}.

這樣一來, 我們就把橢圓方程的可解性化為一個(gè)抽象的和緊算子有關(guān)的可解性問題. 為此我們需要以下定理:

定理8. 設(shè)H是實(shí)的Hilbert空間, K:H\rightarrow H是緊的, 那么:
(1)\operatorname{im}(I-K)=\ker(I-K^*)^\perp
(2)\dim\ker(I-K)=\dim\ker(I-K^*)<\infty.

由這個(gè)定理, 把-\lambda M視為K, 我們可以得知u+\lambda Mu=Mf有解u, 當(dāng)且僅當(dāng)\forall v滿足v+\lambda M^*v=0, 都有(Mf,v)_{L^2}=0. 我們把這句話翻譯一下, 就是Lu=f有解u, 當(dāng)且僅當(dāng)對(duì)任何v滿足L^*v=0, 均有0=(Mf,v)_{L^2}=(f,M^*v)_{L^2}=(f,-v/\lambda)_{L^2}, 或者說(f,v)_{L^2}=0.

總結(jié)起來, 我們有

定理9(第二存在性定理). 設(shè)\partial\Omega\in C^1(或者至少H_0^1(\Omega)\hookrightarrow L^2(\Omega)是緊的), f\in L^2, 則以下兩條等價(jià):
(1)Lu=f有解u\in H_0^1;
(2)對(duì)任何v\in H_0^1L^*v=0的弱解, 均有(f,v)_{L^2}=0.
并且如果Lu=0只有零解時(shí), 對(duì)每個(gè)f\in L^2, Lu=f均有唯一解.

最后我們來討論特征值問題, 即方程Lu=\lambda u, 邊界條件依然是u|_{\partial\Omega}=0. 先做一個(gè)簡(jiǎn)單的觀察: 設(shè)\lambda_0=C(L)如G?rding不等式中所定義, 這樣L-\lambda_0H_0^1H^{-1}的線性同胚,并且對(duì)任何v\in H_0^1, 有\langle -Lv+\lambda_0v,v\rangle\ge\frac{\theta}{2}\|v\|_{H_0^1}^2. 如果Lu=\lambda u, 那么(L-\lambda_0)u=(\lambda-\lambda_0)u, 從而-(\lambda-\lambda_0)\|u\|_{L^2}^2=-\langle(\lambda-\lambda_0)u,u\rangle\ge\frac{\theta}{2}\|u\|_{H_0^1}^2\ge0, 所以我們得知\lambda\le\lambda_0.

為了更進(jìn)一步, 我們還是得把方程改寫為和緊算子有關(guān)的形式: Lu=\lambda u\Leftrightarrow (L-\lambda_0)u=(\lambda-\lambda_0)u\Leftrightarrow u=(\lambda-\lambda_0)Mu, 其中M=(L-\lambda_0)^{-1}L^2上的緊算子. 從此可以看出, Lu=\lambda u有除了零解以外的解, 當(dāng)且僅當(dāng)1/(\lambda-\lambda_0)M的點(diǎn)譜.

那么我們就需要知道緊算子的譜是什么樣的. 事實(shí)上我們有:

定理10.設(shè)H是實(shí)Hilbert空間, \dim H=\infty, K\in B(H)是緊的. 那么
(1)0\in\sigma(K);
(2)\sigma(K)\setminus\{0\}\subset\sigma_p(K);
(3)\sigma(K)\setminus\{0\}要么是有限集, 要么是僅以0為聚點(diǎn)的可數(shù)集.

從這個(gè)定理我們知道, 如果記\Sigma=\{\lambda\in\mathbb{R}|\exists u\ne0, Lu-\lambda u=0\}, 那么\Sigma要么是有限集, 要么\{1/(\lambda-\lambda_0)|\lambda\in\Sigma\}是一列趨于0的實(shí)數(shù)列, 這說明\Sigma僅以-\infty為其聚點(diǎn)(注意\sup\Sigma\le\lambda_0!), 故\Sigma中的元素可以按順序排列為\{\lambda_i\}_{i=1}^\infty, 其中-\infty<\cdots<\lambda_2<\lambda_1\le\lambda_0.

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末锥咸,一起剝皮案震驚了整個(gè)濱河市驯遇,隨后出現(xiàn)的幾起案子绪囱,更是在濱河造成了極大的恐慌跋涣,老刑警劉巖讽营,帶你破解...
    沈念sama閱讀 217,734評(píng)論 6 505
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異阐污,居然都是意外死亡融欧,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,931評(píng)論 3 394
  • 文/潘曉璐 我一進(jìn)店門盯孙,熙熙樓的掌柜王于貴愁眉苦臉地迎上來鲁森,“玉大人,你說我怎么就攤上這事振惰「韪龋” “怎么了?”我有些...
    開封第一講書人閱讀 164,133評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵骑晶,是天一觀的道長(zhǎng)痛垛。 經(jīng)常有香客問我,道長(zhǎng)桶蛔,這世上最難降的妖魔是什么匙头? 我笑而不...
    開封第一講書人閱讀 58,532評(píng)論 1 293
  • 正文 為了忘掉前任,我火速辦了婚禮羽圃,結(jié)果婚禮上乾胶,老公的妹妹穿的比我還像新娘。我一直安慰自己朽寞,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,585評(píng)論 6 392
  • 文/花漫 我一把揭開白布斩郎。 她就那樣靜靜地躺著脑融,像睡著了一般。 火紅的嫁衣襯著肌膚如雪缩宜。 梳的紋絲不亂的頭發(fā)上肘迎,一...
    開封第一講書人閱讀 51,462評(píng)論 1 302
  • 那天甥温,我揣著相機(jī)與錄音,去河邊找鬼妓布。 笑死姻蚓,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的匣沼。 我是一名探鬼主播狰挡,決...
    沈念sama閱讀 40,262評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼释涛!你這毒婦竟也來了加叁?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,153評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤唇撬,失蹤者是張志新(化名)和其女友劉穎它匕,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體窖认,經(jīng)...
    沈念sama閱讀 45,587評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡豫柬,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,792評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了扑浸。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片轮傍。...
    茶點(diǎn)故事閱讀 39,919評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖首装,靈堂內(nèi)的尸體忽然破棺而出创夜,到底是詐尸還是另有隱情,我是刑警寧澤仙逻,帶...
    沈念sama閱讀 35,635評(píng)論 5 345
  • 正文 年R本政府宣布驰吓,位于F島的核電站,受9級(jí)特大地震影響系奉,放射性物質(zhì)發(fā)生泄漏檬贰。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,237評(píng)論 3 329
  • 文/蒙蒙 一缺亮、第九天 我趴在偏房一處隱蔽的房頂上張望翁涤。 院中可真熱鬧,春花似錦萌踱、人聲如沸葵礼。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,855評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)鸳粉。三九已至,卻和暖如春园担,著一層夾襖步出監(jiān)牢的瞬間届谈,已是汗流浹背枯夜。 一陣腳步聲響...
    開封第一講書人閱讀 32,983評(píng)論 1 269
  • 我被黑心中介騙來泰國(guó)打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留艰山,地道東北人湖雹。 一個(gè)月前我還...
    沈念sama閱讀 48,048評(píng)論 3 370
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像曙搬,于是被迫代替她去往敵國(guó)和親摔吏。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,864評(píng)論 2 354