橢圓型方程的差分解法

1、精選優(yōu)質(zhì)文檔-傾情為你奉上橢圓型方程的差分解法1.引言考慮問題二維Poisson方程， 其中為中的一個有界區(qū)域，其邊界為分段光滑曲線。在上滿足下列邊界條件之一：（第一邊值條件），（第二邊值條件），（第三邊值條件），都是連續(xù)函數(shù)，.2.差分格式將區(qū)間作m等分，記為；將區(qū)間作n等分，記為.稱為方向的步長，為方向的步長。2.1 Poisson方程五點差分格式參考單如圖所示：j+1jj-1i+1i-1i以為中心沿方向Taylor展開：由 + 得：整理得到：同理，以為中心沿方向Taylor展開：代入原方程得：得到截斷誤差：其中是原方程光滑解，舍去截斷誤差得到逼近Poisson方程的五點差分方程：考慮到邊

2、值條件，構(gòu)成差分格式：2.2 Poisson方程九點差分格式由上式 + 得：又則得到：舍去截斷誤差得到逼近Poisson方程的九點差分方程：考慮到邊值條件，構(gòu)成差分格式：3. 格式求解3.1 Poisson方程五點差分格式記，矩陣格式改寫為：，其中，可進一步寫為：3.2 Poisson方程九點差分格式記，矩陣格式改寫為：，其中，可進一步寫為：4.數(shù)值例子4.1 Poisson方程五點差分格式計算如下問題：其精確解為：考慮到本例中h1=h2，則有利用Gauss-Seidel迭代方法對k=0,1,2,計算表1 部分結(jié)點處的精確解和取不同步長時所得的數(shù)值解(h1,h2)(x,y)(1/4,1/4)(

3、1/2,1/4)(3/4,1/4)(1/4,3/4)(1/2,3/4)(3/4,3/4)(1/4,1/4)1.2.2.1.2.2.(1/8,1/8)1.2.2.1.2.2.(1/16,1/16)1.2.2.1.2.2.(1/32,1/32)1.2.2.1.2.2.(1/64,1/64)1.2.2.1.2.2.精確解1.2.2.1.2.2.表2 取不同步長時部分結(jié)點處數(shù)值解的誤差絕對值(h1,h2)(x,y)(1/4,1/4)(1/2,1/4)(3/4,1/4)(1/4,3/4)(1/2,3/4)(3/4,3/4)(1/4,1/4)8.4094e-040.00120.00119.0500e-04

4、0.00130.0012(1/8,1/8)2.1769e-043.1498e-042.8202e-042.3451e-043.3891e-043.0635e-04(1/16,1/16)5.4874e-057.9388e-057.1220e-055.9144e-058.5448e-057.7433e-05(1/32,1/32)1.3096e-051.9001e-051.7247e-051.4215e-052.0586e-051.8853e-05(1/64,1/64)6.1867e-068.3780e-064.7959e-065.5548e-067.4932e-064.0549e-06表3 取不同

5、步長時部分結(jié)點處數(shù)值解的最大誤差(h1,h2)(1/4,1/4)0.0016*(1/8,1/8)4.2732e-043.(1/16,1/16)1.0885e-043.(1/32,1/32)2.6148e-054.(1/64,1/64)1.1928e-052.圖1 取h=1/4時所得的數(shù)值解曲線圖2 取h=1/4時所得的誤差曲線圖3 取h=1/16時所得的數(shù)值解曲線圖4 取h=1/16時所得的誤差曲線圖5 取h=1/64時所得的數(shù)值解曲線圖6 精確解曲線圖7 取h=1/64時所得的誤差曲線4.2 Poisson方程九點差分格式計算如下問題：其精確解為考慮到本例中h1=h2，則有利用Gauss-S

6、eidel迭代方法對k=0,1,2,計算表1 部分結(jié)點處的精確解和取不同步長時所得的數(shù)值解(h1,h2)(x,y)(1/4,1/4)(1/2,1/4)(3/4,1/4)(1/4,3/4)(1/2,3/4)(3/4,3/4)(1/4,1/4)1.2.2.1.2.2.(1/8,1/8)1.2.2.1.2.2.(1/16,1/16)1.2.2.1.2.2.(1/32,1/32)1.2.2.1.2.2.(1/64,1/64)1.2.2.1.2.2.精確解1.2.2.1.2.2.表2 取不同步長時部分結(jié)點處數(shù)值解的誤差絕對值(h1,h2)(x,y)(1/4,1/4)(1/2,1/4)(3/4,1/4)(

7、1/4,3/4)(1/2,3/4)(3/4,3/4)(1/4,1/4)8.4094e-040.00120.00119.0500e-040.00130.0012(1/8,1/8)2.1769e-043.1498e-042.8202e-042.3451e-043.3891e-043.0635e-04(1/16,1/16)4.0638e-085.3998e-083.5875e-083.3749e-084.4845e-082.9794e-08(1/32,1/32)1.3096e-051.9001e-051.7247e-051.4215e-052.0586e-051.8853e-05(1/64,1/64

8、)6.1867e-068.3780e-064.7959e-065.5548e-067.4932e-064.0549e-06表3 取不同步長時部分結(jié)點處數(shù)值解的最大誤差(h1,h2)(1/4,1/4)5.1568e-09*(1/8,1/8)5.1568e-091.(1/16,1/16)6.9859e-080.(1/32,1/32)1.0268e-060.(1/64,1/64)1.5435e-050.圖1 取h=1/4時所得的數(shù)值解曲線圖2 取h=1/16時所得的數(shù)值解曲線圖3 取h=1/64時所得的數(shù)值解曲線圖4 取h=1/4時所得的誤差曲線圖5 取h=1/16時所得的誤差曲線圖6 取h=1/6

9、4時所得的誤差曲線5.結(jié)論觀察Poisson方程五點格式，方程以較快速度迭代收縮。對于Poisson方程九點格式，一開始有較快速度和較好結(jié)果，但是與五點格式比起來，迭代速度因步長增加而擴散，結(jié)果很不理想，且運行速度與五點格式相比慢了許多，與預(yù)計結(jié)果有失偏頗。參考文獻：1 孫志忠. 偏微分方程數(shù)值解法M. 北京：科學(xué)出版社， 20052 李榮華. 偏微分方程數(shù)值解法M. 北京：高等教育出版社， 20043 易昆南. 基于數(shù)學(xué)建模的數(shù)學(xué)實驗M. 長沙：中國鐵道出版社， 2014源代碼：%結(jié)點取值，結(jié)點誤差及最大誤差輸出函數(shù)functionp,e,u,x,y,k=five(h,m,n,kmax,ep

10、)%p代表數(shù)值解，h代表步長,m,n,分別代表x、y軸結(jié)點數(shù),kmax代表最大迭代次數(shù)，ep代表最大誤差,e用于存儲精確解syms temp;u=zeros(n+1,m+1);x=0+(0:m)*h;y=0+(0:n)*h;for(j=1:n+1) %邊值條件 u(j,1)=sin(y(j)+cos(y(j); u(j,m+1)=exp(1)*(sin(y(j)+cos(y(j);endfor(i=1:m+1) u(1,i)=exp(x(i); u(n+1,i)=exp(x(i)*(sin(1)+cos(1);endt=zeros(n-1,m-1);for(k=1:kmax) %存儲迭代結(jié)果for(j=2:n)for(i=2:m) temp=(u(j,i+1)+u(j,i-1)+u(j+1,i)+u(j-1,i)/4; t(j,i)=(temp-u(j,i)*(temp-u(j,i); u(j,i)=temp;endend t(j,i)=sqrt(t(j,i);if(kkmax) %結(jié)果滿足kmax次停止迭代break;endif(max(max(t)e); %輸出最大誤差 e=abs(u(j,i)-p(j,i);endif(x(i)=1/4&y(j)=1/4)|(x(i)=1/2&y(j)=1/4)|(x(i)=3/4&y(j)=1/4)|(x(i)=1/4&y(j)=3/