有限差分法matlab两点边值代码,两点边值问题的有限差分法.doc

两点边值问题的有限差分法.doc

学生实验报告实验课程名称偏微分方程数值解开课实验室数统学院学院数统年级2013专业班信计2班学生姓名学号开课时间2015至2016学年第2学期总成绩教师签名数学与统计学院制开课学院、实验室数统学院实验时间2016年月日实验项目类型实验项目名称两点边值问题的有限差分法验证演示综合设计其他指导教师曾芳成绩是一．实验目的通过该实验，要求学生掌握求解两点问题的有限差分法，并能通过计算机语言编程实现。二．实验内容考虑如下的初值问题1,,DUXDUXLUPRQXFABX2,AB其中，，，，是1,PXCB,,RXQFXCABMIN0PXQX,给定常数。将区间等分，设，网点。NBHN,1,IXHIN1．在第三部分写出问题(1)和(2)的差分格式，并给出该格式的局部截断误差。2．根据你写出的差分格式，编写一个有限差分法程序。将所写程序放到第四部分。3．给定参数，，，，问题1的精确解0,1AB3,12PRQ01，其中将及带入方程1可得。分别取21XUE2XUE,3FX，用所编写的程序计算问题1和2。将数值解记为，0,4,86NIU，网点处精确解记为，。然后计算相应的误差1IIU1,N，及收敛阶，将计算结果填0MAXNICIEU120NIIEHU2LNNE入第五部分的表格，并对表格中的结果进行解释4将数值解和精确解画图显示，每种网格上的解画在一张图。三．实验原理、方法(算法)、步骤1差分格式1/H2/2H𝐿ℎ𝑢𝑖𝑝𝑖1/2𝑢𝑖1𝑝𝑖1/2𝑝𝑖‒1/2𝑢𝑖𝑝𝑖‒1/2𝑢𝑖‒1𝑟𝑖𝑢𝑖1‒𝑢𝑖‒1𝑞𝑖𝑢𝑖𝑓𝑖A,𝑢0𝑢𝑁𝐵2局部阶段误差UOH2𝑅𝑖3程序CLEARALLN10A0B1PX1RX2QX3ALPHA0BETA1FX4X22EXPX1HBA/NHZEROSN1,N1GZEROSN1,1FORI1HI,I2PAI1/2HPAI1/2H/H2HQAIHHI,I12PAI1/2H/HRAIHGI2HFAIH2PAI1/2H/HRAIHALPHAENDFORI2N2HI,I12PAI1/2H/HRAIHHI,I2PAI1/2HPAI1/2H/H2HQAIHHI,I12PAI1/2H/HRAIHGI2HFAIHENDFORIN1HI,I12PAI1/2H/HRAIHHI,I2PAI1/2HPAI1/2H/H2HQAIHGI2HFAIH2PAI1/2H/HRAIHBETAENDUH\GUALPHAUBETAXAHBYX2EXPX1PLOTX,UHOLDONPLOTX,YYY ZYU四．实验环境(所用软件、硬件等)及实验数据文件MATLAB五．实验结果及实例分析NNCE收敛阶0NE收敛阶100001042560000735242000002616819341000018348145304000000654120001000004585200008000000163619993000001146200001600000004092000000000028720000N越大只会使绝对误差变小，方法没变，所以收敛阶一致。图示为绿线为解析解，蓝线为计算解N10N20N40N80N160教师签名年月日