二维波动方程的有限差分法.docx

资源ID：1671395 资源大小：98.97KB 全文页数：13页
资源格式： DOCX 下载积分：5金币

快捷下载

会员登录下载

三方登录下载：

下载资源需要5金币

邮箱/手机：
温馨提示：	用户名和密码都是您填写的邮箱或者手机号，方便查询和重复下载（系统自动生成）
支付方式：
验证码：	换一换

加入VIP免费专享

账号：
密码：
验证码：	换一换
当日自动登录忘记密码？

友情提示

1、下载资料失败解决办法

2、PDF文件下载后，可能会被浏览器默认打开，此种情况可以点击浏览器菜单，保存网页到桌面，就可以正常下载了。

3、本站不支持迅雷下载，请使用电脑自带的IE浏览器，或者360浏览器、谷歌浏览器下载即可。

4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩，下载后原文更清晰。

5、试题试卷类文档，如果标题没有明确说明有答案则都视为没有答案，请知晓。

网站客服

侵权投诉

二维波动方程的有限差分法.docx

学生实验报告试疆课程名称偏徽分方程数值解开课试验室数统学院学院数统年级2013专业班信计02班学生姓名学号开课时间2015至2016学年第2学期总成绩老师签名数学与统计学院制开课学院、试验室：数统学院试验时间:2016年6月20日试验项目名称二维波动方程的有限差分法试验项目类型验证演示综合设计其他指导老师曾芳成绩是一.试验目的通过该试验，要求学生驾驭求解二维波动方程的有限差分法，并能通过计算机语言编程实现。二,试验内容考虑如下的初值问题：票=粤+斐,(XMeC=(OjfJW(0,1.4)drx'y'u(x.y.0)=sinsin11y,w(x.y.)=(-t.y)(0.1)*(1)(x,y,r)=0.(x.y)ece,1.41 .在第三部分写出问题(1)三层显格式。2 .依据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。3 .¾=0.1.r=0.I,分别将f=051.0.1.4时刻的数值解画图显示。4 .该问题的解析解为“(x.y.t)=cosVmsin.rsin11y»将四个时刻的数值解的误差画图显示，对数值结果进行简洁的探讨。三.试验原理、方法(算法)、步骤网格划分=0.1.r=().1,故；V=1=10,/=140,hX1=Hi,yf=jh,i,j-0,1,IO>ti-Ar,=0,1.,14()0在内网点(内,)也)，利用二阶中心差商，对(1)建立差分格式："：''-2m*+>m.1.,z-2哈+"1吟”-2“：,+j.1.(2)整理得到：口=/(哈J+哈J+H+%)+(2-4r2M哈(3)其中，i,j=12,9,=1,2,.139,网比r=>0.1,局部截断误差为“(，+/)。考虑边界条件“(.v,/)=0,(X.y)w次Zr,1.4,差分格式为:脸>=$n=Mq=心,n=-=0,140(4)考虑初始条件“(X,1v.0)=sinsin11y.差分格式为：it-i=sin(11xi)sin(11yi)=sin(sin(cjh),/,y=OJ.10(5)若虑初始条件(A-,y,O)=O,(.r,y)(O,1.)：»利用二阶差商近似：u>.,½=ojj=10(6)设K=O时刻的点为内点，则满意差分格式(2),代入上式得到：“：/=/("Mj+心./+“+%)+(2-4-用-吗(7)将(6)得到的结果*=0代入中，整理得到：=9W1.J+1C-U+-+心i)+(J2产)（8）综上（2）、（4）、（5）、（8）得到三层显格式的差分格式为:uZ'='(心j+<.+*T)+(2-4/)-<,/J=1.2.9,=1.2.139<嗑=&z="；.O=心.n=0/=01./40(9)it'j=sin(.r.)sin(11,y7)=sin(-)sin(11jh),i,j=0.1,.1()=2(w.1.'u+吟“+<-1)+(>-2,)<u',7=O1.,10其中亏=0,局部截断误差为。任+的。四.试验环境(所用软件、硬件等)与试验数据文件Mat1.ab%二维波动方程数值计算(关键：怎么运用i,j,k三个指标建立循环)c1.c;%可以将代码换成函数m文件h=0.1.;tau=0.1.*h;%定义步长r=1.auh;/网比x,y1.t=meshgrid(0:h:1.,0:h:1.,0:tau:1.4);%空间网格剖分UU=Cos(sqrt(2)*pi*t).*sin(pi*x).*sin(pi*y);%精确解计算%第一层网点计算u=sin(pi*x).*sin(pi*y);%初始条件UI=U(:,:,D;%因为此时得到的u为11x11x141,故只取第一层%其次层网点计算fori=2:10forj=2:10u(i,j,2)=0.5*r2*(u(i+1.,j,1.)+u(i-1.,j,1.)+u(i,j+1.,1.)+u(i,j-1.,1)+(12*r-2)*u(i,j,1);u(1.1.,:,2)=0;u(:,11,2)=0;endendu2=u(:,:,2);%第3T41层网点计算fork=2:140fori=2:10forj=2:10u(i,j,k+1.)=r'2*(u(i+1.,j,k)+u(i-1.,j,k)+u(i,j+1.,k)+u(i,j-1.,k)÷(2-4*r2)*u(i,j,k)-u(i,jfk-1.);u(1.1.,:,k+1.)=0;u(:»1.1.,k+1.)=0;endendend%结果分析与作图犒哪嬲胱(三1.W嬲黑略用绑%概蹴%wucha=abs(u-uu),求肯定误差矩阵UXUXI41WUChaI=WUCha计算t=0.1时刻的肯定误差矩阵I1.XUwucha2=wucha(:,:,51);%计算t=0.5时刻的肯定误差矩阵I1.x1.1.wucha3=wucha(:,:,101);%计算t=1.0时刻的肯定误差矩阵I1.x1.1.WUCha4=wucha(:,:,141);%计算t=1.4时刻的肯定误差矩阵UXUx0=0:h:1.;y0=0:h:1.;%误差分J%作t=0.1时刻的肯定误差图subp1.ot(2,2,1):mesh(x,y,wucha1.);titIe(,t=0.1时刻的肯定误差')；x1.abe1.(,X变量');y1.abe1.(,y变量');Z1.abe1.C肯定误差值');%作t=0.5时刻的肯定误差图subp1.ot(2,2,2);mesh(x,y,wucha2);tit1.e('t=0.5时刻的肯定误差');XIabe1('X变量)；y1.abe1.('y变量)；Z1.abe1('肯定误差值)；先作t=1.0时刻的肯定误差图subp1.ot(2,2,3)：mesh(x,y,wucha3);Iit1.eCt=1.O时刻的肯定误差')；x1.abe1.CX变量')；y1.abe1.('y变量');Z1.abe1('肯定误差值');学作t=1.4时刻的肯定误差图subp1.ot(2,2,4);mesh(x,y,WUCha4);tit1.e('t=1.4时刻的肯定误差')；x1.abe1.(,X变量'Ay1.abe1.Cy变量');Z1.abe1.C肯定误差值');%四个时刻数值解、精确燃作t=0.k0.5时刻的数值解与精确解subp1.ot(2,2,1);mesh(x,y,u(:,:,11);%作t=0.1时刻的数值解tit1.e('t=O.1时刻的数值解')；x1.abe1.(,X变量')；y1.abe1.('y变量')：z1.abe1.(,U值');subp1.ot(2,2,2)：mesh(x,y,uu(:,:,11);%作t=0.1时刻的精确解tit1.e('t=O.1时刻的精确解x1.abe1.(,X变量');y1.abe1.(,y变量')；ZIabeIeU值')；嬲作t=0.5时刻的数值解与精确解subp1.ot(2,2,3);nesh(x,y,u(:,:,51);%作t=0.5时刻的数值解tit1.e(,t=0.5时刻的数值解')；x1.abe1.CX变量');y1.abe1.(,y变量')jz1.abe1.Cu值')；subp1.ot(2,2,4);nesh(x,y,uu(:,:,51)；%作t=0.5时刻的精确解tit1.e('t=0.5时刻的精确解')；x1.abe1.CX变量');y1.abe1.(,y变量');ZIabe1.cU值');%分别复制粘贴运行嬲嬲WW1.%嬲%w%w三wm%嬲作t=1.O、1.4时刻的数值解与精确解subp1.ot（2,2,1）：mesh（x,y,u（:,:,IOD），作t=1.0时刻的数值解tit1.e（'I=1.O时刻的数值解'）；x1.abe1.（,X变量'）；y1.abe1.（'y变量）;Z1.abe1（'U值'）；subp1.ot（2,2,2）;mcsh（x,y,uu（:,:,101）;%作t=1.0时刻的精确解tit1.e（'t=1.0时刻的精确解'）;x1.abe1.CX变量'）：y1.abe1.（,y变量'）;ZIabeICU值'）；般作t=1.4时刻的数值解与精确解subp1.ot（2,2,3）;mesh（x,y,u（:,:,MD）;%作t=1.4时刻的数值解tit1.e（'t=1.4时刻的数值解x1.abe1.CX变量'）：y1.abe1.（,y变量'）;ZIabe1.（'u值'）;subp1.ot（2,2,4）;mesh（x,y,uu（:,:,14I）作t=1.4时刻的精确解tit1.e（,t=1.4时刻的精确解'）;x1.abe1.Cx变量'）;y1.abe1.（'y变量'）;ZIabe1.（'u值'）;五.试验结果与实例分析1、/=0.1Q5、1.OJ.4时刻的数值解与精确解图图1t0.1、0.5时刻的数值解、精确解图2t=1.0,1.4时刻的数值解、精确解注：上两图为四个时刻的数值解与精确解，r=（=（）.1.vj；（P代表维数），本文/,=2，三层显格式达二阶收敛，不难看出，收敛效果很好，符合理论。下图是四个时刻的肯定误差图像，从图中看出，肯定误差较小，且经过计算得到，收敛阶近似于2,正好符合理论值。图3四个时刻的肯定误差3、四个时刻（t=0.10.5、1.0、1.4）的肯定误差表t=0.1时刻的肯定误差0.00.00.00.00.00.00.00.00.00.00.00000000000000000000000000000000000.00.00.00.00.00.00.00.00.00.00.00000010010020020020020020010010000000010030040010050040040030010000.00.00.00.00.00.00.00.00.00.00.00000020010050060060060050040020000.00.00.00.00.00.00.00.00.00.00.00000020040060070070070060040020000.00.00.00.00.00.00.00.00.00.00.00000020050060070080070060050020000.00.00.00.00.00.00.00.00.00.00.00000020040060070070070060040020000.00.00.00.00.00.00.00.00.00.00.00000020040050060060060050040020000.00.00.00.00.00.00.00.00.00.00.00000010030040040050040040030010000.00.00.00.00.00.00.00.00.00.00.00000010010020020020020020010010000.00.00.00.00.00.00.00.00.00.00.0000000000000000000000000000000000t=0.5时刻的肯定误差0.00.00.00.00.00.00.00.00.00.00.00000000000000000000000000000000000.00.00.00.00.00.00.00.00.00.00.00000070130180210220210180130070000000130250340400420400340250130000.00.00.00.00.00.00.00.00.00.00.00000180340470550580550470340180000.00.00.00.00.00.00.00.00.00.00.00000210400550650680650550400210000.00.00.00.00.00.00.00.00.00.00.00000220420580680710680580420220000.00.00.00.00.00.00.00.00.00.00.00000210400550650680650550400210000.00.00.00.00.00.00.00.00.00.00.00000180340470550580550470340180000.00.00.00.00.00.00.00.00.00.00.00000130250340400420400340250130000.00.00.00.00.00.00.00.00.00.00.00000070130180210220210180130070000.00.00.00.00.00.00.00.00.00.00.0000000000000000000000000000000000t=1.0时刻的肯定误差0.00.00.00.00.00.00.00.00.00.00.00000000000000000000000000000000000.00.00.00.00.00.00.00.00.00.00.00000160310430510530510430310160000000310590820961010960820590310000.00.00.00.00.00.00.00.00.00.00.00000430821131321391321130820430000.00.00.00.00.00.00.00.00.00.00.00000510961321561641561320960510000.00.00.00.00.00.00.00.00.00.00.00000531011391641721641391010530000.00.00.00.00.00.00.00.00.00.00.00000510961321561641561320960510000.00.00.00.00.00.00.00.00.00.00.00000430821131321391321130820430000.00.00.00.00.00.00.00.00.00.00.00000310590820961010960820590310000.00.00.00.00.00.00.00.00.00.00.00000160310430510530510430310160000.00.00.00.00.00.00.00.00.00.00.0000000000000000000000000000000000t=1.4时刻的肯定误差0.00.00.00.00.00.00.00.00.00.00.00000000000000000000000000000000000.00.00.00.00.00.00.00.00.00.00.0000002003005006006006005003002000OOO003007009O1.1.O1.1.O1.1.0090070030000.00.00.00.00.00.00.00.00.00.00.00000050090120150150150120090050000.00.00.00.00.00.00.00.00.00.00.0000006O1.1.015017018017015O1.1.0060000.00.00.00.00.00.00.00.00.00.00.0000006O1.1.015018019018015O1.1.0060000.00.00.00.00.00.00.00.00.00.00.0000006O1.1.015017018017015O1.1.0060000.00.00.00.00.00.00.00.00.00.00.00000050090120150150150120090050000.00.00.00.00.00.00.00.00.00.00.0000003007009O1.1.O1.1.O1.1.0090070030000.00.00.00.00.00.00.00.00.00.00.00000020030050060060060050030020000.00.00.00.00.00.00.00.00.00.00.0000000000000000000000000000000000老师签名

注意事项

本文（二维波动方程的有限差分法.docx）为本站会员（夺命阿水）主动上传，课桌文档仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对上载内容本身不做任何修改或编辑。若此文所含内容侵犯了您的版权或隐私，请立即通知课桌文档（点击联系客服），我们立即给予删除！

温馨提示：如果因为网速或其他原因下载失败请重新下载，重复下载不扣分。