偏微分方程的數(shù)值解法課后習題答案概要

第二章習題答案1.試討論il近對流方程—=0的差分方程。

dt班此廣1―熄〃”一般尸1)/ L+」 -*1)T h2)3 L+u——=0口r 2b的截斷誤差。~

1）解：設點為（犬JJMl）中則吟=認、&)="(,&+])+駕也(―切+0(?).dt口;1；二漢（力一口;1；二漢（力一1，/必）=笈（%.//J+加巧?應十1）dx所以截斷誤差為：裂我，】一口;所以截斷誤差為：裂我，】一口;我;“-我E= L+_i ：t h—(把+包”

dtdxdx「h—dx「h—。圖）-尸…+di=0（r+h）72）?!設點為：（卻，￡葉I）“則嗎=以勺&)=〃(,&+1)+ (f)+0(T>（我2）+口dt（我2）+口算"1="Oj+iJ/J=〃(兀jJ/i)+ ; (卜)+以4；：；="（％_1，4+1）=打（仃，44；：；="（％_1，4+1）=打（仃，4十］）+法（f%gdxM）+”/）*.-.就新誤差為；"qip〃*一尸 也出E=上——+ 」一（―+—）~r h didx也（彳「與心）._ 雙燈&+1）一雙弓?/+】）+—耳—丁+0口）T

加f￥.￡八/川X.f.,、2.試用積分插值法推導畀上。逼近的差分格式”71—般;或71工——+__^=o0~T hMg四、丁口「， 、7c解：11(—+—)^^=|(^,+OTj^=Op?3優(yōu)dr,U JU-附近+玩2匯+笈3｝一亮4&=0E(j-1,n)F(j」n)G(j,n+l)H(廣1通十l)v，4-ftjn2=〃； “3=”： "4=〃；-ih=hr=r,一成h+K"+LL；*h一笈:=Ouu JJ J,P784J.如果如'（0）=0,則稱而是J（0）的駐點（或穩(wěn)定點）。設矩陣A對稱（不必正定），求證而是J。）的駐點的充要條件是！大。是方程Ax=b的解.dj2證；充分性；①⑷=JQo）+a（③oT㈤+二（4,力2①⑷=(刃/一瓦工)+2(工Zx)①'(0)=0而一九X)8若&（0）=0,即5仆一次力=0txw及鼠■飛一8=0即Ax=b^則x0是方程Ax=b的解。必要性：若為是方程Ax=b的解“則j4x0-i=0Q4xq—瓦1)=On$（0）二（念0-瓦力=01所以所是J（x）的駐點，P93p3=證明非齊次兩點邊值問題，+夕歸=/a<x與下面的變分問題等階：求以EH19%.（以）=G使J@.)=mmJ(w)其中u“（2）T，⑷=；以3,切—5%-「⑻短⑻，乙而以3"）如（2.13）（提示，先把邊值條件齊次化）證明；令=vp(x)+v(z)其中印(兀)=a+(1一以)^^(a)=av(a)=0v(b)=0^所以〃TOC\o"1-5"\h\z上以二一;0號)+"=」歐dx 口dr,dwdp、、/ 、r=_~rlp{—十丁)]+Q+M=jax dx ax人「 出, 如、 ，，及 處 、上令 =-—(^―)+^v=/-(--7?--+^)=f^dx ax ax ax所以(1)的等價的形式。rd,孫、， /Lu=--(p-)^^=f1daxaxu(a)=a/(?)=0-^其中左二，一(一字+qw)，axax則由定理2.2知，y?是邊值問題(2)的解的充要條件是JE明且滿足變分方程”以口㈤-5J)=0WE*又9/)一(工⑷u=|(Zvt-fy)tdx+p@)v；(班? (3)y?2①⑷=，&)=J&.+丸〉=-a(u^+為”+龍)-(丁必+上)—p(b)規(guī)以?⑻+龍⑼]21=3)+孫3，/)-(/~@)龍(砌+彳。3乙日3?/）一（〃）- d=’6字小+q以/_fiW-p?恿3）*'，°axax=「（￡氏-/,小+「@）〃；（句“勾一戶（4）4（3）=>（4）所以可證得。3必要性：若電是邊值問題（1）的解.則卬?7=0“:3）二產所以a3.x）一（/X）—「。）#@）二0 ①'（6=0且而用>二中@〃解：由上一題知:(2.28)等價于口rd/小八，/ddwv￡2^="-Qi?)+^v=/-(--^―+^)=/^以ax axaxv(a)=0 v'(^)=0^v|(Zv一f^idx=|p +qvt一jQdx=0*,g gdxdx' 」即a(yQ-Gi?=0"又因為以①⑷一(加)=「加一九紜k+9(如'⑸我如*e所以v(b)=/3)-M?=-?-g卜一衣=Lu-f所 以①/)—(力?=『?-力以天+M句以1(g@)一7。)#3)二o,?2又因為血￡)-(/,￡)“也dudt「一=[p---+qui-ft]dx^、。dx西=p@)t/@)￡@)+|{Lu-所以得⑶#@)=0"即以是(2.28)的解的充分必要條件為，4@〃)一(九0一》@)就@)=?？赑93y5：試建立與邊值問題尸.d%. ,La=-j-+玄=/ a<x<匕》加u(a)=u'(b)=0 %0)=/(5)二0的等價變分問題J解：設Vve/fi且v(a)=v\a)=Ov(b)=v=0^『3 r$fO?RdI(Zw-fy^dx=i——vd?x+Imv-?2 皿辦'「號曲 dv,—rvL_—r—dx.'dx?y治d必必/八一,、 ■,八，八、■‘、、、f"d"ad^v.(b)v(^)-u@>@)+打⑷y⑷+「加戒d"=2⑶ t/@"⑻+a'(a)S⑷+(^~^udx3qd￡又v(a)=v(a)=v(i)=v(h)=0^所以f整小小余介-4廣力 r?V-d所以I{Lu—f)vdx=Iu——-+Iwv-Jvdx^?2 皿靖s_2.4*y 一其中次譏》=|u—r-+uv)dx是一個雙線性泛函。J?。d?x4所以迫值問題的變分問題為“求u&Hg使得da@”)-CAM=0，爪e*且v(a)=v(a)=0找出)=/(b)=UpJ⑸=3(&X>—(/M〃乙P103^3；試就pazsso^方程。.3)的非齊次定值條件(3.31)——+anI=5a:>=0^

沏lr導出等價的變分問題。解:取一特定函數(shù)入2C%)+”■+叫|二戶令v=ti-u^則dndv.門—+^vL="&?則得(3.3) (3.31)的等價問題會—Av=j+A%二尸v加 I p——+dtvl=0~灰lz所以7(v)=1(-Av,V)-(F,p)^2=；①,y)-(F,r)"乙Ti【(少+(第收的+；Iavds-f|Fvdxdy^2-dx砂 2?y-1nr/^\2/次」tfduM也也。2JQ勺)+(—)dxdy-I[-___+—--]^279zdy7也所在中常數(shù)算+—|ailds—|aiitiQds-|\fitdxdy-(\h.ii^udxdy+TOC\o"1-5"\h\z21 )—p?aja」般/a?j ?Ajy ??=Jg)R俘等+=等收砂+J(粵-灼加T他遇心辦工力Bxdydy ;dn -+常如其中J(u)-［f［史史1+絲_曳￡］歐砂+|au"而一(jjudxdy^:以辦辦8； -1；又由格林第一公式知道〃［JuIT仔答+*筌兩辦一］爭功7 7班班?? :獨原問題的變分問題的為〃求li,eH\(u)使得UJQJ=minSa..--wjj|J(y)=—|j[- jau2ds-||judxdy-|肉ds2』.小小 辦辦2； ?』 ；6.P104/4/試就腌圓方程第一邊值問題，口(3.32)-V(r7%)+g=」,(x,y)eG,w|r=g建立等價的極小位能原理和虛功原理，其中左=左(i,y)e/(G),噢>0,"beC@),b20Je￡“◎geC(T),而〃「「「、d,.du.ddie、V(歸V%)二 (k )+ (k )."dxdxdydy解；（1）極小位能原理；d設劭wc"。）為一特定的數(shù)，w0\T=g令解；（1）極小位能原理；d設劭wc"。）為一特定的數(shù)，w0\T=g令v=u-u^則得（3.32）的等價問題；~? a%,-V（A：Vv）+ov=尸=」+—（紀一^）—明勿沙 8JIt=0J(y)=1(―v07切+ —(F,y)2=-||-V(Wv)v^+-||加七二心,|[四心心乙G 乙G G*J||-V(^Vv)vcfx^y加、1)1—)]vdxdy=- v+尢udxdxd2vv]dxdyi v+ V飛dxdxdydykv)dxdy]Gre近%一公式N，dk:5vdk:111 v+ v股o+d皤等等岫砂一年"oxoxcyoy j,o）2=腦少+（即M星砂班dyTOC\o"1-5"\h\z\o"CurrentDocument"：i?dv3udv *.,?令4=H^(——+——f)dxdy+Hcuvdxdy^礪赤方力 。點則九w=Le?）—（艮y）“2下面回到原問題?,八 1J⑺=5以憶切一（凡日乙=；U上噂一筌）2+儀當_等）2]小砂+1。6以一%）2石右2飛?k uy 2飛一口L/+3（左等）+（（上等）-g°]3-%）心功》dxdxdydy1rrri/也、全,,班、）、,,ft-du面介du獨力、,,{t>n,二K伏%-）+上（77）]以的一武丁二一十kK"）"不分十?&xdy2^oxdy 3&x&x&ydy力-II -||JudTdy-|j[^-（^^）* * %浙班第二章第三章第四章第五章第六章pl99“1.試從二次泛函〃出發(fā)導出差分方程。，解：J（〃）=?!4 |（%_-,匚）2小一々（"匕+七2工士無H2?等（%-*，）|*2￡a= 2J（〃）對小走偏導等于0*-,/-%-i/-%-i_%」—%_2r-h--h-2J??T17fdx--|/dx=CM2xu“i-2%+%]但一/-2=/1=12…,N“i%其中0i=77T（I/dx+l了dx）〃2力；_ ；1.用積分插值法導出逼近微分方程（2.1）的差分方程。

JUi?(2.1)Lu=——(p—)+rdxdxJUi?a<x<bp,q,r」eC（l＞可以直接積分"a<x<b在［a,對內任一小區(qū)間［x⑴⑵］上積分有一w(x⑴)-w(x⑶)+J尸竺dx+［q葭dx=J/x("dx* 人其中w(x)=p也在［a,3］上連續(xù).,dxr—ax+

dxIqudx=?。踃⑴，五⑵］為對偶單元［五1，五r—ax+

dxIqudx=W(X］)=—W(X1)吃 啊/ 、du?.?W(X)=p——dx也=皿dxp(x)du,cw(x),——dx=-dxpdx/p(x)R/T%%逢W(x1)f—!—dx(中矩形公式)y… 4/。⑶w（x,1）領/%-i）（j -（；J-y-，三 二P（x）K瓦二。⑶令名=（—f」-dx）（XquXqudxR々+々包％4其中同=——［2 %十如X!q(x)dx。p200.4^4.構造逼近?（p〃"）"+‘打"+廣w=/于（4,b）"（a歸以'（a）=0，〃（5）=〃'（/＞）=03的中心差分格式。6解：取少+1個節(jié)點,以=而＜xi＜……＜號＜……人= 1J=l,2, ,N, h=max無”X1=:(為-1+號)7=12 N")2則有，d〃du %一%-I2 L—J.1L—J1—； 一一： d2udx入dx力如加"3 %+/ %+/~2~ -Fd2d2u“必"[以加一4 %一%-i1 22)[烏p嗎卜/d/2rdd，,rddM[7Pyr〕i-丘。"7^].i信.l+'idxdx1+]axdxr? -信.l+'i[。/+l-[p曲[p^-[pR>\ ]. ： -Pi[-陶+2+自.1)4.12__(公]2__(公]+%)公1 垢1%一%-; 2 -% %+%)%-2（%+如）如d2d2d2d2u1 1[pi+i(/a-2%.i+%〉-2Pi(%.「2%+%.iP210d1,用積分插值法構造逼近方程~（3.21）-V（kVu）="[W.蟲）+2（先生）卜/的第-dxdx分格式，這里k=發(fā)（五，y）> >0^-'解：1)正則內點解：1)正則內點d于G0.上積分(3.21)式，一-J[V(kVu)dxdy=JJfdxdy^由Green第一公式得：小-[^kds

京dn=\\Jdxdy^J+J+LZ、+J4)?上ds=]JJdxdy/du.?du.—kds=———kdn dy也7Jdu? J—"，J?——kds=——ki%=ki—— 飽，dndx，?襯 ，\Jh[r du 以3 —%3I—kds=—k1況i=k ~~士h!加dy國3 % ’f—kds=—k]h2=k]上山_—%~；dndx ，"J力]；.綜上有：“以；，?[—tii,, 以丘1,，i tlqiliq,,-[<i岑2+ki「：J±+k-””+k2方 力2 ’V "1 2、 %="/,其中做/=3"/dxdy?力也G1)非正則內點"

2.用積分插值法構造逼近方程（3。21）的第二邊值問題的五點差分格式。解：1）正則內點，同第一題中1）標2）非正則內點2）非正則內點9同第一題中2）?在界點口處于曲邊三角形ABC上對(3.21)式積分,-jjV(kVu)dxdy=JJfdxdy^hABC IlABC~ff^-kds=JJjdxdy^~(~(.U)—kdsJJ/dxdy/

aabc，B(x,y)kds=3ok，B(x,y)kds=3ok於鑫*一 ?1 ?0— —AB