橢圓形方程的差分法

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)學(xué)與計(jì)算科學(xué)學(xué)院實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)項(xiàng)目名稱 橢圓方程的差分格式 所屬課程名稱 微分方程數(shù)值解 實(shí) 驗(yàn) 類 型 綜合型 實(shí) 驗(yàn) 日 期 2016-4-23 班 級(jí) 信計(jì)1301 學(xué) 號(hào) 1 姓 名 向 溶 成 績(jī) 一、實(shí)驗(yàn)概述：【實(shí)驗(yàn)?zāi)康摹?、建立內(nèi)點(diǎn)的五點(diǎn)差分格式。(取h=1) 2、建立包括邊界點(diǎn)在內(nèi)的五點(diǎn)差分格式方程組。 3、用雅克比迭代迭代法求解方程組。 4、計(jì)算結(jié)果（保留至小數(shù)點(diǎn)后4位）。 5、由計(jì)算結(jié)果，寫(xiě)出結(jié)論【實(shí)驗(yàn)原理】橢圓形方程的差分格式：正方形區(qū)域中的Ladplace方程Dirichlet邊值問(wèn)題的差分模擬：【實(shí)驗(yàn)環(huán)境】Matlab R2010a二、

2、實(shí)驗(yàn)內(nèi)容：【實(shí)驗(yàn)方案】在Matlab環(huán)境下分別用五點(diǎn)差分法編寫(xiě)解橢圓方程的算法，然后調(diào)試運(yùn)行，表述結(jié)果。【實(shí)驗(yàn)過(guò)程】（實(shí)驗(yàn)步驟、記錄、數(shù)據(jù)、分析）【實(shí)驗(yàn)結(jié)論】（結(jié)果）y = 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 3.8000 4.0000【實(shí)驗(yàn)小結(jié)】（收獲體會(huì)）通過(guò)本次實(shí)驗(yàn)，我們不僅對(duì)橢圓方程的差分法有了進(jìn)一步的了解，而且熟悉了雅克比求解這個(gè)矩陣形式，這對(duì)我們以后更好地學(xué)習(xí)它們奠

3、定了一定的實(shí)踐基礎(chǔ)，提高了我們一定的分析問(wèn)題和解決問(wèn)題的能力，也讓我們對(duì)所學(xué)知識(shí)有了更深刻的了解，同時(shí)也提高了我們的動(dòng)手能力。三、指導(dǎo)教師評(píng)語(yǔ)及成績(jī)：評(píng) 語(yǔ)評(píng)語(yǔ)等級(jí)優(yōu)良中及格不及格1.實(shí)驗(yàn)報(bào)告按時(shí)完成,字跡清楚,文字?jǐn)⑹隽鲿?邏輯性強(qiáng)2.實(shí)驗(yàn)方案設(shè)計(jì)合理3.實(shí)驗(yàn)過(guò)程（實(shí)驗(yàn)步驟詳細(xì),記錄完整,數(shù)據(jù)合理,分析透徹）4實(shí)驗(yàn)結(jié)論正確. 成 績(jī)： 指導(dǎo)教師簽名： 批閱日期：附錄1：源 程 序Helmtz.m：function u,x,y= helmtz(f,g,bx0,bxf,by0,byf,D,M,N,tol,maxiter)x0=D(1);xf=D(2);y0=D(3);yf=D(4);dx=(xf

4、-x0)/M;x=x0+0:M*dx;dy=(yf-y0)/N;y=y0+0:N'*dy;%邊界條件for m=1:N+1 u(m,1,M+1)=feval(bx0,y(m),feval(bxf,y(m);endfor n=1:M+1 u(1,N+1,n)=feval(by0,x(n),feval(byf,x(n);end%邊界的平均值作為初始值bvaver=sum(sum(u(2:N,1,M+1),sum(u(1,N+1,2:M);u(2:N,2:M)=bvaver/(2*(M+N-2);for i=1:N for j=1:M F(i,j)=feval(f,x(j),y(i); G(

5、i,j)=feval(g,x(j),y(i); endenddx2=dx*dx;dy2=dy*dy;dxy2=2*(dx2+dy2);rx=dx2/dxy2;ry=dy2/dxy2;rxy=rx*dy2;for itr=1:maxiter for j=2:M for i=2:N u(i,j)=ry*(u(i,j+1)+u(i,j-1)+rx*(u(i+1,j)+u(i-1,j)+rxy*(G(i,j)*u(i,j)-F(i,j); end end if itr>1&max(max(abs(u-u0)<tol break; end u0=u;end ex11.m:f=inli

6、ne('0','x','y');g=inline('0','x','y');x0=0;xf=4;M=20;y0=0;yf=4;N=20;bx0=inline('exp(y)-cos(y)','y');bxf=inline('exp(y)*cos(4)-exp(4)*cos(y)','y');by0=inline('cos(x)-exp(x)','x');byf=inline('exp(4)*cos(x)-exp(x)*cos(4)','x');D=x0 xf y0 yf;maxiter=500;tol=1e-