二维波动方程的有限差分法
9页1、学 生 实 验 报 告实验课程名称 偏微分方程数值解 开课实验室 数统学院 学 院 数 统 年级 2013 专业班 信计02班 学 生 姓 名 学 号 开 课 时 间 2015 至 2016学年第 2 学期总 成 绩教师签名数学与统计学院制开课学院、实验室: 数统学院 实验时间 : 2016年 6月20日实验项目名 称二维波动方程的有限差分法实验项目类型验证演示综合设计其他指导教师曾芳成 绩是一实验目的通过该实验,要求学生掌握求解二维波动方程的有限差分法,并能通过计算机语言编程实现。二实验内容考虑如下的初值问题: (1)1在第三部分写出问题(1)三层显格式。2根据你写出的差分格式,编写有限差分法程序。将所写程序放到第四部分。3取,分别将时刻的数值解画图显示。4. 该问题的解析解为,将四个时刻的数值解的误差画图显示,对数值结果进行简单的讨论。三实验原理、方法(算法)、步骤网格划分,故,。在内网点,利用二阶中心差商,对(1)建立差分格式: (2)整理得到: (3)其中,网比,局部截断误差为。考虑边界条件,差分格式为: (4)考虑初始条件,差分格式为: (5)考虑初始条件,利用二阶差商近似:
2、 (6)设时刻的点为内点,则满足差分格式(2),代入上式得到: (7)将(6)得到的结果代入(7)中,整理得到: (8)综上(2)、(4)、(5)、(8)得到三层显格式的差分格式为: (9)其中,局部截断误差为。四实验环境(所用软件、硬件等)及实验数据文件Matlab%二维波动方程数值计算(关键:怎么运用i,j,k三个指标建立循环)clc;%可以将代码换成函数m文件h=0.1;tau=0.1*h;%定义步长r=tau/h;%网比x,y,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);%初始条件u1=u(:,:,1);%因为此时得到的u为11x11x141,故只取第一层%第二层网点计算for i=2:10 for j=2:10 u(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)+(1-2*r2)*u(i,j,1); u(11,:
3、,2)=0;u(:,11,2)=0; endendu2=u(:,:,2);%第3-141层网点计算for k=2:140 for i=2:10 for j=2:10 u(i,j,k+1)=r2*(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,j,k-1); u(11,:,k+1)=0;u(:,11,k+1)=0; end endend%结果分析与作图%wucha=abs(u-uu);%求绝对误差矩阵11x11x141wucha1=wucha(:,:,11);%计算t=0.1时刻的绝对误差矩阵11x11wucha2=wucha(:,:,51);%计算t=0.5时刻的绝对误差矩阵11x11wucha3=wucha(:,:,101);%计算t=1.0时刻的绝对误差矩阵11x11wucha4=wucha(:,:,141);%计算t=1.4时刻的绝对误差矩阵11x11x0=0:h:1;y0=0:h:1;%误差分析%作t=0.1时刻的绝对误差图subplot(2,2,1);mesh(x0,y0,wucha1);t
《二维波动方程的有限差分法》由会员suns****4568分享,可在线阅读,更多相关《二维波动方程的有限差分法》请在金锄头文库上搜索。
土地管理与地籍测量---第八章界址点测量
人机工程学案例分析(2)
工程安全培训_201303
第9章房地产投资决策分析
第2章房地产经纪制度
ACM程序设计-东北林业大学acm05
《亲爱的汉修先生》读书交流会
中原_深圳新世界尖岗山项目市场汇报_40P_2012年_别墅_项目分析_量价走势
五年级数学质量分析演示文稿
人工智能小镇-智慧小镇建设20180525
景观基本知识及发展历程
建设工程信息管理(2)
机电驱动技术第二章步进驱动技术
工程力学-第9章圆轴扭转时的应力变形分析与强度刚度设计
第一章第二节幼儿园文化环境建设的原则
第一章检测技术的基础知识
第一章__现代表面工程技术
第六章钢结构工程
第9节项目试运行管理
班主任工作经验交流课件(4)
2023-05-26 28页
2023-04-20 24页
2024-03-05 18页
2023-06-19 117页
2023-06-25 9页
2023-08-29 9页
2023-05-12 24页
2024-02-27 19页
2023-12-05 47页
2023-08-11 2页