没有合适的资源？快使用搜索试试~ 我知道了~

文库首页行业研究行业报告二维波动方程的有限差分法.zip_二维波动_二维波动方程_差分方程_波动方程_波动方程差分

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

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

0 下载量 92 浏览量 2024-04-30 21:56:27 上传评论收藏 1.91MB DOC 举报

温馨提示

试读

11页

二维波动方程的有限差分法，与解析解进行了误差比对

资源推荐

资源详情

资源评论

学生实验报告

实验课程名称偏微分方程数值解

开课实验室数统学院

学院数统年级 2013 专业班信计 02 班

学生姓名学号

开课时间 2015 至 2016 学年第 2 学期

总成绩

教师签名

数学与统计学院制

开课学院、实验室：数统学院实验时间： 2016 年 6 月 20 日

实验项目类型

实验项目

名称

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

验证

演示

综合

设计

其他

指导教师

曾芳

成绩

是

一．实验目的

通过该实验，要求学生掌握求解二维波动方程的有限差分法，并能通过计算机语言编程

实现。

二．实验内容

考虑如下的初值问题：

( ) ( ) ( )

( ) ( ) ( ) ( )

( ) ( )

[ ]

2 2 2

, , 0,1 , 0,1.4

, , 0 sin sin , , , 0 0, , 0,1

, , 0, , , 0,1.4

u u u

x y t

t x y

u x y x y u x y x y

u x y t x y t

p p

¶ ¶ ¶

= + Î W = Î

¶ ¶ ¶

= = Î

= Î ¶W Î

(1)

1．在第三部分写出问题（1）三层显格式。

2．根据你写出的差分格式，编写有限差分法程序。将所写程序放到第四部分。

3．取

0.1, 0.1h h

= =

，分别将

0.5,1.0,1.4t =

时刻的数值解画图显示。

4. 该问题的解析解为

( )

, , t cos 2 sin sinu x y t x y

p p p

，将四个时刻的数值解的误差画图显示，

对数值结果进行简单的讨论。

三．实验原理、方法（算法）、步骤

网格划分

0.1, 0.1h h

= =

，故

1 1.4

10, 140N M

= = = =

，

, , , 0,1, ,10

i j

x ih y jh i j= = = L

，

t k

，

0,1, ,140k = L

。在内网点

( )

, ,

i j k

x y t

，利用二阶中心差商，对（1）建立差分格式：

1 1

, , , 1, , 1, , 1 , , 1

2 2 2

k k k k k k k k k

i j i j i j i j i j i j i j i j i j

u u u u u u u u u

h h

+ -

+ - + -

- + - + - +

= +

（2）

整理得到：

( ) ( )

1 2 2 1

, 1, 1, , 1 , 1 , ,

2 4

k k k k k k k

i j i j i j i j i j i j i j

u r u u u u r u u

+ -

+ - + -

= + + + + - -

（3）

其中，

, 1, 2, ,9, 1, 2, ,139i j k= =L L

，网比

0.1r

= =

，局部截断误差为

( )

2 2

o h

。

考虑边界条件

( ) ( )

[ ]

, , 0, , , 0,1.4u x y t x y t= Î ¶W Î

，差分格式为：

0,0 0,N ,0 ,

0, 0,1, ,140

k k k k

N N N

u u u u k= = = = = L

（4）

考虑初始条件

( )

, , 0 sin sinu x y x y

p p

，差分格式为：

( )

( ) ( )

sin sin sin sin , , 0,1, ,10

i j i j

u x y ih jh i j

p p p p

= = = L

（5）

考虑初始条件

( ) ( ) ( )

, , 0 0, , 0,1

u x y x y= Î

，利用二阶差商近似：

1 1

, ,

0, , 0,1, ,10

i j i j

u u

i j

= = L

（6）

设

0k =

时刻的点为内点，则满足差分格式（2），代入上式得到：

( ) ( )

1 2 0 0 0 0 2 0 1

, 1, 1, , 1 , 1 , ,

2 4

i j i j i j i j i j i j i j

u r u u u u r u u

+ - + -

= + + + + - -

（7）

将（6）得到的结果

1 1

, ,i j i j

u u

代入（7）中，整理得到：

( ) ( )

1 2 0 0 0 0 2 0

, 1, 1, , 1 , 1 ,

1 2

i j i j i j i j i j i j

u r u u u u r u

+ - + -

= + + + + -

（8）

综上（2）、（4）、（5）、（8）得到三层显格式的差分格式为：

( ) ( )

( )

( ) ( )

1 2 2 1

, 1, 1, , 1 , 1 , ,

0,0 0,N ,0 ,

1 2 0 0 0

, 1, 1, , 1

2 4

, 1, 2, ,9, 1, 2, ,139

0, 0,1, ,140

sin sin sin sin , , 0,1, ,10

k k k k k k k

i j i j i j i j i j i j i j

k k k k

N N N

i j i j

i j i j i j i j i

u r u u u u r u u

i j k

u u u u k

u x y ih jh i j

u r u u u u

p p p p

+ -

+ - + -

+ - +

= + + + + - -

= =

= = = = =

= = =

= + + +

L L

( ) ( )

0 2 0

, 1 ,

1 2 , , 0,1, ,10

j i j

r u i j

+ - =

（9）

其中

0.1r

= =

，局部截断误差为

( )

2 2

o h

。

四．实验环境（所用软件、硬件等）及实验数据文件

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);%空间网格剖分

剩余10页未读，继续阅读

评论收藏

内容反馈

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

资源评论

资源反馈

评论星级较低，若资源使用遇到问题可联系上传者，3个工作日内问题未解决可申请退款~

手把手教你学AI

粉丝: 8168
资源: 4710

上传资源快速赚钱

我的内容管理展开

我的资源快来上传第一个资源

我的收益

登录查看自己的收益

我的积分登录查看自己的积分

我的C币登录后查看C币余额

我的收藏

我的下载

下载帮助

前往需求广场，查看用户热搜

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

Two dimensional explosion wave simulation.zip_wave_二维波动方程_二维波动方程

五点差分法(matlab)解椭圆型偏微分方程.zip_wudianchafenfa_五点差分_五点差分例题_偏微分_椭圆方程

matlab leinuo.zip_matlab雷诺方程_雷诺_雷诺方程_雷诺方程matlab_雷诺方程求解

基于matlab的热传导多介质有限差分法.zip

使用有限差分法和稳态解的二维热方程附matlab代码.zip

相关实用应用程序（Windows可用）

李飞飞自传 我看见的世界 The World I see

ChatGPT使用总结：150个ChatGPT提示词模板（完整版）

chromedriver-win64.zip

全国计算机二级WPSoffice精选350道选择题题库（含答案）.pdf

哈尔滨工业大学-ChatGPT调研报告-2023.3.6-94页.pdf

智联招聘：2024年大学生就业力调研报告.pdf

第十九届研电赛-技术论文模板

4个亲测好用的ChatGPT4渠道

学术海报模板+论文科研+研究生

北森能力测评题库.zip

农村公交与异构无人机协同配送优化

车载毫米波雷达DOA估计综述博文仿真代码

ST-LINK Utility 4.6.0

认知智能技术与产业研究报告2023

python大作业 含爬虫、数据可视化、地图、报告、及源码（2016-2021全国各地区粮食产量）.rar

2023泛娱乐社交出海手册-ZEGO即构科技

研究生电子设计竞赛论文集

1000份ppt模版，PPT模板优秀PPT

由于找不到iUtils.dll,无法继续执行代码

249个 ChatGPT 关键词汇总 中文版

ST语言规则编程手册全面讲解ST语言

【R213】The Worlds I See 我看见的世界【Fei-Fei_Li 李飞飞】.pdf

GJB 3206B-2022 《技术状态管理》

最新资源

李飞飞自传我看见的世界 The World I see

python大作业含爬虫、数据可视化、地图、报告、及源码（2016-2021全国各地区粮食产量）.rar

249个 ChatGPT 关键词汇总中文版