椭圆形方程的差分法_第1页
椭圆形方程的差分法_第2页
椭圆形方程的差分法_第3页
椭圆形方程的差分法_第4页
椭圆形方程的差分法_第5页
已阅读5页,还剩4页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上数学与计算科学学院实 验 报 告实验项目名称 椭圆方程的差分格式 所属课程名称 微分方程数值解 实 验 类 型 综合型 实 验 日 期 2016-4-23 班 级 信计1301 学 号 1 姓 名 向 溶 成 绩 一、实验概述:【实验目的】1、建立内点的五点差分格式。(取h=1) 2、建立包括边界点在内的五点差分格式方程组。 3、用雅克比迭代迭代法求解方程组。 4、计算结果(保留至小数点后4位)。 5、由计算结果,写出结论【实验原理】椭圆形方程的差分格式:正方形区域中的Ladplace方程Dirichlet边值问题的差分模拟:【实验环境】Matlab R2010a二、

2、实验内容:【实验方案】在Matlab环境下分别用五点差分法编写解椭圆方程的算法,然后调试运行,表述结果。【实验过程】(实验步骤、记录、数据、分析)【实验结论】(结果)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【实验小结】(收获体会)通过本次实验,我们不仅对椭圆方程的差分法有了进一步的了解,而且熟悉了雅克比求解这个矩阵形式,这对我们以后更好地学习它们奠

3、定了一定的实践基础,提高了我们一定的分析问题和解决问题的能力,也让我们对所学知识有了更深刻的了解,同时也提高了我们的动手能力。三、指导教师评语及成绩:评 语评语等级优良中及格不及格1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强2.实验方案设计合理3.实验过程(实验步骤详细,记录完整,数据合理,分析透彻)4实验结论正确. 成 绩: 指导教师签名: 批阅日期:附录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-

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论