




已阅读5页,还剩4页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
学 生 实 验 报 告实验课程名称 偏微分方程数值解 开课实验室 数统学院 学 院 数 统 年级 2013 专业班 信计02班 学 生 姓 名 学 号 开 课 时 间 2015 至 2016学年第 2 学期总 成 绩教师签名数学与统计学院制开课学院、实验室: 数统学院 实验时间 : 2016年 6月20日实验项目名 称二维波动方程的有限差分法实验项目类型验证演示综合设计其他指导教师曾芳成 绩是一实验目的通过该实验,要求学生掌握求解二维波动方程的有限差分法,并能通过计算机语言编程实现。二实验内容考虑如下的初值问题: (1)1在第三部分写出问题(1)三层显格式。2根据你写出的差分格式,编写有限差分法程序。将所写程序放到第四部分。3取,分别将时刻的数值解画图显示。4. 该问题的解析解为,将四个时刻的数值解的误差画图显示,对数值结果进行简单的讨论。三实验原理、方法(算法)、步骤网格划分,故,。在内网点,利用二阶中心差商,对(1)建立差分格式: (2)整理得到: (3)其中,网比,局部截断误差为。考虑边界条件,差分格式为: (4)考虑初始条件,差分格式为: (5)考虑初始条件,利用二阶差商近似: (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,:,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);title(t=0.1时刻的绝对误差);xlabel(x变量);ylabel(y变量);zlabel(绝对误差值);%作t=0.5时刻的绝对误差图subplot(2,2,2);mesh(x0,y0,wucha2);title(t=0.5时刻的绝对误差);xlabel(x变量);ylabel(y变量);zlabel(绝对误差值);%作t=1.0时刻的绝对误差图subplot(2,2,3);mesh(x0,y0,wucha3);title(t=1.0时刻的绝对误差);xlabel(x变量);ylabel(y变量);zlabel(绝对误差值);%作t=1.4时刻的绝对误差图subplot(2,2,4);mesh(x0,y0,wucha4);title(t=1.4时刻的绝对误差);xlabel(x变量);ylabel(y变量);zlabel(绝对误差值);%四个时刻数值解、精确解%作t=0.1、0.5时刻的数值解与精确解subplot(2,2,1);mesh(x0,y0,u(:,:,11);%作t=0.1时刻的数值解title(t=0.1时刻的数值解);xlabel(x变量);ylabel(y变量);zlabel(u值);subplot(2,2,2);mesh(x0,y0,uu(:,:,11);%作t=0.1时刻的精确解title(t=0.1时刻的精确解);xlabel(x变量);ylabel(y变量);zlabel(u值);%作t=0.5时刻的数值解与精确解subplot(2,2,3);mesh(x0,y0,u(:,:,51);%作t=0.5时刻的数值解title(t=0.5时刻的数值解);xlabel(x变量);ylabel(y变量);zlabel(u值);subplot(2,2,4);mesh(x0,y0,uu(:,:,51);%作t=0.5时刻的精确解title(t=0.5时刻的精确解);xlabel(x变量);ylabel(y变量);zlabel(u值);%分别复制粘贴运行%作t=1.0、1.4时刻的数值解与精确解subplot(2,2,1);mesh(x0,y0,u(:,:,101);%作t=1.0时刻的数值解title(t=1.0时刻的数值解);xlabel(x变量);ylabel(y变量);zlabel(u值);subplot(2,2,2);mesh(x0,y0,uu(:,:,101);%作t=1.0时刻的精确解title(t=1.0时刻的精确解);xlabel(x变量);ylabel(y变量);zlabel(u值);%作t=1.4时刻的数值解与精确解subplot(2,2,3);mesh(x0,y0,u(:,:,141);%作t=1.4时刻的数值解title(t=1.4时刻的数值解);xlabel(x变量);ylabel(y变量);zlabel(u值);subplot(2,2,4);mesh(x0,y0,uu(:,:,141);%作t=1.4时刻的精确解title(t=1.4时刻的精确解);xlabel(x变量);ylabel(y变量);zlabel(u值);五实验结果及实例分析1、时刻的数值解与精确解图 图1 t=0.1、0.5时刻的数值解、精确解 图2 t=1.0、1.4时刻的数值解、精确解注:上两图为四个时刻的数值解与精确解,三层显格式达二阶收敛,不难看出,收敛效果很好,符合理论。下图是四个时刻的绝对误差图像,从图中看出,绝对误差较小,且经过计算得到,收敛阶近似于2,正好符合理论值。2、时刻的绝对误差图图3 四个时刻的绝对误差3、四个时刻(t=0.1、0.5、1.0、1.4)的绝对误差表t=0.1时刻的绝对误差0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002 0.0002 0.0002 0.0002 0.0001 0.0001 0.0000 0.0000 0.0001 0.0003 0.0004 0.0004 0.0005 0.0004 0.0004 0.0003 0.0001 0.0000 0.0000 0.0002 0.0004 0.0005 0.0006 0.0006 0.0006 0.0005 0.0004 0.0002 0.0000 0.0000 0.0002 0.0004 0.0006 0.0007 0.0007 0.0007 0.0006 0.0004 0.0002 0.0000 0.0000 0.0002 0.0005 0.0006 0.0007 0.0008 0.0007 0.0006 0.0005 0.0002 0.0000 0.0000 0.0002 0.0004 0.0006 0.0007 0.0007 0.0007 0.0006 0.0004 0.0002 0.0000 0.0000 0.0002 0.0004 0.0005 0.0006 0.0006 0.0006 0.0005 0.0004 0.0002 0.0000 0.0000 0.0001 0.0003 0.0004 0.0004 0.0005 0.0004 0.0004 0.0003 0.0001 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002 0.0002 0.0002 0.0002 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 t=0.5时刻的绝对误差0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0007 0.0013 0.0018 0.0021 0.0022 0.0021 0.0018 0.0013 0.0007 0.0000 0.0000 0.0013 0.0025 0.0034 0.0040 0.0042 0.0040 0.0034 0.0025 0.0013 0.0000 0.0000 0.0018 0.0034 0.0047 0.0055 0.0058 0.0055 0.0047 0.0034 0.0018 0.0000 0.0000 0.0021 0.0040 0.0055 0.0065 0.0068 0.0065 0.0055 0.0040 0.0021 0.0000 0.0000 0.0022 0.0042 0.0058 0.0068 0.0071 0.0068 0.0058 0.0042 0.0022 0.0000 0.0000 0.0021 0.0040 0.0055 0.0065 0.0068 0.0065 0.0055 0.0040 0.0021 0.0000 0.0000 0.0018 0.0034 0.0047 0.0055 0.0058 0.0055 0.0047 0.0034 0.0018 0.0000 0.0000 0.0013 0.0025 0.0034 0.0040 0.0042 0.0040 0.0034 0.0025 0.0013 0.0000 0.0000 0.0007 0.0013 0.0018 0.0021 0.0022 0.0021 0.0018 0.0013 0.0007 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 t=1.0时刻的绝对误差0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0016 0.0031 0.0043 0.0051 0.0053 0.0051 0.0043 0.0031 0.0016 0.0000 0.0000 0.0031 0.0059 0.0082 0.0096 0.0101 0.0096 0.0082 0.0059 0.0031 0.0000 0.0000 0.0043 0.0082 0.0113 0.0132 0.0139 0.0132 0.0113 0.0082 0.0043 0.0000 0.0000 0.0051 0.0096 0.0132 0.0156 0.0164 0.0156 0.0132 0.0096 0.0051 0.0000 0.0000 0.0053 0.0101 0.0139 0.0164 0.0172 0.0164 0.0139 0.0101 0.0053 0.0000 0.0000 0.0051 0.0096 0.0132 0.0156 0.0164 0.0156 0.0132 0.0096 0.0051 0.0000 0.0000 0.0043 0.0082 0.0113 0.0132 0.0139 0.0132 0.0113 0.0082 0.0043 0.0000 0.0000 0.0031 0.0059 0.0082 0.0096 0.0101 0.0096 0.0082 0.0059 0.0031 0.0000 0.0000 0.0016 0.0031 0.0043 0.0051 0.0053 0.0051 0.0043 0.0031 0.0016 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 t=1.4时刻的绝对误差0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0002 0.0003 0.0005 0.0006 0.0006 0.0006 0.0005 0.0003 0.0002 0.0000 0.0000 0.0003 0.0007 0.0009 0.0011 0.0011 0.0011 0.0009
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 难点详解人教版八年级上册物理声现象《声音的特性声的利用》必考点解析试卷(解析版)
- 2025年中国微晶玻璃板行业市场分析及投资价值评估前景预测报告
- 2025年中国椭圆桌子行业市场分析及投资价值评估前景预测报告
- 2024年中山市检察机关招聘司法辅助人员考试真题
- 考点解析人教版九年级物理《生活用电》综合练习试题(含详解)
- 达标测试人教版八年级上册物理声现象《声音的特性声的利用》定向测试试卷(含答案详解版)
- 2025国考白山市党务工作岗位申论题库含答案
- 2025国考天津市阿拉伯语翻译岗位申论高频考点及答案
- 2025国考荆州长航公安行测资料分析题库含答案
- 2025国考白山市政务公开岗位申论模拟题及答案
- 2025-2026学年(人教版)初中数学七年级上册第一次月考 (1-2章)(含答案)
- 2025天津泰达热电能源管理有限公司及所属企业招聘12人考试参考题库及答案解析
- 有机半导体课件
- 城乡燃气管道维护保养技术方案
- 第二单元第四课实现灯光警报呼叫器-课件
- 水土保持治理工应急处置考核试卷及答案
- 初中学生心理健康辅导手册
- 仓库搬运安全培训内容课件
- 气管切开的患者如何护理
- 血透高钾患者个案护理
- 【基于财务状况质量分析理论的京东财务报表分析案例18000字(论文)】
评论
0/150
提交评论