聪明文档网

    聪明文档网

    最新最全的文档下载
    当前位置： 首页> 本科毕业设计--求解热传导方程的高精度隐式差分格式

    本科毕业设计--求解热传导方程的高精度隐式差分格式

    时间：    下载该word文档




    新疆大学毕业论文(设计


    题目:求解热传导方程的高精度隐式差分格式所属院系：数学与系统科学学院专业：信息与计算科学

    I




    声明

    本人郑重声明该毕业论文（设计）是本人在开依沙尔老师指导下独立完成的，本人拥有自主知识产权，没有抄袭、剽窃他人成果，由此造成的知识产权纠纷由本人负责。

    声明人（签名）：
    年月日

    亚库甫江.买买提同学在指导老师的指导下，按照任务书的内容，独立完成了该毕业论文（设计），指导教师已经详细审阅该毕业论文（设计）。

    指导教师（签名）：
    年月日


    II




    新疆大学
    毕业论文（设计）任务书

    班级：信计07-2姓名：亚库甫江.买买提论文（设计）题目：求解热传导方程的高精度隐式差分格式专题：毕业设计论文（设计）来源：教师自拟要求完成的内容：学习和掌握一维热传导方程已有的各种差分
    格式的基础上，扩散方程对空间变量应用紧致格式离散，对时间变量应用梯形方法，构造热传导方程的精度为O2h4数值格式，讨论格式的稳定性，最后数值例子来验证。
    发题日期：2012年12月25日完成日期：2012年5月28日实习实训单位：数学学院地点：数学学院论文页数：19页；图纸张数：4指导教师：开依沙尔老师教研室主任
    院长（系主任）

    III




    摘要

    本文首先对热传导方程经典差分格式进行复习和讨论，然后热传导方程对空间变量四阶紧致格式进行离散,时间变量保持不变，把一维热传导方程转化为常微分方程组的初值问题,再利用梯形方法构造热传导方程方程的时间二阶空间四阶精度的一种差分格式，并稳定性进行分析，数值结果与Crank-Nicholson格式进行比较，数值结果表明,该方法是有效求解热传导方程的数值计算.
    关键词:热传导方程,高精度紧致格式;梯形方法；两层隐格式;Crank-Nicolson格式

    ABSTRACT

    Thispaperfirststudyonsomeclassicalfinitedifferencefortheheatconductionequation,secondelysecondelyweapplycompactfinitedifferenceapproximationoffourthorderfordiscretizingspatialderivativesbutleavethetimevariableContinuous.ThisapproachresultsinasystemofODEs,whichcanthenbeusedtrapezodialformuladerivedfourthorderinspaceandsecondorderintimeunconditionallystableimplicitscheme.thestabilityandlocaltruncationerroroftheobtainedmethodareanalysied.NumericalexperimentsshowsthatthismethodUseful,efficientmethodforsolvingdiffusionequation

    Keywords:Heatconductioneqution;Higher-odercompactscheme;Trapezodialformula;Two-levelimplictscheme;Crank-Nicolsonscheme



    IV





    目录

    引言...................................................................................................................................................1预备知识...........................................................................................................................................21.扩散方程的经典差分格式............................................................................................................31.1显式差分格..........................................................................................................................31.1.1显式的截断误差........................................................................................................41.1.2显式差分格式的稳定性.............................................................................................41.2隐式差分格式........................................................................................................................51.2.1隐式差分格式的截断误差.........................................................................................51.2.2隐式差分格式的稳定性….........................................................................................61.3Crank-Nicolson格式….....................................................................................................61.3.1Crank-Nicolson差分格式的截断误差....................................................................71.3.2Crank-Nicolson差分格式的稳定性…....................................................................82.高精度格式的构造…....................................................................................................................92.1梯形方法…….........................................................................................................................92.2本文格式的构造....................................................................................................................102.3稳定性分析….......................................................................................................................113.数值实验…...................................................................................................................................13结论……..........................................................................................................................................17致谢……..........................................................................................................................................18参考文献…......................................................................................................................................19

    V




    引言
    热传导方程是一类描述物理量随时间的扩散和衰减规律的抛物型微分方程．自然环境、工程设备及生物机体中的许多物理现象，诸如气体的扩散、液体的渗透、热的传导、以及半导体材料中杂质的扩散等都可用热传导方程来描述．由于物理问题本身的复杂性，其精确解往往不容易求得，因此研究其数值求解方法无疑具有非常重要的理论意义和工程应用价值【1】．
    求解该问题的数值方法主要有差分法、有限元法、边界元法等，其中有限差分方法数值求解扩散方程的应用广泛的有效地方法之一。目前求解该问题的主要的差分格式有显式格式，隐式格式，Crank-Nicolson格式等[1,2,4]。虽然显式格式计算简单，但是稳定性有所限制，一般隐式格式和Crank-Nicolson格式分别为一阶和二阶精度的绝对稳定的隐式格式，还显得误差阶不够高，得到的结果也往往不能令人满意,考虑到这些不足文[7]中半离散方法构造O2h2格式结果Crank-Nicolson格式进行比较，在文[10]待定参数法构造精度O3h6的显式格式但是稳定性条件比较苛刻，它文的稳定性条件为r
    a1
    ，本文热2h6
    传导方程对空间变量应用紧致格式离散，对时间变量应用梯形方法，构造热传导方程的精度为O2h4的绝对稳定的隐式差分格式，并讨论了稳定性，数值值结果与经典Crank-Nicholson格式进行比较，数值结果表明,该方法是有效求解扩散方程的数值计算.
    本文分为三大部分，第一部分简单介绍热传导方程的经典差分格式，第二部分主要介绍热传导方程的高精度格式的构造和稳定性，第三部分给出具体的数值算例，结果与Crank-Nicolson格式，准确值进行比较，最后给出结论。


    1




    预备知识
    利用下面的各种数值微分公式得到不同的差分格式
    uxj,tn1uxj,tn

    2
    uxj1,tnuxj,tn
    h
    uxj,tnuxj1,tn
    h2h
    h2
    uOtju2Otj
    nn
    n
    uxj,tn1uxj,tn1
    uOhxjuOhxju2Ohxj
    2u2
    Oh2xj
    n
    nn

    uxj1,tnuxj1,tn
    uxj1,tn2uxj,tnuxj1,tn
    截断误差：一般说来，微分方程的解不会精确地满足差分方程。将差分方程中的各个项同时用微分方程的解在相应点的值代入，利用泰勒展开，就会得到一个误差项，这个误差项就是截断误差。
    相容性：若时间步长以及空间步长h同时趋于0，截断误差Rnj0，就说差分格式与微分方程是相容的。一个差分格式与一个微分方程相容，则表明当
    ,h0时，差分算子与微分算子对任一光滑函数的作用是相同的，所以可用相
    容的差分格式近似相应的微分方程，而截断误差则是对这一近似程度的一个度量。
    收敛性：考察差分格式在理论上的准确解能否任意逼近微分方程的解。如
    n
    果当时间步长以及空间步长h趋于0时，enju(xj,tnuj0，我们称差分
    格式是收敛的，即时间步长以及空间步长h趋于0时，差分格式的解逼近于微分方程的解。
    1
    稳定性：差分格式的计算是逐层计算的，计算第n1层上的unj时，要用
    n1
    到第n层上计算出来的结果。计算unj时的舍入误差，必然会影响到uj的值，从
    而就要分析这种误差传播的情况。因此，一个有实用价值的数值方法应该具有能够控制这种误差影响的性能，这就是数值方法的稳定性。
    精度：如果一个差分格式的截断误差EO(qhp，就说差分格式对时间t是q阶精度的，对空间x是p阶精度的。
    Lax等价定理[5]：给定一个适定的线性初值问题以及与其相容的差分式，则差分格式的稳定性是差分格式收敛性的充分必要条件。
    2




    定理1（vonNeumann条件）微分方程的差分格式稳定的必要条件是当0，
    nT,对所有kR有
    j(G(,k1M，j1,2,,p
    其中G(,k为增长因子（或增长矩阵），j(G(,k表示G(,k的特征值，M为常数。
    定理2如果差分格式的增长矩阵G(,k是正规矩阵，则vonNeumann条件是差分格式稳定的必要且充分条件。
    推论2.1当G(,k为实对称矩阵，酉矩阵，Hermite矩阵时，vonNeumann条件是差分格式稳定的充分必要条件。
    推论2.2当p1时，即G(,k只有一个元素，则vonNeumann条件是差分格式稳定的充要条件。
    定理3如果存在常数K,0使得
    G(,k1K，00，则差分格式是稳定的。
    1.热传导方程的经典差分格式考虑一维热传导方程的初边界问题：
    u2u
    0,2tx
    u(x,0f(x,u(a,t(t,u(b,t(t,
    axb,t0
    axb
    t0t0

    1.1显式差分格式
    uxj,tn1uxj,tnunu我们可以对用向前差分Ottj2u
    用二阶差商x2
    uxj1,tn2uxj,tnuxj1,tn
    h2
    2u
    2Oh2xj
    n
    得到差分格式为
    1ununjj



    nn
    unj12ujuj1
    h2
    0
    (1.1.1
    3




    1.1.1显式差分格式的截断误差
    证：（用taylor展开）
    un22un1
    u(x,ttu(xj,tn[]j[2]j,(011
    t2!t
    uu(xj,tnu
    n
    jn1j
    u
    nj1
    unh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j2
    x2!x3!x4!x
    (021unh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j3
    x2!x3!x4!x
    (031
    u
    nj1
    把上述代入差分格式中，得截断误差为：
    un2un12unh24unh24un
    k(x,t[]j[2]j{[2]j[4]j2[4]j3}
    t2tx24x24x
    n
    j
    un2un2un1h24un4un
    {[]j[2]j}{[2]j{[4]j2[4]j3}}
    tx2t24xxo(h2.(01,2,31
    2从上述可知，截断误差为kn(x,to(h，它对空间方向为一阶截断误差而对j
    时间方向为二阶截断误差。
    1.1.2显式差分格式的稳定性
    证：先把差分格式公式(1.1.1改写为：
    1nnnn
    unu(u2uujjj1jj1其中
    
    h2

    nikjh
    利用稳定性的Fourier方法，令un，并将它代入上式就得到jve
    vn1eikjhvneikjh(eikj2eikjvn
    消去共因子有
    vn1[1(eikj2eikj]vn
    由此得到增长因子
    G(,k1(eikj2eikj(12(eikjeikj

    (122coskh12(1coskh
    kh|12

    hkkh
    114sin2104sin22
    22G(,k||12(1coskh||14sin2
    4




    因为k>0,h>0且0，所以必然左边成立，则右边为
    2sin2
    khkh
    1且sin2122
    21
    1
    2
    显然这个格式是相容的。它在它是条件收敛的（收敛条件1）。
    2
    1
    时稳定的，因为按照Lax定理可知；2
    1.2隐式差分格式
    2uu
    我们可以对用向后差分，2用二阶差商，得到差分格式为：
    tx
    n1
    unujj


    nn
    un2*uuj1jj1
    h
    2
    0(1.2.1
    1.2.1隐式差分格式的截断误差
    证：（用taylor展开）
    un22un1
    u(x,ttu(xj,tn[]j[2]j,(011
    t2!t
    uu(xj,tnu
    n
    jn1j
    u
    nj1
    unj1
    unh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j2
    x2!x3!x4!x
    (021
    unh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j3
    x2!x3!x4!x
    (031
    把上述代入差分格式中，得截断误差为：
    un2un12unh24unh24un
    k(x,t[]j[2]j{[2]j[]j2[]j3}
    t2tx24x424x4
    n
    j
    un2un2un1h24un4un
    {[]j[2]j}{[2]j{[4]j2[4]j3}}
    tx2t24xxo(h2.(01,2,31

    h2，从上述可知，截断误差为kn它对时间方向为一阶截断误差，j(x,to(
    而对空间为二阶截断误差。

    5




    1.2.2隐式差分格式的稳定性
    证:先用差分格式（1.2.1）为：
    1nnnununjj(uj12ujuj1其中
    
    h
    2

    nikjh
    利用稳定性的Fourier方法令，un，并将它代入上式就得到jve
    vn1eikjhvneikjh(eikj2eikjvn
    消去公因子有
    vn1[1(eikj2eikj]vnvn
    由此得到增长因子
    1(eikj
    1
    vn1ikj
    2e
    G(,k
    11

    1(eikj2eikj(12(eikjeikj

    11

    (122coskh12(1coskh
    11
    1
    |12(1coskh|14sin2kh
    2kh
    h0,k0,0且0sin21
    2

    khkh
    4sin2014sin21
    221
    1
    kh
    14sin2
    2
    显然这个格式是相容的。它是无条件稳定，因为按照Lax定理可知，该格式收敛的。
    |G(,k|
    1.3、Crank-Nicolson格式
    我们在在前面讨论的显格式和隐格式，即：
    n1
    unujj
    1n1n1
    un2uuj1jj1

    n1
    unjuj

    h
    2
    0（1.3.1）


    nn
    unj12ujuj1
    h2
    0（1.3.2）
    用乘（1.3.2），用(1乘（1.3.1），把其结果相加就得到一个差分格式
    6


    n1
    unjuj



    nn1n11
    unununj12ujuj1j12ujj1(10（1.3.3）22
    hh其中01，我们乘差分格式公式(1.3.3）为加权隐式格式。
    从上述可以看到，当
    1ununjj
    1
    时的情况，此时我们把它单独写出2



    2h
    n1n1n1nnn[(u2uu(u2uuj1jj1j1jj1]0（1.3.4）2
    此格式一般称作Crank-Nicolson格式。此外我们注意到，当1时，公式
    (1.3.3）为向后差分格式（隐式格式）；当0时，公式(1.3.3）为向前差分格式（显式格式）。
    1.3.1Crank-Nicolson差分格式的截断误差
    证：（用taylor展开）
    n1
    unu(x,t,uu(xj,tn1jjnj
    u
    u
    n1
    j
    un22un33un1
    u(x,ttu(xj,tn[]j[2]j[3]j
    t2!t3!t
    unh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j2
    x2!x3!x4!xunh22unh33unh44un
    u(xx,tu(xj,tnh[]j[2]j[3]j[4]j3
    x2!x3!x4!x
    (011
    (021(031
    n
    j1
    unj1
    把上述代入差分格式中，得截断误差为：
    un2un23un12un1h24un1
    k(x,t[]j[2]j[]j{([2]j[]j4
    t2t6t32x24x4
    h24un12unh24unh24un[4]j5([2]j[4]j2[4]j3}24xx24x24x
    nj
    un2un12un2un23un1
    {[]j([2]j[2]j}{{[2]j[3]j}
    t2xx2t6th24un14un14un4un
    {([4]j4[4]j5([4]j2[4]j3}}
    24xxxxo(2h2.(01,2,3,4,51
    22
    从上述可知，截断误差为kn(x,to(h，它对空间方向为二阶截断误j
    差，而对时间方向为二阶截断误差，则此隐格式的精度为2。
    7




    1.3.2Crank-Nicolson差分格式的稳定性
    证：由差分格式公式（1.3.4）可以写成如下形式

    1111nnn1n11un(1uuu(1uunj1jj1j1jj12222
    其中

    h2消去公因子有

    11
    eikj(1eikj
    2vn2vn1
    11
    eikj(1eikj
    22
    由此得到增长因子
    111
    eikj(1eikj(1(eikjeikj
    22G(,k2
    111eikj(1eikj(1(eikjeikj222
    kh
    12sin2
    1(1coskh2
    1(1coskh12sin2kh
    2
    h0,k0,0且12sin2kh
    21|G(,k|
    kh
    12sin2
    2
    12sin2
    khkh
    12sin222

    显然这个格式是相容的。它是无条件稳定，因为按照Lax定理可知，
    Crank-Nicolson差分格式收敛的。

    8




    2．高精度格式的构造
    本文热传导方程对空间变量应用四阶紧致格式离散，对时间变量应用梯形方
    24
    法，构造扩散方程的精度为Oh的绝对稳定的隐式差分格式

    2.1梯形方法[4]
    求解常微分方程初值问题
    dy
    f(x,ydx
    y(x0y0

    对方程
    dy
    f(x,y两边从xn到xn1积分，得dx
    xn1xn
    y(xn1yxn
    f(x,y(xdx
    (2.1
    用左矩形公式计算上式右侧积分，即

    xn1
    xn
    f(x,y(xdxhf(xn,y(xn
    并用yn作为的近似值y(xn，得
    yn1ynhfxn,yn(2.2）故欧拉法也称为矩形法。
    欧拉法形式简单,但精度低，为了达到较高精度的计算公式，对欧拉法进行改进，用梯形公式计算(2.1式右侧积分，即

    xn1
    xn
    f(x,y(xdx
    h
    fxn,yxnfxn1,yxn12
    并用yn作为y(xn的近似值，得到梯形公式:
    yn1yn
    h
    fxn,ynfxn1,yn12(2.3
    (2.3式称为梯形公式,此方法称为梯形法。


    9




    2.2本文格式的构造
    下面我们考虑非齐次边界条件的热传导方程
    u2u
    0,axb,t02tx
    axbu(x,0f(x,
    u(a,tg1(tt0
    t0u(b,tg2(t
    （2.2.1）的右边对x变量四阶紧致格式离散
    (2.2.1
    x22u
    2uio(h4
    12xh2(1x12(2.2.2
    把(2.2.2代入(2.2.1后得到下面的常微分方程组
    dui(t121222(1xxui，令A2(11x,Bx
    dt1212h

    dui(t
    A1BuiA1C(t(2.2.3dt
    T
    U(t[u1(t,u2(t,...,um2(t,um1(t]T,U(0[g(x1,g(x2,...,g(xm2,g(xm1].
    C(t[（-

    1'1'u0(t+(2u0(t），0，,0，（-uN(t+(2uN(t）]'12h12h

    561
    其中：A12

    11256
    112112

    21
    12B，

    56
    11
    2
    dU(t
    A1BU(t+A1C(t,
    应用梯形方法得到本文格式dt
    U(0U0
    un1(I

    2
    A1B1(I

    1
    A1Bun（IA1B1(A1C(tnA1C(tn12222
    （2.2.4）
    10




    2.3稳定性分析
    定理：本文差分格式（2.2.4）是绝对稳定的；
    引理1假设i(i1,2,...,m1为矩阵A1B的特征值，xi(i12,...,m1为其相应的是特征值向量，则特征值为实数且满足i0证明：
    561
    令A12

    1
    11256
    112112

    21
    12B

    56
    11
    2
    下面来看AB的特征值则：
    令ABXiXiAXBX两边乘X写成矩阵形式xAxxBx
    T
    T
    1

    T
    521521152x1x1x2x2x2x3xn1xnxn612612126521252121252222x1(x1x2x2(x2x3(xnxxn1n624624246192929292192x1x2x3xnxn012412121224xAx
    T
    22xBx2x1x1x22x2x1x222x1
    T
    2
    x1x22xn
    12121222222(x1x22x2(x2x3(x2x32xn2225252222x13x23x33xnxn1
    225252222
    (x13x23x33xnxn01
    22
    由此推出s0
    下面来看un1
    2u的特征值：
    11nIAB2
    11
    I

    A1B




    I
    令C

    2
    A1B

    I
    2
    A1B
    写成矩阵形式：un1Cun
    Ii
    2则C的特征值为：cun1Ii
    2它的谱半径(cmaxc1
    1in

    I
    因此un1


    2u是绝对稳定的。11nIAB2
    A1B
    12




    3数值实验
    数值例子1
    给出下面的常系数一维扩散方程初边值问题
    u2u
    t00x1tx2
    x
    0x1u(x,0e
    u(0,tet;u(1,te1t
    t0

    该方程的准确解为u(x,text

    表一：h=0.1;tao=0.1;T=1，绝对误差比较x隐式格式Crank-Nicolson格式本文格式
    8.1239E-032.7358E-041.3474E-040.1
    1.4847E-025.0077E-042.4818E-040.2
    2.0079E-026.7709E-043.3615E-040.3
    2.3696E-027.9890E-043.9680E-040.4
    2.5545E-028.6123E-044.2767E-040.5
    2.5433E-028.5770E-044.2560E-040.6
    2.3127E-027.8039E-043.8661E-040.7
    1.8348E-026.1935E-043.0559E-040.8
    1.0768E-023.6229E-041.7651E-040.9
    2.5545E-028.6123E-044.2767E-04最大误差

    表二：当=0.001,T=1;时不同空间步长的收敛界的比较隐式格式Crank-Nicolson格式本文格式
    h最大误差收敛阶最大误差收敛阶最大误差收敛阶1/42.9321e-0032.6757e-0038.3652e-0061/89.3104e-0041.65506.7288e-0041.99154.8370e-0074.11221/164.3004e-0041.11441.6966e-0042.00481.0174e-0084.5712

    表三:当h=0.005,w=1T=1;是的不同空间步长的收敛界的比较隐式格式Crank-Nicolson格式本文格式
    最大误差收敛阶最大误差收敛阶最大误差收敛阶1/102.5355e-0024.3558e-0044.3448e-0041/201.2863e-0020.97901.0980e-0041.98811.0871e-0041.99881/406.4786e-0030.98952.8270e-0051.95752.7182e-0051.9998

    13







    h=0.1;tao=0.1;T=1

    h=0.1;tao=0.1;T=1

    14






    h=0.1;tao=0.1;T=1
    数值例子2
    给出下面的常系数一维扩散方程初边值问题
    u2u
    t00x1tx2

    u(x,0sin(x0x1
    u(0,tu(1,t0t0
    该方程的准确解为u(x,tsin(xet
    表一：h=0.1;tao=0.01;T=1，绝对误差比较
    x0.10.20.30.40.50.60.70.80.9最大误差
    隐式格式1.1195E-052.1294E-052.9308E-053.4454E-053.6227E-053.4454E-052.9308E-052.1294E-051.1195E-053.6227E-05
    Crank-Nicolson格式
    1.2118E-062.3050E-063.1726E-063.7296E-063.9215E-063.7296E-063.1726E-062.3050E-061.2118E-063.9215E-06
    15
    2
    本文格式
    1.2133E-072.3079E-073.1765E-073.7343E-073.9264E-073.7343E-073.1765E-072.3079E-071.2133E-073.9264E-07




    h1/41/81/16

    1/101/201/40

    表二：当=0.001,T=1;时不同空间步长的收敛界的比较隐式格式Crank-Nicolson格式本文格式
    最大误差收敛阶最大误差收敛阶最大误差收敛阶3.7094e-0053.3294e-0058.3072e-0079.7900e-0061.92186.9518e-0062.25984.6769e-0084.15074.2936e-0061.18911.6599e-0062.06639.7749e-0105.5803

    表三:当h=0.005,w=1T=1;是的不同空间步长的收敛界的比较隐式格式Crank-Nicolson格式本文格式最大误差收敛阶最大误差收敛阶最大误差收敛阶9.9096e-0043.1582e-0053.1587e-0052.7642e-0041.84209.7016e-0061.70289.7107e-0061.70179.5859e-0051.52792.5387e-0061.93412.5489e-0061.9297

    图一：h=0.1;=0.01;T=1


    16





    结论

    由表1-3和图1可以看出本文针热传导方程构造出的一种两层绝对稳定的隐
    24
    式差分格式截断误差为Oh的，这个理论和数值试验相符合；
    从由表1-3数值试验结果说明本文所建立格式针对于，隐式差分格式和经典的Crank-Nicolson格式比，更为准确的逼近近似解。

    17





    致谢
    本毕业论文是在开依沙尔老师的指导下完成的，我首先要感谢开依沙尔老师，他对我进行了热心的指导并提出严格的要求，提供很多课外的资料，论文研究工作中提出了很多宝贵的意见，使我的研究工作有了目标和方向。在这近几个月的时间里，他对我进行了悉心的指导和教育，使我能够不断地学习提高，还要感谢同组论文设计人员，在他们熟悉部分对我的帮助，感谢父母养育之恩，感谢所有给我关心和帮助的同学和亲朋好友。

    18





    参考文献
    [1]陆金甫,关治.偏微分方程数值解法[M]，第二版.北京：清华大学出版社,2004
    [2]G.D.smithNumericalsolutionofpartialdifferentialequations(finitedifferencemethods[M]，Thirdedition,Oxford:OxfordUniversityPress.1996
    [3]李庆扬,王能超,易大义.数值分析[M],第四版,北京:清华大学出版社,2001年8月
    [4]李瑞霞，何志庆.微分方程数值方法[M],广州：华南理工大学出版社，2005年12月，38-43.[5]R.A.UsmaniandR.P.Agarwal,AnA-stableextendedtrapezoidalruleforthenumericalintegrationofordinarydifferentialequations[J].ComputersMath.Applic.1995，11(12：1183-1191.
    [6]JOHN.H.Mathews,著.陈渝等译.数值方法（MATLAB版）[M]，第三版.北京：电子工业出版社,2000
    [7]马明书,抛物型方程的一个新的显格式[J].纺织高校基础科学学报,2001年6月第14卷第2期,133-135.
    [8]周顺兴.解抛物型偏微分方程的高精度差分格式[J].计算数学,1982,4(2:204～213.
    [9]田振夫.非齐次热传导方程的高精度隐式格式[J].宁夏大学学报,1996,17(3:34～38.
    [10]马明书.一维抛物型方程的一个新的高精度显式差分格式[J].数值计
    算与计算机应用,2001,22(2:156～160.

    19





    新疆大学本科生毕业论文(设计评议书
    学院：数学与系统科学学院
    论文(设计题目:求解热传导方程的高精度隐式差分格式
    学生姓名:亚库甫江.买买提专业:信息与计算科学班级：2007-2指导教师姓名:开依沙尔职称:讲师
    评价内容
    具体要求
    能独立查阅文献和课题调研，能提出较科学、合理、可
    方案论证（15分）
    行得实施方案。
    坚持实事求是科学态度，没有造假和抄袭行为。观点、
    论文(设计内容（30
    结论正确、论证充分、设计合理。内容与专业要求相吻
    分）
    合，理论与实际联系紧密。
    工作量和难度（20分）
    遵守毕业论文（设计）管理制度，按期完成任务书规定的内容，工作量饱满，有一定难度。
    结构合理、条理清楚、文理通顺、用语符合专业要求；
    论文(设计质量（20
    文体格式规范、图表清楚。图样绘制与技术要求符合国
    分）
    家标准，图面质量符合要求。
    创新性与应用价值
    （15分）总分（100分）指导教师评语:
    具有一定的创新性和应用价值。
    11
    85
    得分
    13
    27
    15
    17
    20




    该论文讨论热传导方程的高精度数值解法，应用空间方向四阶紧致格式离散，时间方向梯形方法，构造热传导方程的精度为O2h4两层绝对稳定的隐式差分格式，讨论了稳定性，做了数值实例。与结果经典的隐式差分格式和Crank-Nicholson差分格式比较，数值验证该解法的可行性
    在写作过程中，该同学查阅了查看参考文献了解这方面的有关研究情况，做数值实例时用计算机编写了程序，锻炼了解决实际问题的能力。
    文章研究了一个较有实际意义的问题，在研究中思路清晰，结构严谨，对问题的讨论有一定的深度，但是还有很多问题有待进一步的解决，希望以后的学习中不断进步。
    指导教师（签名）：年月日

    新疆大学本科生毕业论文(设计评议书
    学院：数学与系统科学学院
    论文(设计题目:求解热传导方程的高精度隐式差分格式
    学生姓名:亚库甫江.买买提专业:信息与计算科学班级：2007-2
    21




    指导教师姓名:开依沙尔职称:讲师
    评价内容
    具体要求
    结构合理、条理清楚、文理通顺、用语符合专业要求，
    规范程度
    文体格式规范，图表清楚。图样绘制与技术要求符合
    （25分）
    国家标准，图面质量符合要求，资料齐全。
    论文(设计内容与质量（60分）创新性与应用价值
    具有一定的创新性和应用价值。
    （15分）总分（100分）
    评阅教师评语:
    亚库甫江.买买提同学的论文《求解热传导方程的高精度隐式差分格式》研究了一个较有理论意义的问题，锻炼了对偏微分方程数值解法的应用，深化了学习。
    从该论文可以看出，该同学比较好的掌握了专业课的基础知识，并能灵活地运用到解决实际问题上，具备了较好的独立思考、查阅文献、比较好的达到了毕业设计的要求。
    22
    得分
    21
    观点、结论正确，论证充分，设计合理。内容与专业要求相吻合，理论与实际联系紧密；查阅文献有一定广泛性；有综合归纳资料的能力，有自己的见解；
    55
    10
    86





    评阅教师（签名）：年月日



    新疆大学本科生毕业论文(设计评议书
    学院：数学与系统科学学院
    论文(设计题目:求解热传导方程的高精度隐式差分格式
    学生姓名:亚库甫江.买买提专业:信息与计算科学班级：2007-2指导教师姓名:开依沙尔职称:讲师
    评价内容论文（设计）水平
    （30分）
    具体要求
    论文（设计）内容正确，撰写规范、有一定的创新

    性和应用价值。
    得分
    论文（设计）介绍思路清晰，表达简明扼要，重点
    论文（设计）报告突出，能全面准确介绍论文（设计）内容，报告时
    间符合要求。
    23




    （25分）论文（设计）答辩
    （45分）总分（100分）
    答辩委员会（小组）评语:
    答辩小组认真听取了亚库甫江.买买提同学的《求解热传导方程的高精度隐式差分格式》的答辩，经过充分的讨论。认为亚库甫江.买买提同学在查阅大量文献并进行充分的课题调研后，独立完成了毕业设计工作。该设计在充分的理论学习和仔细阅读参考书和相关文献的基础上，构造出扩散方程的两层隐式格式，最后进行了数值实验。其特点如下：
    1.能独立查阅资料和课题调研，能提出较科学、合理、可行的实施方案。2.内容叙述较充分，工作量较饱满，具有一定的难度。
    3.文体格式规范，图标清楚，图样绘制与技术要求基本符合国家标准。
    答辩小组综合了指导教师、评阅教师的评分，最终作出如上成绩评定结果。
    答辩委员会（小组）签名：年月日

    回答问题正确，有理论依据，基本概念清楚，逻辑性较强。


    24




    指导教师成评阅教师成答辩指导小综合成绩
    论文(设计
    绩30%
    综合成绩


    五级分制成绩
    绩20%
    组成绩50%


    25

    免费下载 Word文档免费下载： 本科毕业设计--求解热传导方程的高精度隐式差分格式

    TOP热门搜索

    相关推荐：

    干部工作材料

    • ￥29.8

      ￥45 每天只需1.0元
      1个月 推荐

    • ￥9.9

      ￥15
      1天

    • ￥59.8

      ￥90
      3个月

    选择支付方式

    • 微信付款
    会员须知 联系在线客服 登录账号
    郑重提醒：支付后，系统自动为您完成注册

    请使用微信扫码支付(元)

    订单号：
    支付后，系统自动为您完成注册
    遇到问题请联系 在线客服

    常用手机号：
    用于找回密码
    图片验证码：
    看不清？点击更换
    短信验证码：
    新密码：
     
    绑定后可用手机号登录
    请不要关闭本页面，支付完成后请点击【支付完成】按钮
    遇到问题请联系 在线客服

    侵权投诉  © 2013-2019 http://www.tcm999.cn   站点地图 |  手机版

    本站文档均来自互联网及网友上传分享，本站只负责收集和整理，有任何问题可通过上访投诉通道进行反馈