• 解决争端根本办法是建立机制 而不是利用社会舆论取得某种优势(原创首发) 2019-03-25
  • 用药晚一步,起效来不及 2019-03-02
  • [新华网]重庆:深化投融资改革 促进企业去杠杆 2019-03-02
  • 他帖子里的明理,就是要人们放弃自己的利益,一切顺从别人的指挥棒转 2019-02-25
  • 运-20完成首次重装空投?为什么说这意义重大 2019-02-14
  • (两会受权发布)最高人民检察院检察长简历 2019-02-14
  • 习近平齐鲁之行的七个感人瞬间 2019-02-05
  • 全国首起公益诉讼调解结案案件:被告同意全部诉讼请求 2019-02-05
  • 大数据,贵在“融”和“用” 2019-01-20
  • 湖州市启动中小学生校外培训机构治理 2019-01-20
  • 牙生·司地克调研我市水资源管理保护工作 2019-01-12
  • 我看“支付宝回收垃圾”这件事不错,应该支持。[微笑][微笑] 首先是提高的回收效率,其次便于集中处理旧物品,防止污染有利。 2019-01-02
  • 楼市下半年或持续降温 房地产长效机制加速推进 2018-12-24
  • 人的本质的演变规律:从原始母系氏族社会的公有者经过父系氏族社会私有和公有双重所有者而演变为私有制阶级社会的私有者,然后经过现代社会公有和私有双重所有者... 2018-12-24
  • 《不起眼女主角培育法》宣布将推出剧场版动画 2018-12-14
  • 分分彩后三做号技巧:抛物型方程有限差分方法的应用- 报告讲诉

    来源:互联网 由 盛开伤花 贡献 责任编辑:鲁倩  
    2015 年秋季学期研究生课程考核

    分分彩软件手机版 www.xpmw.net (读书报告、研究报告)

    考核科目:偏微分方程数值解法

    学生所在院(系):理学院数学系

    学生所在学科:数学

    学生姓名:H i t e r

    学号:1X S012000

    学生类别:

    考核结果阅卷人

    抛物型方程有限差分方法的应用

    摘要

    抛物型偏微分方程是一类比较重要的偏微分方程。热传导方程是最简单的一种抛物型方程。热传导方程研究的是热传导过程的一个简单数学模型。根据热量守恒定律和傅里叶热传导实验定律可以导出热传导方程。在本篇论文中,将先详述抛物型偏微分方程的有限差分法的相关知识,然后给出抛物型方程的两个具体的应用实例。

    关键字:抛物型方程,差分格式,应用

    Abstract

    Parabolic partial differential equation is a kind of important partial differential equation. The heat conduction equation is one of the simplest parabolic equations. The heat conduction equation is a simple mathematical model of the heat conduction process. Heat conduction equation is derived based on the law of conservation of heat and Friyege's law of conduction. In this thesis, we first give a detailed knowledge of the finite difference method for parabolic partial differential equations, and then give two specific examples of the application of the parabolic equation.

    Keywords: parabolic equation, difference scheme, application

    0 前言

    抛物型方程是偏微分方程中的三大方程(另两种为双曲型方程和椭圆型方程)之一,如何去研究抛物型方程的性质在《偏微分方程数值解法》的课程中占有很大的比例。如果要研究抛物型方程,我们一般从以下几个方面来研究,分别是:定界问题、格林函数、极值原理、解的正则性、抛物方程、拟线性蜕化和反应扩散方程。但是由于所学知识有限,所以我们在此只简单的介绍抛物型方程的有限差分法并给出两个应用实例。

    1 抛物型方程有限差分法

    1.1 简单差分法

    考虑一维模型热传导方程

    )(22x f x

    u

    a t u +??=??,T t ≤<0 (1.1) 其中a 为常数。)(x f 是给定的连续函数。(1.1)的定解问题分两类:

    第一,初值问题(Cauchy 问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初始条件:

    ()()x x u ?=0,,∞<<∞-x (1.2)

    第二,初边值问题(也称混合问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初

    始条件:

    ()()x x u ?=0,,l x l <<- ()13.1

    及边值条件

    ()()0,,0==t l u t u ,T t ≤≤0 ()23.1

    假定()x f 和()x ?在相应的区域光滑,并且于()0,0,()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。

    现在考虑边值问题(1.1),(1.3)的差分逼近。取N l h =为空间步长,M

    T

    =τ为时间步长,其中N ,M 是自然数,

    jh x x j ==, ()N j ,,1,0 =; τk y y k ==, ()M k ,,1,0 =

    将矩形域G {}T t l x ≤≤≤≤=0;0分割成矩形网格。其中 ()

    j i y x ,表示网格节点;

    h G 表示网格内点(位于开矩形G 中的网格节点)的集合;

    h G 表示位于闭矩形G 中的网格节点的集合; h Γ表示h G -h G 网格边界点的集合;

    k j u 表示定义在网点()k i t x ,处的待求近似解,N j ≤≤0,M k ≤≤0。

    注意到在节点()k i t x ,处的微商和差商之间的下列关系((,)k

    j k j

    u u x t t t ????

    ≡ ?????): ()()

    ()ττ

    O t u t x u t x u k

    j k j k j +???

    ????=-+,,1 ()()

    ()

    2112,,ττ

    O t u t x u t x u k

    j

    k j k j +???

    ????=--+ ()()()h O x u h t x u t x u k

    j k j k j +??? ????=-+,,1

    ()()

    ()h O x u h

    t x u t x u k

    j k j k j +???

    ????=--,,1 ()()

    ()

    2112,,h O x u h

    t x u t x u k j

    k j k j +???

    ????=--+

    ()()()

    ()

    2

    222

    11,,2,h O x u h

    t x u t x u t x u k

    j

    k j k j k j +???? ????=+--+ 可得到以下几种最简差分格式

    (一)向前差分格式

    =-+τ

    k j

    k j u u 1j k

    j k j k j f h

    u u u a

    ++--+2

    1

    12,()()j j x f f = ()14.1

    ()j j j x u ??==0,k u 0=k

    N u =0 (

    )24.1 其中1,,1,0-=N j ,1,,1,0-=M k 。取2

    h a r τ

    =

    为网比,则进一步有 1+k j u =k j ru 1++()r 21-k j u +k

    1-+j f τ

    ()14.1' 此差分格式是按层计算:首先,令0=k ,得到

    1j u =01+j ru +(

    )r 21-0j u +0

    u +j f τ 于是,利用初值()

    j j j x u ??==0

    和边值k u 0=k N u =0,可算出第一层的1

    1,,1,0-=N j 。

    再由()14.1'

    取1=k ,可利用1j u 和k u 0=k N u =0算出2

    1,,1,0-=N j 。如此下去,即可

    逐层算出所有k

    1,,1,0-=N j ,1,,1,0-=M k )。

    由于第()1+k 层值可以通过第()k 层值直接得到,如此的格式称为显格式。并视k

    j u 为

    ()k j t x u ,的近似值。

    若记

    ()

    T

    k N k k k u u u 1

    -= u ,()()()()T N x x x 1-=???? ,()()()()T N x f x f x f 121,,,-=τττ f 则显格式()14.1'

    可写成向量形式

    ???=-=+=+?

    11

    ,,1,0,u f Au u M k k k 其中

    ???

    ???

    ?

    ? ??----=r r r r r r r r

    r r

    21002100210021

    A

    若记

    22x

    u a t u Lu ??-??=

    ()

    --=

    k j

    k j k

    j

    h u u u L 112

    1

    12h

    u u u a

    k

    j k j k j -++-

    那么截断误差为

    ()=u R k

    j

    ()

    ()[]k j

    Lu t x u L -,1=()ττO t x t u r k j +?

    ??

    ? ??????? ??--)~,~(2112122=()2

    h O +τ (1.5) 其中(,)j k x t 是矩形11+-<u R k

    j k

    j

    x u ????

    ????222τ+()

    2

    τO k

    j

    x u

    h a ???? ?????-442?12 =

    k

    j

    x u ????

    ????222τ+()

    2

    τO ()

    2

    222

    2?1

    12τO t u a h a k

    j

    +???? ??????

    ?- =??????-?-211212ττa h ()2

    22~τO t u k

    j +????

    ???? =?

    ????

    ?--21121r τ()222~τO t u k

    j

    +???? ????=(

    )

    2

    h O +τ。 这里

    ???

    ? ??????=??2222

    44x u x

    a x u a ??? ???????=t u a x a 122???

    ??????=t u x 22t x u ???=23 22t u ????? ??????=t u t 2???

    ? ??????=22x u a t t x u ???=23 故22t u ??44244x u a x u a a ??=???

    ? ?????=,从而=??44x u 221t u a ??? (二)向后差分格式

    =-+τ

    k j

    k j u u 1

    j k j k j k j f h

    u u u a

    ++-+-+++2

    1

    1

    1112 ()()j j x f f = ()16.1

    ()j

    j j x u ??==0, k u 0=k

    N u =0 ()26.1 其中 1,,1,0-=N j ,1,,1,0-=M k 。取2h

    a r τ

    =

    为网比,则进一步有

    抛物方程的有限差分法

    抛物型方程有限差分法

    抛物型方程有限差分法_数学_自然科学_专业资料。抛物型方程有限差分法抛物方程差分法的构造在空间方向上与椭圆方程类似, 在时间方向上用一阶 差商代替代替一阶微...

    有限差分法求解抛物型方程说明

    有限差分法求解抛物型方程说明_数学_自然科学_专业资料。有限差分法求解抛物型方程偏微分方程只是在一些特殊情况下, 才能求得定解问题解的解析式, 对比较复杂的...

    抛物形扩散方程的有限差分法及数值实例_图文

    :抛物形扩散方程的有限差分法及数值实例 学生姓名: 向聘 抛物形扩散方程的有限差分法及数值实例 1.1 抛物型扩散方程抛物型偏微分方程是一类重要的偏微分方程。...

    大连理工大学 高等数值分析 抛物型方程有限差分法...

    大连理工大学 高等数值分析 抛物型方程有限差分法_物理_自然科学_专业资料。抛物型方程有限差分法 抛物方程差分法的构造在空间方向上与椭圆方程类似, 在时间方向上...

    二维抛物方程的有限差分法

    有限差分法求解偏微分方程

    有限差分法求解偏微分方程摘要: 本文主要使用有限差...通过对机械系统应用自然的基本定律和原 理得到带有...4 = 0,则称为抛物型方程; 如果Δ = 2 ? 4 ...

  • 与《抛物型方程有限差分方法的应用- 报告讲诉》相关:
  • 抓时间补弱科不熬夜抓效率
  • 如何利用假期提高中考作文水平?(图)
  • 2016高考生必看:每月高考大事备忘录
  • 备考2016高考:只有1%的高三学生知道这样做
  • 新高三开学了,艺术生该如何开始各科复习?
  • 作为一个过来人写给高三党的建议
  • 高三生的学习和生活应该是这样的!
  • 166中学高考文科状元张子秋:高三生要努力、自信
  • 高考微问答81期:新高三生怎么开始复习?
  • 高三考生培养十大学习习惯 可上名校
  • 新高三复习计划 助你快速提分
  • 2016年高考复习时的六点注意事项
  • 2016年高考改革下各科目备考攻略
  • 新高三必读 成绩提升应循序渐进
  • 高三新学期刚开始 学生别急着打疲劳战
  • 本站网站首页首页教育资格全部考试考试首页首页考试首页职业资格考试最近更新儿童教育综合综合文库22文库2文库作文总结建筑资料库考研建筑专业资料考试首页范文大全公务员考试首页英语首页首页教案模拟考考试pclist爱学首页日记语文古诗赏析教育教育资讯1高考资讯教育头条幼教育儿知识库教育职场育儿留学教育高考公务员考研考试教育资讯1问答教育索引资讯综合学习网站地图学习考试学习方法首页14托福知道备考心经冲刺宝典机经真题名师点睛托??纬?/a>雅思GREGMATSAT留学首页首页作文
    免责声明 - 关于我们 - 联系我们 - 广告联系 - 友情链接 - 分分彩软件手机版 - 分分彩软件手机版
    Copyright © 2017 分分彩软件手机版 www.xpmw.net All Rights Reserved
  • 解决争端根本办法是建立机制 而不是利用社会舆论取得某种优势(原创首发) 2019-03-25
  • 用药晚一步,起效来不及 2019-03-02
  • [新华网]重庆:深化投融资改革 促进企业去杠杆 2019-03-02
  • 他帖子里的明理,就是要人们放弃自己的利益,一切顺从别人的指挥棒转 2019-02-25
  • 运-20完成首次重装空投?为什么说这意义重大 2019-02-14
  • (两会受权发布)最高人民检察院检察长简历 2019-02-14
  • 习近平齐鲁之行的七个感人瞬间 2019-02-05
  • 全国首起公益诉讼调解结案案件:被告同意全部诉讼请求 2019-02-05
  • 大数据,贵在“融”和“用” 2019-01-20
  • 湖州市启动中小学生校外培训机构治理 2019-01-20
  • 牙生·司地克调研我市水资源管理保护工作 2019-01-12
  • 我看“支付宝回收垃圾”这件事不错,应该支持。[微笑][微笑] 首先是提高的回收效率,其次便于集中处理旧物品,防止污染有利。 2019-01-02
  • 楼市下半年或持续降温 房地产长效机制加速推进 2018-12-24
  • 人的本质的演变规律:从原始母系氏族社会的公有者经过父系氏族社会私有和公有双重所有者而演变为私有制阶级社会的私有者,然后经过现代社会公有和私有双重所有者... 2018-12-24
  • 《不起眼女主角培育法》宣布将推出剧场版动画 2018-12-14